哈工大结构动力学作业威尔逊法.docx

《哈工大结构动力学作业威尔逊法.docx》由会员分享，可在线阅读，更多相关《哈工大结构动力学作业威尔逊法.docx（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、构造动力学大作业 威尔逊-法 姓名： 学号： 班级： 专业： 第 6 页威尔逊-法原理及应用【摘要】在求解单自由度体系振动方程时我们用了常加速度法及线加速度法等数值分析方法。在多自由度体系中，也有类似求解方法，即中心差分法及威尔逊-法。实际上后两种方法也能求解单自由度体系振动方程。对于数值方法，有三个重要要求：收敛性、稳定性及精度。本文推导了威尔逊-法的公式，并利用MATLAB编程来研究单自由度体系的动力特性。【关键词】威尔逊-法 冲击荷载 阻尼比【正文】威尔逊-法可以很方便的求解任意荷载作用下单自由度体系振动问题。实际上，当时，威尔逊-法是无条件收敛的。一、威尔逊-法的原理 威尔逊-法是线性

2、加速度法的一种拓展当时，两者一样，其根本思路与实现方法是求出在时间段时刻的运动，其中，然后通过内插得到时刻的运动见图 1.1。1、公式推导推导由时刻的状态求时刻的状态的递推公式：对积分2、MATLAB源程序：clc;clear;K=input(请输入构造刚度k(N/m);M=input(请输入质量(kg);C=input(请输入阻尼(N*s/m);t=sym(t);%产生符号对象tPt=input(请输入荷载);Tp=input(请输入荷载加载时长(s);Tu=input(请输入需要计算的时间长度(s) );dt=input(请输入积分步长(s);Sita=input(请输入);uds=0:d

3、t:Tu;%确定各积分步时刻pds=0:dt:Tp;Lu=length(uds);Lp=length(pds);if isa(Pt,sym)%荷载为函数 P=subs(Pt,t,uds); %将荷载在各时间步离散 if LuLp P(Lp+1:Lu)=0; endelseif isnumeric(Pt)%荷载为散点 if Lu=Lp P=Pt(1:Lu); else P(1:Lp)=Pt; P(Lp+1:Lu)=0; endendy=zeros(1,Lu);%给位移矩阵分配空间y1=zeros(1,Lu);%给速度矩阵分配空间y2=zeros(1,Lu);%给加速度矩阵分配空间pp=zeros

4、(1,Lu-1);%给广义力矩阵分配空间yy=zeros(1,Lu-1);%给y(t+theta*t)矩阵分配FF=zeros(1,Lu);%给内力矩阵分配空间y(1)=input(请输入初位移(m);y1(1)=input(请输入初速度(m/s);%-初始计算-y2(1)=(P(1)-C*y1(1)-K*y(1)/M;%初始加速度FF(1)=P(1)-M*y2(1);l=6/(Sita*dt)2;q=3/(Sita*dt);r=6/(Sita*dt);s=Sita*dt/2; for z=1:Lu-1kk=K+l*M+q*C;pp(z)=P(z)+Sita*(P(z+1)-P(z)+(l*y

5、(z)+r*y1(z)+2*y2(z)*M+(q*y(z)+2*y1(z)+s*y2(z)*C;yy(z)=pp(z)/kk;y2(z+1)=l/Sita*(yy(z)-y(z)-l*dt*y1(z)+(1-3/Sita)*y2(z);y1(z+1)=y1(z)+dt/2*(y2(z+1)+y2(zp);y(z+1)=y(z)+y1(z)*dt+dt*dt/6*(y2(z+1)+2*y2(z);FF(z+1)=P(z+1)-M*y2(z+1);endplot(uds,y,r),xlabel(时间 t),ylabel(位移 y),title(位移图形)二、利用威尔逊-法求冲击荷载下的构造反响1、

6、矩形脉冲研究不同时长脉冲作用下，体系振动位移。取单自由度刚度为1N/m，质量为1/(4*pi2)kg，频率为2*pi，周期为1s，阻尼c=0，荷载为1N，积分步长为0.1，2，初位移为0，初速度为0时的质点位移时间图如下：图2.1 图2.2 图2.3 图2.4 图2.5 由图形可看出：当时，最大位移发生在荷载离开前； 当时，最大位移发生在荷载离开后； 当时，最大位移发生在荷载离开时。特别的，当时，后没有位移。2、其他脉冲图 负斜率直线图2.7 正斜率直线图2.8 二次抛物线图2.9 5次抛物线三、利用威尔逊-法求不同阻尼下构造自振反响本体系度刚度为1N/m，质量为1/(4*pi2)kg，频率为2*pi。故其临界阻尼pi。分别取构造阻尼c为：0.05/pi，0.1/pi，1/pi，1.5/pi，进展计算。计算结果见下列图：图 图 图 图 由图形比照可知，当时，阻尼越大，构造运动衰减越快，此时构造处于小阻尼状态；当体系不能振动，此时的c为临界阻尼；当时，为超阻尼系统，不发生自由振动。四、利用威尔逊-法求实例：如下列图，求：体系位移反响。图4.1 构造 图4.2 荷载由MATLAB计算出的结果见图 。图4.3 位移曲线手算结果：手算结果与电算结果近似相等，说明当冲击荷载作用时间远小于构造自振周期是，可近似认为最大位移S为荷载与坐标轴所围成的面积。