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

t工程力学》增刊 1999年 ·975· Matl ab在地震工程研究中的应用实例 葛映宇 叶献国 (舍罡工业大学建工系t台肥230009) 提要 lIatlab是解决科学工程计算问题的高效工具．本文具体论述了如何用htlab语言实现求地震 反应谱，以及由加速度记录求积分速度，积分位穆，由此说明lIatlah语言在地震工程研究中具有良好的 实用性，值得在土木工程领域推广使用． 关键词 妇tlab语言．地震反应谱，积分速度，积分位移． 1．前言 c¨等语言比较，Matlab的语法规则简单，更加贴近人的思维方式，易学易用。用№tlab写 程序，犹如在一张演算纸上排列公式和求解问题，因而被称为“演算纸方式的”科学工程 的算法语言‘“。众所周知，F0rtr∞、Basic等语言在土木工程计算编程方面具有广泛的 c++语言编程较为繁琐。 应用，但涉及矩阵理论、数值分析等问题时用F0rtran或Visual 而№tlab语言提供了一批功能强大的库函数，在进行矩阵运算等方面特别简洁、高效， 更由于它易于调试，不需要高深的编程技巧和数学知识，可以很方便的解决上述问题。本 文利用Matlab语言实现了直接积分法求积分位移，以及采取线性加速度法求地震反应谱， 说明了它在工程计算与科学研究中的实用价值。 2．用_atlab求积分位移 振动台试验中，由于仪器的误差，加速度纪录的零线会发生跳跃现象，虽然这种变动 对加速度纪录影响很小，但会对积分速度产生一定的影响，而积分位移会产生较大的零线 漂移例。采用普通方法对加速度纪录直接积分，有时会发现积分位移有明显零线漂移的趋 势。可运用最小二乘法原理，对所得数据进行趋势项排除，但有时处理效果并不理想。本 文采用Ma吐ab编程计算取得较好效果。M文件如下： a：＼1122．txt’，’工’) ％打开外部文件％ fid-f。pen(t 葛映宇，1972年】0片出生．硕士研究生 ·976· 《工程力学》增刊 1999年 ％将外部数据加速度记录读入数组a％ [a，count]=fscanf(fid，’。f。)； 8tatus=fclose(fid) {关闭外部文件{ delt=O．02； t=delt+(1：length(a))‘；{生成一个卜数组a的维数单位增量的 向量t} edot=cumsum(a)+de工t； edot=edot—mean(edot)；舌用Matlab函数cumsum对加速度记 录a积分一次，对其进行校正，求得 g 积分速度ed。t epos=cumsumfedot)+delt； epos=epos—mean(ep。s)； ％同理，求积分位移ep。s％ out2=fopen(’七yyl’，’w’) $打开外部文件％ outl=fopen(‘七yy2‘，‘w’) 。ut0=fopen(。tyy3’，’w’) out=fopen(1tyy4’，’w’) ％将时间、加速度记录、积分速度、积 countl=fprintf(out2，’if＼n’，t)j 分位移分别写入外部文件％ fprin七f(outl，’}f＼n。，a)； fprintf(out0，’岳f＼n‘，edot)； fprintf(out，。}f＼n’，ep05)； statusl；fclose(out2) 从程序可以看出，在MaⅡab环境下数组(向量或矩阵)的操作同数的操作一样简单， 这是因为Ma吐ab语言的基本数据单元是