[修改版]中心差分法在单自由度中的应用

《[修改版]中心差分法在单自由度中的应用》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、中心差分法求解单自由度体系的自由振动问题刖吞吋域逐步积分法是根据运动方程，引进某些假设，建立ilitW刻状态向量均、；仃、ui到t+Ar时刻的状态向最终屮、u/+1、“心1的递推关系，从而从t=o时刻的初始状态向最“°、；（）、；；（）出发，逐步求出各时刻的状态向量，由于引进的假设条件不同，可以有各种不同的方法，下面主耍介绍一种吋域逐步积分方法一中心差分法。中心差分法（centraldifferencemethod）原理山中心差分法的基木思路:是将运动方程中的速度向量和加速度向量用位移的某种组介来表示，将微分方程组的求解问题转化为代数方程组的求解问题，并在时间区间内求得每

2、个微小时间区间的递推公式，进而求得整个时程的反应。中心差分法只在相隔一些离散的时间区间内满足运动方程，其基于有限差分代替位移对时间的求导（即速度和加速度），如果采用等时间步长，Arf.=t,则速度与加速度的中心差分近似为:2Ar由式(e)就nJ以根据厶及厶以询时刻的运动，求得如吋刻的运动，如果需要可以用式(a)和式(b)求得体系的速度和加速度。冷+i—2%+给]Ar(b)而离散时间点的运动为冷仃="(—)(i=0」2,3,)由体系运动方程为：muQJ+cw（r）+ku{t}=0(C)将速度和加速度的差分近似公式（a）和式（b）代入式（c）nJ-以得到心时刻的运动方程:w

3、；.i一2色+u-!ifm———十c——+gAr(d)在（d）式中，假设％和是已知的，即在―及-以前时刻的运动已知，则可以把已知项移到方程的右边，整理得到：zmC门2加加C、V+肓叽7一計7丽一百）"I(e)假设给定的初始条件为”()="(0),・・(g)wo=w(0),由式(g)确定—］。在零时刻速度和加速度的中心差分公式为：Wo2Ar(h)••u,-2wn+u将式(i)消去旳得:(j)而零时刻的加速度值uo町以用t=()时的运动方程muo+cuo+ku()=0确定(k)这样就可以根据初始条件知⑷和初始荷载弘就川•以根据上式确定心的值。下面给岀采用中心差分法分析时的具

4、体计算步骤:(1)基本数据准备和初始条件计算••1•Wo=—(―cuo—ku())m人•Ar2••=uo一u()2(2)计算等效刚度和屮心并分计算公式屮的和关系数r2M2mCl—k—aAr2,mcb=Ar22Ar(3)根据匚及耳以前时刻的运动，计算/冲时刻的运动P=-aui—bu_p/k如果需要，可计算妁+1—组2Ar£+1-2心+%•_]Ar2(4)下一步计算用i+1代替i,对于线弹性结构体系，重复第3步，对于非线性结构体系，重复第2步和第3步。以上为中心差分法逐步计算公式，其具冇2阶精度，即误差£*0(&2)；并且为有条件稳定，稳定条件为:271上式小，7；为结构

5、的自振周期，对于多自由度结构体系则为结构的最小白振周期。算例对于一个单层框架结构，假设楼板刚度无限大，且结构质量集中于楼层，其质量M=2000kg、刚度K=50KN/m、阻尼系数C=3KNs/m,假设结构处于线弹性状态，川中心差分法计算结构的白由振动反应。采用MATLAB语言编程，并以单自由度体系为例进行计算，设初位移uO=O和初速度v0=0,T2由稳定条件Ar<^=—=5rad/s,取不同的步长分别计算，以验证屮心差分法的稳定条件Ar<^o先计算卜，则厶0.4兀55所以木次计算取△/=0.1,0.3,0.4,0.41,0.42,0.45分别进行计算MATLAB程序清单f

6、unction[u,v,ac]=centraldifferent(M,C,K,uO,vO,time,dt)%本程序采用中心差分法计算结构的动力响应%本程序是既可以计算单自由度体系又可以计算多口由度体系，且均假设结构体系处于线弹性状态；%%%%%%输入参数%%%%%%%%M质量矩阵%C肌尼矩阵%K刚度矩阵%uO初始位移%vO初始速度%time模拟时间%dt时间步长%・・・％%%%%%输出值％%%%%%%%——%u位移%V速度%ac加速度%%%%%%%%%中心差分法主要公式及原理%%%%%%%%%%%MX“+CX'+KX=0%M*(X(t+dt)-2*X(t)+X(t-dt)

7、)/(dtA2)4-C*(X(t+dt)-X(t-dt))/(2*dt)+K*X(t)=0%(M/dtA2+C/2*dt)*(X(t+dt))=-(K-2*M/dtA2)*X(t+dt)-(M/dtA2-C/2*dt)*X(t-dt)%等效刚度Ke等效荷载Pe和相关系数a,b%Kc=M/dtA2+C/2*dt%a=K-2*M/dtA2%b=M/dtA2-C/2*dt%Pe=-a*X(t)-b*X(t-dt)%X(t+dt)=Pe/Ke%X(t)=(X(t+dt)-X(t-dt))/(2*dt)%X(t)H=(X(t+dt)-2