杜哈梅积分.docVIP

杜哈梅积分程序 （1）任务：1、试编写单自由度杜哈梅数值积分的计算程序。 （2）要求：①可以使用至少2种积分方法； ②可以处理的荷载种类：a.荷载数组；（？）b.脉冲荷载；c.函数荷载。中的两种。 （3）绘出系统响应时程曲线，求出最大响应（位移、速度、加速度）。 （4）算例分析： 图示质量m受爆炸荷载作用，P=100kN，m=6.4kg，结构刚度k=34847.77N/m，阻尼比为5%，荷载如下图示： 所做MATLAB程序有6个文件（真正包含5个有用文件，还有1个为组合梯形积分函数，共参考，未利用）。文件名：lFunction.m，FuncLoad.m，simpson.m，trapezoid.m，stepFuncTion.m。 各文件所实现的功能： example.m ：运行程序，调用FuncLoad函数，传入lFunction函数句柄及时间向量、结构刚度、集中质量、阻尼比，然后由FuncLoad函数生成的位移、速度、加速度向量，绘制“位移—时间”、“速度—时间”、“加速度—时间”曲线，同时找到最大响应值，并做标注。 lFunction.m ：荷载的函数，荷载关于时间t的函数，作为句柄导入FuncLoad函数中。在A(t)、B(t与及的乘积形成新的被积函数。simpson.m ：辛普森积分函数，输入参数为函数句柄、积分上下限（为标量），输出为对应的积分值。 trapezoid.m ：梯形组合积分函数，输入、输出同辛普森积分函数。（功能与FuncLoad.m ：连续荷载处理函数，输出为对应时间向量的输入参数为荷载函数句柄、时间向量、结构刚度、集中质量、阻尼比 stepFuncTion.m ：脉冲荷载处理函数，输出为对应时间向量的输入参数为荷载函数句柄、时间向量、荷载终止时刻、结构刚度、集中质量、阻尼比。调用它，可计算脉冲荷载作用下结构的响应。 lFunction.m中已写出，计算不同荷载的响应，只需将要计算的荷载函数取消注释，运行辛普森积分，速度、加速度的数值算法为边界节点为向前或荷载函数作为函数句柄导入FuncTion、stepFuncTion对于位移响应计算，杜哈梅积分的数值法在理论上应该对于任意荷载都成立，但是本程序在载入跳跃荷载函数后，位移荷载响应图形貌似符合要求，但根据位移响应生成的速度、加速度响应出现问题（出现无数次突变），因此可判断位移响应有问题。 程序文件如下 simpson.m function fv =simpson(f,a,b) %辛普森积分法函数 %input - @f 为被积函数的句柄 % - a 为积分下限 % - b 为积分上限 %output - fv 为积分值 A=a:(b-a)/20000:b; %A为(a,b)的20000等分横坐标向量 B=f(A); %B为A对应的函数值向量 fv=0; h=(b-a)/20000; for j=1:10000; s=h/3*(B(2*j-1)+4*B(2*j)+B(2*j+1)); fv=fv+s; end trapezoid.m function fv=trapezoid(f,a,b) %组合梯形积分法函数 %input - f 为被积函数 % - a 为积分下限 % - b 为积分上限 %output - fv 为积分值 A=a:(b-a)/10000:b; %A为(a,b)的10000等分横坐标向量 B=f(A); %B为A对应的函数值 fv=0; h=(b-a)/10000; for n=1:10000; s=h/2*(B(n)+B(n+1)); %s 为小梯形面积 fv=fv+s; end stepFuncLoad function [u,v,a]=stepFuncLoad(lFunction,t,t1,k,m,Xi) %函数荷载作用下的单自由度体系动力响应——位移u、速度v、加速度a %input - ff2 为何在函数 % - t 为时间向量（单位：s） % - k 为结构刚度（单位：N/m） % - m 为集中质量（单位：kg） % - Xi 为阻尼比ξ %output - u 为位移列向量 % - v 为速度列向量 %