整理得到 A1{un+1}=A2+A3|{un}+A4l{n 其中, A1=M△t+B/2△t+K/3 Dynamic Matrix l=1/3{Pn+1+Pn+Pn-1 Applied force A2l=[2M/△t-K/3 Initial Conditions from Previous M/△t+B/2△t-K/3 Time Step (2)暻态响应分析中的阻尼 +1/W ∑KE 其中 B1=阻尼单元(vsc,DAMP)+B2GG B2=B2PP直接输入矩阵+传递函数 G=整体结构阻尼系数( PARAM.G) W3=感兴趣的整体结构阻尼转化为频率弧度秒( PARAMW3) K1=整体刚度矩阵 GE=单元结构阻尼系数(GE在MAT卡中定义 W4=感兴趣的单元结构阻尼转化为频率弧度秒 PARAMW4 KE=单元刚度矩阵 瞬态响应分析中的不允许复系数,因此结构阻尼转化为等效粘性阻尼进行计算 W3,W4的缺省为0,这时不计阻尼 7.3模态瞬态响应分析 (1)过程 物理坐标与模态坐标变化 无阻尼的动力学方程整理得到 其中， （2） 瞬态响应分析中的阻尼 其中， 瞬态响应分析中的不允许复系数，因此结构阻尼转化为等效粘性阻尼进行计算 W3,W4 的缺省为 0，这时不计阻尼 7.3 模态瞬态响应分析 （1） 过程 物理坐标与模态坐标变化 无阻尼的动力学方程 B1 = 阻尼单元(VISC,DAMP) + B2GG B2 = B2PP 直接输入矩阵+传递函数 G = 整体结构阻尼系数 (PARAM,G) W3 = 感兴趣的整体结构阻尼转化为频率-弧度/秒 (PARAM,W3) K1 = 整体刚度矩阵 GE = 单元结构阻尼系数 (GE 在MATi卡中定义) W4 =感兴趣的单元结构阻尼转化为频率-弧度/秒(PARAM,W4) KE = 单元刚度矩阵