
Bistable composite laminates are smart composites that have been employed for engineering structures due to their superlative offering of features like ability to change shape and low densities. Because of the embedded geometrical nonlinearity factor, a small variation of input parameters leads to significant changes in the response of the bistable composite laminates. In other words, Uncertainty Quantification (UQ) makes a change in the bistability characteristics. As a result, bistability behavior is extremely reliant on geometrical dimensions and elastic material properties as design parameters. Reliability analysis deals with the quantitative assessment of the occurrence probability due to UQ. In this regard, the reliability and sensitivity analysis of bistable composite plate are investigated through the Monte Carlo Simulation (MCS) and multiple types of uncertain parameters, geometry and material properties, are assumed as random variables. The results indicate bistable composite plates have a high probability to be bistability behavior with the assumed statistical properties. Moreover, the sensitivity reliability analysis illustrates that the thickness and coefficient of thermal expansion have more effect on the bistability behavior in comparison to other input parameters. The results are confirmed by comparing them with those determined by the Finite Element Method (FEM).
Citation: Saeid Saberi, Azam Abdollahi, Fawad Inam. Reliability analysis of bistable composite laminates[J]. AIMS Materials Science, 2021, 8(1): 29-41. doi: 10.3934/matersci.2021003
[1] | Mateo Duarte, Johan A. Oquendo, Sebastián Vallejo, Johnattan Vargas, Yamile Cardona-Maya, Cesar A. Isaza . Fabrication of FRP/CNT hybrid laminate composites and their effect on interlaminar and mechanical properties. AIMS Materials Science, 2024, 11(6): 1125-1144. doi: 10.3934/matersci.2024054 |
[2] | Venkatarajan Subbarayalu, Subbu Chinnaraman, Athijayamani Ayyanar, Jayaseelan Chinnapalanisamy . Mechanical properties of vinyl ester hybrid composite laminates reinforced with screw pine and glass fiber. AIMS Materials Science, 2024, 11(1): 114-128. doi: 10.3934/matersci.2024007 |
[3] | Abhishek Sharma, Ashok Kumar Bagha, Dinesh Kumar Shukla, Shashi Bahl . Finite element model updating of metallic and composite structures-A state of the art review. AIMS Materials Science, 2021, 8(3): 390-415. doi: 10.3934/matersci.2021025 |
[4] | Marek Konieczny . Mechanical properties and failure analysis of laminated magnesium-intermetallic composites. AIMS Materials Science, 2022, 9(4): 572-583. doi: 10.3934/matersci.2022034 |
[5] | Mica Grujicic, Jennifer S. Snipes, S. Ramaswami . Penetration resistance and ballistic-impact behavior of Ti/TiAl3 metal/intermetallic laminated composites (MILCs): A computational investigation. AIMS Materials Science, 2016, 3(3): 686-721. doi: 10.3934/matersci.2016.3.686 |
[6] | Elisa Padovano, Francesco Trevisan, Sara Biamino, Claudio Badini . Processing of hybrid laminates integrating ZrB2/SiC and SiC layers. AIMS Materials Science, 2020, 7(5): 552-564. doi: 10.3934/matersci.2020.5.552 |
[7] | Kotaro Mori, Fumio Narita, Yasuhide Shindo . Detection and response characteristics of clamped-free giant magnetostrictive/piezoelectric laminates under concentrated loading. AIMS Materials Science, 2015, 2(4): 401-413. doi: 10.3934/matersci.2015.4.401 |
[8] | Marek Konieczny . Transformation superplasticity of laminated CuAl10Fe3Mn2 bronze-intermetallics composites. AIMS Materials Science, 2020, 7(3): 312-322. doi: 10.3934/matersci.2020.3.312 |
[9] | Adel M Bash, Sulaiman E. Mnawe, Salim A. Salah . Numerical buckling analysis of carbon fibre-epoxy composite plates with different cutouts number by finite element method. AIMS Materials Science, 2020, 7(1): 46-59. doi: 10.3934/matersci.2020.1.46 |
[10] | M. Rajanish, N. V. Nanjundaradhya, Ramesh S. Sharma, H. K. Shivananda, Alok Hegde . Directional Interlaminar Shear Strength (ILSS) of nano-modified epoxy/unidirectional glass fibre composite. AIMS Materials Science, 2018, 5(4): 603-613. doi: 10.3934/matersci.2018.4.603 |
Bistable composite laminates are smart composites that have been employed for engineering structures due to their superlative offering of features like ability to change shape and low densities. Because of the embedded geometrical nonlinearity factor, a small variation of input parameters leads to significant changes in the response of the bistable composite laminates. In other words, Uncertainty Quantification (UQ) makes a change in the bistability characteristics. As a result, bistability behavior is extremely reliant on geometrical dimensions and elastic material properties as design parameters. Reliability analysis deals with the quantitative assessment of the occurrence probability due to UQ. In this regard, the reliability and sensitivity analysis of bistable composite plate are investigated through the Monte Carlo Simulation (MCS) and multiple types of uncertain parameters, geometry and material properties, are assumed as random variables. The results indicate bistable composite plates have a high probability to be bistability behavior with the assumed statistical properties. Moreover, the sensitivity reliability analysis illustrates that the thickness and coefficient of thermal expansion have more effect on the bistability behavior in comparison to other input parameters. The results are confirmed by comparing them with those determined by the Finite Element Method (FEM).
Advanced composites are being employed for many industrial applications [1,2]. Bistable composite laminates are excellent candidates for their applications in morphing structures due to their unique features. They can change shapes and remain in natural equilibrium configuration without any external force. One of their most important applications is in aerospace structures such as airfoil section, air inlet demonstrator, morphing aerofoil and winglet and solar tracking device [3,4,5,6,7]. Moreover, they have considerable applications in energy harvesting fields, robotic gripper [8,9,10,11]. In order to use them effectively, different factors that affect their two stable states must be examined.
However, there are a wide range of undesirable uncertainties related to the design parameters in the composite laminate structures. They are random in nature due to their manufacturing or measuring errors. A probabilistic analysis can determine safer and reliable predictor of their structural response.
A large number of previous researches relating to bistable composite laminates have been published over the past 30 years [12]. It can be said that if the thermal load is applied to the unsymmetric composite laminates, they will have two cylindrical stable states in addition to an unstable saddle shape. This cannot be predicted with the classical lamination theory with these two configurations [13]. The analytical model, used in various reports, was based on the total potential energy and the Rayleigh–Ritz method using a geometrical nonlinearity [14,15]. In the theoretical approach, the displacement fields have a notable effect on the accuracy of the achieved results. A wide range of functions have been used for this aim, so far. Pirrera et al. studied the behavior of bistable composite shells by proposing the high-order polynomials and using a combination of the Ritz model with the path-following algorithm [16]. Lamacchia et al. estimated the transverse displacements by using Legendre polynomials and presented a new method based on the Rayleigh–Ritz approach for analyzing static and thermal behaviors of the multistable composite shells [17]. Material properties, geometry and environmental conditions have a considerable effect on the behavior of bistability of composite laminates. Vidoli and Maurini analyzed the tristability of orthotropic shallow shells with uniform initial curvature and investigated the range of initial curvatures and material parameters leading to tristability [18]. Brampton et al. used the total potential energy and the Rayleigh–Ritz methods and determined the influences of material properties, geometry and environmental conditions uncertainties on the curvature of the bistable composite laminate. The sensitivity of the stable shape was determined by considering a ±5% change in these parameters [19]. Che et al. studied the effect of the side length and stacking sequence on the stable configurations of the fiber metal laminates with different initial curvatures under thermal load using the analytical model, finite element and experimental methods [20]. Carey at el. studied the bistable helical structures by introducing an analytical method and investigated the variation of curvature respect to change in length for two different lay-ups [21].
The character of bistable laminates depends on the variation in temperature as well. Giddings et al. studied the impact of variation of laminate angle on the cured shape of bistable laminates under a range of temperature situations [22]. Eckstein et al. presented a model to estimate the cured shape of multistable plate subject to linear through-thickness temperature variation and for temperature-dependent material properties [23]. Zhang et al. presented an experimental study on the bistable behaviour of anti-symmetric laminated cylindrical shell at various temperatures. Based on the results, as the temperature rises, the snap-back load and curvatures increase but the snap-through force falls [24]. Fujioka et al. made bistable composite rods by the woven carbon fabrics composites and studied the influence of temperature on the deployment force of these structures. According to the result, the increase of temperature can significantly fall the deployment forces of bi-stable rods [25]. Chai et al. systematically analyzed the impact of moisture and temperature on bistable anti-symmetric composite shells and variable stiffness composite structures. The results showed that the temperature and moisture of bistable anti-symmetric composite shells are sensitive for different plies [26].
The aim of this study is to conduct reliability and sensitivity analysis with respect to random variables, involving material properties, geometrical features, which are studied by the combination of minimum potential energy, Rayleigh–Ritz and Monte Carlo methods. The equations of bistable composite laminates are derived by using the principle of minimizing potential energy and the Rayleigh–Ritz method. The limit state function is expressed via the context of bistability. The appropriate distribution is applied to variables and then the bistability probability is evaluated for cross-ply bistable composite by MCS method. Also, the uncertainty effect of the input parameters on out-off-plane displacement at the corner is studied using sensitivity analysis. In this work, FEM is employed to confirm these results as well.
The principle of minimum potential energy is used to establish the equations of the bistable composite plate. The total potential energy Π of the composite laminate is expressed by Eq 1.
Π=∫v12σTϵdv | (1) |
where v is the volume of the plate, σ is the stress tensor, the superscript "T" illustrates the transpose operator and the ϵ demonstrates the strain terms. According to the Classic Laminate Analysis (CLA), the stress tensor is given in Eq 2.
σ=−Q(ϵ−ϵt) | (2) |
where −Q is the reduced stiffness matrix, the superscript "t" illustrates the contribution of thermal deformations and ϵt is the thermal strain terms (ϵt=αΔT). By substituting Eq 2 into Eq 1, total potential energy is defined in Eq 3.
Π=∫v(12−Qijϵiϵj−−QijαiϵjΔT)dv | (3) |
with i,j=(1,2,6), where αi thermal expansion coefficients, ΔΤ is temperature difference between the cure and ambient temperatures.
By considering the assumption of small strains and moderate rotations, the strain field ∈ of bistable composite plates are described by the Von-Karman formulation in Eq 4.
∈=ϵ+zκ | (4) |
where z is the distance from the mid-plane of laminate. Also, ϵ is the mid-plane strains as ϵ={ϵxxϵyyϵxy}T and κ is the mid-plane curvatures as κ={κxxκyyκxy}T, which are expressed by Eq 5.
[ϵxxϵxxϵxx]=[∂U0∂x+12(∂W0∂x)2∂V0∂y+12(∂W0∂y)212(∂U0∂y+∂V0∂x+∂W0∂x∂W0∂y)],[κxxκyyκxy]=[−∂2W0∂x2−∂2W0∂y2−2∂2W0∂x∂y] | (5) |
where W0,U0,V0 illustrate the out-of-plane and in-plane displacements in the z, x and y and directions. In Eq 6, these displacement fields are provided.
U0(x,y)=m∑i=0n∑j=0tijy2j+1x2iV0(x,y)=m∑i=0n∑j=0sijy2jx2i+1W0(x,y)=m∑i=0n∑j=0rijy2jx2i | (6) |
where tij, sij and rij are unknown coefficients. Displacement fields proposed in Eq 3 can satisfy boundary conditions of the plate, which is fixed at center and has free edges.
By substituting Eq 6 in Eq 4 and then in Eq 3, the potential energy will turn into Eq 7.
Π(t)=Ly/2∫−Ly/2Lx/2∫−Lx/2h/2∫−h/2(12−Q11ϵ2xx+−Q12ϵxxϵyy+−Q16ϵxxϵxy+12−Q22ϵ2yy+−Q26ϵyyϵxy+12−Q66ϵ2xy−(−Q11αxx+−Q12αyy+−Q16αxy)ϵxxΔT−(−Q21αxx+−Q22αyy+−Q26αxy)ϵyyΔT−(−Q61αxx+−Q62αyy+−Q66αxy)ϵxyΔT)dzdxdy | (7) |
where Lx and Ly are the side lengths and h is the thickness of plate.
By integrating Eq 7, the total potential energy is obtained as a function of the coefficients of displacement fields. According to the principle of minimizing potential energy, to determine the equilibrium states, it is enough that its variations with respect to the unknown parameters obtained and the resulting expressions are set to zero, which is stated in Eq 8.
δΠ=∂Π∂rijδrij+∂Π∂sijδsij+∂Π∂tijδtij=0 | (8) |
Eq 8 is a system of nonlinear equations that are solved by using the Newton-Raphson numerical method in MATLAB software. There three types of solutions for this equation, which two of them are stable equilibrium states and other is an unstable equilibrium state. The equilibrium configuration of a plate will stable if the Jacobin matrix J of the potential energy of it be definite positive, according to Eq 9.
J=∂2Π∂X2 | (9) |
where X={tij,sij,rij}T is the unknown coefficients vector.
In the next section, the reliability and sensitivity analysis using Monte Carlo Simulation (MCS) are presented.
The uncertainty sources may lead to inappropriate behaviour of the engineering systems. The probabilistic evolution of engineering system may include a significant number of uncertainties in their behavior [27]. To implement probabilistic assessment for an engineering system, the number of identified issues such as correlation between random variables (RVs), high dimension problems, information about rare event scenarios and many interactive responses in the statement of performance criteria could be addressed [28].
Reliability analysis has been extensively used in engineering applications to evaluate the risk caused by such uncertainties [29]. In addition, it has a specific place in the process design of composite materials which leads to design more reliable laminated structures [30,31,32,33,34].
Generally, there are three approaches to estimate the failure probability, namely approximation (first and second order reliability method), simulation (MCS, subset simulation, important sampling, etc.) and adaptive surrogate-modeling-based (Adoptive kriging-MCS, Adaptive polynomial chaos kriging-MCS, etc. [35]. MCS is a robust tool for analyzing complex system. In the reliability problems with small failure probability, the computational cost is high but it is more accurate. The MCS failure probability Pf is defined in Eq 10.
Pf=∫g(x)≤0f(x)dx | (10) |
where f(x) is the joint density function of x1,…,xn. If the limit state function g(x)<0, then the failure is occurred. The integral failure probability Pf of Eq 10 using indicator function Ig(x)≤0(x) can be express by Eq 11.
Pf=E[Ig(x)≤0(x)≤0]=1NN∑k=1Ig(x)≤0(x)≤0=NfN | (11) |
where Ig(x)≤0(x)=1 if g(x)≤0, otherwise Ig(x)≤0(x)=0. The parameter N is the total number of simulations and Nf is number of simulation when g(x)<0. To enhance the efficiency of MCS analysis, some state-of-the-art methodologies can be applied as an alternative to estimate Pf such as stochastic finite element method (SFEM) based on stochastic perturbation [36] which is not under the scope of this study.
Figure 1 shows the schematic exhibition of the safe and failure domains in 2D reliability problem.
To investigation the accuracy of Pf can be used the coefficient of variation (COV) as Eq 12.
COV=√(1−Pf)NPf | (12) |
Sensitivity is applied to explore the rate of change in a model output due to changes in the model inputs in deterministic design, which is frequently calculated by partial derivative analytically or numerically. The sensitivity analysis provides a different meaning for random uncertainty.
The reliability sensitivity analysis is to investigate the rate of change in the probability characteristics of the response, especially Pf, due to the probability characteristics of a basic RV xi, such as ∂Pf∂μxi and ∂Pf∂σxi [37]. Eq 13 defines the MCS sensitivity analysis.
∂Pf∂ξi=∂∂ξi∫DFIg(x)≤0(x)f(x)dx | (13) |
where ξi indicated the probability characteristics of RVs, μxi and σxi, and DF is the domain of failure. The Eq 13 could be written as Eq 14 using the Lebesgue dominated convergence theorem [38].
∂Pf∂ξi=∂∂ξi∫DFIg(x)≤0(x)∂f(x)∂ξidx=∫DFIg(x)≤0(x)∂lnf(x)∂ξif(x)dx=E[Ig(x)≤0(x)∂lnf(x)∂ξi] | (14) |
where ∂lnf(x)∂ξi is introduced as a score function Sξi. Hence, the Eq 14 be converted into Eq 15.
∂Pf∂ξi=E[Ig(x)≤0(x)Sξi] | (15) |
In the current research, the bistability probability of cross-ply bistable composites plate is studied. So, the limit state function g(x) is presented by Eq 16.
g(x)=J=∂2Π∂Xi2 | (16) |
The limit state function in Eq 16 illustrates that g(x)≤0 is the unstable configuration or failure while g(x) > 0 is the stable configuration or non-failure.
The bistable composite plate is composed of Graphite/Epoxy (T300/5208) in the unsymmetric stacking sequences [0/90] with 150 mm × 150 mm side length. The means and the Coefficients of Variation (COV) of the geometry and material properties are listed in Table 1. In this work, nine independent random variables are considered. The uncertain material properties, laminate dimensions are considered as random variables which the distribution type is normal [39,40].
RV index | Property | Mean value | COV |
1 | E11 (GPa) | 146.95 | 0.1 |
2 | E22(GPa) | 10.702 | 0.1 |
3 | G12 (GPa) | 6.977 | 0.1 |
4 | ν12 | 0.3 | 0.1 |
5 | α11 (1/℃) | 5.028 × 10−7 | 0.1 |
6 | α22 (1/℃) | 2.65 × 10−5 | 0.1 |
7 | tply (m) | 365 × 10−6 | 0.1 |
8 | Lx (m) | 0.150 | 0.1 |
9 | Ly (m) | 0.150 | 0.1 |
To study the different behaviors of bistable composite laminates, the first step is to determine the static stable states. There are two cylindrical stable configurations for [0/90] bistable composite laminates so that they are identical and perpendicular to each other and have opposite directions. In other words, each state has two transverse curvatures which are the opposite. For cross-ply laminate one of them is much more than other and the twisting curvature is close to zero. These stable states obtained from analytical and finite element methods are presented in Figure 2.
The numerical results for the corner displacement obtained from analytical and finite element methods are listed in Table 2.
Type of laminate | FE | Analytical | Error (%) |
[0/90] | 11.35 | 11.54 | 1.6 |
*Note: Error=Analytical−FEMFEM×100. |
These results show that the analytical method has acceptable accuracy. The ABAQUS software is employed to study the bistable composite plate by finite element method, the simulated plate is meshed with the S4R element, which consists of 2500 elements and 2601 nodes. The "Static, General" step is used and "Nlgeom" is selected for static analysis due to there is geometric nonlinearity. Figure 3 illustrates the stable states of the desired cross-ply bistable plate after the curing process in the room temperature. The bistability characteristic depends on the geometry dimensions and material properties. In the next section, the effect of these parameters on bistability behavior will be investigated.
According to Table 1, the bistability probability is estimated about 0.0152 through MCS. The MCS is an accurate method which helps the experts to recognize the different factors on the reliability analysis and future behavior of complex problems. The bistability probability is estimated by 10000 samples and the COV is achieved about 0.080.
The sensitivities of bistability probability are shown in Figure 4. It is observed from Figure 4a that the mean values of the basic random variables α22, G12, E22, Lx and Ly have positive influences on reliability, i.e., the reliability improves with increase in mean values of these parameters. In fact, by increasing the mean value of these parameters the bistability probability is increased. However, the reliability reduces with the increase as mean values of tply, α11, E11 and ν12. The mean values of tply, α11 have the most effect on decreasing as well as α22 has the most impact on the increasing bistability probability. According to Figure 4b, standard deviations have the same effects as the mean values on the probability bistability, of course, their values are different.
The sensitivity of the corner displacement in the stable configuration to a ±10% change in each parameter is shown in Figure 5. According to that, α22 has the most influence upon stable shape. In other words, ±10% variation in α22 makes a 23.3% increase in the amount of maximum displacement at the corner. ±10% variation in tply and α11, however, lead to decrease in the amount of maximum displacement at the corner by 18.98% and 0.43% in maximum displacement, respectively.
The bistable composite laminates structure is affected by the design parameters such as their material properties and dimensions in practical application. In this paper, the reliability and sensitivity of the cross-ply composite plate were analyzed. The principle conclusions of this work are:
1.Reliability analysis estimates the bistable probability about 0.0152.
2.Sensitivity analysis indicates that thickness and coefficient of thermal expansions have the most important effect on bistability character of the composite laminates.
3.The transverse displacement at the corner increases with the rise of the transverse coefficient of thermal expansion. However, it decreases with the rise of the longitudinal coefficient of thermal expansion and the thickness of the plate.
Moreover, the finite element model of the bistable composite plate was established using ABAQUS software. It could be seen that the finite element analysis results were consistent with the theoretical analysis results. The impact of inter-dependency between RVs can be investigated in future studies. In addition, it's strongly recommended that probability analysis of multi-stable composite structures will be studied in the next work.
There is no conflict of interest.
[1] |
Emam SA, Inman DJ (2015) A review on bistable composite laminates for morphing and energy harvesting. Appl Mech Rev 67: 060803. doi: 10.1115/1.4032037
![]() |
[2] |
Chillara VSC, Dapino MJ (2020) Review of morphing laminated composites. Appl Mech Rev 72: 010801. doi: 10.1115/1.4044269
![]() |
[3] |
Daynes S, Weaver P, Trevarthen J (2011) A morphing composite air inlet with multiple stable shapes. J Intel Mat Syst Str 22: 961–973. doi: 10.1177/1045389X11399943
![]() |
[4] |
Nicassio S, Scarselli G, Pinto F, et al. (2018) Low energy actuation technique of bistable composites for aircraft morphing. Aerospac Sci Technol 75: 35–46. doi: 10.1016/j.ast.2017.12.040
![]() |
[5] |
Kuder IK, Fasel U, Ermanni P, et al. (2016) Concurrent design of a morphing aerofoil with variable stiffness bi-stable laminates. Smart Mater Struct 25: 115001. doi: 10.1088/0964-1726/25/11/115001
![]() |
[6] |
Arrieta AF, Bilgen O, Friswell MI, et al. (2013) Modelling and configuration control of wing-shaped bi-stable piezoelectric composites under aerodynamic loads. Aerospa Sci Technol 29: 453–461. doi: 10.1016/j.ast.2013.05.004
![]() |
[7] |
Zhang Z, Pei K, Sun M, et al. (2020) A novel solar tracking model integrated with bistable composite structures and bimetallic strips. Compos Struct 148: 112506. doi: 10.1016/j.compstruct.2020.112506
![]() |
[8] | Kim SW, Koh J, Cho M, et al. (2011) Design & analysis a flytrap robot using bi-stable composite, 2011 IEEE International Conference on Robotics and Automation, IEEE, 215–220. |
[9] |
Zhang Z, Li X, Yu X, et al. (2019) Magnetic actuation bionic robotic gripper with bistable morphing structure. Compos Struct 229: 111422. doi: 10.1016/j.compstruct.2019.111422
![]() |
[10] |
Lee AJ, Inman DJ (2018) A multifunctional bistable laminate: Snap-through morphing enabled by broadband energy harvesting. J Intel Mat Syst Str 29: 2528–2543. doi: 10.1177/1045389X18770895
![]() |
[11] |
Rehman TU, Qaiser Z, Johnson S (2019) Tuning bifurcation loads in bistable composites with tunable stiffness mechanisms. Mech Mach Theory 142: 103585. doi: 10.1016/j.mechmachtheory.2019.103585
![]() |
[12] |
Zhang Z, Li Y, Yu X, et al. (2019) Bistable morphing composite structures: A review. Thin Wall Struct 142: 74–97. doi: 10.1016/j.tws.2019.04.040
![]() |
[13] |
Hyer MW (1981) Some observations on the cured shape of thin unsymmetric laminates. J Compos Mater 15: 175–194. doi: 10.1177/002199838101500207
![]() |
[14] |
Hamamoto A, Hyer MW (1987) Non-linear temperature-curvature relationships for unsymmetric graphite-epoxy laminates. Int J Solids Struct 23: 919–935. doi: 10.1016/0020-7683(87)90087-4
![]() |
[15] |
Peeters L, Powell P, Warnet L (1996) Thermally-induced shapes of unsymmetric laminates. J Compos Mater 30: 603–626. doi: 10.1177/002199839603000504
![]() |
[16] |
Pirrera A, Avitabile D, Weaver PM (2010) Bistable plates for morphing structures: a refined analytical approach with high-order polynomials. Int Jo Solids Struct 47: 3412–3425. doi: 10.1016/j.ijsolstr.2010.08.019
![]() |
[17] |
Lamacchia E, Pirrera A, Chenchiah IV, et al. (2015) Morphing shell structures: a generalised modelling approach. Compos Struct 131: 1017–1027. doi: 10.1016/j.compstruct.2015.06.051
![]() |
[18] | Vidoli S, Maurini C (2008) Tristability of thin orthotropic shells with uniform initial curvature. P Roy Soc A-Math Phy 464: 2949–2966. |
[19] |
Brampton CJ, Betts DN, Bowen CR, et al. (2013) Sensitivity of bistable laminates to uncertainties in material properties, geometry and environmental conditions. Compos Struct 102: 276–286. doi: 10.1016/j.compstruct.2013.03.005
![]() |
[20] |
Che L, Fang G, Wu Z, et al. (2020) Investigation of curing deformation behavior of curved fiber metal laminates. Compos Struct 232: 111570. doi: 10.1016/j.compstruct.2019.111570
![]() |
[21] | Carey S, Robert T, Vincenzo O, et al. (2019) Bistable composite helices with thermal effects. P Roy Soc A-Math Phy 475: 20190295. |
[22] |
Giddings PF, Bowen CR, Salo AIT, et al. (2010) Bistable composite laminates: effects of laminate composition on cured shape and response to thermal load. Compos Struct 92: 2220–2225. doi: 10.1016/j.compstruct.2009.08.043
![]() |
[23] |
Eckstein E, Pirrera A, Weaver P (2013) Morphing high-temperature composite plates utilizing thermal gradients. Compos Struct 100: 363–372. doi: 10.1016/j.compstruct.2012.12.049
![]() |
[24] |
Zhang Z, Wu H, Ye G, et al. (2016) Experimental study on bistable behaviour of anti-symmetric laminated cylindrical shells in thermal environments. Compos Struct 144: 24–32. doi: 10.1016/j.compstruct.2016.02.062
![]() |
[25] |
Fujioka E, Yokozeki T, Watanabe A, et al. (2019) Analysis on temperature-dependent deployment behavior of bi-stable composite rods. Adv Compos Mater 28: 245–257. doi: 10.1080/09243046.2018.1499411
![]() |
[26] |
Chai H, Li Y, Zhang Z, et al. (2020) Systematic analysis of bistable anti-symmetric composite cylindrical shells and variable stiffness composite structures in hygrothermal environment. Int J Adv Manuf Tech 108: 1091–1107. doi: 10.1007/s00170-019-04887-5
![]() |
[27] |
Xiao S, Reuschen S, Köse G, et al. (2019) Estimation of small failure probabilities based on thermodynamic integration and parallel tempering. Mech Syst Signal Pr 133: 106248. doi: 10.1016/j.ymssp.2019.106248
![]() |
[28] |
Tee KF, Khan LR, Li H (2014) Application of subset simulation in reliability estimation of underground pipelines. Reliab Eng Syst Safe 130: 125–131. doi: 10.1016/j.ress.2014.05.006
![]() |
[29] | Liu B (2010) Uncertain risk analysis and uncertain reliability analysis. J Uncertain Syst 4: 163–170. |
[30] | Yoo K, Bacarreza O, Aliabadi MF (2020) A novel multi-fidelity modelling-based framework for reliability-based design optimisation of composite structures. Eng Comput. Available from: https://link.springer.com/article/10.1007/s00366-020-01084-x. |
[31] |
Su H, Kam T (2020) Reliability analysis of composite wind turbine blades considering material degradation of blades. Compos Struct 234: 111663. doi: 10.1016/j.compstruct.2019.111663
![]() |
[32] | Shi B, Den Z, Tan LS, et al. (2020) Micromechanics-based reliability analysis method for laminated composite structures. Mech Adv Mater Struct. Available from: https://doi.org/10.1080/15376494.2020.1717022. |
[33] |
Sadeghian E, Toosi S (2019) Reliability analysis of a composite laminate using estimation theory. J Cent South Univ 26: 665–672. doi: 10.1007/s11771-019-4037-3
![]() |
[34] |
Peng X, Yuliang G, Chan Q, et al. (2020) Reliability optimization design for composite laminated plate considering multiple types of uncertain parameters. Eng Optimiz 53: 221–236. doi: 10.1080/0305215X.2019.1705289
![]() |
[35] | Marelli S, Schöbi R, Sudret B (2017) Uqlab user manual-structural reliability. Report UQLab-V1, 1–107. Available from: https://usermanual.wiki/Document/UserManualReliability.1348874550. |
[36] |
Kamiński M, Świta P (2015) Structural stability and reliability of the underground steel tanks with the stochastic finite element method. Arch Civ Mech Eng 15: 593–602. doi: 10.1016/j.acme.2014.04.010
![]() |
[37] |
Song S, Lu Z, Qiao H (2009) Subset simulation for structural reliability sensitivity analysis. Reliab Eng Syst Safe 94: 658–665. doi: 10.1016/j.ress.2008.07.006
![]() |
[38] |
MiarNaeimi F, Azizyan G, Rashki M (2019) Reliability sensitivity analysis method based on subset simulation hybrid techniques. Appl Math Model 75: 607–626. doi: 10.1016/j.apm.2019.05.038
![]() |
[39] |
Blanchard J, Sobey A (2019) Comparative design of E-glass and flax structures based on reliability. Compos Struct 225: 111037. doi: 10.1016/j.compstruct.2019.111037
![]() |
[40] |
Li X, Lv Z, Qiu Z (2018) A novel univariate method for mixed reliability evaluation of composite laminate with random and interval parameters. Compos Struct 203: 153–163. doi: 10.1016/j.compstruct.2018.06.097
![]() |
1. | S. Saberi, M. Ghayour, H. R. Mirdamadi, M. Ghamami, Free vibration analysis and mode management of bistable composite laminates using deep learning, 2021, 0939-1533, 10.1007/s00419-021-01930-4 | |
2. | Saeid Saberi, Alireza Sadat Hosseini, Fatemeh Yazdanifar, Saullo G. P. Castro, Developing Equations for Free Vibration Parameters of Bistable Composite Plates Using Multi-Objective Genetic Programming, 2022, 14, 2073-4360, 1559, 10.3390/polym14081559 | |
3. | Vishrut Deshpande, Oliver Myers, Suyi Li, Jae-Hung Han, Shima Shahab, Jinkyu Yang, 2022, Switchable structures using asymmetric fiber composite laminates: two case studies, 9781510649613, 33, 10.1117/12.2612988 | |
4. | Rebai Billel, Contribution to study the effect of (Reuss, LRVE, Tamura) models on the axial and shear stress of sandwich FGM plate (Ti–6A1–4V/ZrO2) subjected on linear and nonlinear thermal loads, 2023, 10, 2372-0484, 26, 10.3934/matersci.2023002 | |
5. | Ammar Elsheikh, Bistable Morphing Composites for Energy-Harvesting Applications, 2022, 14, 2073-4360, 1893, 10.3390/polym14091893 | |
6. | Saeid Saberi, Azam Abdollahi, Michael I. Friswell, Probability analysis of bistable composite laminates using the subset simulation method, 2021, 271, 02638223, 114120, 10.1016/j.compstruct.2021.114120 | |
7. | Petro Pukach, Roman Kvit, Tetyana Salo, Myroslava Vovk, A Probable Approach to Reliability Assessment of Reinforced Plates, 2023, 6, 2571-5577, 73, 10.3390/asi6040073 | |
8. | Saeid Saberi, Hamid Nasiri, Omid Ghorbani, Michael I. Friswell, Saullo G. P. Castro, Explainable Artificial Intelligence to Investigate the Contribution of Design Variables to the Static Characteristics of Bistable Composite Laminates, 2023, 16, 1996-1944, 5381, 10.3390/ma16155381 | |
9. | K.S. Suraj, P.M. Anilkumar, C.G. Krishnanunni, B.N. Rao, Uncertainty quantification of bistable variable stiffness laminate using machine learning assisted perturbation approach, 2023, 319, 02638223, 117072, 10.1016/j.compstruct.2023.117072 | |
10. | K. S. Suraj, P. M. Anilkumar, C. G. Krishnanunni, B. N. Rao, 2023, Chapter 12, 978-981-99-5048-5, 141, 10.1007/978-981-99-5049-2_12 | |
11. | Paulomi Mukherjee, Banamali Sahu, A. Arockiarajan, Shaikh Faruque Ali, Design optimization of hybrid bistable symmetric laminate using machine learning, 2025, 1537-6494, 1, 10.1080/15376494.2025.2488512 |
RV index | Property | Mean value | COV |
1 | E11 (GPa) | 146.95 | 0.1 |
2 | E22(GPa) | 10.702 | 0.1 |
3 | G12 (GPa) | 6.977 | 0.1 |
4 | ν12 | 0.3 | 0.1 |
5 | α11 (1/℃) | 5.028 × 10−7 | 0.1 |
6 | α22 (1/℃) | 2.65 × 10−5 | 0.1 |
7 | tply (m) | 365 × 10−6 | 0.1 |
8 | Lx (m) | 0.150 | 0.1 |
9 | Ly (m) | 0.150 | 0.1 |
Type of laminate | FE | Analytical | Error (%) |
[0/90] | 11.35 | 11.54 | 1.6 |
*Note: Error=Analytical−FEMFEM×100. |
RV index | Property | Mean value | COV |
1 | E11 (GPa) | 146.95 | 0.1 |
2 | E22(GPa) | 10.702 | 0.1 |
3 | G12 (GPa) | 6.977 | 0.1 |
4 | ν12 | 0.3 | 0.1 |
5 | α11 (1/℃) | 5.028 × 10−7 | 0.1 |
6 | α22 (1/℃) | 2.65 × 10−5 | 0.1 |
7 | tply (m) | 365 × 10−6 | 0.1 |
8 | Lx (m) | 0.150 | 0.1 |
9 | Ly (m) | 0.150 | 0.1 |
Type of laminate | FE | Analytical | Error (%) |
[0/90] | 11.35 | 11.54 | 1.6 |
*Note: Error=Analytical−FEMFEM×100. |