几何非线性弹性梁坠撞动力学响应研究.pdfVIP

l993年第 振 动 与 冲 击 总第 48期 ———— 一 撩磊 蜀 锄 几何非线性弹性梁坠撞动力学响应研究 2／，(cf) ，一 ／ 逯兆乾 朱世晋 (南京航空学院) T U 3ff ， — ／ 摘 要 本文对计凡几侗非线性 的弹性梁坠落撞地响应问题提 出一种简便实际的算法 ，并作了算倒计算 作 者还用应变片技术对算例的钢粱做了幕列实验测量。文中对计算和实验结果进行丁分析讨论 。 一 引 言 结构坠落碰撞是工程实际中常会遇到的一类问题，譬如直升机失事或粗暴着陆使飞机结构受 到冲击等。因为碰撞物理过程短暂、作用力大而且变化剧烈复杂 斓E由于结构碰撞时大变形而必须 考虑其几何非线性，有时还有材料非线性 问题；此外物体形状类型及碰撞情况不同都可引起碰撞响 应有很大的不同 这些都使坠撞响应问题变得很复杂 ，除一些经典问题外 ，至今尚缺少系统探入的 研究 。 本文作为一种基础性研究，对几何非线性弹性梁的坠落撞地响应 问题进行了理论分析计算和 实验研究，以期能对这类结构物坠撞响应的基本特性有所认识。 二 理论分析 1．有限元方程推导 对于几何非线性梁的动响应问题 ，为了使推导过程简洁、直观 ，本文采用 自然刚度法的有限元 方法 一。该法推导平面粱问题时，首先将平面粱单元二节点处共3×2=6个物理位移转化为与梁单 元的三个刚体模态及三个 自然变形模态相对应的共六个广义位移 其次，由于计入梁有大变形的几何非线性，采用小步增量法 ．即考虑单元节点力矢量增量和位 移矢量增量间关系．可导出为 {dp)= [K$-j{出 ) (1) 其中单元的切线刚度矩阵 -KT：= [K ]1-CK ] (2) EKE]是单元线弹性刚度矩阵和寻常的相同；而[Ko]是单元几何刚度矩阵，用 自然刚度法可写出为 EKo：==EKo：]+ Ore] (3) 其 中[ 为 由于正应力 引起的几何刚度矩阵，可导 出形如Pn ·EB]l[K ]为由于切应力 r引 起的几何刚度矩阵，可导出形如s ·ED] 以上P 是梁单元截面拉力，是粱单元截面切力 ，EB]、 [D]是依赖于单元相对于总坐标系倾角 的矩阵，其具体推导见文[2]。周P 、S、a等都随单元变值 而改变，所以[ 、[K ]为单元刚度矩阵的非线性部分。 ． 振 动 与 冲 击 由粱单元刚度矩阵、质量矩阵、阻尼矩阵可以组装成整个梁的刚度、质量、阻尼矩阵。从而可 以 得到整个梁的增量形式的动力学微分方程式 [ ](At3)+ [c](△ + [K —Kc]{AU)= {Af(t)) (4) 最终的方程可用某种逐步积分法．本文用 Newmank法求解其动态响应 问题 。因为系统是非线性 的，所以积分的每一步中还须进行若干次迭代 ，达到一定精度后再做下一步。 2．碰撞的处理 解决结构碰撞 问题常用方法是 ：将结构作有限元离散，用若干集中质量 ，再用若干线性或非线 性弹簧来代替碰撞处结构的作用 用这类模型求解整结构坠落