声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2493|回复: 13

[编程技巧] matlab在建筑抗震工程中的应用(帮忙看看)

[复制链接]
发表于 2010-11-16 21:00 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
一个弹性时程分析程序设计例题:某三层钢筋混凝土结构地震波采用200gal el centre波,采样周期为0.02s.用弹性时程分析发求解结构地震反应.程序如下:
clear;
h=[4000,3300,3300]
m=[2.762 2.760 2.300]*1e+3;
k0=[2.485 1.921 1.522]*1e+5;
cn=length(m)  %总层数
%earthquake parameter
ct=1.4  %Wilson-法的值
dt=0.02
load ('dzhbo1')
xs=200/max(abs(dzhbo1))  %调整地震输入加速度幅值200gal
ag=dzhbo1*0.01*xs;         %地震波,单位m/s2

ndzh=400
ag1=ag(1:ndzh)
ag2=ag(2:ndzh+1)
agtao=ct*(ag2-ag1) %initial value
chsh=zero(cn,1)
wyi1=chsh
sdu1=chsh
jsdu1=chsh
wyimt=chsh
sdumt=chsh
jsdumt=chsh
unit=one(cn,1)
m=diag(m)
[ik]=matrixju(k0,cn)
[x,d]=eig(ik,m)
d=sqrt(d)
w=sort(diag(d))
a=2*w(1)*w(2)*(0.05*w(2)-0.07*w(1))/(w(2)^2-w(1)^2);
b=2*(0.07*w(2)-0.05*w(1))/(w(2)^2-w(1)^2);
c0=a*m+b*ik
运行后说红字部分错误?那个地震波不是个.m文件吗?是不是就是地震的数据呢?

评分

1

查看全部评分

回复
分享到:

使用道具 举报

发表于 2010-11-16 21:41 | 显示全部楼层
你把错误的提示信息贴出来
 楼主| 发表于 2010-11-16 21:45 | 显示全部楼层
本帖最后由 499543921 于 2010-11-16 21:49 编辑

我觉得可能是地震波引用的不太对!可是不太会处理.
 楼主| 发表于 2010-11-16 21:47 | 显示全部楼层
回复 2 # qibbxxt 的帖子

??? Attempt to execute SCRIPT load as a function.

Error in ==> load at 12
load ('dzhbo1')
谢谢
发表于 2010-11-16 23:16 | 显示全部楼层
本帖最后由 Happy99 于 2010-11-16 23:21 编辑

dzhbo1是什麽档案? dzhbo1.mat?
"...那个地震波不是个.m文件吗?是不是就是地震的数据呢?"
LZ不知道吗? 程序不是LZ写的?:@)
发表于 2010-11-16 23:48 | 显示全部楼层
ref
常见的程序出错问题整理 (eight), 6F
http://forum.vibunion.com/forum/thread-46001-1-1.html
发表于 2010-11-17 15:47 | 显示全部楼层
兰州都说是个例题 当然不是他自己写的啦
这应该是要导入一个原数据
要不你重命名了 要不还在你的压缩包里面放着  
你搜索一下 找到后 放到work下面呗
 楼主| 发表于 2010-11-17 21:50 | 显示全部楼层
回复 6 # ChaChing 的帖子

拜读中!先谢谢啦!
 楼主| 发表于 2010-11-17 21:53 | 显示全部楼层
回复 7 # moon555sss 的帖子

试了下,不行呢...唉我是直接把那个dzhbo.txt文件剪切到.m文件里的!不行吗这样?是些这样的数据:埃而森特罗波(EI Centro)地震 时间间隔0.02秒,持续时间取30s 适合于2类场地土
EL CENTRO  NS  1940  CALIF.                 2688   0.020 341.700(10F8.3)  
南北方向,  
(高振世等著,建筑结构抗震设计,中国建筑工业出版社,1995年)

-1.400
-10.800
-10.100  
-8.800  
-9.500
-12.000
-14.200
-12.800
-11.000  
-8.500
-8.500
-13.100
-17.600
-19.400
-16.200
-14.400
-10.800  
-8.200  
-4.200  
-6.600
-13.100
-19.000
-19.600  
-6.600   
3.000  
14.100  
-4.900
-12.800
-14.400
-20.300
-26.000
-32.500
-30.600
-17.200
-19.700
-16.300
-16.400  
-6.700   
2.500  
15.000
23.600  
25.200  
33.600  
46.300  
49.200  
41.900  
35.900  
27.100  
23.500  
33.900
41.200  
53.000  
63.900  
73.200  
65.200  
59.900  
40.000  
40.000   
6.300
-51.500
-78.700
-60.300
-48.400
-25.000  
-5.900  
13.400  
30.800  
49.900  
71.000  
99.500
121.900
152.900
144.900
115.500  
93.500  
89.200  
92.600  
83.900  
90.100  
99.300
120.900  
32.800
-147.500
-206.600
-198.900
-203.400
-181.600
-172.500
-175.200
-175.300
-180.500
-163.000
-134.700
-108.700
-78.200
-42.900  
-1.700  
36.000  
78.500
116.400
159.800
196.000
241.200
272.900
303.600
320.000
341.700
282.100
232.400
-119.800
-237.300
-164.000
-186.500
-109.500
-75.300
-17.300  
11.300  
53.300  
89.500
118.600
175.700  
57.600
-263.100
-154.700
-172.900
-101.200
-57.900  
23.700
-67.000
-198.000
-164.100
-168.500
-148.100
-123.100
-100.100
-75.100
-52.300
-27.100  
-4.400  
18.800
-9.500
-43.300
-83.800
-95.100
-71.600
-59.900
-33.400
-10.800  
18.500  
42.000
67.300  
-9.700
-37.200  
-4.000   
1.100  
34.400  
56.500  
 楼主| 发表于 2010-11-17 21:54 | 显示全部楼层
回复 5 # Happy99 的帖子

哦?我查一下.呵呵,新手..
发表于 2010-11-18 00:03 | 显示全部楼层
资料剪切到.m文件, 也是可以, 但要命名为dzhbo矩阵
一般直接使用load即可, 但dzhbo.txt文件的开头注解要删除

还有求助完整格式:出错代码和出错提示
 楼主| 发表于 2010-11-18 17:21 | 显示全部楼层
十分感谢大家的帮忙!终于运行出来啦.我直接把那个地震波数据.txt格式用load把它引用过来的!
dzhbo=load ('elcentero.txt')
呵呵  还是集思广益好解决问题!

评分

1

查看全部评分

发表于 2010-11-19 00:06 | 显示全部楼层
回复 12 # 499543921 的帖子

其实这个很基本, 建议LZ看下基本书并多练习!:@)
发表于 2011-12-22 11:38 | 显示全部楼层
你好,我也正在做这方面的编程,很多问题搞不懂。。。需要向你请教
我的扣扣 383611339
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-13 14:00 , Processed in 0.072998 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表