数值计算求解非线性动力方程.pdfVIP

科学之友 FriendofScienceAmateurs 2011年04月 数值计算求解非线性动力方程 郭 冲 (山西省公路局长治分局勘测设计所 ，山西 长治 046000) 摘 要：由于非线性动力方程中的外荷载或结构体系随时间发生改变，其解析解往往很 难获得 ，需要采用数值算法进行求解。文章首先介绍 了激励线性插值法、中心差分法、 Newmark的平均加速度法和线性加速度法的基本原理，随后给 出了在Matlab中编制的函 数范例及各方法的动力响应求解结果。经与理论解比较，验证 了所编制函数的正确性．可 作为工程设计人员求解一般非线性动力响应问题的一条简便途径。 关键词：非线性动力方程；激励线性插值法；中心差分法；平均加速度法；线性加速度法 中图分类号：TB122 文献标识码：A 文章编号：1000—8136(2011)11—0053—03 1 引言 一 n(寿 ) 如果外荷载或者结构的刚度、质量 、阻尼等任何一项随时间 变化 ，那么这样的结构体系就是非线性的，这种结构体系的动力 -一 Afcosset-寿 方程也是非线性的。对于此类非线性动力方程，其解析解往往较 难由理论推导得到，这类问题可以通过数值计算方法对动力响 c=}{- l／ 一△ J＼ 应方程进行求解。常用的数值计算方法包括激励线性插值法、中 心差分法、Newmark的平均加速度法和线性加速度法等。 sin∞△t一古c。s∞△t 2 各数值算法的基本原理 。= 出( 刚 )} 2．1 激励线性插值法 (E，：外荷载频率； 对于线性体系，通过在每个时间间隔里对激励进行插值，并 ：结构固有频率。 利用线性体系响应解析解，能够推导出一种有效的数值计算方法。 给定初始条件，代入式 (1)即可得到结构在任意时刻的响应。 只要保证插值的时间间隔较短，就可以得到令人满意的求解结果。 2．2 中心差分法 对于欠阻尼体系(1)，有非线性动力方程 该方法是基于对位移时间导数的有限差分近似进行的。取固 定的时间步长△ ，则时刻 i的速度和加速度的中心差分格式为： m 在时间间隔0≤ ≤Ati内，反应u(T)由3部分组成： 警 (1)r=0时刻的初位移 。和初速度如引起的自由振动响应； 将上式的速度加速度代人动力方程中得： (2)零初始条件下对阶跃力P的反应； k +1 (3)零初始条件下对斜坡力(Ap]At)r的反应。 ． 经推导可得出第 i+1个时间步的位移与速度公式如下： 其中：后盎 IAi+I=A + i+c i+上 。+1 Ui+l=AMi+日U_+c 。+D (1) _一|I【A()一2Atl—一r寺(△￡)J] 为了计算 u 需要知道位移u和u。 式 =e ( si+c～ ) 一 。。 若已知初始位移