基于回传射线矩阵法的非均匀ig梁多轴响应.docxVIP

基于回传射线矩阵法的非均匀ig梁多轴响应 回传波束矩阵法是一种新型发展趋势的分析方法，用于解决空间和时间变量的线性偏差方程和复杂边界条件。在分析平面结构的瞬态响应时，ho等人提出了该方法，并在结构和层的瞬态响应分析中成功应用于结构和层的瞬态响应分析。tian等人讨论了梁缝的检测问题，并进一步介绍了该方法在实际工程中的应用价值。pao等人进一步改进了回传波束矩阵法，并与传统的传递矩阵法进行了详细比较，在估计静态和初始瞬态响应的情况下，决定了其准确性、数值稳定性和列式统一等优点。这是一种理想的弹性动力学分析。 回传射线矩阵法的主要思路是:采用连续单元模型和一套对偶局部坐标系统,应用积分变换手段将控制偏微分方程组转换为变换域内的常微分方程组;以波幅为基本未知量,并根据物理意义将弹性波分成离开波(出射波)和到达波(入射波),引入节点处力平衡条件和位移协调条件,建立散射关系,从数学上清晰地表示出节点处的出射波和入射波波幅之间的关系;联合根据对偶局部坐标系的特点建立的相位关系,求解波幅未知量;最后利用积分逆变换方法得到所关心的真实物理空间的解.迄今为止,回传射线矩阵法列式仅能考虑有限个集中力荷载作用情况,无法处理其他荷载形式,如分布荷载和移动荷载. 梁结构在分布荷载作用下的动力响应是实际工程中经常出现的问题,具有很强的现实意义.Timoshenko等人以Euler梁为对象,利用模态叠加法求解了均匀简支梁在不同分布荷载下的动力响应.Tong等人基于传递矩阵法对简谐分布荷载作用下的非均匀Timoshenko梁进行了分析.Doyle提出了频谱单元法,将分布荷载等效为节点荷载.本文从回传射线矩阵法的基本思路出发,以非均匀Timoshenko梁为例,对任意分布荷载作用下结构的动力响应进行分析.本文不仅在理论上对回传射线矩阵法进行了深化,而且拓展了其应用范围,使其能更好地为实际工程服务. 1 d4vd4+2密度7结构参数 考虑一根受竖向分布荷载q(x,t)作用的非均匀梁,长度为L,如图1(a)所示.梁的材料密度ρ、弹性模量E和横截面面积A可沿梁的长度方向任意变化.对于这样一根非均匀梁,直接求解往往非常困难,很多情况下难以求得解析解.为此,首先将梁沿长度方向划分为N段,当每段长度足够小时,可近似认为各梁段内的材料参数均为常数,由此就建立了分段均匀 梁模型,如图1(b)所示.视两端边界或相邻梁段的交界点为节点,每个梁段用其两端的节点编号表示,与该梁段相关的各个物理量采用节点编号作上标以示区别,如lIJ表示梁段IJ的长度. 在回传射线矩阵法中,对每个梁段IJ引入一对对偶局部坐标系(xIJ,yIJ)和(xJI,yJI),坐标原点分别位于两端节点处,如图1(c)所示.采用Timoshenko梁理论,考虑剪切和转动惯量的影响,局部坐标系下每个梁段(忽略上标IJ或JI)的运动方程为 EΙz?2??x2+κAG(?v?x-?)=ρΙz?2??t2,κAG(?2v?x2-???x)=ρA?2v?t2-q(x,t).}(1) 式中:E为弹性模量,G为剪切模量,ρ为密度,A为横截面面积,Iz为横截面惯性矩,κ为剪切系数,v为挠度,?为转角. 假设L0、A0、E0、ρ0为相应的长度、面积、模量和密度的参考量,令ξ=x/L0,τ=t/t0,ˉv=v/L0,ˉq(ξ,τ)=qL0/(E0A0),其中t0=L0/c1,c1=√E0/ρ0为纵波波速,则方程(1)可无量纲化为 ?2ˉ??ξ2+mn[?ˉv?ξ-ˉ?]=mEmρ?2ˉ??τ2,?2ˉv?ξ2-?ˉ??ξ=nmEmρ?2ˉv?τ2-nmAmEˉq.}(2) 式中:符号上一横表示无量纲量;n=E/(κG)为弹性模量与修正剪切模量之比;m=AL20/I是与横截面尺寸有关的无量纲系数;mΑ=A0/A,mΕ=E0/E,mρ=ρ0/ρ分别梁段面积、弹性模量、密度的比值. 假设梁受到简谐荷载ˉq(ξ,τ)=ˉq0(ξ,ˉω)eiˉωτ的作用,式中ˉω=ωt0为无量纲频率.不失一般性,考虑梁段材料参数相同(mΕ=mρ=1)情形.设ˉv=V(ξ)eiˉωτ,ˉ?=Φ(ξ)eiˉωτ,方程(2)可简化为 d2Φdξ2+mn(dVdξ-Φ)=-ˉω2Φ,d2Vdξ2-dΦdξ=-nˉω2V-nmAˉq0(ξ).}(3) 从式(3)中消去变量Φ可得 d4Vdξ4+ˉω2(n+1)d2Vdξ2+(nˉω2-m)ˉω2V=(m-nˉω2)mAˉq0-nmAd2ˉq0dξ2. 易得上式的齐次解为 V(ξ)=a2eik2ξ+d2e-ik2ξ+a3eik3ξ+d3e-ik3ξ. (4) 式中:a2、a3为待定入射波波幅;d2、d3为待定出射波波幅;k2、k3为波数,按下式确定: k2,3=√[(n+1)ˉω2±√(n+1)2ˉω4-4ˉω2