基于横向poisson效应的梯度功能梁的波动方程.docxVIP

基于横向poisson效应的梯度功能梁的波动方程 随着科学研究的发展，新材料也层出不穷。作为一种新型的功能性材料，梯度功能材料在航天工业、能源工业、电子工业、化学工程和生物医学工程等重要领域得到了广泛应用和开发。梯度功能材料中波的传播性质的研究已备受国内外学者的关注。目前对梯度功能材料中波的传播进行的研究,大多数是采用分层理论和有限元方法,对功能梯度材料中声表面波的传播进行研究,并取得了一定的研究成果。由于数学上的复杂性,对梯度功能材料中波的传播问题的解析研究,只有很少的研究成果见诸报道。Han X对弹性模量E(x)、密度ρ(x)以幂函数形式变化的梯度功能材料中波的传播问题进行了研究,得出了贝赛尔函数解;Z. Wesolowski研究了在两种均匀介质间存在密度连续变化介质的情况下的波的传播问题;张立刚等研究了梯度功能材料中SH波的传播问题,利用WKBJ理论求出了梯度功能材料位移的近似解,还对均匀覆层梯度功能半空间中的Love波频散问题进行了研究;刘睫和王子昆研究了梯度材料层状结构中的Love波和梯度介质半空间Rayleigh面波的传播性能。这些研究基本上都是对功能梯度材料中表面波的传播进行研究,而对于梯度功能材料中非线性波传播的研究则很少涉及。 本文在考虑有限变形和横向Poisson效应的情况下,借助Hamilton 变分原理建立了弹性模量和质量密度沿厚度连续变化时的一维非线性波动方程,根据非线性波动方程具有孤波解的性质,从基本的sech(x)函数出发构造出满足方程的形函数解,并利用待定法解得待定系数从而求得了非线性波动方程的位移孤波解,并讨论了当材料属性沿厚度方向指数形式变化和抛物线形式变化时,不同梯度功能材料参数和波传播时的波速对孤波变化的影响。 1 功能梯度材料参数 无限长均质矩形等厚度的梯度功能梁厚度为h,宽度为b,波沿x轴正向传播,如图1所示。材料的泊松比μ为常量,弹性模量为E,质量密度为ρ,且E和ρ均为厚度z的函数,即: E=E(z),ρ=ρ(z) (1) 当功能梯度材料的参数沿厚度呈指数变化时,有: E(z)=E0eα(z+h2),ρ(z)=ρ0eβ(z+h2)E(z)=E0eα(z+h2),ρ(z)=ρ0eβ(z+h2), (2) 当功能梯度材料的参数沿厚度呈抛物线变化时,有: E(z)=E0[1+α1(z+h2)]2ρ(z)=ρ0[1+β1(z+h2)]2(3)E(z)=E0[1+α1(z+h2)]2ρ(z)=ρ0[1+β1(z+h2)]2(3) 其中,E0和ρ0为常数,即为梁底面(z=-h/2)处的材料常数;α,β,α1,β10,称为功能梯度材料参数。 设平面波沿x轴的正方向传播,梁只受单轴应力作用,不发生弯曲和剪切变形,即梁中σx≠0,σy=σz=0;变形前为平面的横截面在变形后始终保持平面;设梁中沿x,y,z方向的位移分别为u,v,w,当考虑有限变形的影响且计入横向Poisson效应时,梁中沿x,y,z方向的应变分别为:εx=?u?x+12(?u?x)2,εy=εz=-μ?u?xεx=?u?x+12(?u?x)2,εy=εz=?μ?u?x;则侧向位移分量可简写为v=-μy?u?x?w=-μz?u?x。 基于以上假定,梁单位长度上的动能为: 其中,梯度功能材料参数沿厚度呈指数变化时: ρ1=∫∫Sρdydz=∫h2-h2∫b2-b2ρ0eβ(z+h2)dzdy=ρ0bβ(eβh-1)(5) κ=μ2∫∫Sρ(y2+z2)dzdy=ρ0bμ21βeβh[(h2-1β)2+1β2+b212]-ρ0bμ21β[(h2+1β)2+1β2+b212](6) 梯度功能材料参数沿厚度呈抛物线变化时: ρ1=∫h2-h2∫b2-b2ρ0[1+β1(z+h2)]2dzdy=ρ0bh(1+β1h+β21h23)(7) κ=μ2∫h2-h2∫b2-b2ρ0[1+β1(z+h2)]2(y2+z2)dzdy=ρ0μ2bh[112(b2+h2)+β1h12(b2+h2)+β21h2(b236+h230)](8) 上述式子中,S为梁的横截面积。 假设材料服从线性的本构关系,梁单位长度上的弹性势能为: 其中,梯度功能材料参数沿厚度呈指数变化时: E1=∫∫SEdzdy=∫h2-h2∫b2-b2E0eα(z+h2)dzdy=E0bα(eαh-1)(10) 梯度功能材料参数沿厚度呈抛物线变化时: E1=∫h2-h2∫b2-b2E0[1+α1(z+h2)]2dzdy=E0bh(1+α1h+α21h23)(11) 利用Hamilton变分原理 其中: 根据Euler方程,若L(u,ut,ux,uu,uxx,uxt,… )=0,则: ?L?u-??t?L?ut-??x?L?ux+?2?t2?L?uu+?2?x2?L?uxx+