动力反应数值分析方法.ppt

* * 动力反应数值分析方法 分段解析法 中心差分法 Newmark--β法 Wilson--θ法 lunge-kutta法 引言： 1、当P(t)为解析函数时 u一般为解析函数 当P(t)为复杂时 寻求数值解 并且，P(t)数值较小时，即满足线弹性时，叠加原理适用； P(t)数值较大时，叠加原理不再适用 几何非线性 材料非线性 2、逐步积分法主要考虑以下几个方面 收敛性，精度，稳定性，效率 一、分段解析法 把解析函数表示的P(t)分成若干小段，在每一小段上认为P(t)直线，即P(t)由小的直线段组成 当 时 在此时间段内，振动方程为： 初始条件： 特解： 齐次通解： 全解： 代入初始条件，确定A,B 其中， 当 此式给出了根据ti时刻运动及外荷载，计算ti+1时刻的运动递推公式 二、中心差分法 如果是等步长 代入运动方程， 此即为，已知ti，ti-1时刻运动，确定ti+1时刻运动----二步法 关于起步问题： 1、对于地震作用，可设： 消去u1，得 2、对于非零初始问题 ： 主要需要确定u-1 其中， 练习题： 3m 3m 已知，阻尼 梁柱均为C30混凝土，截面300*300 三、Newmark—β法 单步法（自起步） 仍然是把时间离散化，运动方程仅需在离散的时间点上满足即可 已知 求 假设： 之间的加速度是介于 的某一常量a 对a 在 时间段上积分 在 时刻 把上页得到的 代入 初始值 然后用 确定 四、wilson—θ法 假定加速度在 内线性变化 按照线性关系写出 时刻的运动 然后内插获得 时刻的运动 按照加速度线性假设： 积分： 当 时 由后式得： 代入前式得： 运动方程 其中， 线性外推获得 把得到的 代入运动方程得 由 解出 代入 得 再把 代入 并令 得 以下开始内插 令 得 再看原来 表达式的两个积分式 此4个式子便是wilson-θ法 练习题： 已知，阻尼 EI1=∞ C m=17.5T 1、准确解曲线 2、中心差分法时程曲线 3、Newmark--β法时程曲线 要求画在同一幅图上 *