哈工大结构动力学作业_威尔逊_θ法.docxVIP

结构动力学大作业 （威尔逊-法） 姓名: 学号: 班级: 专业: 威尔逊-法原理及应用 【摘要】在求解单自由度体系振动方程时我们用了常加速度法及线加速度法等数值 分析方法。在多自由度体系中，也有类似求解方法，即中心差分法及威尔逊 -法。 实际上后两种方法也能求解单自由度体系振动方程。对于数值方法，有三个重要要 求：收敛性、稳定性及精度。本文推导了威尔逊 -法的公式，并利用MATLABS程 来研究单自由度体系的动力特性。 【关键词】威尔逊-法冲击荷载阻尼比 【正文】威尔逊-法可以很方便的求解任意荷载作用下单自由度体系振动问题。 实 际上，当 1.37时，威尔逊-法是无条件收敛的。 一、威尔逊-法的原理 威尔逊-法是线性加速度法的一种拓展（当 1时，两者相同），其基本思路 和实现方法是求出在时间段 t,t t时刻的运动，其中 1，然后通过内插得到 ti t时刻的运动（见图1.1 ）。 1、公式推导 推导由t时刻的状态求t t时刻的状态的递推公式: y t y t —^( y t t y t) 对积分 -(y t t yt) 3 n(yt t y t) t( y t t y t) t)2( y t) 2 yt y t) Rt 6 t? Pt (Pt t Pt) 6 (t)2 7 y t 2 y t) c(^yt 2 y t 2、 MATLAB源程序: clc;clear; K=in put(请输入结构刚度k(N/m)); M=input(请输入质量(kg)); C=input(请输入阻尼(N*s/m)); t=sym(t);% 产生符号对象t Pt=input(请输入荷载); Tp=i nput(请输入荷载加载时长(s)); Tu=in put(请输入需要计算的时间长度(s)); dt=input( 请输入积分步长(s)); Sita=i nput( 请输入)； uds=O:dt:Tu;%确定各积分步时刻 pds=O:dt:Tp; Lu=le ngth(uds); Lp=le ngth(pds); if isa(Pt,sym)% 荷载为函数 P=subs(Pt,t,uds); % 将荷载在各时间步离散 if LuLp P(Lp+1: Lu)=0; end elseif isnu meric(Pt)% 荷载为散点 if Lu=Lp P=Pt(1:L u); else P(1: Lp)=Pt; P(Lp+1:L u)=0; end end y=zeros(1,Lu);% 给位移矩阵分配空间 y1=zeros(1,Lu);% 给速度矩阵分配空间 y2=zeros(1,Lu);% 给加速度矩阵分配空间 pp=zeros(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⑴-C*y1(1)-K*y(1))/M;% 初始加速度 FF(1)=P(1)-M*y2(1); l=6/(Sita*dt)A2; q=3/(Sita*dt); r=6/(Sita*dt); s=Sita*dt/2; for z=1:Lu-1 kk=K+l*M+q*C; pp(z)=P(z)+Sita*(P(z+1)-P(z))+(l*y(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); end plot（uds,y,丫）,xlabel（时间 t）,ylabel（ 位移 y）,title（ plot（uds,y,丫）,xlabel（ 、利用威尔逊-法求冲击荷载下的结构反应 1、矩形脉冲 研究不同时长脉冲作用下，体系振动位移。取单自由度刚度为 iN/m质量为 1/(4*piA2)kg ，频率为2*pi s 1，周期为is，阻尼c=0,荷载为1N,积分步长为0.1， =1.42，初位移为0,初速度为0时的质点位移时间图如下： 图 2.1 td 1/ 4s 1/ 4Tn 图 2.2 td 1/ 2s 1/ 2Tn 图 2.3 td 1s Tn 2 位移團形 1 1 / \ 1.5 I / \ 「 1 j /