球面网壳失稳运动特性研究.pdfVIP

球面网壳失稳的运动特性研究 夏开全姚卫星 董石麟 (南京航空航天大学．南京21C016)(浙江走学，杭州310027) 提要 本文从运动学的角度．基十非线性运动平衡方程，对球面网壳结构发生跳跃失稳时的运动特 性进行了研究。在结构达到临界失稳状态的前后分别采用基于静力的几何非线性分析方法和基于动力响应 的几l可非线性分析方法，描述结构的实际行为特性。对厨壳结构先稳时盼加速度、速度与位移等结构响应 的研究，揭示了失稳的运动本质，为结构稳定性分析、失稳机制和结构失稳的传播特性的研究提供了另一 条途径。 关键词球面罔壳非线性稳定性运动特性 一、引 言 单层球面网壳结构几何上为正高斯曲率网壳，具有刚性而不可展开的几何形状，是一 种几何软化结构，其节点失稳将引起局部曲面的反转，伴以强烈的运动特性。因此，由网壳 结构节点引起的局部失稳多为跳跃失稳，具有很大的危险性。这种失稳往往由于理想壳的几 何初始缺陷或杆件的失稳所引起。由于跳跃失稳实际上是一个运动过程，局部失稳可能会逐 渐扩展传播而形成大范围的捷体失稳{L2】。 失稳传播的研究20世纪，O年代始于并主要限于对海底管道屈曲传播的研究nq。失稳 的传播是19 10年美国ButtlcColumnbus实验室在模型实验中首次发现的。而对网壳结构失 稳传播的研究却起步较晚。Oioncu和Lenza等学者对铰接单层柱面网壳的失稳传播有一定 研究口一。失稳运动特性的研究能揭示失稳的运动本质．对正确认识失稳的机制以及失稳的 传播有着重要的意义。 以往的失稳分析，往往是基于静力的失稳分析。当结构达到临界失稳状态时，其后临 界状态的跟踪分析采取卸载方式。这样所得到的后临界分析是一种人为的假设，与失稳的实 际情况不符。事实上，结构处于临界状态时．微小的外力扰动将导致结构的不稳定平衡，经 历一个运动过程而达到新的稳定平衡状态。本文从临界点处开始采用动力响应分析方法用描 述结构的不稳定平衡过程，以揭示结构失稳的本质。首先根据D’Alcmbcrt原理建立结构的 几何非线性增量有限元运动平衡方程，然后采用Newmark-B隐式积分法嗍求解该方程，得 到结构的后临界状态的运动响应，以研究球面网壳是否存在跳跃失稳的传播特性。 二、非线性增量有限元运动平衡方程的建立 考虑几何非线性的有限元结构动力分析的基本方程为 阮】+k。船}+[c玲}+阻玲}-妒O)} (1) 式中t[丘]为结构的弹性刚度矩阵，[《]为结构的几何剐度矩阵，[q为结构的阻尼矩阵，【^扣 为结构的质量矩阵，{F(0}为随时间变化的等效节点载荷，忙)为位移矢量，口}为速度矢 量t f田为加速度矢量。 和结构的弹性HⅡ度矩阵一样，结构的阻尼矩阵和质量矩阵分别由单元的阻尼矩阵和质 量矩阵组集而成，组集的方法与结构弹性刚度矩阵的组集方法完全相同。在本文中，忽略阻 尼的影响，质量矩阵采用一致质量矩阵，可由下式推导得 川=』Nr,oNd (2) Ⅳ为插值函数。 平面杆单元的质量矩阵为 2 0 l 0 O 2 0 l M=譬l 0 2 0 0 l O 2 空间杆单元的质量矩阵为 2 0 0 l 0 0 0 2 0 O l o 0 0 2 0 o I ㈨=警l 0 0 2 o o 0 1 0 0 2 o 0 0 l O 0 2 平面梁单元的质量矩阵为