duhamel积分matlab,动力学方程数值解法：振型叠加法与Duhamel积分

Duhamel积分matlab实例

%%%%%%%%%%%%%%%%%

%单自由度Duhamel积分

%内部采用Simpson积分

%初始状态为静止状态

%可计算无阻尼和阻尼强迫震动

%输入可为数组或函数荷载

%只能得出位移时程图

%%%%%%%%%%%%%%%%%

clear all;

%w=30

w=input('输入自振频率 ω:');

%n=10

n=input('输入荷载步数n ：');

%m=96.6/32.3

m=input('输入单元质量m ：');

%T=0.05

T=input('输入外荷载持续时间T：');

%si=0.0

si=input('输入阻尼系数ξ：');

%k=2700

k=input('输入刚度系数k：');

deta=T/n;

wD=w*(1-si^2)^0.5;

G=deta/(m*wD)/3;

t1=[0:deta:T];

reply = input('输入荷载数组直接回车或敲击Y，输入函数荷载输入N： [Y]:','s');

if isempty(reply)

reply = 'Y';

end

if reply=='Y&#