哈工大结构动力学作业威尔逊θ法.docx
- 文档编号:6349527
- 上传时间:2023-01-05
- 格式:DOCX
- 页数:9
- 大小:156.14KB
哈工大结构动力学作业威尔逊θ法.docx
《哈工大结构动力学作业威尔逊θ法.docx》由会员分享,可在线阅读,更多相关《哈工大结构动力学作业威尔逊θ法.docx(9页珍藏版)》请在冰豆网上搜索。
哈工大结构动力学作业威尔逊θ法
结构动力学大作业
(威尔逊“法)
姓名:
学号:
班级:
专业:
威尔逊-法原理及应用
【摘要】在求解单自由度体系振动方程时我们用了常加速度法及线加速度法等数值分析方法。
在多自由度体系中,也有类似求解方法,即中心差分法及威尔逊「法。
实际上后两种方法也能求解单自由度体系振动方程。
对于数值方法,有三个重要要求:
收敛性、稳定性及精度。
本文推导了威尔逊-法的公式,并利用MATLAB编程来研究单自由度体系的动力特性。
【关键词】威尔逊-法冲击荷载阻尼比
【正文】威尔逊7法可以很方便的求解任意荷载作用下单自由度体系振动问题。
实际上,当1.37时,威尔逊-法是无条件收敛的。
一、威尔逊…法的原理
威尔逊“法是线性加速度法的一种拓展(当二=1时,两者相同),其基本思路和实现方法是求出在时间段H,r^.:
t1时刻的运动,其中二_1,然后通过内插得到1迸时刻的运动(见图1.1)。
1、公式推导
推导由t时刻的状态求tt时刻的状态的递推公式:
2
23
冷L」y「'yJ石$厂6t("'曰一S)
.-n't
:
y:
tr=:
y「八t、y厂扌八tcy〔宀-2t)
2
yt.t.-'ytt"yl(◎(:
y;.-12:
y[)
6
:
y;=(r:
fCy[Tt_$•〔)_2:
y;
:
y】Tt=^Cy;i虫一:
y】)一2:
yi-予:
y]
m,'y.;:
「Cty;亠k‘yfP:
mW-t'CVt-t'k-yA■-1二
k'7i■t=Rtut
"kJmi^ci
(仙)29At
R>pL*Pi.t」P:
t)m旨6yMt吕;St・2如)cl(走St2^t号⑺t)
2、MATLAB源程序:
clc;clear;
K=input('请输入结构刚度k(N/m)');
M=input('请输入质量(kg)');
C=input('请输入阻尼(N*s/m)');
t=sym('t');%产生符号对象t
Pt=input('请输入荷载);
Tp=input('请输入荷载加载时长(s)');
Tu=input('请输入需要计算的时间长度(s)');
dt=input('请输入积分步长⑸');
Sita=input('请输入二');
uds=O:
dt:
Tu;%确定各积分步时刻
pds=0:
dt:
Tp;
Lu=length(uds);
Lp=length(pds);
ifisa(Pt,'sym')%荷载为函数
P=subs(Pt,t,uds);%将荷载在各时间步离散ifLu>Lp
P(Lp+1:
Lu)=0;
end
elseifisnumeric(Pt)%荷载为散点ifLu<=Lp
P=Pt(1:
Lu);
else
P(1:
Lp)=Pt;
P(Lp+1:
Lu)=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
(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;
forz=1:
Lu-1kk=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('位移图形')
二、利用威尔逊7法求冲击荷载下的结构反应
1、矩形脉冲
研究不同时长脉冲作用下,体系振动位移。
取单自由度刚度为1N/m,质量为
1/(4*piA2)kg,频率为2*pisJ,周期为1s,阻尼c=0,荷载为1N,积分步长为0.1,
二=1.42,初位移为0,初速度为0时的质点位移时间图如下:
位移图形
图2.1td=1/4s=1/4Tn
图2.2td-1/2s=1/2Tn
位移團形
图2.3td=1s=Tn
位移團形
图2.4td=1.25s=1.25Tn
图2.5td=1.5s".5Tn
由图形可看出:
当tdTn/2时,最大位移发生在荷载离开前;
当td 当td二Tn/2时,最大位移发生在荷载离开时。 特别的,当td=人时,t丸后没有位移。 2、其他脉冲 图2.6负斜率直线 图2.7正斜率直线 图2.8二次抛物线 位移團形 时间t 图2.95次抛物线 三、利用威尔逊宀法求不同阻尼下结构自振反应 本体系度刚度为1N/m,质量为1/(4*pW2)kg,频率为2*pis'。 故其临界阻尼 Cr=2m•=1/pi。 分别取结构阻尼c为: 0.05/pi,0.1/pi,1/pi,1.5/pi,进行计算。 计算结果见下图: 图3.1--0.05 图3.2『: =0.1 图3.3=1 图3.4=1.5 由图形对比可知,当;: : 1时,阻尼越大,结构运动衰减越快,此时结构处于小 阻尼状态;当“体系不能振动,此时的c为临界阻尼;当•・1时,为超阻尼系统,不发生自由振动。 四、利用威尔逊—法求实例 已知: 如下图,W=438.18kN;k=40181.1kN/m 求: 体系位移反应。 由MATLAB计算出的结果见图4.3 图4.3位移曲线手算结果: 3 =43.8810kN T=—=0.208s co 133 S0.02543010=5.375103 2 3 5.375103 4.05710m 43.881030.19 说明当冲击荷载作用时间远小于结构自振周期 手算结果与电算结果近似相等, 是,可近似认为最大位移^― m® (S为荷载与坐标轴所围成的面积)。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 哈工大 结构 动力学 作业 威尔逊
![提示](https://static.bdocx.com/images/bang_tan.gif)