newmark法程序法计算多自由度体系动力响应.pdfVIP

. Newmark - β 用 matlab 编程实现 法 计算多自由度体系的动力响应 Newmark - β 用 matlab 编程实现 法 计算多自由度体系的动力响应 一、 Newmark - β法的基本原理 Newmark- 法是一种逐步积分的方法，避免了任何叠加的应用，能很好的适 应非线性的反应分析。 Newmark- 法假定： { u} t t { u} t [(1 ){ u} t { u} t t ] t (1-1) . . 1 2 { u} t t { u} t { u} t [( ){ u} t { u} t t ] t (1-2) 2 式中， 和 是按积分的精度和稳定性要求进行调整的参数。 当 =0.5 ， =0.25 时，为常平均加速度法，即假定从 t 到 t + t 时刻的速度不变，取为常数 1 ({ u} t { u} t t ) 。研究表明， 当 ≥0.5 ， ≥0.25(0.5+ )2 时，Newmark- 2 法是一种无条件稳定的格式。 由式(2-141) 和式 (2-142) 可得到用 { u} t t 及 { u} t ，{ u} t ，{ u} t 表示的 { u} t t ，{ u} t t 表达式，即有 1 1 1 { u} t t 2 ({ u} t t { u} t ) { u} t ( 1){ u} t (1-3) t t 2 { u} t t ({ u} t t { u} t ) (1 ){ u} t (1 ) t{ u} t (1-4) t 2 考虑 t + t 时刻的振动微分方程为： [ M ]{ u} t t [C]{ u} t t [ K ]{ u} t t { R} t t (1-5) 将式 (2-143) 、式 (2-144) 代入 (2-145) ，得到关于 ut + t 的方程 [ K ]{ u} t t { R} t t (1-6) 式中