Matlab在地震工程研究中的应用实例研究.pdfVIP

《工程力学》增刊 1999年 .975 。 Matlab在地震工程研究中的应用实例 葛映宇 叶献国 (合肥工业大学建工系，合肥230009) 提要 Matlab是解决科学工程计算问瓜的高效工具.本文具体论述了如何用Matlab语言实现求地震 反应谱，以及由加速度记录求积分速度，积分位移，由此说明Hatlab语言在地展工程研究中具有良好的 实用性，值得在土木工程领域推广使用。 关健词 Hatlab语言，地震反应谱，积分速度，积分位移。 1.前言 Matlab语言是一种高效率的用于科学工程计算的高级语言。与Basic,Fortran,Visual C+十等语言比较，Matlab的语法规则简单，更加贴近人的思维方式，易学易用。用Matlab写 程序，犹如在一张演算纸上排列公式和求解问题，因而被称为 “演算纸方式的”科学工程 的算法语言u〔。众所周知，Fortran,Basi。等语言在土木工程计算编程方面具有广泛的 应用，但涉及矩阵理论、数值分析等问题时用Fortran或VisualC++语言编程较为繁琐。 而Matlab语言提供了一批功能强大的库函数，在进行矩阵运算等方面特别简洁、高效， 更由于它易于调试，不需要高深的编程技巧和数学知识，可以很方便的解决上述问题。本 文利用Matlab语言实现了直接积分法求积分位移，以及采取线性加速度法求地震反应谱， 说明了它在工程计算与科学研究中的实用价值。 2.用Matlab求积分位移 振动台试验中，由于仪器的误差，加速度纪录的零线会发生跳跃现象，虽然这种变动 对加速度纪录影响很小，但会对积分速度产生一定的影响，而积分位移会产生较大的零线 漂移3(]采用普通方法对加速度纪录直接积分，有时会发现积分位移有明显零线漂移的趋 势。可运用最小二乘法原理，对所得数据进行趋势项排除，但有时处理效果并不理想。本 文采用Matlab编程计算取得较好效果。M文件如下: fid=fopen(a八1122.txtI，r’ %打开外部文件% 蕊映字，1972年12月出生，硕士研究生 工《程力学》增刊 1999年 [a,count]=fscanf(fid,’一f %将外部数据加速度记录读入数组a% status=fclose(fid) 务关闭外部文件毛 deft=0.02; t=delt*(l:length(a)) 令生成一个1-数组a的维数单位增量的 向量七% edot=cumsum(a)*delt; edot=edot-mean(edot); %用Matlab函数cumsum对加速度记 录a积分一次。对其进行校正，求得 积分谏度edot3 epos=cumsum(edot)*delt; epos=epos-mean(epos); 店同理，求积分位移epos% out2=fopen(tyyl’，’w’ 云打开外部文件% outl=fopen(’tyy2’，w out0=fopen(tyy3’，’w’ out=fopen(’tyy4，w) countl=fprintf(out2,’3f\n,t); %将时间、加速度记录、积分速度、积 分位移分别写入外部文件% fprintf(outl,’-f\n’,a); fprintf(out0,’}-f\n,edot); fprintf(out,’tf\n,epos); statusl=fclose(out2) 从程序可以看出，在Madab环境下数组 (向量或矩阵)的操作同数的操作一样简单，