单自由度非线性体系线加速度法程序编写

单自由度非线性体系线加速度法程序编写 本程序主要应用线加速度法数值求解单自由度非线性体系的位移曲线。整个程序通过Matlab编写。 1.线加速法原理 1.1增量方程 t时刻：    时刻：    其中：    则可得增量方程    1.2线加速度法 设          并消去      得        其中    2.程序编写图    N Y Y y(i+1) N Y y(i+1)∈[A,B] N N N v(i+1)<0 v(i+1)>0 Y Y 3.具体应用算例 通过所编写的程序，得到以下结果： 当Δt=0.1s时，ymax=0.1423,fsmax=3 fs-t曲线 y-t曲线 当Δt=0.05s时，ymax=0.1563,fsmax=3 当Δt=0.025s时，ymax=0.1583,fsmax=3 从以上图像中可以看到，当Δt越小时，所得到的结构精度越高，但可以看到当Δt=0.1s时的位移曲线与Δt=0.05s和Δt=0.025s有区别，说明所编程序存在小误差，仍需改进。