Next Article in Journal
Hydrodynamics Analysis of an Underwater Foldable Arm
Previous Article in Journal
Sequence Stratigraphy and Implications for Shale Gas Exploration in the Southern Sichuan Basin, South China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Progressive Trigonometric Mixed Response Surface Method for Double-Loop Interval Optimization

by
Guanchao Zheng
1,2,
Wen Xu
3,4,
Mohd Khairol Anuar Bin Mohd Ariffin
2,*,
Nuraini Binti Abdul Aziz
2 and
Siti Azfanizam Binti Ahmad
2
1
Department of Marine Engineering Equipment Technology, Faculty of Marine Engineering, Tianjin Maritime College, Tianjin 300350, China
2
Department of Mechanical and Manufacturing Engineering, Faculty of Engineering, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
3
Department of Economics and Management, Tianjin Bohai Vocational Technology College, Tianjin 300400, China
4
Department of Mechanical and Manufacturing Engineering, School of Business and Economics, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2023, 11(7), 1394; https://doi.org/10.3390/jmse11071394
Submission received: 10 May 2023 / Revised: 22 June 2023 / Accepted: 26 June 2023 / Published: 10 July 2023

Abstract

:
For some highly nonlinear problems, the general second-order response surface method (RSM) cannot satisfy the accuracy requirement. To improve accuracy, the highest order number has to be determined in advance. Thus, a progressive trigonometric mixed response surface method (PTMRSM) was proposed to enhance the approximation accuracy and define the highest order number, rather than determining it in advance. After that, a double-loop interval optimization process could be constructed using this PTMRSM to save time while maintaining accuracy when compared to other experimental or computational methods. Unfortunately, the traditional double-loop interval optimization method had issues with the probability of reliable constraints. Then, for the construction of this double-loop interval optimization process, the modified reliable constraints were introduced. A more reliable and effective double-loop interval optimization was introduced for addressing practical engineering problems using the effective approximate method of the PTMRSM and the amended reliable constraints. Two numerical test functions and a composite submersible hull were performed to verify the accuracy and effectiveness of the PTMRSM and the double-loop interval optimization.

1. Introduction

Regarding engineering design optimization, it is usual to assess structural behavior using deterministic parameters. In practice, however, the uncertainty associated with actual parameter values is ubiquitous, ranging from simple models to complex systems. Geometric dimensions, material properties, loads, boundary conditions, manufacturing tolerance, and so on are all inherently uncertain in real-world problems [1,2,3]. Recently, a probabilistic method based on a probability density function, or distribution function, has become widely used for dealing with uncertainties. Unfortunately, it seems inevitable that there will be precise probabilistic density functions or a wealth of probabilistic data for some complicated engineering problems. Furthermore, even a small, inaccurate deviation may lead to a large error in structural problems [4,5].
To resolve this problem, a method named interval analysis was proposed by Moore [6] in the 1960s. In this book, a new interval analysis is described as neglecting the premise of a probabilistic density function or huge amounts of probabilistic data. Since then, as an efficient method, interval analysis has been used to cope with a large number of structural problems. Elishakoff and Ben-Haim [7,8] may first have applied this interval analysis as a bounded, uncertain approach for structural engineering. Qiu [5,9] employed the anti-optimization technique to solve linear interval equations for interval static displacement bounds of structures.
Nonetheless, investigating only linear programming problems in all of these reviews was an incomplete analysis since most engineering problems are nonlinear. Therefore, a nonlinear interval number programming (NINP) method is introduced to transform the uncertain optimization into deterministic multi-objective optimization using the penalty function method and the first-order Taylor expansion method introduced by Jiang [10,11,12]. It is worth noting that a large perturbation deviation of variables may cause inaccurate results in highly nonlinear systems when using the first-order Taylor expansion method for interval analysis (See Figure 1). In Figure 1, it is shown that the variation ( Δ y ˜ L or Δ y ˜ U ) of the first-order Taylor expansion method is obviously larger than the variation ( Δ y L or Δ y U ) of the actual objective function.
To precisely obtain more reliable results, Li et al. [3,13] formulated a nested-loop optimization procedure for engineering design optimization. Zhao [14] proposed a double-loop optimization and rebuilt the objective and constraints at each iteration step. However, the computational cost of the nested-loop optimization algorithm is always prohibitive, whether using the finite element method (FEM) or the experimental method; this may cause several weeks or months to be spent on only one optimization design [15]. To overcome this limitation of expensive computing, a surrogate model that can replace FEM or an experimental method can be expressed. Wu [1,16] developed a novel interval method for nonlinear systems based on Chebyshev polynomial series to achieve tighter and sharper bounds. Cheng [17] and Li [13] integrated the Kriging approximate model and the nested genetic algorithm to improve computational efficiency. Zhao [14] and Cheng [15] formed radial basis functions (RBF) to solve the nested interval optimization for mechanical performance. All these surrogate-based optimization methods can reduce the computational time to a few hours for one optimization design [18,19,20].
Although these surrogate models can achieve a good simulation effect, more sample points should be added to get more accurate simulation results. Basudhar [21] built a precise approximated explicit decision function based on an adaptive sample process to update this decision function. Another new sequential improvement criterion was implemented to simultaneously get the robust optimization solution and advance the suitability of the RBF proposed by Havinga [22]. Cheng [15,23] employed an adaptive Kriging model for more effective computing results by adding more sample points. Wang [24] constructed an adaptive response surface method (ARSM) that required only a reduced number of design experiments for high-dimensional engineering design problems. Chakraborty [25] performed an adaptive RSM in an iterative finite element model updating procedure that proved efficient. It is important to note that more sample points mean not only a more accurate response but also a higher computational cost. For signal classification problems, a weighted fuzzy process neural network model is proposed, and the accuracy of this method is higher than other existing automatic identification methods by S. Xu et al. [26].
Lots of methods for obtaining more precise approximations with fewer sample points have been developed and implemented. In order to construct a more efficient and accurate approximating model, Kim [27] and Youn [28] combined sensitivity information and the moving least squares method (MLSM). Kaymaz [29] and Taflanidis [30] implement a response surface model using the moving least squares method (MLSM) for the reduction of the computational burden. Li [31] proposed a doubly weighted moving least squares method, including the normal weight factor of MLSM and the distance between the most probable failure points. S. Lee [32] proposed a new approach toward reliability-based design optimization using the response surface augmented moment method (RSMM) for the progressive update of the response surface. Even though the aforementioned scholars did their best to establish an accurate surrogate model with a small number of sample points, sensitive data or other information is still needed, which can still increase computing costs.
Consequently, the trigonometric mixed response surface (TMRSM) is proposed as an updated response surface method (RSM) for obtaining perfect approximate responses in the same sample points [33]. Specifically, most of the surrogate model was constructed by the second-order polynomial [24,27,28,29,31]. Only the quadratic RSM model may be appropriate for engineering performance in some highly nonlinear cases. The conventional approach is to artificially select the order of the terms in advance, then establish this RSM function. Few papers have conducted more extensive research on the determination of high nonlinear programming for RSM. Gavin [34] introduced an algorithm process that can iteratively decide the fittest order polynomial. Along this path, this paper improved the trigonometric mixed response surface (TMRSM) to the progressive trigonometric mixed response surface (PTMRSM), which can automatically determine the order number of the surrogate model using the t-test method. Following that, a double-loop interval optimization can be built using the PTMRSM and reliable constraints. Although the possibility degree may be first introduced by C. Jiang [10] in the treatment of the interval inequality constraints, this paper fixed a minor problem with respect to this possibility degree and proposed reliable constraints in accordance with the amended possibility degree for the double-loop interval optimization.
This paper will be organized as follows. In Section 2, a novel progressive response surface method (PTMRSM) is introduced to replace the more time-consuming finite element method (FEM) and determine the order number automatically. The accuracy of the PTMRSM is verified by two general numerical test functions and one practical submersible example in Section 3. Then, a new double-loop interval optimization is proposed to improve the accuracy of interval boundaries using reliable constraints in Section 4. Finally, in Section 5, a practical submersible composite hull design is carried out to verify the efficiency of the double-loop interval optimization. The conclusion is in Section 6.

2. Progressive Trigonometric Mixed RSM (PTMSRM)

2.1. Identification of Response Surface Method

The response surface method (RSM) was originally proposed as a simple mathematical expression to express the relationship between the input and output of an experiment. Then, this effective method may be first introduced by Box and Hunter [35] and extended to other fields, especially engineering problems. The RSM function can be presented as:
y x = y ^ x + ε
in which y x expresses the output function; y ^ x illustrates the polynomial form of x = x 1 , x 2 , , x n T ; ε indicates the experimental residual error, assuming that it is a normal distribution with E ε = 0 and V a r ε = σ 2 .
In general, approximate polynomial functions can be modeled as linear, second-order, or higher-order polynomials. Then, the common formula of the RSM model can be indicated as:
y = β 0 + i = 1 n β i x i + i = 1 n j > i n β i j x i x j + i = 1 n k = 2 K β i k x i k + ε
where β = β 0 , β 1 , , β l T is the actual vector of the coefficient.

2.2. Moving Least Squares Method (MLSM)

In the traditional least squares method (LSM), the response surface function is formed as:
y ^ = p x T β ^
where p x = 1 , p 1 x , p 2 x , , p l x T is a vector of l polynomial basis functions, and β ^ = β ^ 0 , β ^ 1 , β ^ 2 , , β ^ l T is an evaluated vector of the coefficient. It is noted that these evaluated coefficients are constant with global regression, while x changes. That is to say, the whole sample points are considered and weighted equally. However, in practical application, a precise approximated model requires greater weight on the sample points closer to the prediction point. Consequently, a weighted regression method called the moving least squares method (MLSM) is proposed to improve the coefficients β ^ . This implies that the coefficients β ^ x become functions of x . And the coefficients β ^ x can be calculated by minimizing the residual error between the output value and the response surface value.
To construct this function, m sample points of variables in the design space must be chosen based on optimal Latin hypercube design (OLHD). The output value can be calculated from these OLHD sample points using computer simulation, such as the finite element method (FEM). Once the required sample data was obtained, the polynomial constructed for all sample points can be shown as:
y = P β ^ + ε
where P can be defined as:
P = p T x 1 p T x 2 p T x m m × l + 1
After that, the experimental residual error between the output value and the response surface value can be expressed as:
L x = ε T W x ε = y p x T β ^ x T W x y p x T β ^ x
where can W x be defined as:
W x = w x x 1 0 0 0 w x x 2 0 0 0 w x x m
In Equation (6), the weight w x x I can be presented as:
w x x I = i f x x I / D I 1 L i n e a r 1 x x I / D I Q u a d r a t i c 1 x x I / D I 2 E x p o n e n t i a l e α x x I / D I G a u s s i a n e α x x I / D I 2 e α 2 1 e α 2 o t h e r w i s e 0
where x x I is the Euclid distance between the prediction point and the sample point, and D I is the influence range whose appropriate size should be selected among a sufficient number of neighboring sample points to avoid singularity in the solution. It is useful to remember that D I should contain at least l sample points, and the weight function should vanish outside the D I influence range. And α is the coefficient of the exponential and the Gaussian weight function. Among all the weight functions, the Gaussian function is considered to be an effective weight by most scholars for MLSM [30,31].
Figure 2 shows the shape of different weight functions using Equation (8). From Figure 2, the exponential weight function and the Gaussian weight function may have a better weight effect both in the domain range close to the sample point and away from the sample point. Nevertheless, using the exponential weight function, the weight value is not zero near the boundary of the domain area. In this regard, the Gaussian weight function may have a more suitable performance for the MLSM compared to the exponential weight function.
Figure 3 demonstrates the shape of the Gaussian weight function under different coefficients α . Specifically, the value of this coefficient suggested α = 2.5 has a nice distribution characteristic from Figure 3 [30].
The minimum of the experimental residual error L x can be gained by the partial derivatives with respect to coefficients, as below:
L x β ^ = 0
Then, the evaluated coefficient can be achieved from Equation (9):
β ^ x = P T W x P 1 P T W x y
Therefore, the contribution of the sample points will be decreased as the distance between the prediction point x and x I is increased. This implies the MLSM can denote both a local and global approximation over the entire region.

2.3. Trigonometric Mixed Response Surface Method

Normally, the base function of the response surface method is polynomial. However, if the surrogate relationship has considerable curvature and is extremely complicated, only high-order nonlinear programming may fail to meet the accuracy requirement. For this condition, Lee and Lin [36,37] proposed a novel type of RSM to solve these complex composite problems using trigonometric functions as the basic functions. In this method, the trigonometric functions ( sin x and cos x ) are chosen as the basic elements to replace the traditional polynomial functions. Then, the vector of the trigonometric base functions can be presented as:
p x = [ 1 , sin x 1 , sin x 2 , , sin x n , cos x 1 , cos x 2 , , cos x n , sin x 1 sin x 2 , sin x 1 sin x 3 , , sin x 1 sin x n , sin x 2 sin x 3 , , sin x n 1 sin x n , sin x 1 cos x 1 , sin x 1 cos x 2 , , sin x 1 cos x n , , sin x 2 cos x 1 , , sin x 2 cos x n , , sin x n cos x n , cos x 1 cos x 2 , cos x 1 cos x 3 , , cos x 1 cos x n , cos x 2 cos x 3 , cos x n 1 cos x n , sin x 1 2 , sin x 2 2 , , sin x n 2 , cos x 1 2 , , cos x n 2 , sin x 1 3 , sin x 2 3 , , sin x n 3 , cos x 1 3 , , cos x n 3 , , sin x 1 k , sin x 2 k , , sin x n k , cos x 1 k , , cos x n k ]
Therefore, a trigonometric mixed response surface method (TMRSM) is proposed to integrate the trigonometric functions and the MLSM based on OLHD [33].

2.4. Progressive Trigonometric Mixed Response Surface Method

Usually, most RSM models in the literature take second-order polynomials as their basic form. TMRSM, though, has successfully improved the accuracy of the surrogate model, and the determination of the higher-order polynomial should be tested in advance for the given sample data. To circumvent this deficiency, a progressive trigonometric mixed response surface method (PTMRSM) is proposed in this paper to approximate these sample points and test the necessity of the higher-order terms using a t-test.
According to Equations (4) and (10), the experimental residual error can be written as:
ε = y P β ^ = y P P T W x P 1 P T W x y = I m P P T W x P 1 P T W x y
If H = P P T W x P 1 P T W x , then the residual error can be abbreviated to:
ε = I H y
Hence, the mathematical expectation and the variance can be expressed as:
E ε = E y P β ^ = y P E β ^ = y P β = 0
V a r ε = V a r I H y = σ 2 I H I H T
where β is the expectation of the evaluated coefficient E β ^ . For the MLSM, the expectation of the estimated coefficient is the actual coefficient for the sample points. If c = c 1 , c 2 , , c n T R n , the variance of c T y is defined as V a r c T y = c T I c σ 2 . In which the variance of y is indicated as V a r y = σ 2 .
Furthermore, the expectation of the sum of the squared residual error is:
E i = 1 m ε i 2 = i = 1 m E ε i 2 0 = i = 1 m E ε i 2 E 2 ε i = i = 1 m V a r ε i = σ 2 t r I H I H T
If the trace is simplified as s = t r I H I H T , there is the variance of y :
E i = 1 m ε i 2 = i = 1 m V a r ε i = σ 2 s
Therefore, the unbiased estimation of σ 2 is:
σ ^ 2 = 1 s i = 1 m ε j 2
The squared residual error can comply with ε i ~ N 0 , σ 2 , and then it can be written as ε i σ 2 ~ N 0 , 1 . Hence, this squared residual error obeys the Chi-square distribution:
ε i 2 σ 2 ~ χ 2 1
According to Equations (18) and (19), the summation of the square residual error should abide by χ 2 s :
i = 1 m ε i 2 σ 2 = s σ ^ 2 σ 2 ~ χ 2 s
From Equation (10), the mathematical expectation of this evaluated coefficient can be expressed as:
E β ^ = E P T W x P 1 P T W x y = P T W x P 1 P T W x E y = P T W x P 1 P T W x P β = β
So, the variance of the evaluated coefficient can be indicated as:
V a r β ^ = C o v β ^ , β ^ = C o v P T W x P 1 P T W x y , P T W x P 1 P T W x y = σ 2 P T W x P 1 P T W x P T W x P 1 P T W x T = σ 2 P T W x P 1 P T W x W x P P T W x P 1
Make C j j = d i a g P T W x P 1 P T W x W x P P T W x P 1 , then we have β ^ j ~ N β j , σ 2 C j j . And a standard normal distribution can be achieved:
β ^ j β j σ C j j ~ N 0 , 1
Combining Equations (20) and (23), the t distribution can be represented as:
β ^ j β j σ C j j i = 1 m ε i 2 s σ 2 = β ^ j β j C j j 1 s i = 1 m ε i 2 ~ t s
where is a t-distribution of the s degree of freedom.
Determining the significance of the higher-order terms means measuring the importance of these evaluated coefficients. This indicates that the insignificant coefficients stand for the unimportant higher-order terms. Equation (10) shows the vector of the coefficients, and it is easy to find the coefficient values β ^ j of the highest order. Then, the focus of this work becomes to check whether the highest coefficients have significant effects on the output response. It is equivalent to whether the highest order coefficients β ^ j are equal to the values of 0 (0 shows the insignificant impact). Assume a null hypothesis H 0 : β ^ j = 0 , and the test can be performed by calculating the t-test value.
t j = β ^ j 0 S d j ~ t s
where t s indicates t distribution of s degree of the freedom (see Equation (24)); S d is the standard deviation of the coefficients; j is the number of the terms of the model, which can be selected from the first number of the highest-order terms to the last number of the highest-order terms.
S d j 2 = C j j 1 s i = 1 m ε i 2
For this t-test, a two-sided test and 90% confidence intervals are usually chosen in this paper [34]. Then, the boundary value t 0.9 can be calculated by t-distribution functions. If t j > t 0.9 , the null hypothesis H 0 can be rejected. This implies these coefficients β ^ j 0 and the relevant high-order terms may have a significant influence. Otherwise, the null hypothesis H 0 cannot be dismissed. That is to say, the coefficients β ^ j = 0 and these high-order terms may not have a substantial effect.
In addition to the t-test criterion, other criteria should exist to ensure the relationships in which the odd or even order terms have a significant effect on the response results. Only testing whether the highest order has an important effect may easily lead to inaccurate approximations of these special functional relationships. To avoid this deficiency, the determination of the coefficient R 2 and the mean relative error M R E should be introduced for the m sample points, respectively.
R 2 = 1 S S E S S T
M R E = 1 m i = 1 m y i y ^ i y i
where S S E represents the sum of squares for error, which is caused by experimental error. And S S T expresses the sum of squares for the total, representing the total variation of values of m sample observations. These two values can be shown as:
S S E = i = 1 m y i y ^ i 2
S S T = i = 1 m y i y ¯ 2
where y ¯ is the mean value of output response y , which can be expressed as:
y ¯ = 1 2 i = 1 m y i
The determination of the coefficient is a statistical indicator that reflects the reliability of the response surface models to account for changes in terms of dependent variables. The closer this value is to 1, the better the approximated effect will be. If the determination of the coefficient R 2 > 0.99 for the m sample points, it means that the approximation degree is sufficient [38].
The mean relative error is another indicator that can illustrate the relative error between the actual response and the approximate response. The lesser the value M R E is, the smaller the error is. For most practical engineering problems, it can be viewed as an appropriate limit when M R E 0.02 [38].
Eventually, another stopping criterion is the highest order limit. Given the possibility of over-fitting for the approximate functions, this model takes a stopping criterion with respect to the sixth highest order. That is due to the fact that the highest terms are really rare to exceed the fifth order for most engineering problems [34].
The proposed PTMRSM can be started with the order k = 3 , and the program can determine whether the t-test values of the third-order terms are sufficient for the 90% confidence intervals. If not, the determination of the coefficient and the mean relative error will be utilized to check the reliability of the fitting performance and the relative error. If these criteria are satisfied, stop this process, and output the response. If not, continue this process until the highest order reaches the sixth order. Then, the response of PTMRSM will be outputted, displaying that “Warning: the highest order has reached the maximum value (6th order)”. Figure 4 shows the step-by-step procedure of this algorithm:
  • Perform the OLHD sample points, and obtain the response with respect to these sample points based on FEM or specific functions;
  • Transform all the sample points and the input parameters x into the trigonometric functions sin x and cos x . Then, choose order k = 3 as the beginning, and construct this third-order TMRSM;
  • Calculate the values of the t-test criterion t j . Meanwhile, the determination of the coefficient R 2 and the mean relative error M R E can also be calculated from Equations (28) and (29);
  • Check the t-test criterion and decide whether this operation should be stopped. If the coefficients t j t 0.9 , it is concluded that the highest-order terms are insignificant. Then, stop the process, and construct the model with (k − 1)-th order TMRSM;
  • If the t-test criterion is not satisfied, check the criteria for the determination of the coefficient R 2 0.99 and the mean relative error M R E 0.02 . If these criteria are fulfilled, it implies the performance and the accuracy of the approximated degree are sufficient. Then, stop this operation and build the model with k-th order TMRSM;
  • Otherwise, select (k + 1)-th order as the next step order, and check whether the highest order is k = 6 satisfied. If not, repeat steps 3 to step 6 until one of these criteria is satisfied;
  • When the highest order k = 6 is fulfilled, establish the model with the sixth order TMRSM, “Warning: the highest order number has reached the maximum value (6th order)”.
The pseudocode for the PTMRSM process is listed in the Algorithm 1.
Algorithm 1: The pseudocode for the PTMRSM process:
1:Read input variable matrix and output response matrix:
2:load Sample.mat;
3:Set the confidence level of t-test criterion, the determination coefficient criterion and
4:the mean relative error criterion.
5:Transform input variables and sample matrix into trigonometric functions sin θ and cos θ.
6:for sn = 1:Number of response column
7:  for k = 3:6
8:   Construct the X matrix with MLSM using k order
9:   Evaluated coefficient: b = (P.′ ∗ W ∗ P)\(P.′ ∗ W ∗ Y);
10:   Calculate errors between actual response values and predicted values
11:   Calculate the value of t-test criterion:
12:    H= P *(P.′ ∗ W ∗ P)\(P.′ ∗ W);
13:    s = trace((eye(m)-H) ∗ (eye(m)-H).’);
14:    SSE = e.’ ∗ e; %e is the error between actual value and evaluated value
15:    C_jj = diag((P.′ ∗ W ∗ P)\(P.′ ∗ W ∗ W ∗ P)/(P.′ ∗ W ∗ P));
16:    s_d = (SSE./s). ∗ C_jj;
17:    tj = b./(s_d).^0.5;
18:   Calculate the determination coefficient criterion value: R2 = 1-SSE./SST;
19:   Calculate the mean relative error criterion value: MSE = mean(abs(e./Y));
20:   if (t < =tinv(confidence level,s) ∗ ones(highest order number,1))
21:    k = k − 1;
22:    break
23:   end
24:   if (R2 > = 0.99)&(MSE < = 0.02)
25:    k = k;
26:    break
27:   end
28:  end
29:  Adopt the new order and reconstruct the TMRSM model using updated order
30:  Output the response value:
31:    f (sn) = p ∗ b;
32:    h_order(sn) = k;
33:    sn = sn + 1;
34:end

3. Accuracy Analysis

3.1. Numerical Rastrigin Function

The first numerical function is the well-known Rastrigin function, which can be stated as:
f x = 20 + x 1 2 + x 2 2 10 cos 2 π x 1 + cos 2 π x 2 , 1 x 1 1 , 1 x 2 1
In this Rastrigin function, 40 sample points were designed using optimal LHD, and 100 random test points were adopted for the verification of each surrogate model. Figure 5 shows the shape function or responses of the original function, second-order RSM with MLSM, second-order TMRSM and PTMRSM under 40 sample points.
It is obvious that the PTMRSM model has exactly the same shape as the original model from Figure 5a,d. Namely, the PTMRSM model has the best approximate performance in comparison with the other surrogate models. And Table 1 represents the performance and the accurate measures of each surrogate model under 100 random test points. From Table 1, it is further confirmed that the PTMRSM model has a better fitting performance and smaller mean relative error compared to other surrogate models.
For the PTMRSM model, Figure 6 denotes the number of its order, which can be decided by the step-by-step procedure for 100 random test points. It can be seen that the number of orders for the PTMRSM model is third and fourth order. That implies that the highest order number is not a fixed value for PTMRSM.

3.2. Numerical Example Function

The next numerical function is an example function, which has been utilized in the literature [34]. This function can be indicated below.
f x = 0.16 x 1 1 3 x 2 + 4 0.04 cos x 1 x 2 5 x 1 5 , 5 x 2 5
In this example function, there are also 40 sample points that were designed based on optimal LHD and 100 random test points that were taken to verify each surrogate model. Figure 7 expresses the shape function or responses of the original function, second-order RSM with MLSM, second-order TMRSM, and PTMRSM under 40 sample points. In this Figure, it is evident that both the PTMRSM model and the original function have almost the same shape of the response.
Table 2 supports the same point of view as the Rastrigin function because the PTMRSM model also has the largest R 2 and the least M R E .
Figure 8 represents the order number of the PTMRSM model through 100 random test points. As seen, for these test points, all the highest order numbers of the constructed PTMRSM models are third order. And this means the higher order terms can improve the performance and accuracy.

3.3. Composite Submersible Hull Example

In this paper, a composite submersible hull was put forward to withstand underwater pressure, whose type is the well-known Myring, as a practical example [39]. For this composite structure optimization, the stacking sequence is one of the non-trivial factors in terms of buckling performance. The conventional stacking sequence can be designed in a symmetrical form. This paper takes the symmetric type stacking sequence with a total of 10 layers, and critical buckling pressure P c r is considered the pivotal design objective. Another design constraint can be the maximum buckling displacement d max for the underwater pressure. Consequently, the orientation angle of each layer can be defined as the designed vector θ and the thickness of the layer can be considered an uncertain factor t because of the manufacturing tolerance. Figure 9 shows the diagram of the composite submersible hull.
Table 3 displays the performance and the accuracy of different RSM models under 100 random test points with respect to the critical buckling pressure P c r and the maximum buckling displacement d max , respectively.
From this table, it is obvious that the PTMRSM model has more precise results and better fitting responses than other RSM models for both the critical buckling pressure and the maximum buckling displacement.
Figure 10 and Figure 11 display the highest order distribution under 100 testing points for the critical buckling pressure and the maximum buckling displacement, respectively. For this practical submersible hull, the third-highest order is adequate in terms of both the buckling performance and the deformation performance.

4. Double-Loop Interval Optimization Method

4.1. Interval Optimization Method

Based on the interval analysis method [11,40,41], the real numbers are bounded by intervals that the interval vector can be expressed in the following form:
α I = α _ , α ¯ = α C α W , α C α W = α C + 1 , 1 α W
α C = α _ + α ¯ 2 , α W = α ¯ α _ 2
where α _ and α ¯ denote the lower and upper bounds, respectively; α C and α W show the center point vector and radius vector, respectively. Considering the interval vectors, the objective and constraint functions will be constructed as lower and upper intervals, respectively:
f _ x , α = min α Γ f x , α , f ¯ x , α = max α Γ f x , α , Γ = α α _ α α ¯
g _ x , α = min α Γ g x , α , g ¯ x , α = max α Γ g x , α , Γ = α α _ α α ¯
where x is the design variables, and α is the uncertain interval parameters.
C. Jiang performed the first-order Taylor expansion method, which is a simple approximated method for the evaluation under relatively small uncertainty [10,11,42,43].
f _ x , α = f x , α C j = 1 q f x , α C α j α j W f ¯ x , α = f x , α C + j = 1 q f x , α C α j α j W
g _ x , α = g x , α C j = 1 q g x , α C α j α j W g ¯ x , α = g x , α C + j = 1 q g x , α C α j α j W
Although this Taylor expansion method is feasible and convenient, inaccurate bounds can be obtained for high nonlinear complex structural problems under large perturbation deviations. Thus, a double-loop interval optimization method can be proposed using sequential quadratic programming (SQP) by J. Wu and F. Li [2,13,44].
F i n d x , α ; M i n i m i z e f I x , α S u b j e c t   t o x _ x x ¯ , α I = α _ , α ¯
F i n d x , α ; M a x i m i z e f I x , α S u b j e c t   t o x _ x x ¯ , α I = α _ , α ¯
In these equations, f _ x , α = M i n i m i z e f I x , α and f ¯ x , α = M a x i m i z e f I x , α . Similarly, the lower bound g _ x , α and upper bound g ¯ x , α of optimization constraints can also be written as:
F i n d x , α ; M i n i m i z e g I x , α S u b j e c t   t o x _ x x ¯ , α I = α _ , α ¯
F i n d x , α ; M a x i m i z e g I x , α S u b j e c t t o x _ x x ¯ , α I = α _ , α ¯
With reference to Equation (35), the interval objective and constraint functions can be described as their center points and radii, respectively:
f C x = f _ x , α + f ¯ x , α 2 f W x = f ¯ x , α f _ x , α 2
g C x = g _ x , α + g ¯ x , α 2 g W x = g ¯ x , α g _ x , α 2
Therefore, the above interval optimal problem can be transformed into the deterministic two-objective problem below:
min x Ω f C x , f W x g x , α b I = b _ , b ¯ Ω = x | x _ x x ¯

4.2. Reliable Constraints

To compare two interval numbers, a certain degree named possibility degree is introduced, which can convert two interval numbers into a deterministic relationship. This possibility degree may be first introduced by C. Jiang to express which interval values are larger or smaller than the other ones [28,45,46]. Figure 12 shows the whole six positional cases with respect to these two interval values A and B . In order to describe this deterministic relationship, the possibility degrees of the cases in Figure 12 can be displayed for ranking these two interval numbers, A and B , as:
P A B = 1 , A ¯ B _ 1 1 2 A ¯ B _ A ¯ A _ A ¯ B _ B ¯ B _ , A _ B _ < A ¯ B ¯ B ¯ A ¯ B ¯ B _ + 1 2 A ¯ A _ B ¯ B _ , B _ < A _ < A ¯ B ¯ 1 2 B ¯ A _ B ¯ B _ B ¯ A _ A ¯ A _ , B _ < A _ B ¯ < A ¯ B _ A _ A ¯ A _ + 1 2 B ¯ B _ A ¯ A _ , A _ B _ < B ¯ < A ¯ 0 B ¯ < A _
For Case 2 (in Figure 12), there exist four events with respect to the two random parameters a and b :
H 1 : a A _ , B _ b B _ , A ¯ ; H 2 : a A _ , B _ b A ¯ , B ¯ ; H 3 : a B _ , A ¯ b B _ , A ¯ ; H 4 : a B _ , A ¯ b A ¯ , B ¯
However, when the event H 3 occurs, three possible events ( A < B , A = B and A > B ) may be encountered. An assumption that all these events have equal chance was made by P. Sevastianov [47]. That means when the event H 3 occurs, the possibility is:
P A < B | H 3 = P A = B | H 3 = P A > B | H 3 = 1 3
Then, the conditional possibilities can be expressed as:
P A B | H 1 = 1 ; P A B | H 2 = 1 ; P A B | H 3 = 2 3 ; P A B | H 4 = 1
Consequently, the possibility of the event A B can be indicated as:
P A B = 1 1 3 A ¯ B _ 2 A ¯ A _ B ¯ B _
In Case 3 of Figure 12, there are also three possible events: A < B , A = B and A > B (see Figure 13). Define three events in terms of the random parameters a and b in given subintervals when Case 3 appears:
H 1 : a A _ , A ¯ b B _ , A _ ; H 2 : a A _ , A ¯ b A _ , A ¯ ; H 3 : a A _ , A ¯ b A ¯ , B ¯ .
Thus, the possibilities of A < B , A = B and A > B under given event H 2 can be expressed as:
P A < B | H 2 = P A = B | H 2 = P A > B | H 2 = 1 3
Then, the whole conditional possibilities of A B given H 1 , H 2 and H 3 can be indicated as:
P A B | H 1 = 0 ; P A B | H 2 = 2 3 ; P A B | H 3 = 1 .
Therefore, the probability degree of the event A B can be attained as:
P A B = B ¯ A ¯ B ¯ B _ + 2 3 A ¯ A _ B ¯ B _
This is a different value from Cases 2 and 3 of Equation (45). It implies C. Jiang only makes the assumption of the events B > A and B = A , rather than the assumption of the event of B < A . This assumption, however, is at odds with reality. Similarly, the events A B in Case 4 and 5 can be illustrated, respectively, as:
P A B = 2 3 B ¯ A _ 2 A ¯ A _ B ¯ B _
P A B = B _ A _ A ¯ A _ + 1 3 B ¯ B _ A ¯ A _
So the amended possibility degree can be improved by integrating Equations (47), (51) and (55)–(57).
On the other hand, the possibility degree under the opposite rank of the two interval numbers A and B can be shown:
P A B = 0 , A ¯ B _ 2 3 A ¯ B _ B ¯ B _ A ¯ B _ A ¯ A _ , A _ < B _ < A ¯ < B ¯ A _ B _ B ¯ B _ + 2 3 A ¯ A _ B ¯ B _ , B _ A _ < A ¯ B ¯ 1 1 3 B ¯ A _ A ¯ A _ B ¯ A _ B ¯ B _ , B _ < A _ < B ¯ < A ¯ A ¯ B ¯ A ¯ A _ + 2 3 B ¯ B _ A ¯ A _ , A _ B _ < B ¯ A ¯ 1 , A _ B ¯
When there is an event A = B , the possibility degrees of Cases 1 and 6 are zero. Additionally, this possibility degree is one only if A _ = B _ and A ¯ = B ¯ . For the condition of Case 2 and 4, the possibility degrees can be displayed, respectively, as:
P A = B = A ¯ B _ B ¯ A _
P A = B = B ¯ A _ A ¯ B _
Analogously, for Cases 3 and 5, the possibility degrees can be written, respectively, as:
P A = B = A ¯ A _ B ¯ B _
P A = B = B ¯ B _ A ¯ A _
Thus, the whole of possibility degree under the rank of event A = B can be expressed as:
P A = B = 1 , A _ = B _ & A ¯ = B ¯ A ¯ B _ / B ¯ A _ , A _ < B _ < A ¯ < B ¯ A ¯ A _ / B ¯ B _ , B _ < A _ < A ¯ < B ¯ B ¯ A _ / A ¯ B _ , B _ < A _ < B ¯ < A ¯ B ¯ B _ / A ¯ A _ , A _ < B _ < B ¯ < A ¯ 0 O t h e r w i s e
It should be mentioned that P A < B = λ means the possibility degree level is λ when A is smaller than B . λ = 1 means that A is always smaller than B , and λ = 0 represents that is A always larger than B . Under this amended possibility degree method, the reliable constraints should fulfill a confidence degree level, and these reliable constraints are transformed into simple deterministic constraints for the double-loop interval optimization. After the premise of the possibility degree has been ensured, the simple deterministic constraints can be expressed as:
P G i I < B i I λ i , i = 1 , 2 , , n g
where G i I = g _ i , g ¯ i means the interval values of constraints; and B i I = B _ i , B ¯ i represents the relationship with regard to the constraints interval; n g is the number of the constraints. And G i I can be calculated by Equations (42) and (43).

4.3. Double-Loop Interval Optimization Method

For the sake of dealing with the two-objective problem (see Equation (46)), a linear combination of the two objectives is introduced, which can be generally defined as
F = β f C x , α φ + 1 β f W x , α ϕ
where β is the weight values of each objective function; φ and ϕ are the normalization factors that can make f C x , α and f W x , α the same magnitude. It has to emphasize that the interval optimization problem may be either a minimization problem or a maximization problem. In this paper, the objective is a maximization issue. Consequently, the linear combination of objectives and reliable constraints can be written as:
min x Ω , α Γ β f W x , α f C x , α + 1 β f W x , α ϕ s . t . P g I x , α < B I λ 1 , P g I x , α < B I λ 2 , B I = B _ , B ¯ , α I = α _ , α ¯ , x _ x x ¯
where ϕ is the normalization factor that can make f W x , α f C x , α and f W x , α as the same magnitude; λ 1 and λ 2 are the lower and upper bounds of the preset possibility degree.
In order to perform this double-loop interval optimization process, the bounds of the interval objectives and constraints should be evaluated in advance. Therefore, an inner-loop optimization should be presented to solve this min-max problem. Note that sequential quadratic programming (SQP) [3,48,49] and particle swarm optimization (PSO) [49,50] can be selected as the inner and outer optimization solvers, respectively. Figure 14 shows the flowchart of the double-loop interval optimization method.
It is clear that this flowchart is a nesting double-loop interval optimization problem. The samples of design and uncertain space are performed using OLHD, and the response of these sample points can be obtained by ANSYS simulation. For the sake of saving computational cost, the aforementioned PTMRSM surrogate model can be formed as the basic objective and constraint functions instead of the more time-consuming ANSYS simulation. In the inner loop interval optimization process, SQP can be executed to compute the bounds of the inner loop objective and constraint with respect to uncertain interval vectors α . Then, during the outer loop interval optimization process, PSO is carried out to optimize the linear combination objective with reliable constraints for the optimal design vector x .

5. Practical Composite Submersible Hull Design

5.1. The Objective Function of the Submersible Hull Example

As mentioned above, the critical buckling pressure P c r can be regarded as the inner loop objective function; meanwhile, the maximum buckling displacement d max should also be another significant effect that can be considered the inner loop constraint. In this submersible optimization design, the larger the buckling pressure is, the better the optimal solution is. Thus, the center point of the critical buckling P c r C θ , t can be treated as the maximization problem. To circumvent this problem, P c r W θ , t / P c r C θ , t can be considered as one of the linear combination functions instead of P c r C θ , t . The following single-objective function with a penalty function can be formed based on the double-loop interval optimization method.
min x Ω , α Γ β P c r W θ , t P c r C θ , t + 1 β P c r W θ , t ϕ s . t . P g I x , α < B I λ 1 , P g I x , α < B I λ 2 , B I = 50 , 60 mm , t I = 0.5 , 0.7 mm , P o u t e r I = 11 , 9 × 10 6 P a , 90 ° θ 90 °
where P c r C θ , t and P c r W θ , t are center point and a radius value of the critical buckling pressure; P g I x , α < B I is the possibility degree of the reliable constraint. β is the weight coefficient, which can be selected as 0.5, since the center point and its radius should be minimized simultaneously. ϕ is the normalization factor with the value of 2.0 × 105. B I is the defined maximum buckling displacement interval. t I is the layer thickness interval. P o u t e r I is the outer pressure interval. λ 1 and λ 2 are the mean preset possibility degree level, which can be selected between 0 and 1. The larger λ is, the stricter the constraints will be and may lead to a narrower feasible region [10]. In this paper, a common preset possibility level can be defined as:
λ 1 = 0.4 , λ 2 = 0.6 N e u t r a l l e v e l λ 1 = 0.6 , λ 2 = 0.8 S l i g h t l e v e l λ 1 = 0.8 , λ 2 = 1.0 S t r o n g l e v e l λ 1 = 1.0 , λ 2 = 1.0 A b s o l u t e l e v e l

5.2. Results and Analysis

Table 4 shows the optimal design results of the double-loop interval optimization for the 3rd order RSM with MLSM and the PTMRSM, and the nonlinear interval number programming (NINP) method for the PTMRSM model under neutral level possibility degree. The NINP method was proposed by C. Jiang [11,42]. And this NINP method can be treated as a comparative item.
From this table, it can be seen that the double-loop interval optimization for PTMRSM has the smallest radius of the critical buckling pressure, which means the deviation of the optimal design solution can be least affected by the uncertain variables. In other words, the interval optimization for PTMRSM can obtain a more robust and reliable optimal design solution. Furthermore, the similarity between the optimal design solutions of the interval optimization for PTMRSM and the NINP method for TMRSM is greater than that of the interval optimization for 3rd-order RSM with MLSM.
Figure 15 shows the lower and the upper bounds of the critical buckling pressure for one variable orientation angle and the other four fixed orientation angles under the neutral possibility degree. Figure 15a indicates the lower and the upper bounds of the critical buckling pressure under one variable ( θ 1 ) and four invariants ( θ 2 , θ 3 , θ 4 and θ 5 ). Figure 15b expresses the bounds of the critical buckling pressure under one variable ( θ 2 ) and four invariants ( θ 1 , θ 3 , θ 4 and θ 5 ). Figure 15c displays the same bounds under one variable ( θ 3 ) and four invariants ( θ 1 , θ 2 , θ 4 and θ 5 ). Figure 15d shows the same bounds under one variable ( θ 4 ) and four invariants ( θ 1 , θ 2 , θ 3 and θ 5 ).
From these figures, it should be noted that the deviation of the bounds of the interval optimization for PTMRSM is shaper than the deviations of the bounds of the other two optimization methods.
Table 5, Table 6 and Table 7 express the optimal design solutions under the slight level, the strong level, and the absolute level possibility degree.
From Table 4, Table 5, Table 6 and Table 7, it is obvious that the center point of the critical buckling pressure becomes smaller with the increment of the possible degree. This phenomenon may be due to the fact that the increasing possibility degree limits the space of the optimal solutions. Among the whole optimal design solutions, it can be seen that the interval optimization for PTMRSM always has the smallest radius of the critical buckling pressure. This implies that the interval optimization for PTMRSM can attain robust and efficient optimal design solutions.
Figure 16, Figure 17 and Figure 18 display the lower and the upper bounds of the critical buckling pressure under the slight level, the strong level, and the absolute level possibility degree. From these whole figures, it is clear that the results of the interval optimization for PTMRSM also have the smallest deviations of the bounds compared to the other two optimization methods. It means that the double-loop interval optimization method with PTMRSM has enough efficiency and robustness in contrast to the same interval optimization method with 3rd-order RSM with MLSM and the NINP optimization method.

6. Conclusions

In this paper, a progressive trigonometric mixed response surface method (PTMRSM) was introduced to increase the accuracy of approximation for some high nonlinear engineering problems. The results of two numerical examples and a composite submersible hull were utilized to verify the efficiency and accuracy. For these highly nonlinear complex problems, the PTMRSM may have a better fitting performance and avoid the previous determination of the highest order number of polynomial terms.
Then, the PTMRSM is utilized as one approximate step of the novel double-loop interval optimization method. With the help of the PTMRSM and the amended reliable constraints, the novel double-loop interval optimization process can be constructed, and a composite submersible hull is also performed to verify this new interval optimization method compared with the NINP method and other traditional RSM with the MLSM model. From the comparison results of the composite submersible hull, it can be seen that the double-loop interval optimization for PTMRSM can obtain a robust and efficient optimal design solution.
Furthermore, the center point of the critical buckling pressure becomes smaller with the increase in the possibility degree. This phenomenon may be due to the fact that the increasing possibility degree limits the space of the optimal solutions.
Among the whole optimal design solutions, it can be seen that the interval optimization for PTMRSM always has the smallest radius and deviations of the critical buckling pressure. This implies the double-loop interval optimization for PTMRSM can attain a robust and efficient optimal design solution.

Author Contributions

Conceptualization, G.Z. and W.X.; methodology, G.Z.; software, G.Z.; validation, W.X., M.K.A.B.M.A., N.B.A.A. and S.A.B.A.; writing—original draft preparation, G.Z.; writing—review and editing, M.K.A.B.M.A., N.B.A.A. and S.A.B.A.; supervision, M.K.A.B.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This paper did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors wish to thank Soren Ebbesen for sharing the particle swarm optimization code.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationship that could have appeared to influence the work reported in this paper.

References

  1. Wu, J.; Zhang, Y.; Chen, L.; Luo, Z. A Chebyshev Interval Method for Nonlinear Dynamic Systems under Uncertainty. Appl. Math. Model. 2013, 37, 4578–4591. [Google Scholar] [CrossRef]
  2. Wu, J.; Luo, Z.; Zhang, N.; Zhang, Y. A New Interval Uncertain Optimization Method for Structures Using Chebyshev Surrogate Models. Comput. Struct. 2015, 146, 185–196. [Google Scholar] [CrossRef]
  3. Li, F.; Luo, Z.; Sun, G.; Zhang, N. An Uncertain Multidisciplinary Design Optimization Method Using Interval Convex Models. Eng. Optim. 2013, 45, 697–718. [Google Scholar] [CrossRef]
  4. Chen, S.H.; Wu, J. Interval Optimization of Dynamic Response for Uncertain Structures with Natural Frequency Constraints. Eng. Struct. 2004, 26, 221–232. [Google Scholar] [CrossRef]
  5. Qiu, Z.; Elishakoff, I. Anti-Optimization Technique—A Generalization of Interval Analysis for Nonprobabilistic Treatment of Uncertainty. Chaos Solitons Fractals 2001, 12, 1747–1759. [Google Scholar] [CrossRef]
  6. Moore, R.E.; Kearfott, R.B.; Cloud, M.J. Introduction to Interval Analysis; SIAM: Philadelphia, PA, USA, 2009; Volume 110, ISBN 089871771X. [Google Scholar]
  7. Elishakoff, I.; Haftka, R.T.; Fang, J. Structural Design under Bounded Uncertainty—Optimization with Anti-Optimization. Comput. Struct. 1994, 53, 1401–1405. [Google Scholar] [CrossRef]
  8. Ben-Haim, Y. Convex Models of Uncertainty: Applications and Implications. Erkenntnis 1994, 41, 139–156. [Google Scholar] [CrossRef]
  9. Qiu, Z.; Elishakoff, I. Antioptimization of Structures with Large Uncertain-but-Non-Random Parameters via Interval Analysis. Comput. Methods Appl. Mech. Eng. 1998, 152, 361–372. [Google Scholar] [CrossRef]
  10. Jiang, C.; Han, X.; Liu, G.R.; Liu, G.P. A Nonlinear Interval Number Programming Method for Uncertain Optimization Problems. Eur. J. Oper. Res. 2008, 188, 1–13. [Google Scholar] [CrossRef]
  11. Jiang, C.; Han, X.; Guan, F.J.; Li, Y.H. An Uncertain Structural Optimization Method Based on Nonlinear Interval Number Programming and Interval Analysis Method. Eng. Struct. 2007, 29, 3168–3177. [Google Scholar] [CrossRef]
  12. Jiang, C.; Zhang, Z.G.; Zhang, Q.F.; Han, X.; Xie, H.C.; Liu, J. A New Nonlinear Interval Programming Method for Uncertain Problems with Dependent Interval Variables. Eur. J. Oper. Res. 2014, 238, 245–253. [Google Scholar] [CrossRef]
  13. Li, F.; Luo, Z.; Rong, J.; Zhang, N. Interval Multi-Objective Optimisation of Structures Using Adaptive Kriging Approximations. Comput. Struct. 2013, 119, 68–84. [Google Scholar] [CrossRef]
  14. Zhao, Z.; Han, X.; Jiang, C.; Zhou, X. A Nonlinear Interval-Based Optimization Method with Local-Densifying Approximation Technique. Struct. Multidiscip. Optim. 2010, 42, 559–573. [Google Scholar] [CrossRef]
  15. Cheng, J.; Duan, G.F.; Liu, Z.Y.; Li, X.G.; Feng, Y.X.; Chen, X.H. Interval Multiobjective Optimization of Structures Based on Radial Basis Function, Interval Analysis, and NSGA-II. J. Zhejiang Univ. Sci. A 2014, 15, 774–788. [Google Scholar] [CrossRef] [Green Version]
  16. Wu, J.; Luo, Z.; Zhang, Y.; Zhang, N.; Chen, L. Interval Uncertain Method for Multibody Mechanical Systems Using Chebyshev Inclusion Functions. Int. J. Numer. Methods Eng. 2013, 95, 608–630. [Google Scholar] [CrossRef]
  17. Cheng, J.; Liu, Z.; Wu, Z.; Tang, M.; Tan, J. Direct Optimization of Uncertain Structures Based on Degree of Interval Constraint Violation. Comput. Struct. 2016, 164, 83–94. [Google Scholar] [CrossRef]
  18. Xia, H.; Wang, L.; Liu, Y. Uncertainty-Oriented Topology Optimization of Interval Parametric Structures with Local Stress and Displacement Reliability Constraints. Comput. Methods Appl. Mech. Eng. 2020, 358, 112644. [Google Scholar] [CrossRef]
  19. Song, X.; Lv, L.; Sun, W.; Zhang, J. A Radial Basis Function-Based Multi-Fidelity Surrogate Model: Exploring Correlation between High-Fidelity and Low-Fidelity Models. Struct. Multidiscip. Optim. 2019, 60, 965–981. [Google Scholar] [CrossRef]
  20. Parnianifard, A.; Azfanizam, A.S.; Ariffin, M.K.A.; Ismail, M.I.S. Crossing Weighted Uncertainty Scenarios Assisted Distribution-Free Metamodel-Based Robust Simulation Optimization. Eng. Comput. 2020, 36, 139–150. [Google Scholar] [CrossRef]
  21. Basudhar, A.; Missoum, S. Adaptive Explicit Decision Functions for Probabilistic Design and Optimization Using Support Vector Machines. Comput. Struct. 2008, 86, 1904–1917. [Google Scholar] [CrossRef]
  22. Havinga, J.; van den Boogaard, A.H.; Klaseboer, G. Sequential Improvement for Robust Optimization Using an Uncertainty Measure for Radial Basis Functions. Struct. Multidiscip. Optim. 2017, 55, 1345–1363. [Google Scholar] [CrossRef] [Green Version]
  23. Cheng, J.; Liu, Z.; Wu, Z.; Li, X.; Tan, J. Robust Optimization of Structural Dynamic Characteristics Based on Adaptive Kriging Model and CNSGA. Struct. Multidiscip. Optim. 2015, 51, 423–437. [Google Scholar] [CrossRef]
  24. Wang, G.G. Adaptive Response Surface Method Using Inherited Latin Hypercube Design Points. J. Mech. Des. Trans. ASME 2003, 125, 210–220. [Google Scholar] [CrossRef] [Green Version]
  25. Chakraborty, S.; Sen, A. Adaptive Response Surface Based Efficient Finite Element Model Updating. Finite Elem. Anal. Des. 2014, 80, 33–40. [Google Scholar] [CrossRef]
  26. Xu, S.; Feng, N.; Liu, K.; Liang, Y.; Liu, X. A Weighted Fuzzy Process Neural Network Model and Its Application in Mixed-Process Signal Classification. Expert Syst. Appl. 2021, 172, 114642. [Google Scholar] [CrossRef]
  27. Kim, C.; Wang, S.; Choi, K.K. Efficient Response Surface Modeling by Using Moving Least-Squares Method and Sensitivity. AIAA J. 2005, 43, 2404–2411. [Google Scholar] [CrossRef]
  28. Youn, B.D.; Choi, K.K. A New Response Surface Methodology for Reliability-Based Design Optimization. Comput. Struct. 2004, 82, 241–256. [Google Scholar] [CrossRef]
  29. Kaymaz, I.; McMahon, C.A. A Response Surface Method Based on Weighted Regression for Structural Reliability Analysis. Probabilistic Eng. Mech. 2005, 20, 11–17. [Google Scholar] [CrossRef]
  30. Taflanidis, A.A.; Cheung, S. Stochastic Sampling Using Moving Least Squares Response Surface Approximations. Probabilistic Eng. Mech. 2012, 28, 216–224. [Google Scholar] [CrossRef]
  31. Li, J.; Wang, H.; Kim, N.H. Doubly Weighted Moving Least Squares and Its Application to Structural Reliability Analysis. Struct. Multidiscip. Optim. 2012, 46, 69–82. [Google Scholar] [CrossRef]
  32. Lee, S. Reliability Based Design Optimization Using Response Surface Augmented Moment Method. J. Mech. Sci. Technol. 2019, 33, 1751–1759. [Google Scholar] [CrossRef]
  33. Zheng, G.; Ariffin, M.K.A.B.M.; Ahmad, S.A.B.; Abdul Aziz, N.B.; Xu, W. Robust Optimization of Composite Cylindrical Shell by A Trigonometric Mixed Response Surface Method. Int. J. Numer. Methods Eng. 2022, 123, 5010–5035. [Google Scholar] [CrossRef]
  34. Gavin, H.P.; Yau, S.C. High-Order Limit State Functions in the Response Surface Method for Structural Reliability Analysis. Struct. Saf. 2008, 30, 162–179. [Google Scholar] [CrossRef]
  35. Box, G.E.P.; Hunter, J.S. Multi-Factor Experimental Designs for Exploring Response Surfaces. Ann. Math. Stat. 1957, 28, 195–241. [Google Scholar] [CrossRef]
  36. Lee, Y.J.; Lin, C.C. Regression of the Response Surface of Laminated Composite Structures. Compos. Struct. 2003, 62, 91–105. [Google Scholar] [CrossRef]
  37. Lin, C.C.; Lee, Y.J. Stacking Sequence Optimization of Laminated Composite Structures Using Genetic Algorithm with Local Improvement. Compos. Struct. 2004, 63, 339–345. [Google Scholar] [CrossRef]
  38. Zheng, G.; Ariffin, M.K.A.B.M.; Ahmad, S.A.B.; Aziz, N.B.A.; Xu, W. A Novel Progressive Response Surface Method for High-Order Polynomial Metamodel. J. Phys. Conf. Ser. 2022, 2224, 12045. [Google Scholar] [CrossRef]
  39. Zheng, G.; Yang, X. Studies of the Resistance Optimization of Underwater Vehicle Based on Multiple-Speed Approximate Model. Proceedings of 2018 2nd International Conference on Functional Materials and Chemical Engineering (ICFMCE 2018), Abu Dhabi, United Arab Emirates, 20–22 November 2018; Volume 272, p. 1029. [Google Scholar]
  40. Muscolino, G.; Sofi, A. Stochastic Analysis of Structures with Uncertain-but-Bounded Parameters via Improved Interval Analysis. Probabilistic Eng. Mech. 2012, 28, 152–163. [Google Scholar] [CrossRef]
  41. Wu, J.; Zhang, Y.; Chen, L.; Chen, P.; Qin, G. Uncertain Analysis of Vehicle Handling Using Interval Method. Int. J. Veh. Des. 2011, 56, 81–105. [Google Scholar] [CrossRef]
  42. Jiang, C.; Han, X.; Liu, G.P. A Sequential Nonlinear Interval Number Programming Method for Uncertain Structures. Comput. Methods Appl. Mech. Eng. 2008, 197, 4250–4265. [Google Scholar] [CrossRef]
  43. Jiang, C.; Han, X.; Liu, G.R. Optimization of Structures with Uncertain Constraints Based on Convex Model and Satisfaction Degree of Interval. Comput. Methods Appl. Mech. Eng. 2007, 196, 4791–4800. [Google Scholar] [CrossRef]
  44. Wu, J.; Luo, Z.; Zhang, Y.; Zhang, N. An Interval Uncertain Optimization Method for Vehicle Suspensions Using Chebyshev Metamodels. Appl. Math. Model. 2014, 38, 3706–3723. [Google Scholar] [CrossRef]
  45. Shin, Y.S.; Grandhi, R.V. A Global Structural Optimization Technique Using an Interval Method. Struct. Multidiscip. Optim. 2001, 22, 351–363. [Google Scholar] [CrossRef]
  46. Wang, L.; Xiong, C.; Wang, X.; Xu, M.; Li, Y. A Dimension-Wise Method and Its Improvement for Multidisciplinary Interval Uncertainty Analysis. Appl. Math. Model. 2018, 59, 680–695. [Google Scholar] [CrossRef]
  47. Sevastianov, P. Numerical Methods for Interval and Fuzzy Number Comparison Based on the Probabilistic Approach and Dempster-Shafer Theory. Inf. Sci. 2007, 177, 4645–4661. [Google Scholar] [CrossRef]
  48. Victoire, T.A.A.; Jeyakumar, A.E. Hybrid PSO-SQP for Economic Dispatch with Valve-Point Effect. Electr. Power Syst. Res. 2004, 71, 51–59. [Google Scholar] [CrossRef]
  49. Zhou, J.; Cheng, S.; Li, M. Sequential Quadratic Programming for Robust Optimization with Interval Uncertainty. J. Mech. Des. 2012, 134, 100913. [Google Scholar] [CrossRef]
  50. Javidrad, F.; Nazari, M.; Javidrad, H.R. Optimum Stacking Sequence Design of Laminates Using a Hybrid PSO-SA Method. Compos. Struct. 2018, 185, 607–618. [Google Scholar] [CrossRef]
Figure 1. The diagram of the first-order Taylor expansion estimate.
Figure 1. The diagram of the first-order Taylor expansion estimate.
Jmse 11 01394 g001
Figure 2. The shape of different weight functions.
Figure 2. The shape of different weight functions.
Jmse 11 01394 g002
Figure 3. The shape of Gaussian weight function under different coefficients.
Figure 3. The shape of Gaussian weight function under different coefficients.
Jmse 11 01394 g003
Figure 4. The flowchart of the PTMRSM.
Figure 4. The flowchart of the PTMRSM.
Jmse 11 01394 g004
Figure 5. The shape of different responses for Rastrigin.
Figure 5. The shape of different responses for Rastrigin.
Jmse 11 01394 g005
Figure 6. The order number of PTMRSM for Rastrigin.
Figure 6. The order number of PTMRSM for Rastrigin.
Jmse 11 01394 g006
Figure 7. The shape of different responses for an example function.
Figure 7. The shape of different responses for an example function.
Jmse 11 01394 g007
Figure 8. The order number of PTMRSM for six hump camelback.
Figure 8. The order number of PTMRSM for six hump camelback.
Jmse 11 01394 g008
Figure 9. The diagram of the composite submersible hull.
Figure 9. The diagram of the composite submersible hull.
Jmse 11 01394 g009
Figure 10. The order number of PTMRSM for buckling pressure.
Figure 10. The order number of PTMRSM for buckling pressure.
Jmse 11 01394 g010
Figure 11. The order number of PTMRSM for buckling displacement.
Figure 11. The order number of PTMRSM for buckling displacement.
Jmse 11 01394 g011
Figure 12. Six positional cases between two interval values.
Figure 12. Six positional cases between two interval values.
Jmse 11 01394 g012
Figure 13. The diagram of possible events for Case 3. (a) A < B , (b) A = B and (c) A > B .
Figure 13. The diagram of possible events for Case 3. (a) A < B , (b) A = B and (c) A > B .
Jmse 11 01394 g013
Figure 14. The flowchart of double-loop interval optimization.
Figure 14. The flowchart of double-loop interval optimization.
Jmse 11 01394 g014
Figure 15. The lower and the upper bounds under neutral possibility degree.
Figure 15. The lower and the upper bounds under neutral possibility degree.
Jmse 11 01394 g015
Figure 16. The lower and the upper bounds under slight possibility degree.
Figure 16. The lower and the upper bounds under slight possibility degree.
Jmse 11 01394 g016
Figure 17. The lower and the upper bounds under strong possibility degree.
Figure 17. The lower and the upper bounds under strong possibility degree.
Jmse 11 01394 g017
Figure 18. The lower and the upper bounds under absolute possibility degree.
Figure 18. The lower and the upper bounds under absolute possibility degree.
Jmse 11 01394 g018
Table 1. The accuracy of different models for Rastrigin.
Table 1. The accuracy of different models for Rastrigin.
ItemModelR2MRE
b2nd-order RSM with MLSM0.10540.7639
c2nd-order TMRSM0.68590.3394
dPTMRSM0.99630.0275
Table 2. The accuracy of different models for an example function.
Table 2. The accuracy of different models for an example function.
ItemModelR2MRE
b2nd-order RSM with MLSM0.97111.6317
c2nd-order TMRSM0.99870.3321
dPTMRSM0.99990.0757
Table 3. The different models for critical buckling pressure and maximum displacement.
Table 3. The different models for critical buckling pressure and maximum displacement.
ItemModelPcrdmax
R2MRER2MRE
b2nd-order RSM with MLSM0.90830.17150.40760.3166
c2nd-order TMRSM0.94580.10040.91920.1326
dPTMRSM0.94640.09880.92870.1313
Table 4. Optimal solution under the neutral possibility degree.
Table 4. Optimal solution under the neutral possibility degree.
Type of FunctionInterval Optimization for 3rd RSM with MLSMInterval Optimization for PTMRSMNINP for PTMRSM
Optimal solutionθ146.571.55−16.55
θ27.7314.745.45
θ326.6632.9515.64
θ433.3786.56−4.52
θ5−27.87−90.00−90.00
PcrC1.9400 × 1061.6820 × 1061.4592 × 106
PcrW7.7297 × 1054.4201 × 1057.0716 × 105
gC0.05610.05400.0538
gW0.00500.00490.0058
P(gI < BI)0.40000.59930.6000
Table 5. Optimal solution under the slight possibility degree.
Table 5. Optimal solution under the slight possibility degree.
Type of FunctionInterval Optimization for 3rd RSM with MLSMInterval Optimization for PTMRSMNINP for PTMRSM
Optimal solutionθ143.63−14.86−21.83
θ25.172.4811.33
θ323.628.94−8.98
θ430.48−4.8419.22
θ5−26.60−90.00−90.00
PcrC1.8943 × 1061.4092 × 1061.4130 × 106
PcrW7.7191 × 1053.5448 × 1057.4074 × 105
gC0.05390.05080.0519
gW0.00510.00640.0075
P(gI < BI)0.60000.79700.7040
Table 6. Optimal solution under the strong possibility degree.
Table 6. Optimal solution under the strong possibility degree.
Type of FunctionInterval Optimization for 3rd RSM with MLSMInterval Optimization for PTMRSMNINP for PTMRSM
Optimal solutionθ138.92−16.13−12.77
θ23.951.31−4.64
θ319.5614.8411.70
θ427.9216.6816.12
θ5−24.45−89.82−90.00
PcrC1.8283 × 1061.3162 × 1061.3674 × 106
PcrW7.6970 × 1053.4002 × 1056.7439 × 105
gC0.80000.04730.0462
gW0.05120.00520.0056
P(gI < BI)0.00530.96990.9853
Table 7. Optimal solution under the absolute possibility degree.
Table 7. Optimal solution under the absolute possibility degree.
Type of FunctionInterval Optimization for 3rd RSM with MLSMInterval Optimization for PTMRSMNINP for PTMRSM
Optimal solutionθ118.77−12.68−15.05
θ211.150.560.02
θ38.7310.309.76
θ421.9821.7135.82
θ5−15.89−90.00−85.33
PcrC1.6591 × 1061.2902 × 1061.3251 × 106
PcrW7.6861 × 1053.3537 × 1056.9319 × 105
gC0.04430.04350.0435
gW0.00570.00560.0065
P(gI < BI)1.00001.00001.0000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zheng, G.; Xu, W.; Mohd Ariffin, M.K.A.B.; Abdul Aziz, N.B.; Ahmad, S.A.B. A Progressive Trigonometric Mixed Response Surface Method for Double-Loop Interval Optimization. J. Mar. Sci. Eng. 2023, 11, 1394. https://doi.org/10.3390/jmse11071394

AMA Style

Zheng G, Xu W, Mohd Ariffin MKAB, Abdul Aziz NB, Ahmad SAB. A Progressive Trigonometric Mixed Response Surface Method for Double-Loop Interval Optimization. Journal of Marine Science and Engineering. 2023; 11(7):1394. https://doi.org/10.3390/jmse11071394

Chicago/Turabian Style

Zheng, Guanchao, Wen Xu, Mohd Khairol Anuar Bin Mohd Ariffin, Nuraini Binti Abdul Aziz, and Siti Azfanizam Binti Ahmad. 2023. "A Progressive Trigonometric Mixed Response Surface Method for Double-Loop Interval Optimization" Journal of Marine Science and Engineering 11, no. 7: 1394. https://doi.org/10.3390/jmse11071394

APA Style

Zheng, G., Xu, W., Mohd Ariffin, M. K. A. B., Abdul Aziz, N. B., & Ahmad, S. A. B. (2023). A Progressive Trigonometric Mixed Response Surface Method for Double-Loop Interval Optimization. Journal of Marine Science and Engineering, 11(7), 1394. https://doi.org/10.3390/jmse11071394

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop