波纹柱壳在均布载荷作用下非线性动态解.docVIP

图 1 波纹柱壳 Fig.1 corrugated cylindrical shell 国内外学者对波纹壳结构进行了不少研究工作，主要为横向波纹壳结构[1-6]； 但对于纵向波纹壳结构的研究较少，研究的成果也很少，因此对纵向波纹壳的研 究更为迫切。以前的研究中，主要将波纹柱壳进行简化，采用简化模型计算，将 其简化为二维的杆结构，把内部的压力载荷视为径向的弹簧系统[7]，以及将其等 效为各向异性的壳结构—波纹杆和光滑的壳体和各向异性结构[8],不过以上简化 模型都很难准确的反映结构的几何非线性特性，应用的也不广泛。而近年来，只 有杨琳等[9-10]精确考虑波纹柱壳几何结构进行建模计算，以 Nayfeh 和 Pai[11]的理 论体系为依据，采用推广的 Hamilton 原理，建立了曲线坐标系下任意曲边柱壳 的非线性控制方程，并采用解析和数值积分的方法得到了该结构在均匀静态内压 作用下的线性解析解和非线性数值解。 对于波纹柱壳相关结构的动态特性的研究也很有限，Tesar[12]采用改进拉格 朗日运动方程、并行处理 FETM 和自适应网格加密技术对波纹薄壳的极限动力 学行为进行了线性和非线性数值分析。Samanta 和 Mukhopadhyay[13]将梯形波纹 板简化为考虑拉伸和弯曲刚度的等效各向正交模型，并对该等效模型进行了静动 态有限元分析。Semenyuk 等[14]研究得到了波纹壳的固有频率和模态，进一步考 察了波纹长度和幅度对其基频的影响。Liang 等[15]研究了受均布冲击压力作用下 的波纹板的动态响应。Liew 等[16]将加强型和非加强型正弦以及梯形波纹板简化 为等效正交板模型，采用无网格 Galerkin 方法进行了自由振动特性分析。本文 在杨琳等[9-10]得到的波纹壳受均布载荷下的非线性静态解的基础上，进一步研究 波纹柱壳结构受均布内压载荷作用下的动态响应，推导了波纹柱壳的非线性动态 控制方程，经过对控制方程的反复研究与分析，提出了一种基于静力变形模式的 动态响应计算方法，采用此方法得到了波纹柱壳的非线性动态解，并与 LS-DYNA 计算结果进行了对比验证，最后分析了加载速率和载荷形式对波纹柱壳动态变形 的影响。 2 波纹柱壳动态膨胀问题的计算方法 2.1 动态控制方程及边界条件 与波纹柱壳的静态变形计算[9-10]有所不同的是，波纹微段的力平衡方程和矩 平衡方程中应含有惯性项即  ?y?????????y  dy??  dy?? pd yi3????hwdyi c????hvdyib?? 0  (1)  ?M 2 ?y  dy?? Q2d y??? 0  (2)  因此，与静态变形推导类似[9-10]，可以得到波纹柱壳的非线性动态控制方程  ?hw????  ??????????????????????? ?y?????????????????????y  [(By?? v y )Q2 ]?? (By?? v y ) p?? 0  (3)  ?hv????  ??????????????????????? ?y????????????????????y  [(C y?? wy )Q2 ]?? (C y?? wy ) p?? 0  (4)  N2?? K1[ 1?? 2(Cy wy?? Byv y )?? w2?? v 2?? 1]  (5) M2?? D1[Byy wy?? Cyyvy?? (By?? vy )wyy?? (Cy?? wy )vyy ]  (6)  Q2??  ?M 2 ?y  (7)  相应的线性化增量控制方程及线性化增量边界条件为  y y 0 0  y y 0 0  ?N2?? K1(Cy?? wy )?wy?? K1(By?? vy )?vy ?M 2?? D1 (Byy?? vyy )?wy?? D1(C yy?? wyy )?v y ?D1(By?? vy )?wyy?? D1(C y?? wy )?vyy  (8) (9) (10) (11)  线性化增量边界条件为 ?(?M 2 ) ?y  ???Q2  (12)  C y??w?? By??v?? 0? ? ?  y?? 0 or yG  (13)  2.2 求解动力学方程