Home Technology Fast analysis approach for instability problems of thin shells utilizing ANNs and a Bayesian regularization back-propagation algorithm
Article Open Access

Fast analysis approach for instability problems of thin shells utilizing ANNs and a Bayesian regularization back-propagation algorithm

  • T. N. Nguyen EMAIL logo , Dongsheng Zhang and Pairod Singhatanadgid EMAIL logo
Published/Copyright: June 24, 2024
Become an author with De Gruyter Brill

Abstract

This research develops a data-driven methodology for structural instability problems with highly nonlinear, difficult, noisy, and small data. A fast analysis and prediction (FAP) approach for instability problems of thin shells is first proposed. This approach contains two phases: the fast numerical analysis and the pure prediction utilizing artificial neural networks (ANNs) incorporated with the Bayesian regularization (B-R) algorithm as follows: (1) in Phase 1 (the fast numerical analysis), post-buckling analysis is conducted utilizing a minor amount of load steps. The load–displacement relation achieved from Phase 1 is not exact because of the small number of load steps utilized; (2) in Phase 2 (the prediction), the loads and deflections achieved from Phase 1 were employed as the data for training ANNs. The trained networks, including the load and displacement networks, were employed to fast predict loads and deflections at any step of the post-buckling analysis. After utilizing Phase 2, a smooth, complete and exact load–displacement curve was achieved. In Phase 1, the available formulation for post-buckling analysis of thin shells in the literature was utilized. Five popular types of instabilities chosen to confirm the effectiveness and exactness of the FAP were snap-through, snap-back, softening–hardening, kink instabilities, and delamination buckling and post-buckling of composites. The high exactness and effectiveness of the FAP were confirmed in the numerical verification section. The present approach saves a huge computation compared to the other ones. It was found that ANNs incorporated with the B-R algorithm have notable advantages compared to numerous neural networks. The proposed approach is applicable to simulations or experiments where data are “expensive”, highly nonlinear, difficult, and limited. Utilizing the proposed approach for these fields can dramatically save time and money.

1 Introduction

A structure under a compressive load can present an instability phenomenon. When the compressive load is large enough or reaches the buckling load, the structure can reduce its stiffness, experience a notable change in geometry, and become unstable. When instability occurs, the structure can reduce its load-carrying capacity and is incapable of maintaining a stable equilibrium configuration. Investigations into structural instability and post-buckling behavior are necessary and play an important role in structural engineering, especially when post-buckling behavior is unpredictable. Load–displacement equilibrium paths are usually employed to evaluate structural instability. Several approaches, such as an analytical approach, a semi-analytical approach, and a numerical approach, can be employed to obtain those equilibrium paths. The numerical approach is convenient and powerful. It is applicable to problems with complex loading, boundaries, and geometries. Therefore, it is popularly employed in both academic and industrial problems. Utilizing the numerical approach for solving an instability problem needs a complete combination of an appropriate formulation, a numerical technique, and a nonlinear solver. One of the difficulties of solving instability problems utilizing the numerical approach is high time-consuming because of a huge computation and numerous computational iterations. This reality is clearly seen in large-scale problems utilizing thousands or millions degrees of freedom. Another difficulty in solving instability problems is divergent in some cases, especially when the equilibrium path passes the limit or special points, the inflection points. To overcome the aforementioned difficulties, a lot of approaches, nonlinear solvers, and techniques were proposed to decrease the number of computational iterations and to successfully pass limit or special points on the equilibrium paths including the Newton–Raphson (N-R) technique [1], the modified N-R [2], the modified Riks (M–R) technique [3], the iterative techniques upon an optimization technique [4] and upon residual areas [5], the Koiter (K) technique [6], the accelerated K technique [7], the Koiter-Newton (K-N) technique [8,9], a new technique upon the discretization of governing equations [10], a general path-following methodology [11], the methodology that improves effectiveness of the Newton technique [12], the dynamic relaxation techniques [13,14], a technique without utilizing the predictor [15], and the multi-point techniques [16]. Although numerous nonlinear solvers, approaches, and techniques have been developed and proposed, high time-consuming is still a difficulty in solving instability problems utilizing the numerical approach. To completely overcome this disadvantage, the FAP for instability problems of thin shells is proposed in this work. The present approach includes two phases: the fast analysis utilizing pure isogeometric analysis (IGA) and the prediction utilizing ANNs. In Phase 1, a post-buckling analysis is conducted utilizing a minor amount of load steps. In Phase 2, the obtained data from Phase 1 are utilized to train ANNs to fast and completely predict post-buckling behavior of the structure for many load steps. Interestingly, the computational cost of the FAP is dramatically lower than that of other ones.

ANNs are computing systems simulated according to the operation of biological neurons in the human brain. An ANN is a group of connected neurons. An ANN can mimic a human brain to fast and simply resolve complex problems that possess nonlinear relationships [17]. Up to now, ANNs have been successfully used in many areas, including bridge deck structure [18], optimal material structure [19], aircraft design [20], prediction of the fracture energy of concrete [21], buckling capacity assessment of steel structures [22], compressive strength prediction of cement [23], an optimization for hydrogen purification performance [24], an optimization for fatigue prediction [25], structural damage identification [26], and damage detection [27]. For nonlinear problems utilizing the numerical methodology, ANNs were employed to improve an algorithm for nonlinear dynamic analysis [28] and the Newton iterative technique [29]. Recently, the FAP for nonlinear bending analysis of structures was proposed [30], in which, N-R technique was applied to solve the nonlinear equation system in Phase 1. In the study by Nguyen et al. [30], the approach was proposed and investigated for pure bending problems of plates/shells. As known, loads are priorly known in the bending problems, while displacements are unknown. Thus, ANNs incorporated with the B-R algorithm were only employed to predict series of displacements in Phase 2 in the study by Nguyen et al. [30]. In this study, the FAP is proposed for instability problems of thin shells. As known, solving instability problems is much difficult than solving nonlinear bending problems because the equilibrium path can pass the limit or special points, the inflection points. Solving instability problems can be divergent in some cases. Thus, in this study, a powerful nonlinear solver named the modified Riks technique is applied to solve the nonlinear equation system in Phase 1. Besides different from the nonlinear bending problems, loads and displacements are both unknown in the instability problems. Thus, ANNs incorporated with the B-R algorithm are employed to predict the series of displacements and loads in Phase 2. It is emphasized that the present approach employs a minor amount of data points in the training phase. Thus, the time for training a network is very small and can be neglected in evaluating the effectiveness of the proposed approach. The idea of this work is simple in its implementation and can dramatically save computation compared to the another approach. The size of the dataset in this work is relatively small. Both deep neural networks (DNNs) including multiple hidden layers and artificial neural networks (ANNs) including one hidden layer can be used for training. However, we chose ANNs combined with Bayesian regularization back-propagation (BPP) algorithm because of their computational efficiency compared to DNNs. Training a DNN is usually longer than an ANN. This can be explained as follows. DNNs use multiple hidden layers and many neurons on each layer. Thus, there are many computations in each layer, and the weight optimization is more complex and longer. Inversely, ANNs use only one hidden layer, and so time for training is shorter than that for DNNs. However, we recommend using DNNs if the size of dataset is large, for example, 10,000 data points or more. For such big data, DNNs are more accurate than ANNs in predictions. Structure of this article is organized as follows. The next section shows post-buckling analysis of isotropic thin shells utilizing the first-order shear deformation shell theory (FSDT). Section 3 presents the isotropic shell formulation for post-buckling analysis utilizing FSDT and IGA. The FAP for instability problems of thin shells is presented in Section 4. The numerical verification is shown in Section 5. Finally, several notable conclusions are drawn.

2 Post-buckling analysis of isotropic shells utilizing FSDT

The aim of this article is to propose a fast analysis and prediction (FAP) approach for instability problems. Without losing the generality of the proposed approach, shell instabilities were chosen to confirm the effectiveness and exactness of the approach. For the analysis employed in Phase 1, the current formulation was reemployed from a published article on post-buckling analysis of shells. Therefore, a brief formulation of the analysis is presented in this article. The nonlinearity of the shell formulation is established upon the von Karman hypothesis. We now consider an isotropic shell depicted in Figure 1. The shell strain vectors upon the FSDT are

(1) ε = { ε x x ε y y γ x y } T = ε 0 + z κ b , γ = { γ x z γ y z } T = ε s ,

where

(2) ε 0 = ε L + ε N ; ε L = u 0 , x + w 0 R v 0 , y u 0 , y + v 0 , x ; ε N = 1 2 w 0 , x 2 w 0 , y 2 2 w 0 , x y ; κ b = β x , x β y , y β x , y + β y , x ; ε s = u 0 R + w 0 , x + β x w 0 , y + β y ,

ε N is defined as the nonlinear strain vector and can be rewritten as

(3) ε N = 1 2 A θ ; and A = w 0 , x 0 0 w 0 , y w 0 , y w 0 , x ; θ = w 0 , x w 0 , y ,

where Ω is the initial shell configuration. In the total Lagrangian methodology, the virtual work equation is expressed as

(4) Ω σ ˆ T δ ε ˆ d Ω = Ω δ u ¯ T f s d Ω ,

where f s = { f x f y f z } T is the external load vector. The vector u ¯ T = { u x u y u z } , while u z = w 0 and w 0 is the radial displacement. If an isotropic shell is only under a radial load f z = λ f 0 , the virtual work equation is re-expressed as

(5) Ω σ ˆ T δ ε ˆ d Ω = λ Ω δ w 0 f 0 d Ω .

λ is known as the load factor. σ ˆ is known as the stress resultant vector determined as

(6) σ ˆ = { σ p σ b σ s } T ,

where the in-plane stress vector is

(7) σ p = { N x N y N x y } T = h 2 h 2 ( σ x σ y τ x y ) d z T .

The bending and shear stress vectors, respectively, are

(8) σ b = { M x M y M x y } T = h 2 h 2 ( σ x σ y τ x y ) z d z T ,

(9) σ s = { Q x Q y } T = h 2 h 2 ( τ x z τ y z ) d z T ,

Figure 1 
               Panel subjected to a radial pressure.
Figure 1

Panel subjected to a radial pressure.

For isotropic shells, the generalized strain vector ε ˆ is related to the stress resultant vector σ ˆ through Hooke’s law as

(10) σ ˆ = D ˆ ε ˆ ; D ˆ = D p 0 0 0 D b 0 0 0 D s ; ε ˆ = ε L κ b ε s + ε N 0 0 ,

where

(11) D p = E h 1 υ 2 D ¯ ; D b = E h 3 12 ( 1 υ 2 ) D ¯ ; D s = κ E h 2 ( 1 + υ ) I ,

(12) D ¯ = 1 v 0 v 1 0 0 0 ( 1 v ) 2 ; I = 1 0 0 1 ,

where κ = 5 6 is the shear correction factor, h is the thickness, E is Young’s modulus, and v is Poisson’s ratio of shells.

3 Isotropic shell formulation for post-buckling analysis utilizing FSDT and IGA

This part briefly presents the post-buckling analysis of isotropic shells upon FSDT and IGA. As introduced earlier, the detailed formulation can be found in the study for the post-buckling analysis of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) shells utilizing FSDT and IGA. As known, an isotropic shell is a separate case of an FG-CNTRC shell. Therefore, the formulation for isotropic shells is the same as for FG-CNTRC shells, except there is a small difference in determining the matrices for material presented in Eqs (10)–(12). At the m th load increment and the i th iteration, the system of linear incremental equations is

(13) K T ( q m ) Δ q m i = F ext , m i F int , m i ,

where the tangent stiffness matrix is

(14) K T = Ω B A L B A b B A s + B A N 0 0 T D p 0 0 0 D b 0 0 0 D s B A L B A b B A s + B A N 0 0 d Ω + Ω ( B A g ) T N x N x y N x y N y B A g d Ω ,

and the load vector

(15) F ext , m i = ( λ m i + Δ λ m i ) Ω f 0 0 0 N A 0 0 T d Ω = ( λ m i + Δ λ m i ) F 0 ,

where F 0 is the referenced load vector, and N A is the non-uniform rational B-spline basic function. Besides, the internal force vector can be computed as

(16) F int , m i = K m i q m i ,

where

(17) K m i = Ω B A L B A b B A s + B A N 0 0 T D p 0 0 0 D b 0 0 0 D s B A L B A b B A s + 0.5 B A N 0 0 d Ω .

In this article, the modified Riks technique [3] is employed to solve the system of linear incremental equations. With an arbitrary load increment, an iterative process is performed until the following convergence criterion is satisfied as

(18) e = λ m i F 0 F int , m i ( λ m i + Δ λ m i ) F 0 < 1 0 6 .

Once the incremental solutions are achieved from Eq. (13), the load factor λ and the displacement vector q of the iteration are updated as [3]

(19) λ m i + 1 = λ m i + Δ λ m i , q m i + 1 = q m i + Δ q m i , Δ q m i = Δ q R , m i + Δ λ m i q F , m ,

where q F , m is the displacement produced by the reference force. Δ q R , m i is another displacement produced by the residual load determined as

(20) Δ q R , m i = [ K T ( q m ) ] 1 ( λ m i F 0 F int , m i ) , q F , m = [ K T ( q m ) ] 1 F 0 .

4 FAP approach for instabilities of thin shells

4.1 Brief introduction to ANNs

ANNs are systems simulated according to the operation of biological neurons in the human brain [21,31]. An ANN is a group of connected nodes as depicted in Figure 2(a). In this network, signals are spread from these neurons to those neurons via connections. Each neuron in the network has a weight adjusted during the learning process. The weight may decrease and increase the signal strength at the connection. During the learning process, the signals are transmitted from the first layer to the last layer numerous times until optimization of all the weights is achieved. The neuron’s output is computed via the function ( f ), which is sum of the neuron’s input, bias, and weight, as depicted in Figure 2(b). The formula is [21]

(21) y = f i = 1 n w i x i + b ,

where w i and x i , respectively, are the weight and input of the i th neuron connected to the considered one and b is a bias.

Figure 2 
                  Typical architecture and output calculation of an ANN: (a) architecture of a typical ANN and (b) calculating output for a neuron with three inputs.
Figure 2

Typical architecture and output calculation of an ANN: (a) architecture of a typical ANN and (b) calculating output for a neuron with three inputs.

ANNs are established upon the feed-forward BPP algorithm. In application, the input data together with the desired output data are first provided for ANNs. Then, the neuron’s weight and bias are modified (once or many times) to decrease the distinction between the actual and desired outputs of the network. The target of training ANNs is to minimize network error (the distinction between the actual output and the desired output of the network) or optimize the performance of the networks in predictions. The training process can end when the stopping criterion upon the mean square error (MSE) is met. The MSE is calculated as [21,31]

(22) MSE = 1 N i = 1 N ( e i ) 2 = 1 N i = 1 N ( t i y i ) 2 ,

where N is known as the number of samples. For the i th sample, y i is the predicted value (the output of the network) and t i is the target value. In this article, networks are trained utilizing two layers, which are the hidden layer and the output layer. The activation functions f are the PURELIN function for the output layer, while the TANSIG function for the hidden one. These activation functions are written as follows [21,31]:

(23) f = tansig ( x ) = 2 1 + e 2 x 1 f = purelin ( x ) = x .

Several optimization techniques can be employed to speed up the convergence of the BPP algorithm including the scaled conjugate gradient (SCG) technique, the B-R technique, and the Levenberg–Marquardt (L-M) technique. The details of these optimization techniques can be found in previous studies [21,31]. For all the networks in this article, the B-R technique is utilized.

4.2 B-R BPP algorithm

The objective function of ANNs utilizing B-R is written as [32]

(24) F = β E D + α E W ,

where E D is the sum of squared errors, while E W is the sum of squared weights. α and β are known as the parameters of the objective function. The network weights are the random variables when the B-R is employed. After collecting data, the density function for the weights is computed as

(25) P ( w D , α , β , M ) = P ( D w , β , M ) P ( w α , M ) P ( D α , β , M ) ,

where w is the weight vector, M is the employed model while D represents the dataset. P ( w α , M ) is the prior density, and P ( D w , β , M ) is the likelihood function. P ( D α , β , M ) is the normalization factor that ensures that the total probability is one. If we assume that both the noise in the training data and the noise in the prior distribution for the weights are Gaussian, the probability densities are written as

(26) P ( D w , β , M ) = 1 Z D ( β ) exp ( β E D ) , P ( w α , M ) = 1 Z W ( α ) exp ( α E W ) ,

where Z D ( β ) = ( π β ) N 2 and Z W ( α ) = ( π α ) N 2 . Substituting Eq. (26) into Eq. (25), we achieve

(27) P ( w D , α , β , M ) = 1 Z F ( α , β ) exp ( F ( w ) ) .

In the B-R BPP algorithm, the optimal weights are found when the posterior probability P ( w D , α , β , M ) is maximized. Maximizing the posterior probability P ( w D , α , β , M ) is equivalent to minimizing the objective function F = β E D + α E W . Now, we apply Bayes’ rule to optimize the parameters α and β of the objective function. We have

(28) P ( α , β D , M ) = P ( D α , β , M ) P ( α , β M ) P ( D M ) .

We assume that the prior density P ( α , β M ) is uniform for the parameters α and β of the objective function. Maximizing the posterior probability can be achieved by maximizing the likelihood function P ( D α , β , M ) . As seen, the likelihood function is the normalization factor in Eq. (25). The normalization factor is calculated from Eq. (25) as

(29) P ( D α , β , M ) = P ( D w , β , M ) P ( w α , M ) P ( w D , α , β , M ) .

Note that all the probabilities have a Gaussian form. Thus, Eq. (29) can be re-expressed as

(30) P ( D α , β , M ) = Z F ( α , β ) Z D ( β ) Z W ( α ) ,

where Z D ( β ) and Z W ( α ) are the constants and defined before. Z F ( α , β ) can be computed utilizing a Taylor series expansion as

(31) Z F ( 2 π ) N 2 ( det ( ( H ) 1 ) ) 1 2 exp ( F ( w ) ) .

H = β 2 E D + α 2 E W is known as the Hessian matrix. Placing the aforementioned result in Eq. (30), we obtain the optimal values for α and β as

(32) α = γ 2 E W ( w ) and β = N γ 2 E D ( w ) ,

where γ = N 2 α tr ( H ) 1 is the effective number of parameters, while N is the total number of parameters in the network. Finally, optimization of the parameters α and β is performed through some of the following steps [32]:

  • (0) Initialize the weights w , the parameters α and β .

  • (1) Employ the L-M algorithm for finding the minimum of the objective function.

  • (2) Calculate γ = N 2 α tr ( H ) 1 . Herein, utilize the Gauss-Newton approximation for the Hessian matrix as follows: H = 2 F ( w ) 2 β J T J + 2 α I N .

  • (3) Compute new values of α and β utilizing Eq. (32).

  • (4) Steps 1–3 are iterated until we obtain the convergence.

In Step 1, the weight is updated after each iteration of the L-M algorithm as

(33) w i + 1 = w i ( J i T J i + μ i I ) 1 J i T e i ,

where I is the unit matrix, e is the error factor, and μ represents the learning parameter. The global minimum of the objective function is reached when its value does not dramatically change in subsequent iterations.

4.3 FAP approach for instability problems of thin shells

Some popular types of thin shell instabilities chosen to confirm the effectiveness and exactness of the FAP were snap-through, snap-back, kink, and softening–hardening instabilities. One of the great difficulties of solving structural instability problems is highly time-consuming because of numerous computational iterations and a huge computation. This reality is clearly seen in large-scale problems utilizing thousands or millions of degrees of freedom. Another difficulty in solving instability problems is divergence in many cases especially when the equilibrium path passes the limit points, the inflection points, the neighborhoods of turning points, and bifurcations. To overcome the aforementioned difficulties, the FAP for instability problems of thin shells is proposed in this article. It is noted that in instability or post-buckling analysis of structures, loads and displacement vectors of analysis steps are unknown. The present methodology contains two phases: fast numerical analysis and pure prediction utilizing ANNs as follows: (1) for Phase 1, post-buckling analysis is conducted utilizing a minor amount of load steps. The obtained load–displacement relation is not exact because of a minor amount of load steps utilized; (2) for Phase 2, the loads and the deflections obtained from Phase 1 are considered as the data for training networks. Then, the trained networks, including the load network and the displacement network, are utilized to fast predict loads and deflections at any step of the post-buckling analysis. Upon utilizing Phase 2, an exact and smooth load–displacement curve is achieved. Figures 3 and 4, respectively present the description and the flowchart of the FAP. The numerical verification section demonstrates the high exactness and effectiveness of the FAP for instability problems. Besides, the high ability of the FAP to detect limit and inflection points is also confirmed in the next section.

Figure 3 
                  FAP for instability problems utilizing ANNs.
Figure 3

FAP for instability problems utilizing ANNs.

Figure 4 
                  Workflow of the present methodology for instability problems.
Figure 4

Workflow of the present methodology for instability problems.

5 Numerical verification

The effectiveness and exactness of the present approach for instability problems of thin shells are verified in this section. Without losing the generality of the proposed approach, four popular types of thin shell instabilities chosen to confirm the effectiveness and exactness of the present approach were snap-through, snap-back, kink and softening–hardening instabilities. For convenience and consistency, the following parameters are fixed for the whole manuscript.

  • For the analysis in Phase 1, the computational models are analyzed utilizing 14 × 14 cubic elements. The hinged boundary is described as u 0 = v 0 = w 0 = β y = 0 , at x = 0 , α R . The clamped boundary is imposed as u 0 = v 0 = w 0 = β x = β y = 0 . Besides, the geometric data of shells are provided as the thickness h , and the length L = B = 508 mm. The radius is R = 2,540 mm. Note that the thickness h differs in the first three problems in this article. The material properties of all shells are given as Young’s modulus E = 3.103 kN/mm 2 and Poisson’s ratio v = 0.3 .

  • For the prediction packages in Phase 2, the structure of all ANNs includes the hidden and output layers. The number of neurons in the first and second layers is 15 and 1, respectively. ANNs utilizing the B-R BPP algorithm did not require a validation dataset to evaluate the models. The algorithm has a built-in validation technique. Therefore, in the training phase, the dataset was separated into two data groups: the training data ( 85 % ) and the testing data ( 15 % ). Utilizing the aforementioned parameters, we proved numerically that all the following networks are well trained and performed in all the problems. Notably, the cost for both training and application of an ANN is very low because of the very small datasets are employed together with the advantages of ANNs. Thus, it can be neglected in the effectiveness computation of the FAP. This fact is mentioned in the first problem.

  • All the networks were trained utilizing the NN-fitting application built in MATLAB R2021a. The ANN’s parameters and architecture are listed in table 1 in the study by Nguyen et al. [30].

For verification study, the results of the FAP are compared with that of pure IGA and reference solutions. The pure IGA is conducted utilizing the same formulation as the FAP but needs much more load steps.

5.1 FAP approach for snap-through instability problems of thin shells

The effectiveness and exactness of the FAP for snap-through instability problems of shells and the performance of ANNs incorporated with the B-R are studied in this subsection. We consider a hinged-hinged isotropic cylindrical panel under a central point load, as in Figure 5. The geometric data and the material properties of the panel are provided at the beginning of this section. Note that in this problem, the panel thickness is h = 12.7 mm. In Phase 1, we perform a fast post-buckling analysis of the thin isotropic panel utilizing 13 equal load steps. The obtained load–displacement relation is shown in Figure 6. This relation is not exact because of a minor amount of load steps utilized. To achieve a smooth, complete, and correct load–displacement curve, we must apply Phase 2 of the FAP utilizing ANNs. We need to build up two ANNs in Phase 2. They are the networks for displacement prediction and load prediction. The central deflections and loads obtained from Phase 1 are considered as target data for training ANNs in Phase 2. Data structures of Phases 1 and 2 for displacement and load predictions are, respectively, presented in Figures 7 and 8. In Phase 2, the training process and the application of the displacement network are the same as that of the load network. These processes are explained in detail below.

  • For the training process of the displacement network, the data sizes of the input, output, and target are equal. The number of load steps in Phase 1 is 13. Thirteen achieved deflections are utilized as the target data for training the displacement network. For training networks, we utilize a typical input data as [ 1 , 2 , , 14 ] in order that the sizes of the input data and target data are equal as in Figure 7 (the left side). Notably, this typical input data is employed for simplicity and convenience without losing of the generality and exactness of the FAP. After that, the displacement network is trained utilizing the network parameters provided at the beginning of Section 5.

  • For the application of the displacement network: if we need to predict full deflections of 31 equal load steps utilizing the trained ANN, the input is still from 1 to 14 but the increment in the data is ( 14 1 ) n = 0.419 . In which, n = 31 corresponds to 31 load steps in order to the sizes of the input data, the output data, and the target data are equal as in Figure 7 (the right side). The obtained output is compared to the target achieved from pure IGA utilizing 31 equal load steps. A satisfactory agreement is reached.

  • The training process and the application of the load network are the same as that of the displacement network. The details can be found in Figure 8.

Figure 5 
                  Hinged-hinged panel under a point load at the center.
Figure 5

Hinged-hinged panel under a point load at the center.

Figure 6 
                  Analysis of the hinged-hinged panel utilizing the FAP: Phase 1 (fast analysis, 13 load steps), 
                        
                           
                           
                              h
                           
                           h
                        
                      = 12.7 mm.
Figure 6

Analysis of the hinged-hinged panel utilizing the FAP: Phase 1 (fast analysis, 13 load steps), h = 12.7 mm.

Figure 7 
                  Data structure of the first and second phases for displacement prediction.
Figure 7

Data structure of the first and second phases for displacement prediction.

Figure 8 
                  Data structure of the first and second phases for load prediction.
Figure 8

Data structure of the first and second phases for load prediction.

Figure 9 presents the training performance of the displacement and load networks. The training processes of the two networks can be, respectively, stopped after 443 and 99 epochs because the MSEs are very small. The results of applications of the trained ANNs are presented in Figure 10. The very good agreements are found between the obtained outputs (from the networks) and the targets (from pure IGA). The networks are well trained. Interestingly, the cost for training and application of a network is very low (approximately 3 seconds) because the sizes of data are very small together with the advantages of ANNs. Thus, this cost is neglected in evaluation of the effectiveness of the FAP. Besides, Figure 11 presents the obtained results of the FAP utilizing 13 fast analysis steps in the comparison with that of pure IGA utilizing 31 load steps and that of Sabir and Djoudi [33] utilizing pure finite element (FE) analysis with 31 load steps. These results match well. A smooth, complete, and exact curve is achieved after utilizing Phase 2. Also, the upper and lower limit points are detected successfully upon ANNs. The present approach saves 58% of the computational effort compared to other approaches. The high effectiveness and exactness of the present approach are verified. Figure 12 describes the effect of the number of fast analysis steps on the solution. The higher number of fast analysis steps, the higher exactness of the FAP solution. Utilizing 13 analysis steps in Phase 1 can produce a great solution that balances exactness and effectiveness. For snap-through instability problems, we recommend to utilize the present approach with 13 fast analysis steps to ensure exactness and effectiveness in solutions.

Figure 9 
                  Training performance of ANNs in Phase 2: (a) for displacement network and (b) for load network.
Figure 9

Training performance of ANNs in Phase 2: (a) for displacement network and (b) for load network.

Figure 10 
                  Applying the trained networks in Phase 2: (a) for displacement (mm) network and (b) for load (kN) network.
Figure 10

Applying the trained networks in Phase 2: (a) for displacement (mm) network and (b) for load (kN) network.

Figure 11 
                  Central deflection of the hinged-hinged panel utilizing various approaches, 
                        
                           
                           
                              h
                           
                           h
                        
                      = 12.7 mm.
Figure 11

Central deflection of the hinged-hinged panel utilizing various approaches, h = 12.7 mm.

Figure 12 
                  Number of fast analysis steps affects the FAP solution for a hinged-hinged panel, 
                        
                           
                           
                              h
                           
                           h
                        
                      = 12.7 mm.
Figure 12

Number of fast analysis steps affects the FAP solution for a hinged-hinged panel, h = 12.7 mm.

5.2 FAP approach for softening–hardening instability problems of shells

The hinged-hinged panel under a central point load depicted in Figure 5 is continued considering. The geometric data and the material data of the panel are unchanged, except the panel thickness in this problem is h = 25.4 mm. Phase 1 is conducted utilizing seven load steps. The obtained load–displacement relation is shown in Figure 13. As seen, this relation is incomplete and not smooth because of a minor amount of load steps utilized. To achieve a smooth, complete, and correct load–displacement response, we must apply Phase 2 of the FAP utilizing ANNs. The result of Phase 2 is presented in Figure 14 in the comparison with that of pure IGA with 14 load steps and that of Sabir and Djoudi [33] utilizing FE analysis with 13 load steps. These results match well. As expected, a smooth, complete, and exact response was achieved after utilizing Phase 2. Besides, the inflection point is detected successfully upon ANNs. The proposed approach can save about 50 % of computational effort compared to other ones. The high effectiveness and exactness of the FAP was confirmed. Figure 15 shows the effect of the number of fast analysis steps on the FAP solution. Again, higher number of fast analysis steps, the higher exactness of the FAP solution. Utilizing seven fast analysis steps produces the great solution that balances exactness and effectiveness. For softening–hardening instability problems, we recommend to utilize the FAP with seven fast analysis steps to ensure exactness and effectiveness in solutions.

Figure 13 
                  Analysis of the hinged-hinged panel utilizing the FAP: Phase 1 (fast analysis, seven load steps), 
                        
                           
                           
                              h
                           
                           h
                        
                      = 25.4 mm.
Figure 13

Analysis of the hinged-hinged panel utilizing the FAP: Phase 1 (fast analysis, seven load steps), h = 25.4 mm.

Figure 14 
                  Central deflection of the hinged-hinged panel utilizing various approaches, 
                        
                           
                           
                              h
                           
                           h
                        
                      = 25.4 mm.
Figure 14

Central deflection of the hinged-hinged panel utilizing various approaches, h = 25.4 mm.

Figure 15 
                  Number of fast analysis steps affects the FAP solution for a hinged-hinged panel, 
                        
                           
                           
                              h
                           
                           h
                        
                      = 25.4 mm.
Figure 15

Number of fast analysis steps affects the FAP solution for a hinged-hinged panel, h = 25.4 mm.

5.3 FAP approach for snap-back instability problems of shells

In the this problem, the clamped-hinged panel under a central point load depicted in Figure 16 is considered. The geometric data and material data of the panel are the same as those in the first problem, except the panel thickness, in this case, is h = 6.35 mm. Phase 1 is performed utilizing 14 load steps. The load–displacement relation achieved from Phase 1 is shown in Figure 17. Similar to the previous observations in the first and second problems, the load–displacement relation is a broken line and incomplete, especially in the areas around the limit points. However after application of Phase 2, the load–displacement response becomes smooth, complete, and exact. This fact can be seen in Figure 18 in a comparison of the response of fully pure IGA with 43 load steps. It is concluded that the limit points and the snap-back instability of the panel are detected successfully upon ANNs. Note that to trace successfully the snap-back response of the fully pure IGA, we need to employ a smaller arc length compared to that in the two previous problems. As a result, the number of load steps (43) in the pure analysis of this problem is higher than that of the two previous problems. Figure 19 shows the effect of the number of fast analysis steps on the FAP solution. Utilizing 14 fast analysis steps can produce the good solution with the balance between exactness and effectiveness. For snap-back instability problems, we recommend to utilize the FAP with 14 fast analysis steps to ensure exactness and effectiveness in solutions.

Figure 16 
                  Clamped-hinged panel under a central load, 
                        
                           
                           
                              h
                           
                           h
                        
                      = 6.35 mm.
Figure 16

Clamped-hinged panel under a central load, h = 6.35 mm.

Figure 17 
                  Analysis of the clamped-hinged panel utilizing the FAP: Phase 1 (fast analysis, 14 load steps), 
                        
                           
                           
                              h
                           
                           h
                        
                      = 6.35 mm.
Figure 17

Analysis of the clamped-hinged panel utilizing the FAP: Phase 1 (fast analysis, 14 load steps), h = 6.35 mm.

Figure 18 
                  Central deflection of the clamped-hinged panel utilizing various approaches, 
                        
                           
                           
                              h
                           
                           h
                        
                      = 6.35 mm.
Figure 18

Central deflection of the clamped-hinged panel utilizing various approaches, h = 6.35 mm.

Figure 19 
                  Number of fast analysis steps affects the FAP solution for a clamped-hinged panel, 
                        
                           
                           
                              h
                           
                           h
                        
                      = 6.35 mm.
Figure 19

Number of fast analysis steps affects the FAP solution for a clamped-hinged panel, h = 6.35 mm.

ANNs incorporated with the B-R algorithm recommended for the FAP have some notable advantages compared to numerous existing networks. To confirm this, we solved this problem utilizing different types of NNs in Phase 2 including ANNs incorporated with the L-M algorithm [21,31], ANNs incorporated with the SCG algorithm [21,31], DNNs [34], and the group method of data handling (GMDH) [35]. For this problem, the data obtained in Phase 1 were re-utilized for training networks in Phase 2. Notably, the ANN’s structure, which includes number of neurons and layers, is unchanged. The change is the SCG algorithm, and the L-M algorithm is utilized in this investigation. All the ANNs were trained utilizing MATLAB R2021a as mentioned earlier. For all types of networks in this investigation and in this article, data splitting is performed as follows: the training set (70%), the testing set (15%), and the validation set (15%). Group method of data handling is utilized in many fields including forecasting, optimization, and pattern recognition. In this article, a GMDH-based NN [35] includes 15 neurons per layer and 4 hidden layers. We choose the selection pressure factor α = 0.6 . The DNN [34,36,37] is a network utilizing many hidden layers. It can be employed to solve complex nonlinear problems and identify complicated nonlinear systems. DNNs have been investigated for numerous fields including speech recognition, translation, and computer vision. Similar to the GMDH-based network, a DNN includes 15 neurons per layer and 4 hidden layers. The rectified-linear unit function is utilized as an activation function. Besides, the adaptive moment estimation is utilized as an optimizer. MSE is chosen as a loss function. We fix the learning rate at 0.001 and the number of epochs at 1,000. The aforementioned optimizer and functions are well known and presented in detail in the study by Kim et al. [34]. Figures 20, 21, 22, 23, 24 present the applications of the mentioned networks. In the aforementioned figures, the ground truth or the targets are the values obtained from pure IGA utilizing many load steps. It is found that utilizing ANNs incorporated with the L-M algorithm, ANNs incorporated with the SCG algorithm, GMDH, and DNNs produces incorrect predictions. It should be noted that for each type of network, we need to build up and train two networks to predict displacements and loads. Figure 25 presents the results of the present approach utilizing different types of NNs. These results are compared to that of pure IGA utilizing many load steps. The solution of the present approach utilizing ANNs incorporated with B-R agrees very well with that of the pure IGA, while the remaining solutions are completely incorrect and unstable. It should be noted that DNNs and GMDH are the improved networks (compared to ANNs) and structured with 15 neurons per layer and 4 hidden layers, while ANNs contain 15 neurons and 1 hidden layer. However, utilizing ANNs incorporated with B-R produces the best solution, as found in Figure 25. This confirms the notable advantages of ANNs incorporated with the B-R algorithm for training with highly nonlinear, difficult, and small data. We recommend to utilize ANNs incorporated with B-R in Phase 2 to achieve exact predictions. Figure 26 shows the effect of the network’s architecture on the FAP solution. It can be concluded that the more neurons, the more exact solution; the more data, the more exact solution. To obtain exact predictions, we recommend to utilize the present approach with the network’s architecture as: at least 85% data for training and 15 neurons in the hidden layer.

Figure 20 
                  Applications of the networks utilizing ANNs and the L-M algorithm: (a) network for displacement (mm) prediction and (b) network for load (kN) prediction.
Figure 20

Applications of the networks utilizing ANNs and the L-M algorithm: (a) network for displacement (mm) prediction and (b) network for load (kN) prediction.

Figure 21 
                  Applications of the networks utilizing ANNs and the SCG algorithm: (a) network for displacement (mm) prediction and (b) network for load (kN) prediction.
Figure 21

Applications of the networks utilizing ANNs and the SCG algorithm: (a) network for displacement (mm) prediction and (b) network for load (kN) prediction.

Figure 22 
                  Applications of the networks utilizing the GMDH: (a) network for displacement (mm) prediction and (b) network for load (kN) prediction.
Figure 22

Applications of the networks utilizing the GMDH: (a) network for displacement (mm) prediction and (b) network for load (kN) prediction.

Figure 23 
                  Training performances of DNNs: (a) network for displacement prediction and (b) network for load prediction.
Figure 23

Training performances of DNNs: (a) network for displacement prediction and (b) network for load prediction.

Figure 24 
                  Applications of the trained networks utilizing the DNNs: (a) network for displacement (mm) prediction and (b) network for load (kN) prediction.
Figure 24

Applications of the trained networks utilizing the DNNs: (a) network for displacement (mm) prediction and (b) network for load (kN) prediction.

Figure 25 
                  Analysis of a clamped-hinged panel utilizing the FAP and different NNs.
Figure 25

Analysis of a clamped-hinged panel utilizing the FAP and different NNs.

Figure 26 
                  Network’s architecture affects the FAP solution: (a) number of neurons, data-splitting ratio: 85%, 15% and (b) data-splitting ratio (training, testing), 15 neurons in the hidden layer.
Figure 26

Network’s architecture affects the FAP solution: (a) number of neurons, data-splitting ratio: 85%, 15% and (b) data-splitting ratio (training, testing), 15 neurons in the hidden layer.

5.4 FAP approach for kink instability problems of shells

The panel in Figure 5 and Section 5.1 is continued considering. In this problem, the panel thickness is h = 6.35 mm and the laminated composite material is employed with [ 0 0 9 0 0 0 0 ] stacking sequence. The material properties are provided as [38]: E 1 = 3.3 kN/mm 2 , E 2 = 1.1 kN/mm 2 , G 12 = G 13 = 0.66 kN/mm 2 , G 12 = G 23 , ν 12 = 0.25 . A total of 18 data points on the equilibrium path, which has a kink, are extracted from the study by Sze et al. [38] as the result of Phase 1. These data points are employed to train networks for load and displacement in Phase 2. The trained networks are then employed to predict loads and displacements for 36 steps. The results of Phases 1 and 2 are presented in Figure 27 in a comparison with that of Sze et al. [38] utilizing FEM and 36 steps. A very good agreement is found. It is seen that the result of Phase 1 is not smooth and not completely exact. This is clearly seen when displacement is from 0 to 20. However, after application of the trained networks, the equilibrium path is smooth and exact. Besides, the FAP saves 50% computational effort compared to FEM.

Figure 27 
                  Analysis of the hinged-hinged composite panel with 
                        
                           
                           
                              h
                           
                           h
                        
                      = 6.35 mm and [
                        
                           
                           
                              
                                 
                                    0
                                 
                                 
                                    0
                                 
                              
                              ∕
                              9
                              
                                 
                                    0
                                 
                                 
                                    0
                                 
                              
                              ∕
                              
                                 
                                    0
                                 
                                 
                                    0
                                 
                              
                           
                           {0}^{0}/9{0}^{0}/{0}^{0}
                        
                     ] stacking sequence utilizing the FAP: (a) Phase 1 (fast analysis) and (b) Phase 2 (prediction) and comparison.
Figure 27

Analysis of the hinged-hinged composite panel with h = 6.35 mm and [ 0 0 9 0 0 0 0 ] stacking sequence utilizing the FAP: (a) Phase 1 (fast analysis) and (b) Phase 2 (prediction) and comparison.

5.5 FAP approach for an experiment

In the last problem, the FAP is applied to predict delamination buckling and post-buckling of composite laminates utilizing data from experiment. The compression test was conducted on HYE-3574 OH graphite/epoxy composites with built-in delamination [39]. The material properties are found in table 1 in the study by Gu and Chattopadhyay [39]. Geometry of specimen and compression test setup are shown in Figure 28 and Figure 4 in the study by Gu and Chattopadhyay [39], respectively. The delamination length is 3.0 in. The ply stacking sequence of the test specimen was [ 4 5 0 4 5 0 ] 9 S . The delamination is near the surface and at the distance 1.7 mm from the surface. The test specimen was clamped at both ends. The axial displacement and compressive load data were digitally recorded. Ten data points on the compressive stress–axial displacement curve are extracted from figure 5d in the study by Gu and Chattopadhyay [39] as the result of Phase 1. These data points are employed to train networks for stress and displacement in Phase 2. The trained networks are then employed to predict stresses and displacements for any load steps. The results of Phases 1 and 2 are presented in Figure 29 in a comparison with the experimental result [39]. A very good agreement is found. As known, the experimental data are affected by many factors including specimen fabrication technique, test setup, and technique of recording data. Thus, training a network utilizing experimental data is more difficult than that utilizing computational data. Interestingly, the present approach is exact for both the computational and experimental data. Although the aforementioned networks were trained utilizing ten experimental data points, they can be employed to predict stresses and displacements for any load steps. This can save labor, time, and money for the experiment.

Figure 28 
                  An experiment of delamination postbuckling of composite laminates: Geometry of specimen [39].
Figure 28

An experiment of delamination postbuckling of composite laminates: Geometry of specimen [39].

Figure 29 
                  Postbuckling path for the test specimen utilizing the FAP: (a) Phase 1 and (b) Phase 2 (prediction) and comparison.
Figure 29

Postbuckling path for the test specimen utilizing the FAP: (a) Phase 1 and (b) Phase 2 (prediction) and comparison.

Finally, the computational effectiveness of the FAP for various types of instabilities of thin shells is summarized and presented in Tables 1 and 2. It is concluded that utilizing the FAP can save a huge amount of computation compared to the other ones (fully pure isogeometric and fully pure FE analyses) and the experiment. Notably, the computational costs for problems of kink instability and experiment were not available in the literature [38,39]. Therefore, these costs are not present in Table 2.

Table 1

Effectiveness of the FAP according to load steps

Instabilities Total steps per analysis Saved computation (%) compared
Pure analysis FAP approach with pure analysis approach
Snap-through 31 13 58
Softening–hardening 14 7 50
Snap–back 43 14 67
Kink 36 18 50
Experiment 19 10 47
Table 2

Effectiveness of the FAP according to computational cost

Instabilities Total cost per analysis (s) Saved cost (%) compared
Pure analysis FAP approach with pure analysis approach
Snap-through 261 126 52
Softening-hardening 116 64 45
Snap-back 403 132 67

6 Conclusions

An FAP approach for instability problems of thin shells has been proposed. The present approach contains two phases: fast analysis and pure prediction. (1) for Phase 1, fast post-buckling analysis is conducted utilizing a minor amount of load steps. The obtained load–displacement relation from Phase 1 is incomplete and not exact because of a minor amount of load steps utilized; (2) for Phase 2, the loads and the deflections obtained from Phase 1 are considered as the data for training ANNs incorporated with Bayesian regularization (B-R). After that, the trained networks, including the load network and the displacement network, are utilized to fast predict load and deflection at any step of the post-buckling analysis. After utilizing Phase 2, an exact, complete, and smooth load–displacement curve is achieved. Interestingly, the cost for training and application of an ANN is very low (approximately 3 s) because the utilized data sizes are very small together with the advantages of ANNs. From the numerical verification, several notable conclusions are drawn as follows:

  • Good agreement between the results of the FAP and the reference results was found for all problems in this work. The present approach saves a huge computation compared to other ones. High effectiveness and exactness of the proposed approach for post-buckling analysis of thin shells were demonstrated.

  • The FAP solutions for instability problems are stable and exact. The more number of fast analysis steps, the more exact solution but high computational cost is required.

  • To produce good solutions for the FAP with a balance between exactness and effectiveness (the lowest computational cost), some conclusions are drawn as follows: For snap-through instability problems, we recommend to utilize the present approach with 13 load steps in Phase 1. For softening–hardening instability problems, we recommend to utilize the present approach with seven load steps in Phase 1. For snap-back instability problems, we recommend to utilize the present approach with 14 load steps in Phase 1.

  • ANNs incorporated with B-R, which are recommended for the present approach, have some notable advantages compared to numerous existing networks. We recommend to utilize ANNs incorporated with B-R in Phase 2 to achieve exact predictions. Besides, we recommend to utilize the present approach with the network’s architecture as: at least 85% data for training and 15 neurons in the hidden layer.

  • The FAP is applicable to experiments and highly nonlinear problems, which can have a kink on their equilibrium paths. This can save labor, time, and money for the experiments.

  1. Funding information: This research project was supported by the Second Century Fund (C2F), Chulalongkorn University.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: The datasets generated during and analyzed during the current study are available from the corresponding author on reasonable request.

References

[1] Zhang LW, Lei ZX, Liew KM, Yu JL. Large deflection geometrically nonlinear analysis of carbon nanotube-reinforced functionally graded cylindrical panels. Comput Meth Appl Mech Eng. 2014;273:1–18. 10.1016/j.cma.2014.01.024Search in Google Scholar

[2] Crisfield MA. A faster modified Newton-Raphson iteration. Comput Meth Appl Mech Eng. 1979;20(3):267–78. 10.1016/0045-7825(79)90002-1Search in Google Scholar

[3] Crisfield MA. A fast incremental/iterative solution procedure that handles ‘snap-through. Comput Struct. 1981;13(1):55–62. 10.1016/0045-7949(81)90108-5Search in Google Scholar

[4] Rezaiee-Pajand M, Naserian R. Geometrical nonlinear analysis based on optimization technique. Appl Math Model. 2018;53:32–48. 10.1016/j.apm.2017.08.003Search in Google Scholar

[5] Rezaiee-Pajand M, Naserian R. Using residual areas for geometrically nonlinear structural analysis. Ocean Eng. 2015;105:327–35. 10.1016/j.oceaneng.2015.06.043Search in Google Scholar

[6] Castro SGP, Jansen EL. Displacement-based formulation of Koiter’s method: Application to multi-modal post-buckling finite element analysis of plates. Thin-Walled Struct. 2021;159:107217. 10.1016/j.tws.2020.107217Search in Google Scholar

[7] Sun Y, Tian K, Li R, Wang B. Accelerated Koiter method for post-buckling analysis of thin-walled shells under axial compression. Thin-Walled Struct. 2020;155:106962. 10.1016/j.tws.2020.106962Search in Google Scholar

[8] Liang K, Ruess M, Abdalla M. The Koiter-Newton approach using von xxxx kinematics for buckling analyses of imperfection sensitive structures. Comput Meth Appl Mech Eng. 2014;279:440–68. 10.1016/j.cma.2014.07.008Search in Google Scholar

[9] Magisano D, Liang K, Garcea G, Leonetti L, Ruess M. An efficient mixed variational reduced-order model formulation for nonlinear analyses of elastic shells. Int J Numer Methods Eng. 2018;113(4):634–55. 10.1002/nme.5629Search in Google Scholar

[10] Mei Y, Hurtado DE, Pant S, Aggarwal A. On improving the numerical convergence of highly nonlinear elasticity problems. Comput Meth Appl Mech Eng. 2018;337:110–27. 10.1016/j.cma.2018.03.033Search in Google Scholar

[11] Groh RMJ, Avitabile D, Pirrera A. Generalised path-following for well-behaved nonlinear structures. Comput Meth Appl Mech Eng. 2018;331:394–426. 10.1016/j.cma.2017.12.001Search in Google Scholar

[12] Magisano D, Leonetti L, Garcea G. How to improve efficiency and robustness of the Newton method in geometrically non-linear structural problem discretized via displacement-based finite elements. Comput Meth Appl Mech Eng. 2017;313:986–1005. 10.1016/j.cma.2016.10.023Search in Google Scholar

[13] Rezaiee-Pajand M, Estiri H. Geometrically nonlinear analysis of shells by various dynamic relaxation methods. World J Eng. 2017;14(5):381–405. 10.1108/WJE-10-2016-0109Search in Google Scholar

[14] Rezaiee-Pajand M, Estiri H. Finding equilibrium paths by minimizing external work in dynamic relaxation method. Appl Math Model. 2016;40(23):10300–22. 10.1016/j.apm.2016.07.017Search in Google Scholar

[15] Rezaiee-Pajand M, Afsharimoghadam H. An incremental iterative solution procedure without predictor step. Eur J Comput Mech. 2018;27(1):58–87. 10.1080/17797179.2018.1455028Search in Google Scholar

[16] Maghami A, Shahabian F, Hosseini SM. Path following techniques for geometrically nonlinear structures based on Multi-point methods. Comput Struct. 2018;208:130–42. 10.1016/j.compstruc.2018.07.005Search in Google Scholar

[17] Chojaczyk AA, Teixeira AP, Neves LC, Cardoso JB, Guedes Soares C. Review and application of Artificial Neural Networks models in reliability analysis of steel structures. Struct Saf. 2015;52:78–89. 10.1016/j.strusafe.2014.09.002Search in Google Scholar

[18] Fahmy AS, El-Madawy MET, Atef Gobran Y. Using artificial neural networks in the design of orthotropic bridge decks. Alex Eng J. 2016;55(4):3195–203. 10.1016/j.aej.2016.06.034Search in Google Scholar

[19] Sebaaly H, Varma S, Maina JW. Optimizing asphalt mix design process using artificial neural network and genetic algorithm. Constr Build Mater. 2018;168:660–70. 10.1016/j.conbuildmat.2018.02.118Search in Google Scholar

[20] Wang S, Sun G, Chen W, Zhong Y. Database self-expansion based on artificial neural network: An approach in aircraft design. Aerosp Sci Technol. 2018;72:77–83. 10.1016/j.ast.2017.10.037Search in Google Scholar

[21] Nikbin IM, Rahimi I, Allahyari RS. A new empirical formula for prediction of fracture energy of concrete based on the artificial neural network. Eng Fract Mech. 2017;186:466–82. 10.1016/j.engfracmech.2017.11.010Search in Google Scholar

[22] Tohidi S, Sharifi Y. Neural networks for inelastic distortional buckling capacity assessment of steel I-beams. Thin-Walled Struct. 2015;94:359–71. 10.1016/j.tws.2015.04.023Search in Google Scholar

[23] Eskandari-Naddaf H, Kazemi R. ANN prediction of cement mortar compressive strength, influence of cement strength class. Constr Build Mater. 2017;138:1–11. 10.1016/j.conbuildmat.2017.01.132Search in Google Scholar

[24] Ye F, Ma S, Tong L, Xiao J, Baaaanard P, Chahine R. Artificial neural network based optimization for hydrogen purification performance of pressure swing adsorption. The 6th International Conference on Energy, Engineering and Environmental Engineering. 2019. Vol. 44. No. 11. p. 5334–44. 10.1016/j.ijhydene.2018.08.104Search in Google Scholar

[25] Kong YS, Abdullah S, Schramm D, Omar MZ, Haris SM. Optimization of spring fatigue life prediction model for vehicle ride using hybrid multi-layer perceptron artificial neural networks. Mech Syst Signal Proc. 2019;122:597–621. 10.1016/j.ymssp.2018.12.046Search in Google Scholar

[26] Pathirage CSN, Li J, Li L, Hao H, Liu W, Ni P. Structural damage identification based on autoencoder neural networks and deep learning. Eng Struct. 2018;172:13–28. 10.1016/j.engstruct.2018.05.109Search in Google Scholar

[27] Tan ZX, Thambiratnam DP, Chan THT, Abdul Razak H. Detecting damage in steel beams using modal strain energy based damage index and Artificial Neural Network. Eng Fail Anal. 2017;79:253–62. 10.1016/j.engfailanal.2017.04.035Search in Google Scholar

[28] Suk JW, Kim JH, Kim YH. A predictor algorithm for fast geometrically-nonlinear dynamic analysis. Comput Meth Appl Mech Eng. 2003;192(22):2521–38. 10.1016/S0045-7825(03)00270-6Search in Google Scholar

[29] Kim JH, Kim YH. A predictor–corrector method for structural nonlinear analysis. Comput Meth Appl Mech Eng. 2001;191(8):959–74. 10.1016/S0045-7825(01)00296-1Search in Google Scholar

[30] Nguyen TN, Zhang D, Mirrashid M, Nguyen DK, Singhatanadgid P. Fast analysis and prediction approach for geometrically nonlinear bending analysis of plates and shells using artificial neural networks. Mechanics of Advanced Materials and Structures. Philadelphia: Taylor & Francis. 2023. p. 1–19. 10.1080/15376494.2023.2286626. Search in Google Scholar

[31] Mouloodi S, Rahmanpanah H, Gohari S, Burvill C, Davies HMS. Feedforward backpropagation artificial neural networks for predicting mechanical responses in complex nonlinear structures: A study on a long bone. J Mech Behav Biomed Mater. 2022;128:105079. 10.1016/j.jmbbm.2022.105079Search in Google Scholar PubMed

[32] Dan Foresee F, Hagan MT. Gauss-Newton approximation to Bayesian learning. In: Proceedings of International Conference on Neural Networks (ICNN’97). vol. 3; 1997. p. 1930–5. 10.1109/ICNN.1997.614194Search in Google Scholar

[33] Sabir AB, Djoudi MS. Shallow shell finite element for the large deflection geometrically nonlinear analysis of shells and plates. Thin-Walled Struct. 1995;21(3):253–67. 10.1016/0263-8231(94)00005-KSearch in Google Scholar

[34] Kim DJ, Kim GW, Baek Jh, Nam B, Kim HS. Prediction of stress-strain behavior of carbon fabric woven composites by deep neural network. Compos Struct. 2023;318:117073. 10.1016/j.compstruct.2023.117073Search in Google Scholar

[35] Ivakhnenko AG. Polynomial theory of complex systems. IEEE Trans Syst Man Cybern. 1971;SMC-1(4):364–78. Number: 4. 10.1109/TSMC.1971.4308320Search in Google Scholar

[36] Samaniego E, Anitescu C, Goswami S, Nguyen-Thanh VM, Guo H, Hamdia K, et al. An energy approach to the solution of partial differential equations in computational mechanics via machine learning: Concepts, implementation and applications. Comput Meth Appl Mech Eng. 2020;362:112790. 10.1016/j.cma.2019.112790Search in Google Scholar

[37] Zhuang X, Guo H, Alajlan N, Zhu H, Rabczuk T. Deep autoencoder based energy method for the bending, vibration, and buckling analysis of Kirchhoff plates with transfer learning. Eur J Mech A-Solids. 2021;87:104225. 10.1016/j.euromechsol.2021.104225Search in Google Scholar

[38] Sze KY, Liu XH, Lo SH. Popular benchmark problems for geometric nonlinear analysis of shells. Finite Elem Anal Des. 2004;40(11):1551–69. 10.1016/j.finel.2003.11.001Search in Google Scholar

[39] Gu H, Chattopadhyay A. An experimental investigation of delamination buckling and postbuckling of composite laminates. Compos Sci Technol. 1999;59(6):903–10. 10.1016/S0266-3538(98)00130-4Search in Google Scholar

Received: 2024-01-22
Revised: 2024-04-17
Accepted: 2024-04-22
Published Online: 2024-06-24

© 2024 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Editorial
  2. Focus on NLENG 2023 Volume 12 Issue 1
  3. Research Articles
  4. Seismic vulnerability signal analysis of low tower cable-stayed bridges method based on convolutional attention network
  5. Robust passivity-based nonlinear controller design for bilateral teleoperation system under variable time delay and variable load disturbance
  6. A physically consistent AI-based SPH emulator for computational fluid dynamics
  7. Asymmetrical novel hyperchaotic system with two exponential functions and an application to image encryption
  8. A novel framework for effective structural vulnerability assessment of tubular structures using machine learning algorithms (GA and ANN) for hybrid simulations
  9. Flow and irreversible mechanism of pure and hybridized non-Newtonian nanofluids through elastic surfaces with melting effects
  10. Stability analysis of the corruption dynamics under fractional-order interventions
  11. Solutions of certain initial-boundary value problems via a new extended Laplace transform
  12. Numerical solution of two-dimensional fractional differential equations using Laplace transform with residual power series method
  13. Fractional-order lead networks to avoid limit cycle in control loops with dead zone and plant servo system
  14. Modeling anomalous transport in fractal porous media: A study of fractional diffusion PDEs using numerical method
  15. Analysis of nonlinear dynamics of RC slabs under blast loads: A hybrid machine learning approach
  16. On theoretical and numerical analysis of fractal--fractional non-linear hybrid differential equations
  17. Traveling wave solutions, numerical solutions, and stability analysis of the (2+1) conformal time-fractional generalized q-deformed sinh-Gordon equation
  18. Influence of damage on large displacement buckling analysis of beams
  19. Approximate numerical procedures for the Navier–Stokes system through the generalized method of lines
  20. Mathematical analysis of a combustible viscoelastic material in a cylindrical channel taking into account induced electric field: A spectral approach
  21. A new operational matrix method to solve nonlinear fractional differential equations
  22. New solutions for the generalized q-deformed wave equation with q-translation symmetry
  23. Optimize the corrosion behaviour and mechanical properties of AISI 316 stainless steel under heat treatment and previous cold working
  24. Soliton dynamics of the KdV–mKdV equation using three distinct exact methods in nonlinear phenomena
  25. Investigation of the lubrication performance of a marine diesel engine crankshaft using a thermo-electrohydrodynamic model
  26. Modeling credit risk with mixed fractional Brownian motion: An application to barrier options
  27. Method of feature extraction of abnormal communication signal in network based on nonlinear technology
  28. An innovative binocular vision-based method for displacement measurement in membrane structures
  29. An analysis of exponential kernel fractional difference operator for delta positivity
  30. Novel analytic solutions of strain wave model in micro-structured solids
  31. Conditions for the existence of soliton solutions: An analysis of coefficients in the generalized Wu–Zhang system and generalized Sawada–Kotera model
  32. Scale-3 Haar wavelet-based method of fractal-fractional differential equations with power law kernel and exponential decay kernel
  33. Non-linear influences of track dynamic irregularities on vertical levelling loss of heavy-haul railway track geometry under cyclic loadings
  34. Fast analysis approach for instability problems of thin shells utilizing ANNs and a Bayesian regularization back-propagation algorithm
  35. Validity and error analysis of calculating matrix exponential function and vector product
  36. Optimizing execution time and cost while scheduling scientific workflow in edge data center with fault tolerance awareness
  37. Estimating the dynamics of the drinking epidemic model with control interventions: A sensitivity analysis
  38. Online and offline physical education quality assessment based on mobile edge computing
  39. Discovering optical solutions to a nonlinear Schrödinger equation and its bifurcation and chaos analysis
  40. New convolved Fibonacci collocation procedure for the Fitzhugh–Nagumo non-linear equation
  41. Study of weakly nonlinear double-diffusive magneto-convection with throughflow under concentration modulation
  42. Variable sampling time discrete sliding mode control for a flapping wing micro air vehicle using flapping frequency as the control input
  43. Error analysis of arbitrarily high-order stepping schemes for fractional integro-differential equations with weakly singular kernels
  44. Solitary and periodic pattern solutions for time-fractional generalized nonlinear Schrödinger equation
  45. An unconditionally stable numerical scheme for solving nonlinear Fisher equation
  46. Effect of modulated boundary on heat and mass transport of Walter-B viscoelastic fluid saturated in porous medium
  47. Analysis of heat mass transfer in a squeezed Carreau nanofluid flow due to a sensor surface with variable thermal conductivity
  48. Navigating waves: Advancing ocean dynamics through the nonlinear Schrödinger equation
  49. Experimental and numerical investigations into torsional-flexural behaviours of railway composite sleepers and bearers
  50. Novel dynamics of the fractional KFG equation through the unified and unified solver schemes with stability and multistability analysis
  51. Analysis of the magnetohydrodynamic effects on non-Newtonian fluid flow in an inclined non-uniform channel under long-wavelength, low-Reynolds number conditions
  52. Convergence analysis of non-matching finite elements for a linear monotone additive Schwarz scheme for semi-linear elliptic problems
  53. Global well-posedness and exponential decay estimates for semilinear Newell–Whitehead–Segel equation
  54. Petrov-Galerkin method for small deflections in fourth-order beam equations in civil engineering
  55. Solution of third-order nonlinear integro-differential equations with parallel computing for intelligent IoT and wireless networks using the Haar wavelet method
  56. Mathematical modeling and computational analysis of hepatitis B virus transmission using the higher-order Galerkin scheme
  57. Mathematical model based on nonlinear differential equations and its control algorithm
  58. Bifurcation and chaos: Unraveling soliton solutions in a couple fractional-order nonlinear evolution equation
  59. Space–time variable-order carbon nanotube model using modified Atangana–Baleanu–Caputo derivative
  60. Minimal universal laser network model: Synchronization, extreme events, and multistability
  61. Valuation of forward start option with mean reverting stock model for uncertain markets
  62. Geometric nonlinear analysis based on the generalized displacement control method and orthogonal iteration
  63. Fuzzy neural network with backpropagation for fuzzy quadratic programming problems and portfolio optimization problems
  64. B-spline curve theory: An overview and applications in real life
  65. Nonlinearity modeling for online estimation of industrial cooling fan speed subject to model uncertainties and state-dependent measurement noise
  66. Quantitative analysis and modeling of ride sharing behavior based on internet of vehicles
  67. Review Article
  68. Bond performance of recycled coarse aggregate concrete with rebar under freeze–thaw environment: A review
  69. Retraction
  70. Retraction of “Convolutional neural network for UAV image processing and navigation in tree plantations based on deep learning”
  71. Special Issue: Dynamic Engineering and Control Methods for the Nonlinear Systems - Part II
  72. Improved nonlinear model predictive control with inequality constraints using particle filtering for nonlinear and highly coupled dynamical systems
  73. Anti-control of Hopf bifurcation for a chaotic system
  74. Special Issue: Decision and Control in Nonlinear Systems - Part I
  75. Addressing target loss and actuator saturation in visual servoing of multirotors: A nonrecursive augmented dynamics control approach
  76. Collaborative control of multi-manipulator systems in intelligent manufacturing based on event-triggered and adaptive strategy
  77. Greenhouse monitoring system integrating NB-IOT technology and a cloud service framework
  78. Special Issue: Unleashing the Power of AI and ML in Dynamical System Research
  79. Computational analysis of the Covid-19 model using the continuous Galerkin–Petrov scheme
  80. Special Issue: Nonlinear Analysis and Design of Communication Networks for IoT Applications - Part I
  81. Research on the role of multi-sensor system information fusion in improving hardware control accuracy of intelligent system
  82. Advanced integration of IoT and AI algorithms for comprehensive smart meter data analysis in smart grids
Downloaded on 8.12.2025 from https://www.degruyterbrill.com/document/doi/10.1515/nleng-2024-0012/html
Scroll to top button