结构动力学电算报告.doc

结构动力学电算报告

摘要：应用matlab软件，采用线加速度差值方法对单自由度体系的弹塑性材料〔两段式恢复力模型〕编写程序，了解弹塑性材料单自由度体系的运动过程与特点。

关键词：弹塑性，线加速度，振动

正文

软件简介

MATLAB是由美国MATHWORKS公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言〔如C、Fortran〕的编辑模式，代表了当今国际科学计算软件的先进水平。

题型简述

f(N)

f(N)

y(m)

y(m)

恢复力与位移关系

恢复力与位移关系

阻尼力与的速度关系

阻尼力与的速度关系

计算公式推导

设加速度方程：

那么：

所以：

故：

又：

且有：

可得：

并且：

假设：，，，，，

可算出：，，，

进而可算出：，，

如此循环即可求出质点各时程的位移

程序流程图

开始

开始

输入质量

输入质量,刚度，进入塑性阶段的位移，时段长，阻尼系数，外荷载，初始位移，初始速度

判断速度是否反向？

判断速度是否反向？

计算

计算

是否

是

否

算出速度反向时的位移和下次进入塑性阶段的位移

算出速度反向时的位移和下次进入塑性阶段的位移

计算

计算，，

判断是否进入塑形阶段？

判断是否进入塑形阶段？

计算

计算

判断进入塑性阶段？

判断进入塑性阶段？

计算

计算，，，，

计算

计算

判断是否继续运算？

判断是否继续运算？

是

是

否

否

绘

绘制位移时程图与恢复力时程图

结束

结束

程序说明

该程序程序运行后，可根据提示先后输入：质量，弹性阶段刚度，进入塑性阶段的位移，时段长，阻尼系数，外荷载，初始位移，初始速度。

其中须注意的是弹性阶段刚度为常数，阻尼系数也为常数。外荷载采用一维矩阵散点输入。

另外，当塑性阶段的位移很大，其振动最大位移不可能到达时，本程序计算的是线弹性体系。

源程序

functionwangxiangjun()

m=input(请输入质量m);

k=input(请输入弹性阶段刚度k);

ff=0;

Y(1)=input(请输入进入塑性阶段的位移y);

w=2*pi/sqrt(k/m)/10;

T=input([请输入时段长ΔT，〔ΔT应小于等于num2str(w))]);

c=input(请输入阻尼系数c);

HH=input(外荷载以函数形式输入请输入1，以矩阵形式输入请输入2);

ifHH==1

symsx

ph=input(请输入外荷载P的随时间x的函数,s);

hh=input(请输入外荷载作用时间);

n=floor(hh/T);

fori=0:n

p(i+1)=subs(ph,i*T);

end

else

p=input(请输入外荷载P);

end

sj=input(请输入外荷载作用结束后，还需绘制出多少倍外荷载作用时间的位移时程图);

n=length(p);

fori=n+1:(sj+1)*n

p(i)=0;

end

y(1)=input(请输入初始位移y(0));

v(1)=input(请输入初始速度v(0));

ifabs(y(1))Y(1)

y1=y(1);

f(1)=-k*y1;

M=1;

elseify(1)=-Y(1)

y1=-Y(1);

f(1)=-k*y1;

M=0;

elseify(1)=Y(1)

y1=Y(1);

f(1)=-k*y1;

M=0;

end

end

end

a(1)=(p(1)-c*v(1)-k*y1)/m;

K(1)=k+6*m/T/T+3*c/T;

PP(1)=(p(2)-p(1))+m*(6/T*v(1)+3*a(1))+c*(3*v(1)+T*a(1)/2);

yy(1)=PP(1)/K(1);

vv(1)=3*yy(1)