非线性有限元法5---7.pptVIP

4 方程解法 4.1 显示方法（Explicit Methods） 在计算力学和计算物理中应用最广泛的显式方法是中心差分法 定义时间步增量为 * 4 方程解法 速度的中心差分公式为 积分公式为 加速度及积分公式 * 4 方程解法 在第n时间步的有限元半离散方程为 广义约束方程 代半离散方程到加速度积分公式中 这样，在任意时间步n，已知位移 ，通过顺序计算应变-位移方程、由 或 形式表示的本构方程和外力，可确定节点力 ，上式右端全部赋值，可获得 ，由位移的积分公式，可确定 。 * 4 方程解法 特点 a) M矩阵是常数阵，在整个求解过程中不变; b) 如果M为对角阵，实现节点速度和节点位移的更新可以不用求解任何方程； c) 因此，在显示方法中，离散的动量方程一般都用团聚质量矩阵。 —→在显示方法中，对离散的动量方程的时间积分不需求解任何方程 * 4 方程解法 流程图 * 4 方程解法 流程图（续） * 4 方程解法 优点－－ 方法简单，不求解方程 显式时间积分是很强健的（Robust），很少因为数值失败而终止 缺点－－ 显式方法是条件稳定的 如果时间步长超过一个临界值，其计算结果会增长至无穷 临界时间步长 对于采用率无关材料的常应变单元的网格，稳定时间步长为 线性系统的最大频率， 单元e的特征长度 ， 单元e的当前波速， 是折减系数。 * 4 方程解法 有限差分的临界步长问题称为Courant条件 1928年Courant，Friedrichs和Lewy首先发表这一结果。 临界步长随网格细划和材料刚度的增加而减小 对于弹－塑性问题，塑性响应会减慢波速，但并不能增加临界时间步长，因为，在数值计算中常常发生卸载 从单元时间步长确定网格时间步长 * 4 方程解法 4.2平衡解法和隐式时间积分（Implicit time integration） 离散的动量方程可写为 Newmark 方程 更新的位移和速度为 的意义为人工粘性或人工阻尼。 * 4 方程解法