Dynamic Buckling of Functionally Graded Timoshenko Beam under Axial Load
-
摘要: 假设功能梯度材料Timoshenko梁各项物性参数只沿厚度方向按幂函数进行连续变化,研究了功能梯度材料Timoshenko梁的动力屈曲。基于一阶剪切理论,采用Hamilton原理推导出轴向载荷作用下,功能梯度材料Timoshenko梁动力屈曲的控制方程。利用里兹法与棣莫弗公式相结合,获得了功能梯度材料Timoshenko梁在夹支-固支边界条件下动力屈曲临界载荷的解析表达式和屈曲解。应用MATLAB编程计算,讨论了功能梯度材料Timoshenko梁的几何尺寸、梯度指数、模态数、材料构成、泊松比以及弹性模量对临界载荷的影响。结果表明:功能梯度材料Timoshenko梁动力屈曲临界载荷随梁长度的增大而减小,随着梯度指数的增大而减小,随模态数的增大而增大,说明冲击载荷越大,高阶模态越容易被激发;随着泊松比和弹性模量的增大而增大,且泊松比的影响较小,而弹性模量的影响较大。由于剪切项的影响,临界载荷-临界长度的关系曲线在加载端变化趋势平缓。随着模态数的增大,梁的屈曲模态越为复杂。Abstract: In this study, we investigated the dynamic buckling of the functionally graded Timoshenko beam whose property parameters continuously change according to the power function along the thickness direction.Based on the first order shear deformation theory, we derived the governing equation of the dynamic buckling of functionally graded material Timoshenko beams under axial step loading by using the Hamilton's principle.Using the Ritz method combining with the de Moivre's formula, we obtained the buckling solution and the expression of the critical load of the dynamic buckling of functionally graded material Timoshenko beam under the clamped-fixed boundary condition.Then, the influence of geometric size, gradient index, modal number, material composition, Poisson's ratio and elastic modulus on the critical load by MATLAB calculation was discussed.The results show that the critical load of the functionally graded material Timoshenko beam decreases with the increase of beam length and the gradient index, and increases with the increase of the modal number, showing that the higher modal number is more easily excited by the increase of impact load.Furthermore, the critical load increases with the increase of the Poisson's ratio and the elastic modulus, and the effect of elastic modulus is greater than Poisson's ratio.The critical load-critical length curve tends to be gentle at the loading end because of the influence of shear term.Buckling mode of beam becomes more complicated when the modal number increases.
-
Key words:
- functionally graded material /
- dynamic buckling /
- critical load /
- Ritz method /
- de Moivre's formula
-
功能梯度材料是一种集各组分于一体的新型材料,其微观结构和综合性能都呈连续性变化,从而降低了材料的应力集中,提高了材料的性能与功能。因此,功能梯度材料在航天航空、船舶工业、汽车工业、能源工业以及生物医学工业等工程领域有广泛的应用前景,国内外众多学者对功能梯度材料进行了深入的分析和探究[1-3]。Kocaturk和Akbaş[4-5]考虑几何非线性,使用有限元法对不同边界条件下功能梯度材料Timoshenko梁在热载作用下的后屈曲进行了分析。Paul[6]对面内热负荷作用下功能梯度材料Timoshenko梁的后屈曲进行了研究。Şimşek[7]基于Timoshenko梁理论,首次对功能梯度材料在双向递变情况下的梁动力屈曲进行了讨论。Eltaher等[8]基于Timoshenko梁理论和最小势能原理,运用有限元法对功能梯度材料Timoshenko梁的动力屈曲进行了分析。Rahimi等[9]基于Timoshenko梁理论,采用Von-Karman非线性应变-位移关系,对功能梯度材料梁的后屈曲进行了研究。Li[10]比较了轴向载荷作用下功能梯度材料Timoshenko梁和欧拉-伯努利梁的临界屈曲载荷之间的关系。Kahya等[11]采用有限元软件对功能梯度梁的动力屈曲以及自由振动进行了分析。Kiani[12]基于欧拉-伯努利梁理论,对不同热载荷作用下功能梯度材料梁的动力屈曲进行了研究。Anandrao[13]运用里兹法和有限元法对均匀细长功能梯度梁热后屈曲进行了讨论。Akbaş等[14]利用Newton-Raphson迭代法结合有限元法研究了热载荷作用下功能梯度梁的后屈曲行为。Rychlewska[15]通过等效梁近似功能梯度梁,对功能梯度梁进行了动力屈曲分析。综上所述,功能梯度梁的屈曲研究主要以数值计算为主。
本研究基于一阶剪切理论,利用Hamilton原理推导出功能梯度材料Timoshenko梁动力屈曲的控制方程;采用里兹法结合棣莫弗公式,导出功能梯度材料Timoshenko梁在夹支-固支边界条件下动力屈曲临界载荷的解析表达式;应用MATLAB软件编程计算,讨论功能梯度材料Timoshenko梁的几何尺寸、梯度指数、模态数、材料构成、泊松比以及弹性模量对临界载荷的影响。
1. 功能梯度材料物性参数
功能梯度材料Timoshenko梁的各项物性参数只沿厚度(z)方向按幂函数连续变化,表示为[16]
弹性模量 Ez=(Ec−Em)(12+zh)k+Em (1) 密度 ρz=(ρc−ρm)(12+zh)k+ρm (2) 泊松比 μz=(μc−μm)(12+zh)k+μm (3) 剪切模量 Gz=Ez2(1−μz) (4) 式中:Ec是陶瓷的弹性模量,Em是金属的弹性模量,ρc是陶瓷的密度,ρm是金属的密度,μc是陶瓷的泊松比,μm是金属的泊松比, k是梯度指数。
2. 动力控制方程的推导
如图 1所示,功能梯度材料Timoshenko梁的长为l,宽为b,高为h,材料各项性质只沿高度方向呈连续性变化,剖面中心线为x轴,设其轴向位移为u(x, t),横向位移为w(x, t),横截面转角为φ。基于平截面假定,由一阶剪切理论,Timoshenko梁的位移及应变[17]可表示为
{u=u0−zφw=w0 (5) {εx=ε0x+zkxγxz=γ0xz (6) 其中 ε0x=∂u0∂x,γ0xz=∂w0∂x−φ,kx=−∂φ∂x (7) 式中:u0是杆的中性面在x轴方向的位移,w0是杆的中性面在z轴方向的位移,εx0是杆的中性面在x方向的正应变,γxz0是杆的中性面在xz方向的切应变,kx是杆件任意点在x轴方向上的曲率。
其内力、剪力及内力矩可表示为
内力 Nx=A11ε0x+B11kx (8) 剪力 Qx=KA55(∂w0∂x−φ) (9) 内力矩 Mx=B11ε0x+D11kx (10) 式中:A11为拉伸刚度系数; B11为耦合刚度系数; D11为扭转刚度系数; A55为剪切刚度系数; K为剪切修正系数,取K=5/6。
A11=∫h/2−h/2Ez1−μ2zdz (11) B11=∫h/2−h/2zEz1−μ2zdz (12) D11=∫h/2−h/2z2Ez1−μ2zdz (13) A55=∫h/2−h/2Ez2(1+μz)dz (14) 应变能 U=12∫l0∫b0∫h/2−h/2(Nxε0x+Mxkx+Qxγxz)dxdydz (15) 动能 T=12∫l0∫b0∫h/2−h/2ρz[(∂u∂t)2+(∂w∂t)2]dxdydz (16) 外力功 W=12∫l0∫b0∫h/2−h/2Nx(∂w∂x)2dxdydz (17) 定义广义变量为
(I0,I1,I2)=∫h/2−h/2ρz(1,z,z2)dz (18) 将(15)式~(18)式代入Hamilton原理,即δ∫t2t1(T−U+W)dt=0,由u0、v0、φ的变分系数取零可得
δu0:−I0∂2u0∂t2+I1∂2φ∂t2+A11∂2u0∂x2−B11∂2φ∂x2=0 (19) δw0:−I0∂2w0∂t2−N(t)∂2w0∂x2+KA55∂2w0∂x2−KA55∂φ∂x=0 (20) δφ:I1∂2u0∂t2+I2∂2φ∂t2−B11∂2u0∂x2+D11∂2φ∂x2−KA55φ+KA55∂w0∂x=0 (21) 不考虑轴向惯性[18],将(21)式对x求导可得
−I2∂3φ∂t2∂x+D11∂3φ∂x3−KA55∂φ∂x+KA55∂2w0∂x2=0 (22) 对(20)式进行偏导计算,得
{KA55∂φ∂x=−I0∂2w0∂t2−N(t)∂2w0∂x2+KA55∂2w0∂x2KA55∂3φ∂x∂t2=−I0∂4w0∂t4−N(t)∂4w0∂x2∂t2+KA55∂4w0∂x2∂t2KA55∂2φ∂x2=−I0∂3w0∂t2∂x−N(t)∂3w0∂x3+KA55∂3w0∂x3KA55∂3φ∂x3=−I0∂4w0∂t2∂x2−N(t)∂4w0∂x4+KA55∂4w0∂x4 (23) 将(23)式代入(22)式,得到控制方程为
D11KA55[−I0∂4w∂x2∂t2−N(t)∂4w∂x4]+D11∂4w∂x4+I2KA55[I0∂4w∂t4+N(t)∂4w∂x2∂t2]−I2∂4w∂x2∂t2+I0∂2w∂t2+N(t)∂2w∂x2=0 (24) 3. 动力屈曲临界载荷
3.1 用里兹法化简控制方程
设一端夹支、一端固支梁的动力屈曲解的形式为[19]
w=T(t)(Asinn1πxl+Bsinn2πxl) (25) 式中:n1、n2为屈曲的模态阶数,A、B为待求系数。动力屈曲解要满足下列一端夹支、一端固支梁的边界条件
w(0)=∂w∂x|x=0=w(l)=∂w∂x|x=l=0 (26) 将(26)式代入(25)式可得
{An1πl+Bn2πl=0An1πlcos(n1π)+Bn2πlcos(n2π)=0 (27) 若(27)式有解则cos(n1π)与cos(n2π)奇偶性相同,即n2=n1+2。又An1πl+Bn2πl=0,即AB=−n2n1。取A=1,B=−n1n2。则梁的动力屈曲解的形式化为
w=T(t)(sinn1πxl−n1n2sinn2πxl) (28) 式中:n1、n2是屈曲的模态阶数,且n1=1, 2, 3,…; n2=n1+2。
对(28)式求偏导可得
{wxx=Tn1π2l2(−n1sinn1πxl+n2sinn2πxl)wtt=Ttt(sinn1πxl−n1n2sinn2πxl)wxxtt=Tttn1π2l2(−n1sinn1πxl+n2sinn2πxl)wxxxx=Tn1π4l4(n31sinn1πxl−n32sinn2πxl)wtttt=Ttttt(sinn1πxl−n1n2sinn2πxl) (29) 3.2 由棣莫弗公式求解临界载荷
由棣莫弗公式(cosθ+isinθ)n=cos(nθ)+isin(nθ)可得
sinnπxlsinπxl=C1n(cosπxl)n−1−C3n(cosπxl)n−3(sinπxl)2+⋯ (30) 将(29)式除以sinπxl,当梁发生屈曲,即x=l时有
{wxxsinπxl=Tn1π2l2(−n21+n22)wttsinπxl=Ttt(n1−n1)=0wxxttsinπxl=Tttn1π2l2(−n21+n22)wxxxxsinπxl=Tn1π4l4(n41−n42)wttttsinπxl=Ttttt(n1−n1)=0 (31) 将(31)式代入控制方程得
Ttt+T(α−β)=0 (32) 其中
α=N(t)D11(n21+n22)π2+KA55l2[−D11I0+I2N(t)−I2KA55],β=D11(n21+n22)π2KA55[−D11I0+I2N(t)−I2KA55]l2 当α-β<0时, 方程的解为指数形式的发散解,表明系统经受微小的扰动后,运动为无界,因此系统呈现不稳定的状态。当α-β>0时,表明系统围绕平衡位置做微小振动,因此系统呈现稳定平衡的状态[19]。当α-β=0时为Timoshenko梁屈曲临界状态,由此推出动力屈曲临界载荷的表达式为
Ncr=D11(n21+n22)π2KA55D11(n21+n22)π2+KA55l2cr (33) (33)式表明:剪切项对动力屈曲临界载荷有一定的影响,当临界长度lcr=0时,临界载荷为KA55。
4. 数值计算
运用MATLAB软件对(33)式进行编程计算,材料选取陶瓷/钛功能梯度材料、陶瓷/铁功能梯度材料和陶瓷/铜功能梯度材料。梁的宽度b=2cm,厚度h=3cm。讨论梁厚度、模态数(n)、梯度指数(k)、泊松比(μ)、弹性模量(E)以及材料构成对临界载荷的影响。材料各参数如表 1所示。
表 1 材料各项参数Table 1. Material parametersMaterial Elasticity modulus
E/GPaDensity
ρ/ (g·cm-3)Poisson's ratio
μCeramic 385 3.96 0.23 Titanium 108.5 4.54 0.41 Iron 155 7.86 0.291 Copper 119 8.96 0.326 图 2所示为n=1、k=1/2时,不同材料构成下临界载荷Ncr与临界长度lcr的关系曲线。当h、n、k、μ以及E一定,材料构成不同时,其Ncr随着lcr的增加而减小。当lcr相同时,陶瓷/铁材料的临界载荷值高于陶瓷-钛材料,且陶瓷/钛材料梁与陶瓷/铜材料梁的Ncr-lcr曲线在lcr超过1m时接近重合。由于剪切项的影响,曲线在加载端变化趋势平缓。
图 3所示为陶瓷/钛功能梯度材料梁在n=1时,不同k下Ncr与lcr的关系曲线。当梁的h、n、μ、E以及材料构成一定,k不同时,其Ncr随着lcr的增大而减小。当lcr相同时,Ncr随着的k的增大而减小。由于剪切项的影响,曲线在加载端变化趋势也平缓,如图 2所示。
图 4所示为陶瓷/钛功能梯度材料梁在k=1/2时,不同n下Ncr与lcr的关系曲线。当梁的h、k、μ、E以及材料构成一定,n不同时,其Ncr随着lcr的增大而减小。当lcr相同时,Ncr随着n的增大而增大。这表明冲击载荷越大,高阶模态越容易被激发。
图 5表示陶瓷/钛功能梯度材料梁在k=1/2、n=1时,不同h下Ncr与lcr的关系曲线。当梁的n、k、μ、E以及材料构成一定,h不同时,其Ncr随着lcr的增大而减小。当lcr相同时,Ncr随着h的增大而增大。
图 6表示k=1/2、n=1、μ=0.41时,不同E下Ncr与lcr的关系曲线。当梁的h、k、μ、n一定,E不同时,其Ncr随着lcr的增大而减小。当lcr相同时,Ncr随着E的增大而增大。
图 7表示k=1/3、n=1、E=100GPa时,不同μ下Ncr与lcr的关系曲线。当梁的h、k、E、n一定,μ不同时,其Ncr随着lcr的增大而减小。当lcr相同时,不同μ下曲线接近重合,表明μ对Ncr的影响不大。
图 8~图 11是陶瓷/钛功能梯度材料梁在k=1/2时,不同模态数下的屈曲模态图。可见,随着模态数n的增大,梁的屈曲模态越发复杂。固定端屈曲模态的幅值较大,对比不同模态下的固支端与夹支端,其幅值比例范围为1.62~2.88。
5. 结论
(1) 基于一阶剪切理论,通过Hamilton原理推导出功能梯度材料Timoshenko梁的动力屈曲控制方程。
(2) 用里兹法与棣莫弗公式相结合,得到了功能梯度材料Timoshenko梁在夹支-固支边界条件下动力屈曲临界载荷的解析表达式及屈曲解。
(3) 应用MATLAB软件,讨论了功能梯度材料Timoshenko梁的几何尺寸、梯度指数、模态数、材料构成、泊松比以及弹性模量对临界载荷的影响。结果表明:功能梯度材料Timoshenko梁的动力屈曲临界载荷随着梁长度、梯度指数的增大而减小;临界载荷随模态数的增大而增大,且冲击载荷越大,高阶模态越容易被激发;临界载荷随泊松比及弹性模量的增大而增大,且泊松比的影响不大,而弹性模量的影响较大。由于剪切项的影响,临界载荷-临界长度的关系曲线在加载端变化趋势平缓;随着模态数的增大,梁的屈曲模态越发复杂。
-
表 1 材料各项参数
Table 1. Material parameters
Material Elasticity modulus
E/GPaDensity
ρ/ (g·cm-3)Poisson's ratio
μCeramic 385 3.96 0.23 Titanium 108.5 4.54 0.41 Iron 155 7.86 0.291 Copper 119 8.96 0.326 -
[1] BIRMAN V, BYRD L W.Modeling and analysis of functionally graded materials and structures[J]. Applied Mechanics Reviews, 2007, 60(5):195-216. doi: 10.1115/1.2777164 [2] GUPTA A, TALHA M.Recent development in modeling and analysis of functionally graded materials and structures[J]. Progress in Aerospace Sciences, 2015, 79:1-14. doi: 10.1016/j.paerosci.2015.07.001 [3] NAEBE M, SHIRVANIMOGHADDAM K.Functionally graded materials:a review of fabrication and properties[J]. Applied Materials Today, 2016, 5:223-245. doi: 10.1016/j.apmt.2016.10.001 [4] KOCATURK T, AKBAŞ Ş D.Post-buckling analysis of Timoshenko beams made of functionally graded material under thermal loading[J]. Structural Engineering & Mechanics, 2012, 41(6):775-789. doi: 10.1007/s10999-010-9132-4 [5] KOCATURK T, AKBAŞ Ş D.Post-buckling analysis of Timoshenko beams with various boundary conditions under non-uniform thermal loading[J]. Structural Engineering & Mechanics, 2011, 40(3):347-371. http://www.koreascience.or.kr/article/ArticleFullRecord.jsp?cn=KJKHB9_2011_v40n3_347 [6] PAUL A.Non-linear thermal post-buckling analysis of FGM Timoshenko beam under non-uniform temperature rise across thickness[J]. Engineering Science & Technology:An International Journal, 2016, 19(3):1608-1625. https://www.sciencedirect.com/science/article/pii/S2215098616302853 [7] ŞIMŞEK M.Buckling of Timoshenko beams composed of two-dimensional functionally graded material (2D-FGM) having different boundary conditions[J]. Composite Structures, 2016, 149:304-314. doi: 10.1016/j.compstruct.2016.04.034 [8] ELTAHER M A, KHAIRY A, SADOUN A M, et al.Static and buckling analysis of functionally graded Timoshenko nanobeams[J]. Applied Mathematics &Computation, 2014, 229(229):283-295. https://www.sciencedirect.com/science/article/pii/S0096300313013556 [9] RAHIMI G H, GAZOR M S, HEMMATNEZHAD M, et al.On the postbuckling and free vibrations of FG Timoshenko beams[J]. Composite Structures, 2013, 95(1):247-253. https://www.sciencedirect.com/science/article/pii/S0263822312003625 [10] LI S R.Relations between buckling loads of functionally graded Timoshenko and homogeneous Euler-Bernoulli beams[J]. Composite Structures, 2013, 95(1):5-9. [11] KAHYA V, TURAN M.Finite element model for vibration and buckling of functionally graded beams based on the first-order shear deformation theory[J]. Composites Part B:Engineering, 2017, 109:108-115. doi: 10.1016/j.compositesb.2016.10.039 [12] KIANI Y.Thermal buckling analysis of functionally graded material beams[J]. International Journal of Mechanics & Materials in Design, 2010, 6(3):229-238. doi: 10.1007/s10999-010-9132-4 [13] ANANDRAO K S.Thermal post-buckling analysis of uniform slender functionally graded material beams[J]. Structural Engineering & Mechanics, 2010, 36(5):545-560. https://www.researchgate.net/publication/264077187_Thermal_post-buckling_analysis_of_uniform_slender_functionally_graded_material_beams [14] AKBAŞ Ş D, KOCATÜRK T.Post-buckling analysis of functionally graded three-dimensional beams under the influence of temperature[J]. Journal of Thermal Stresses, 2013, 36(12):1233-1254. doi: 10.1080/01495739.2013.788397 [15] RYCHLEWSKA J.A new approach for buckling analysis of axially functionally graded beams[J]. Journal of Applied Mathematics & Computational Mechanics, 2015, 14(2):95-102. [16] FEREZQI H Z, TAHANI M, TOUSSI H E.Analytical approach to free vibrations of cracked Timoshenko beams made of functionally graded materials[J]. Mechanics of Advanced Materials & Structures, 2010, 17(5):353-365. doi: 10.1080/15376494.2010.488608?scroll=top&needAccess=true [17] 徐华, 李世荣.一阶剪切理论下功能梯度梁与均匀梁静态解之间的相似关系[J].工程力学, 2012, 29(4):161-167. http://www.cnki.com.cn/Article/CJFDTotal-GTLX201403013.htmXU H, LI S R.Analogous relationship between the static solutions of functionally graded beams and homogeneous beams based on the first-order shear deformation theory[J]. Engineering Mechanics, 2012, 29(4):161-167. http://www.cnki.com.cn/Article/CJFDTotal-GTLX201403013.htm [18] KIRCHHOFF G.Ueber das Gleichgewicht und die Bewegung eines unendlich dünnen elastischen Stabes[J]. Journal FüR Die Reine Und Angewandte Mathematik, 2009, 1859(56):285-313. https://eudml.org/doc/147766 [19] 李楠, 韩志军, 路国运.基于里兹法研究复合材料层合板的动力屈曲问题[J].振动与冲击, 2016, 35(10):180-184. http://industry.wanfangdata.com.cn/jt/Detail/Periodical?id=Periodical_zdycj201610029LI N, HAN Z J, LU G Y.Research on dynamic buckling of laminated composite plates using Ritz method[J]. Journal of Vibration and Shock, 2016, 35(10):180-184. http://industry.wanfangdata.com.cn/jt/Detail/Periodical?id=Periodical_zdycj201610029 -