物理非线性杆动力分析的有限元方法.pdfVIP

·980· ‘工程力学’增刊 1999年 物理非线性杆动力分析 的有限元方法 高经武 蔡冰 蔡中民 (太原理工大学应用力学研究所030∞4) (太原市建筑设计院) (太原理工大掌应用力学研究所030024) 提基 本文利用有限元方法对物理非线性粘弹性等直杆受轴向简谐外力激励的动态响应进行了数值 计算，分析了其发生分叉和混沌的参数特性，发现当作为参数的外激励振幅小于某一临界值时，系统的 响应是规则的周期振动，若大于该值，系统将出现分叉和混沌，文中给出了特定参数下的奇怪吸引子图 形以及参数空间上的分叉图． 关t词动力系统，有限元，～致质量矩阵，分叉，混沌． 1前言 众所周知，杆的轴向受迫振动是动力机械、抗震结构等设备中的常见情形，因而对杆 在轴向激励下的动态响应的研究也变得至关重要，而且随着材料科学的进步。各种新型复 合材料以及其它一些高分子材料如聚合物等已经获得了广泛的应用，由于材料的物理非线 性所导致的复杂的动态响应也逐渐引起人们广泛的关注。目前固体力学中对于分叉和混沌 的研究主要集中于几何非线性性情彤下的梁和板的横向振动11l口lⅢ，并且大都采用Galerkjn 法t选取某三角函数作为其正交函数系，而且在结果中一般只取其首项或者说一个模态来 进行分析，但所取函数系是否完备，所得级数是否收敛，都必须进行考虑，尤其当系统的 运动表现为混沌的时候。由于非线性问题的复杂性，目前尚无一个统一有效的解决办法， 每一个具体问题似乎都要求发明特殊的算法，运用新颖的技巧。本文利用有限元方法对物 理非线性杆进行分析，其本质上是将位移、速度等在空间的分布用有限个节点信息来代替． 这样做必然要以损失一部分空间关联和演化信息为代价；但却减少了系统的自由度，简化 了计算，而且避免了Gal删n法主观上引入的所谓模态概念。有限元方法虽已广泛应用于 工程、机械以及连续介质的运动分析，并取得了良好的效果。但将其应用于物理非线性材 料分叉与混沌运动的研究，还未见报道。本文就是针对此方法在固体力学混沌研究中的可 行性作出的一些偿试。 膏最武．勇．1螂．1l出生．工拳曩士，讲师 山西省直嚣科学基金量璐疆茸 ‘工程力学)增刊 1999年 ‘98l· 2数学模型 设有一长为L，横截面积为A，密度为p的非线性粘弹性圆柱杆，初始时刻处于自然状 态，左端固定，自由端受到简谐外力激励，则其所满足的基本方程如下： 运动方程： ∥西=昙(一盯) (1) 苏 本构方程； 盯=Es+甜2+肛3+刁占 (2) 几何方程： 占=罢 叠：芒兰 (3) OX c!n宵 初值条件： 如缈=0；西(x，0)=O 边值条件：“rD，u：D：盯(f，f)：曼磐 (41 把(2)和(3)式代入(1)式，则有： 窘=扣胁筹+邛r耖，窘峙急 ㈤ 式中：如=掣加窘，㈣㈨J是位移，口为应力，￡为纵向应变，E和 口，∥是材料常数，T1为粘性系数。由(4)式和(5)式组成杆的基本运动方程式。 3结构离散 对以上偏微分方程直接进行分析求解是很困难的，我们先通过有限元离散的方法，将 问题转化为一个二阶非线性常微分方程组，然后再进行数值计算。离散形式如下： 我们把长度为L的杆分割为n个相等的单元．D=％工』…x．…工。=工，x。 内采用线性位移模式，在惯性力的计算中采用～致质量矩阵，并令：，；z加，引入无量刚 量： 并在最后的结果中去掉+号。则可得以常微分方程组形式表示的结构动力方挥： ·982·