Next Article in Journal
Prediction of Structural Behavior of Continuous Reinforced Concrete Beams with Hybrid CFRP-Steel Bars
Next Article in Special Issue
The Effect of Different Morphologies of WO3/GO Nanocomposite on Photocatalytic Performance
Previous Article in Journal
Nanostructured Hybrid Hydrogels for Solar-Driven Clean Water Harvesting from the Atmosphere
Previous Article in Special Issue
Nonlocal Free Vibration of Embedded Short-Fiber-Reinforced Nano-/Micro-Rods with Deformable Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bimetallic Thin-Walled Box Beam Thermal Buckling Response

Department of Engineering Mechanics, Faculty of Engineering, University of Rijeka, 51000 Rijeka, Croatia
*
Author to whom correspondence should be addressed.
Submission received: 28 September 2022 / Revised: 11 October 2022 / Accepted: 24 October 2022 / Published: 27 October 2022
(This article belongs to the Special Issue Feature Papers in Materials Simulation and Design)

Abstract

:
A beam model for thermal buckling analysis of a bimetallic box beam is presented. The Euler–Bernoulli–Vlasov beam theory is employed considering large rotations but small strains. The nonlinear stability analysis is performed using an updated Lagrangian formulation. In order to account for the thermal effects of temperature-dependent (TD) and temperature-independent (TID) materials, a uniform temperature rise through beam wall thickness is considered. The numerical results for thin-walled box beams are presented to investigate the effects of different boundary conditions, beam lengths and material thickness ratios on the critical buckling temperature and post-buckling responses. The effectiveness and accuracy of the proposed model are verified by means of comparison with a shell model. It is revealed that all of the abovementioned effects are invaluable for buckling analysis of thin-walled beams under thermal load. Moreover, it is shown that the TD solutions give lower values than the TID one, emphasizing the importance of TD materials in beams.

1. Introduction

Thin-walled beams and structures are increasingly used in engineering branches, in standalone forms and as a stiffeners for plate- and shell-like structures, due to their high strength and light weight. However, these structures show susceptibility to local buckling and buckling failure [1,2]. Buckling analysis and the post-buckling response of such weight-optimized structures have been the topic of many research papers, such as [3,4,5,6], especially in the field of composite materials [7,8,9,10,11].
If the thermal environment is considered, the stability of structures has received significant attention in recent years: Duan et al. [12] performed thermal analysis of a beam element, Saha and Ali [13] presented a post-buckling mathematical model of a slender road under uniform temperature rise, while Cui and Hu [14] analyzed the thermal buckling and vibration of a beam. Jeyaraj et al. [15] investigated experimental and theoretical non-uniform heating of an isotropic beam. Burgreen and Mannit [16] and Burgreen and Regal [17] analysed the thermal buckling of bimetallic beams. In the case of composite beams, Aydogdu [18] obtained critical buckling temperatures of composite beams, Luan et al. [19] presented an analytical solution for buckling and vibration of FG beams, Kiani and Eslami [20,21] investigated buckling analysis under different types of thermal loads, while Giunta [22] analyzed FG beams under thermal/mechanical load using the Carrera unified formulation. However, there are not many papers about thermal buckling analyses of thin-walled structures: Libresceu [23] studied stability problems in a high-temperature environment and Ziane et al. [24] studied analytical methods for buckling and vibration responses of porous beams under thermomechanical loads.
In the present work, thermal buckling analysis of a thin-walled bimetallic box beam and frame structures is presented. The material is assumed to be linear, elastic and isotropic. The model is based on Euler–Bernoulli–Vlasov theory and on assumptions of large rotations and small strains. It is also assumed that the cross-section is not deformed in its own plane and that there are no shear strains in the middle surface. The nonlinear displacement field, which includes nonlinear displacement terms due to large rotation effects, is implemented. Using the UL description, the element geometric stiffness is derived. As an incremental iterative solution scheme, the Newton–Raphson method is used. Furthermore, this paper is a continuation of the research in which thermal buckling analysis of temperature-independent materials was conducted [25], which has now been further expanded with temperature-dependent materials’ properties. As far as the authors are aware, there is no beam model solution for thermal buckling analysis of thin-walled beam-type structures with temperature-dependent materials’ properties.
The numerical results for thin-walled box beams are presented to investigate the effects of different boundary conditions, namely clamped–clamped, simply supported and clamped–simply supported, beam lengths and material thickness ratios on the critical buckling temperature and post buckling responses. In order to demonstrate the accuracy of the numerical algorithm, benchmark examples using shell FEM code were developed. Numerical results show that the abovementioned effects have a huge impact on the buckling analysis.

2. Materials and Methods

2.1. Kinematics

Two sets of coordinate systems related to the angle of orientation β are considered. The first one is a Cartesian (z, x, y) coordinate system where the z-axis coincides with the longitudinal beam that passes through the centroid O of each cross section, while the x- and y-axes are principal axes. The second one is a contour coordinate system where the s-axis is tangential to the middle surface directed along the contour line of the cross-section while the n-axis is perpendicular to the s-axis.
The field of incremental displacement measures of a cross section are defined as [8]:
w 0 = w 0 ( z ) ;   u s = u s ( z ) ;   v s = v s ( z ) ;   φ z = φ z ( z ) ; φ x = φ x ( z ) = d v s d z ; φ y = φ y ( z ) = d u s d z ; θ = θ ( z ) = d φ z d z ( z ) ,
where w 0 ,   u s and v s are the rigid-body translations of the cross-section centroid in the z-, x- and y-direction, respectively, while φ z , φ x and φ y are the rigid-body rotations about the aforementioned axis; θ is a warping parameter of the cross-section.
In the case of small rotations, the incremental displacement field consists of the first-order displacement values:
u z ( z , x , y ) = w 0 ( z ) y d v s d z ( z ) x d u s d z ( z ) ω ( x , y ) d φ z d z ( z ) , u x ( z , x , y ) = u s ( z ) ( y y s ) φ z ( z ) u y ( z , x , y ) = v s ( z ) + ( x x s ) φ z ( z ) ,
where u z , u x and u y are the linear displacement increments of an arbitrary point on the cross-section defined by the x and y coordinates and the warping function ω ( x , y ) . When the large rotations are considered, nonlinear displacement increments are expressed as follows:
u ˜ z ( z , x , y ) = 0.5 [ ( x x s ) φ z φ x + ( y y s ) φ z φ y ] , u ˜ x ( z , x , y ) = 0.5 { φ x φ y y [ φ z 2 + φ y 2 ] x + x s φ z 2 } , u ˜ y ( z , x , y ) = 0.5 { φ x φ y x [ φ z 2 + φ x 2 ] y + y s φ z 2 } ,  
and should be added to those from Equation (2).
Considering the nonlinear displacement field, the Green–Lagrange strain tensor components can be written as:
ε i j = 1 2 [ ( u i + u ˜ i ) , j + ( u j + u ˜ j ) , i + ( u k + u ˜ k ) i + ( u k + u ˜ k ) , j ] e i j + η i j + e ˜ i j , 2 e i j = u i , j + u j , i 2 η i j = u k , i + u k , j 2 e ˜ i j = u ˜ i , j + u ˜ j , i  
where e i j and η i j are the linear and nonlinear strain components corresponding to the linear displacement, while e ˜ i j is the linear strain component corresponding to the nonlinear displacement due to the large rotations.
The contour mid-line displacement w ¯ , and v ¯ can be seen more detail in [26].
Due to the in-plane rigidity hypothesis of the cross-section, the non-zero strain components are [7]:
e z z = w z   ,     e z s = w s + v z   ,  
η z z = 1 2 [ ( w z ) 2 + ( u z ) 2 + ( v z ) 2 ] , η z s = w z w s + u z u s + v z v s ,
e ˜ z z = w z   ,     e ˜ z s = w s + v z .  
The stress resultants of the beam can be defined as:
F z = A σ z d n d s ,   M x = A σ z ( y n cos β ) d n d s ,     M y = A σ z ( x + n sin β ) d n d s ,     M t = A τ s z ( n + F s t ) d n d s ,     M ω = A σ z ( ω n q ) d n d s ,
where F z represents axial force, M x and M y are bending moments with respect to x- and y-axis, respectively, M t is the torsion moment and M ω is the warping moment (bimoment). t is the thickness of the closed section contour and F s is the St. Venant circuit flow [26].

2.2. Constitutive Equations

Consider a bimetallic beam made of two different metals. A beam wall with a core thickness of Ti–6Al–4V λ t on the outer surface and SUS304 in the inner part of cross-section beam wall is shown in Figure 1. It is assumed that the layers of materials are perfectly bonded.
If the thermo-elastic material properties are considered as a function of temperature T, they can be calculated for each material, as described in [19,27]:
P ( T ) = P 0 ( 1 + P 1 T 1 + P 1 T + P 2 T 2 + P 3 T 3 ) ,
where P represents Young’s modulus E and thermal expansion coefficient α , while P 0 , P 1 ,   P 1 , P 2 and P 3 are temperature-dependent coefficients listed in Table 1 for different metals [28]. For simplicity, Poisson’s ratio ν is assumed to be constant, ν = 0.3 .
It is assumed that the temperature of the whole beam is uniform and increased from the current ambient temperature T 0 to the critical value in incremental steps of 1 °C. If the axial displacements are prevented, the temperature at a point T ( n , z ) may be raised to T + Δ T , in the way that the beam buckles. Δ T is the temperature rise. The temperature that is read as the critical buckling temperature is the temperature difference compared to the ambient temperature. The process can be described as quasi-adiabatic since the heat exchange between the environment and the beam is neglected.
The stress–strain relations of the bimetallic beam can be written as:
( σ z τ s z ) = ( E ( n , z , T ) 0 0 G ( n , z , T ) ) · ( e z z α ( n , z , T ) Δ T γ s z ) ,
where
G ( n , z ) = E ( n , z , T ) 2 [ 1 + ν ( n , z , T ) ] .
Using Equations (5), (10) and (11), the beam forces can be expressed in a matrix form as:
{ F z M y M x M ω M t } = [ R 11 R 12 R 13 R 14 0 R 21 R 22 R 23 R 24 0 R 31 R 32 R 33 R 34 0 R 41 R 42 R 43 R 44 0 0 0 0 0 R 55 ] { d w 0 / d z d 2 u s / d z 2 d 2 v s / d z 2 d 2 φ z / d z 2 2 d φ z / d z } { N z T M y T M x T M ω T 0 } ,
where R i j represents the thin-walled beam stiffness, as shown in Appendix A. N z T ,   M x T ,   M y T and M ω T are thermal force and thermal moments, respectively:
N z T = A E ( n , z , T ) α ( n , z , T ) Δ T d n d s ; M x T = A E ( n , z , T ) α ( n , z , T ) Δ T ( y n cos β ) d n d s M y T = A E ( n , z , T ) α ( n , z , T ) Δ T ( x + n sin β ) d n d s M ω T = A E ( n , z , T ) α ( n , z , T ) Δ T ( ω n q ) d n d s

2.3. Finite Element Formulation

A two-nodded beam element with 14 degrees of freedom is shown in Figure 2. The nodal displacements and nodal force vectors are as follows:
( u e ) T = { w A   u A   v A   φ z A   φ x A   φ y A   w B   u B   v B   φ z B   φ x B   φ y B   θ A   θ B }  
( f e ) T = { F z A   F x A   F y A   M z A   M x A   M y A   F z B   F x B   F y B   M z B   M x B   M y B   M ω A   M ω B }
where the superscript e denotes the eth finite element. It should be noted that nodal displacement w and nodal forces F z ,   M x and M y are defined in the centroid O , while other nodal components are defined in the shear center.
Applying the principle of the virtual work, the incremental equilibrium equations of a beam element in linearized form are:
δ U E + δ U G = δ 2 W δ 1 W ,  
where the equations from the left side consist of incremental virtual elastic strain energy:
δ U E = 1 V C 1 i j k l   1   e k l   δ 1   e i j 1   d V ,  
and the incremental virtual geometric potential:
δ U G = 1 V S 1 i j   δ 1   η i j 1   d V + V S 1 i j   δ 1   e ˜ i j 1   d V A 1 σ t 1 i   δ u ˜ i 1   d A σ .  
On the right side of the equations, the terms represent the virtual work completed by external forces at the end and at the beginning of the present increment:
δ 2 W = A 1 σ t 2 i   δ u i 1   d A σ , δ 1 W = V 1 S 1 i j   δ 1   e i j 1   d V = A 1 σ t 1 i   δ u i 1   d A σ .  
In these equations, S i j is the second Piola–Kirchoff stress tensor, t i denotes the surface tractions, C i j k l presents the stress–strain tensor and the symbol δ indicates virtual quantities. By applying the linear interpolation functions for w 0 displacement and cubic interpolations for w s , u s and v s , one can obtain:
δ U E = V S ij δ e ij d V = ( δ u e ) T k E e u e ,  
δ U G = V S ij ( δ η ij + δ e ˜ ij ) d V A σ t i δ u ˜ i d A σ = ( δ u e ) T k G e u e
δ W = A σ t i δ u i d A σ = ( δ u e ) T f e  
where f e is the nodal force vector, k E e is the elastic stiffness matrix and k G e is the geometric stiffness matrix of the beam element. Nonlinear equilibrium equations are solved using the Newton–Raphson method as incremental iterative approach [29,30], and the explicit form of the terms given in nonlinear components were described previously in [31].

3. Results and Discussion

In numerical examples, the thin-walled box beam with height h = 100   mm , width b = 150   mm and thickness t = 10   mm is considered (Figure 3). For verification purposes, the critical buckling temperatures were obtained by shell FEM commercial code [32]. In order to simulate the bimetallic material, the beam walls were divided into two layers of different metals with a variable thickness ratio λ. Note that for λ = 0 , the beam wall is fully SUS304, while as the index thickness ratio λ increases, the beam wall becomes fully Ti-6Al-4V.

3.1. Box Beam

In the first example, the eigenvalue results of the box beam for different boundary conditions, which are clamped–clamped (C-C), clamped–simply supported (C-S) and simply supported (S-S), beam lengths of L 1 = 6   m , L 2 = 8   m and L 3 = 10   m and different material thickness ratios λ are given in Table 2, Table 3 and Table 4. Critical buckling temperatures are given for the first two flexural buckling modes. As the thickness ratio λ increases, the critical buckling temperatures increase as well due to the material properties of TI-6Al-4V. Furthermore, longer beams obtain lower critical buckling temperatures. As expected, clamped–clamped beams exhibit the highest eigenvalues for every beam length. Good agreement of the present results and solutions derived from the 2D model for both flexural buckling modes is achieved.
In order to further perform the nonlinear stability analysis of the box beam, a perturbation force of intensity Δ F = 500   N is introduced in the y-axis direction at the midpoint of the clamped–clamped beam. Temperature–displacement curves of the shortest beam with the comparison of temperature-dependent (TD) and temperature-independent (TID) materials are shown in Figure 4, Figure 5 and Figure 6. The results are given for different thickness ratios: λ = 0.2 , λ = 0.5 and λ = 0.8 . It can be seen that curves match very well with the critical buckling temperatures achieved in the eigenvalue manner. As expected, TD solutions obtained lower critical buckling temperatures. In the case of λ = 0.2 , the difference in critical temperature is around 4%, for λ = 0.5 the difference is 3.3%, and for λ = 0.8 it is 1.5%.

3.2. L-Frame

Furthermore, the model is tested for thermal buckling analysis of an L-frame with the length of both legs being L = 5   m and cross-section described in the previous chapter (Figure 7). The frame is fixed at points A and C, while at point B, in-plane translations are prevented.
To verify the results, the critical buckling temperatures for the full SUS304,   Δ T c r ( S U S 304 ) = 16.09   , and full Ti-6Al-4V Δ T c r ( Ti 6 Al 4 V ) = 35.52   , sections are solved by a shell commercial code. To perform nonlinear analysis, a small perturbation force Δ F = 50   N acting in the z-axis direction at point B is applied. The results are shown in Figure 8 for pure metals and for λ = 0.2 t , λ = 0.5 t and λ = 0.8 t . The good recognition of the critical values can be noted. It can be observed that with an increase in the proportion of Ti-6Al-4V material, higher critical temperatures are achieved.

4. Conclusions

A thin-walled beam model capable of thermal buckling analysis has been presented. By means of the updated Lagrangian formulation, the incremental equilibrium equations have been developed using the nonlinear displacement field of the cross-section, taking into account the effects of large rotations. The reliability of the present model was verified by studying the benchmark examples, and the values obtained with the proposed model are in good agreement with those of the shell model. The effects of boundary conditions, the length of the beam and material thickness ratio on the critical buckling temperature and post-buckling response are of great importance. Additionally, it is shown that the TD solutions provide lower values than those of the TID solutions. The model was found to be efficient in predicting eigenvalues and nonlinear buckling behavior.

Author Contributions

Conceptualization, S.K.S. and D.L.; methodology, S.K.S., G.T. and D.L.; software, S.K.S. and D.B.; validation, S.K.S. and D.B.; formal analysis, S.K.S., G.T. and D.L.; investigation, S.K.S.; resources, D.L. and G.T.; data curation, S.K.S.; writing—original draft preparation, S.K.S.; writing—review and editing, S.K.S. and D.L.; visualization, S.K.S.; supervision D.L. and G.T.; project administration, D.L.; funding acquisition, D.L. and G.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Croatian Science Foundation, grant number: IP-2019-04-8615 and University of Rijeka, grant numbers: uniri-tehnic-18-107 and uniri-tehnic-18-139.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The research presented in this paper was made possible by the financial support of Croatian Science Foundation.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

R 11 = A E ( n , z , T ) d n d s , R 12 = R 21 = A E ( n , z , T ) ( x + n sin β ) d n d s , R 13 = R 31 = A E ( n , z , T ) ( y n cos β ) d n d s , R 14 = R 41 = A E ( n , z , T ) ( ω n q ) d n d s , R 22 = A E ( n , z , T ) ( x + n sin β ) 2 d n d s , R 23 = R 32 = A E ( n , z , T ) ( x + n sin β ) ( y n cos β ) d n d s , R 23 = R 32 = A E ( n , z , T ) ( x + n sin β ) ( y n cos β ) d n d s , R 24 = R 42 = A E ( n , z , T ) ( x + n sin β ) ( ω n q ) d n d s , R 33 = A E ( n , z , T ) ( y n cos β ) 2 d n d s , R 34 = R 43 = A E ( n , z , T ) ( y n cos β ) ( ω n q ) d n d s , R 44 = A E ( n , z , T ) ( ω n q ) 2 d n d s , R 55 = A G ( n , z , T ) ( n + F s t ) 2 d n d s .

References

  1. Trahair, N.S. Flexural-Torsional Buckling of Structures; CRC Press: Boca Raton, FL, USA, 1993. [Google Scholar]
  2. Alfutov, N.A. Stability of Elastic Structures; Springer: Berlin/Heidelberg, Germany, 2000. [Google Scholar]
  3. Bin Kamarudin, M.N.; Mohamed Ali, J.S.; Aabid, A.; Ibrahim, Y.E. Buckling Analysis of a Thin-Walled Structure Using Finite Element Method and Design of Experiments. Aerospace 2022, 9, 541. [Google Scholar] [CrossRef]
  4. Turkalj, G.; Brnić, J. Nonlinear Stability Analysis of Thin-Walled Frames Using Ul–Esa Formulation. Int. J. Struct. Stab. Dyn. 2004, 4, 45–67. [Google Scholar] [CrossRef]
  5. Nguyen, N.D.; Vo, T.P.; Nguyen, T.K. An Improved Shear Deformable Theory for Bending and Buckling Response of Thin-Walled FG Sandwich I-Beams Resting on the Elastic Foundation. Compos. Struct. 2020, 254, 112823. [Google Scholar] [CrossRef]
  6. Pagani, A.; Augello, R.; Carrera, E. Frequency and Mode Change in the Large Deflection and Post-Buckling of Compact and Thin-Walled Beams. J. Sound Vib. 2018, 432, 88–104. [Google Scholar] [CrossRef]
  7. Lanc, D.; Turkalj, G.; Vo, T.P.; Brnić, J. Nonlinear Buckling Behaviours of Thin-Walled Functionally Graded Open Section Beams. Compos. Struct. 2016, 152, 829–839. [Google Scholar] [CrossRef] [Green Version]
  8. Lanc, D.; Vo, T.P.; Turkalj, G.; Lee, J. Buckling Analysis of Thin-Walled Functionally Graded Sandwich Box Beams. Thin-Walled Struct. 2015, 86, 148–156. [Google Scholar] [CrossRef] [Green Version]
  9. Banat, D.; Mania, R.J.; Degenhardt, R. Stress State Failure Analysis of Thin-Walled GLARE Composite Members Subjected to Axial Loading in the Post-Buckling Range. Compos. Struct. 2022, 289, 115468. [Google Scholar] [CrossRef]
  10. Nguyen, T.-T.; Kim, N.-I.; Lee, J. Analysis of Thin-Walled Open-Section Beams with Functionally Graded Materials. Compos. Struct. 2016, 138, 75–83. [Google Scholar] [CrossRef]
  11. Nguyen, T.; Lee, J. Optimal Design of Thin-Walled Functionally Graded Beams for Buckling Problems. Compos. Struct. 2017, 179, 459–467. [Google Scholar] [CrossRef]
  12. Jin, M.X.; Xiang, D.Z. Thermally-Induced Bending-Torsion Coupling Vibration of Large Scale Space Structures. Comput. Mech. 2007, 40, 707–723. [Google Scholar] [CrossRef]
  13. Saha, S. Thermal Buckling and Postbuckling Characteristics of Extensional Slender Elastic Rods. J. Mech. Eng. 2009, 40, 1–8. [Google Scholar] [CrossRef] [Green Version]
  14. Cui, D.F.; Hu, H.Y. Thermal Buckling and Natural Vibration of the Beam with an Axial Stick—Slip—Stop Boundary. J. Sound Vib. 2014, 333, 2271–2282. [Google Scholar] [CrossRef]
  15. George, N.; Jeyaraj, P.; Murigendrappa, S.M. Thin-Walled Structures Buckling of Non-Uniformly Heated Isotropic Beam: Experimental and Theoretical Investigations. Thin Walled Struct. 2016, 108, 245–255. [Google Scholar] [CrossRef]
  16. Burgreen, D.; Manitt, P.J. Thermal Buckling of a Bimetallic Beams. J. Eng. Mech. Div. 1969, 95, 421–432. [Google Scholar] [CrossRef]
  17. Burgreen, D.; Regal, D. Higher Mode Buckling of Bimetallic Beam. J. Eng. Mech. Div. 1971, 97, 1045–1056. [Google Scholar] [CrossRef]
  18. Aydogdu, M. Thermal Buckling Analysis of Cross-Ply Laminated Composite Beams with General Boundary Conditions. Compos. Sci. Technol. 2007, 67, 1096–1104. [Google Scholar] [CrossRef]
  19. Trinh, L.C.; Vo, T.P.; Thai, H.T.; Nguyen, T.K. An Analytical Method for the Vibration and Buckling of Functionally Graded Beams under Mechanical and Thermal Loads. Compos. Part B Eng. 2016, 100, 152–163. [Google Scholar] [CrossRef]
  20. Kiani, Y.; Eslami, M.R. Thermal Buckling Analysis of Functionally Graded Material Beams. Int. J. Mech. Mater. Des. 2010, 6, 229–238. [Google Scholar] [CrossRef]
  21. Kiani, Y.; Eslami, M. Thermomechanical Buckling Dependent FGM Beams Of. Lat. Am. J. Solids Struct. 2013, 10, 223–246. [Google Scholar] [CrossRef]
  22. Giunta, G.; Crisafulli, D.; Belouettar, S.; Carrera, E. A Thermo-Mechanical Analysis of Functionally Graded Beams via Hierarchical Modelling. Compos. Struct. 2013, 95, 676–690. [Google Scholar] [CrossRef]
  23. Librescu, L.; Oh, S.Y.; Song, O. Thin-Walled Beams Made of Functionally Graded Materials and Operating in a High Temperature Environment: Vibration and Stability. J. Therm. Stress. 2005, 28, 649–712. [Google Scholar] [CrossRef]
  24. Ziane, N.; Meftah, S.A.; Ruta, G.; Tounsi, A. Thermal Effects on the Instabilities of Porous FGM Box Beams. Eng. Struct. 2017, 134, 150–158. [Google Scholar] [CrossRef]
  25. Kvaternik Simonetti, S.; Turkalj, G.; Lanc, D. Thin-Walled Structures Thermal Buckling Analysis of Thin-Walled Closed Section FG Beam-Type Structures. Thin-Walled Struct. 2022, 181, 110075. [Google Scholar] [CrossRef]
  26. Gjelsvik, A. The Theory of Thin-Walled Bars; Wiley: New York, NY, USA, 1981. [Google Scholar]
  27. Shen, H.S.; Wang, Z.X. Nonlinear Analysis of Shear Deformable FGM Beams Resting on Elastic Foundations in Thermal Environments. Int. J. Mech. Sci. 2014, 81, 195–206. [Google Scholar] [CrossRef]
  28. Reddy, J.N.; Chin, C.D. Thermomechanical Analysis of Functionally Graded Cylinders and Plates. J. Therm. Stress. 1998, 21, 593–626. [Google Scholar] [CrossRef]
  29. Turkalj, G.; Brnic, J.; Kravanja, S. A Beam Model for Large Displacement Analysis of Flexibly Connected Thin-Walled Beam-Type Structures. Thin-Walled Struct. 2011, 49, 1007–1016. [Google Scholar] [CrossRef]
  30. Bathe, K.J. Finite Element Procedure; Klaus-Jurgen Bathe, 2007. Available online: https://www.amazon.com/Finite-Element-Procedures-K-J-Bathe/dp/097900490X (accessed on 27 September 2022).
  31. Turkalj, G.; Brnic, J.; Prpic-Orsic, J. Large Rotation Analysis of Elastic Thin-Walled Beam-Type Structures Using ESA Approach. Comput. Struct. 2003, 81, 1851–1864. [Google Scholar] [CrossRef]
  32. Ansys Mechanical; Release 18.1; Help System. ANSYS, Inc.: Canonsburg, PA, USA. Available online: https://www.ansys.com/academic/terms-and-conditions (accessed on 27 September 2022).
Figure 1. Bimetallic beam wall.
Figure 1. Bimetallic beam wall.
Materials 15 07537 g001
Figure 2. Nodal force vectors and displacements.
Figure 2. Nodal force vectors and displacements.
Materials 15 07537 g002
Figure 3. Box beam cross-section.
Figure 3. Box beam cross-section.
Materials 15 07537 g003
Figure 4. Critical temperatures vs. displacement for a clamped–clamped beam L 1 = 6   m with λ = 0.2 .
Figure 4. Critical temperatures vs. displacement for a clamped–clamped beam L 1 = 6   m with λ = 0.2 .
Materials 15 07537 g004
Figure 5. Critical temperatures vs. displacement for a clamped–clamped beam L 1 = 6   m with λ = 0.5 .
Figure 5. Critical temperatures vs. displacement for a clamped–clamped beam L 1 = 6   m with λ = 0.5 .
Materials 15 07537 g005
Figure 6. Critical temperatures vs. displacement for a clamped–clamped beam L 1 = 6   m with λ = 0.8 .
Figure 6. Critical temperatures vs. displacement for a clamped–clamped beam L 1 = 6   m with λ = 0.8 .
Materials 15 07537 g006
Figure 7. L-frame.
Figure 7. L-frame.
Materials 15 07537 g007
Figure 8. Displacements of the L-frame for several values of the material thickness ratio.
Figure 8. Displacements of the L-frame for several values of the material thickness ratio.
Materials 15 07537 g008
Table 1. Temperature dependent coefficients [19].
Table 1. Temperature dependent coefficients [19].
MaterialPropertiesP0P1P1P2P3
Ti–6Al–4VE (Pa)122.56 × 1090.0−4.586 × 10−4 0.00.0
α (1/K)7.5788 × 10−6 0.06.638 × 10−4–3.147 × 10−60.0
SUS304E (Pa)201.04 × 1090.03.079 × 10−4 –6.534 × 10−7 0.0
α (1/K)12.330 × 10−6 0.08.086 × 10−4 0.00.0
Table 2. Critical buckling temperatures of the beam L 1 = 6   m for different boundary conditions and material thickness ratios.
Table 2. Critical buckling temperatures of the beam L 1 = 6   m for different boundary conditions and material thickness ratios.
λ
BCModeMethod00.20.40.50.60.81
C-CYPresent128.87135.26145.83153.64164.09199.19284.51
Shell130.22138.77151.18159.72170.71205.83287.97
XPresent239.24252.56273.89288.9308.62373.43528.16
Shell238.95254.69277.45293.11313.24377.6527.52
C-SYPresent66.5869.9175.3979.4484.84102.97146.98
Shell65.9870.31376.60480.93686.504104.29145.66
XPresent122.87130.02140.09148.58158.73192.03271.48
Shell120.46128.39139.88147.78157.93190.37265.92
S-SYPresent32.7234.36837.0739.0641.7250.6272.24
Shell32.19634.31437.3939.50742.22650.90471.077
XPresent60.2860.2969.0772.8677.8394.15133.08
Shell58.72562.668.20872.06377.01592.827129.64
Table 3. Critical buckling temperatures of the beam L 2 = 8   m for different boundary conditions and material thickness ratios.
Table 3. Critical buckling temperatures of the beam L 2 = 8   m for different boundary conditions and material thickness ratios.
λ
BCModeMethod00.20.40.50.60.81
C-CYPresent73.1776.8382.8587.393.24113.16161.54
Shell73.54278.36885.37590.296.405116.24162.35
XPresent135.52142.96154.88163.38174.53211.15298.53
Shell135.02143.91156.78165.62177213.37298.08
C-SYPresent37.6339.5242.6344.9147.9758.2183.07
Shell37.32339.77443.3345.78348.93358.99782.396
XPresent69.3573.3379.4583.8189.53108.31153.09
Shell68.25672.75479.26183.73689.487107.87150.68
S-SYPresent18.4519.3820.922.0323.5328.5540.73
Shell18.21119.40921.14922.34623.88528.79340.204
XPresent33.9535.938.941.0443.8453.0374.95
Shell33.27835.47338.65140.83643.64252.60373.466
Table 4. Critical buckling temperatures of the beam L 3 = 10   m for different boundary conditions and material thickness ratios.
Table 4. Critical buckling temperatures of the beam L 3 = 10   m for different boundary conditions and material thickness ratios.
λ
BCModeMethod00.20.40.50.60.81
C-CYPresent47.0349.3953.2756.1359.9572.75103.83
Shell47.15850.25254.74557.83961.81874.536104.11
XPresent86.7491.7299.37104.82111.98135.47191.49
Shell86.60592.308100.56106.23113.53136.86191.19
C-SYPresent24.1425.3527.3528.8130.7737.3453.28
Shell23.96525.53827.82329.39731.41937.88152.906
XPresent44.4346.9950.9153.7157.3869.498.09
Shell43.87246.76350.94553.82257.51869.33596.854
S-SYPresent11.8212.4213.3914.1115.0718.2926.09
Shell11.69512.46413.58214.3515.33818.49125.818
XPresent21.7422.9924.9126.2828.0833.9647.99
Shell21.39622.80724.8526.25528.05933.8247.234
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Simonetti, S.K.; Turkalj, G.; Banić, D.; Lanc, D. Bimetallic Thin-Walled Box Beam Thermal Buckling Response. Materials 2022, 15, 7537. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15217537

AMA Style

Simonetti SK, Turkalj G, Banić D, Lanc D. Bimetallic Thin-Walled Box Beam Thermal Buckling Response. Materials. 2022; 15(21):7537. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15217537

Chicago/Turabian Style

Simonetti, Sandra Kvaternik, Goran Turkalj, Damjan Banić, and Domagoj Lanc. 2022. "Bimetallic Thin-Walled Box Beam Thermal Buckling Response" Materials 15, no. 21: 7537. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15217537

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

Article Metrics

Back to TopTop