基于有限元的复合梁的建模方法.docxVIP

基于有限元的复合梁的建模方法 在现代工程中，管道屋顶被广泛应用，通常由包裹着合成材料的金属杆制成。为了研究复合材料梁的振动特性,通常采用刚体离散元法或有限段法将复合材料梁当作多刚体机构,建立刚体-弹簧振动模型,但精度较低,无法模拟材料的变形,而且该模型无法进行应力计算和疲劳强度计算。也可以采用求解析解和有限元的方法建立复合材料梁动力学模型,如Dancila和Armanios给出复合材料梁固有振动问题的精确解。Vanderpool提出了一种单斜晶体厚壁圆柱壳建模方法。Ghoneim等与申智春分别在Vanderpoo的基础上改进了圆截面空心复合梁的建模方法,并通过实验验证。向宇等采用一种齐次扩容精细积分技术建立了一种新的矩阵方法,研究了层合圆柱壳的振动特性和阻尼特性。Lin提出一种新型的一维两节点层叠梁单元。但对于大量采用这类复合构件的大型结构,不便求解解析公式,并且有限单元数过多,从而给建模和计算带来巨大的困难。 作者根据圆柱坐标系下管状梁的本构关系,采用复合梁单元建立管状复合材料梁的动力学方程,并以某型号变压器连线杆为例,建立复合材料连线杆模型。通过模态实验辨识模型阻尼系数,对振动模型进行修正,建立有阻尼模型。仿真结果表明该文建立的管状复合材料梁振动模型是有效和准确的。 1 材料梁的应力应变关系 考虑如图1(a)所示的管状复合材料梁,由n层不同种类的各向同性材料构成。第i层管状梁的弹性模量为Ei(i=1,…,n),密度为ρi,内径为Ri-1,外径为Ri(i=1,…,n),当基梁内径R0=0时,基杆为实心杆。 在主要考虑梁的横向和纵向振动的前提下,该管状复合材料梁服从Timoshenko梁的假设。建立如图1所示的坐标系,x轴与梁轴线重合,梁的任意横截面上有3个位移分量ux、uy和uz,以及绕x轴、y轴和z轴的转角α、β和?。 假定各层有相同的轴向位移ux和横向位移uy、uz。则σrr、σθθ和τθr可忽略不计,如图1(b)所示。由此得到空间柱坐标系下第i层应力应变关系为 (σixxτixθτirx)=(Ei000κGi000κGi)(εixxγixθγirx)(1) μi为第i层管状梁泊松比,Gi=Ei/(2+2μi),κ为剪切系数,对于圆形及环形截面κ=0.9。由几何关系可知,式(1)中应变分量为 εixx=?uix?x-?βi?xrcosθi+?αi?xrsinθiγixθ=(?uiz?x+αi)cosθi-(?uiy?x-βi)sinθi+r??i?xγirx=(?uiz?x+αi)sinθi-(?uiy?x-βi)cosθi(2) 2 系统动力学方程 2.1 《目标函数》记化 长度为L的单层管状梁单元的应变能可表示为 Ui=12∫L0∫2π0∫RiRi-1(σixxεixx+τixθγixθ+τirxγirx)rdrdθdx(i=1,?,n)(3) 记Ai1=2(1+μi),Ai3=Ei/(2+2μi),将式(1)和式(2)代入式(3),整理可得 Ui=12∫L0[bppi(?uix?x)2+bmmi((?αi?x)2+(?βi?x)2)+bvvi((?uiz?x+αi)2+(?uiy?x-βi)2)+btti(??i?x)2]dx(4) 式中,bppi=π(R2i-R2i-1)Ai3Ai1;bvvi=κπ(R2i-R2i-1)Ai3;bmmi=π(R4i-R4i-1)Ai3Ai1/4;btti=κπ(R4i-R4i-1)Ai3/2。 2.2 vix属性的测定 上述单层管状梁单元的动能可表示为 Τi=12∫L0∫2π0∫RiRi-1ρi(v2ix+v2iy+v2iz)rdrdθdx(i=1,??n)(5) 式中,vix=˙uix+r˙αisinθi-r˙βicosθi;viy=˙uiy-r˙?isinθi;viz=˙uiz-r˙?icosθi。代入式(5)可得 Τi=12∫L0[mi(˙u2ix+˙u2iy+˙u2iz)+Ιi2(˙α2i+˙β2i)+Ιi˙?2i]dx(6) 式中,mi=ρiπ(R2i-R2i-1);Ii=ρiπ(R4i-R4i-1)/2。 2.3 面临的刚度矩阵 当振动幅值很小时,可认为同一段复合梁单元内各层管状梁单元具有相同的节点位移,即uix=ujx,uiy=ujy,uiz=ujz,αi=αj,βi=βj,?i=?j(i≠j)。由此可视包含多层管状材料的复合材料梁单元在两端各有1个节点,每个节点有6个自由度,即(uxuyuzαβ?)Τ。则管状复合材料梁单元的自由度可表示为 δe=(ux1uy1uz1α1β1?1u2xu2yu2zα2β2?2)Τ(7) 各个变量采用线性形函数进行独立插值。将式(7)和形函数代入式(3),可得到以节点位移描述的管状复合材料梁单元的总应变能 U=n∑i=1Ui=n∑i