有限单元法(第十章).ppt

河北工业大学土木工程学院 这种方法的优点是计算简单，缺点则在于它是有条件稳定的，即当时间步长Δt过大时，积分是不稳定的。 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 因at+Δt是利用t时刻运动方程得到的，K不出现在刚度矩阵 中，因此中心差分法称为显示积分法。 式（10.37）是用相邻时刻的位移表示的代数方程组，由此可解出at+Δt。 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 河北工业大学土木工程学院 于是，t+τ时刻的速度和位移可积分求得 （2）线性加速度法 （a） 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 假定在时刻[t,t+Δt]内加速度呈线性变化，即 （b） （c） 将式（a）代入上式， τ=Δt得到 河北工业大学土木工程学院 由式（c）得到 （d） 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 将式（d）代入式（b）得到 （e） t+Δt时刻的位移解答at+Δt 可通过满足t+Δt时刻动力学方程 （10.41） 将式（d）（e）代入式（10.41）到由 河北工业大学土木工程学院 其中 （10.42） 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 将式（d）和（e）代入式（10.41））得到由 计算 的公式 （10.43） （10.44） 在线性加速度法中， at+Δt是利用t+Δt时刻的时刻运动方程得到的，K出现在刚度矩阵 中。因此称为隐式积分法。可以证明，线性加速度法也是有条件稳定的。 河北工业大学土木工程学院 其中，γ和β是按积分稳定性要求而决定的参数。当γ=1/2,β=1/6时，式(a)和(b)相应于线性加速度法。 （3）Newmark法 （a） 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 Newmark（1959）假设 （b） 在Newmark法中， t+Δt时刻的位移解答at+Δt 也是通过满足t+Δt时刻动力学方程得到 （10.45） 河北工业大学土木工程学院 代入式（a），再一并代入式（10.45），从而得由 计算 的公式， 首先从式（b）中解得 （10.46） 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 其中 （10.47） （10.48） 河北工业大学土木工程学院 （4）Willson法 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 积分上式可得 （a） Willson法假定加速度 在[t,t+θΔt]内呈线性变化，参数θ≥1.0。以τ表示时间增量，0≤τ≤ θΔt，则在[t,t+θΔt]内有 （b） （c） 河北工业大学土木工程学院 在式（b）和（c）中令τ=θΔt，得 10.3.2 直接积分法 第十章 动力分析 10.3 结构动力响应 由此可得到用at+θΔt表示的 （d） （f） （e） （g） * * 结构物可能受到强烈地震、机械振动或海浪冲击等动荷载作用，使结构发生不允许的变形甚至破坏。因此，结构设计必须考虑动荷载效应，进而要求进行动力学分析。 河北工业大学土木工程学院 第十章 动力分析 本章以动力学问题为例，阐明如何应用有限元进行结构动力分析，内容包括： （1）动力有限元方程及系数矩阵； （2）结构固有特性计算； （3）结构动力响应计算； （4）解的稳定性分析。 河北工业大学土木工程学院 10.1 动力有限元方程 10.1.1有限元方程 第十章 动力分析 振动系统中由质量、阻尼和刚度所引起的三种力以及外荷载，将振动系统离散后，便可得到有限元动力平衡方程组。 （1）一般方程 （10.1） 其中，M、C、K分别为整体质量矩阵，阻尼矩阵； 分别为整体节点加速度，速度和动荷载向量。 河北工业大学土木工程学院 10.1 动力有限元方程 10.1.1有限元方程 第十章 动力分析 对于地震反应问题，直接作用在结构上的荷载P=0,且阻尼力和弹性恢复力分别只与相对速度和相对位移有关，而惯性力则与绝对加速度有关。现仍以 分别表示相对加速度、相对速度和相对位移， 表示地震牵连加速度（地面加速度）则式（10.1）变为