声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2101|回复: 7

[综合讨论] 地震过程中结构的动态仿真问题

[复制链接]
发表于 2007-12-7 20:59 | 显示全部楼层 |阅读模式

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

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

x
qdzb
k=[1,1,1,1,1]*4.2e+7;
m0=[1,1,1,1,1]*2.0e+5;
g1=1.5e+7;
g2=2.01e+7;
nn=2;
cn=5;
as=3e-2;
t=1.3e-2;
O=zeros(cn);
I=eye(cn);
k0=z(k cn);
m=diag(mo);
[x,d]=eig(k0,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);
co=a*m+b*k0;
nd=[1,1,1,1,1];
kd1=nn*g1*as/t;                                                                                                                                                                                       
cd1=nn*g2*as/w(1)/t;
kd=matrixju(kd1*nd,cn);
cd=matrixju(cd1*nd,cn);
A=[O I;-inv(m)*k0-inv(m)*c0];
B=[O;I];
R=inv(m)*[kd cd];
II=diag(I);
Ca=zeros(cn,2*cn);
Ca(cn,2*cn)=1;
Cw=zeros(cn,2*cn);
Cw(cn,cn)=1;
D=diag(O);

[kcju]=matrixju(korc,cn)
kcju=zeros(cn);
for i=1:cn-1
     kcju(i,i)=korc(i)+korc(i+1);
     kcju(i,i+1)=-korc(i+1);
     kcju(i+1,i)=-korc(i+1);
end
kcju(cn,cn)=korc(cn);
qdzb怎么定义呀,怎样进行dzb模块的定义?请各位帮忙。我这个输进去,运行就是出错。

[ 本帖最后由 eight 于 2007-12-10 15:53 编辑 ]

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2007-12-8 10:51 | 显示全部楼层
原帖由 冬哲 于 2007-12-7 20:59 发表
qdzb
k=[1,1,1,1,1]*4.2e+7;
m0=[1,1,1,1,1]*2.0e+5;
g1=1.5e+7;
g2=2.01e+7;
nn=2;
cn=5;
as=3e-2;
t=1.3e-2;
O=zeros(cn);
I=eye(cn);
k0=z(k cn);
m=diag(mo);
[x,d]=eig(k0,m);
d=sqrt(d);
...

不要拿别人的残缺程序让我们修改。

这种事情好像从程序的出处找答案比较有效

[ 本帖最后由 eight 于 2007-12-10 15:52 编辑 ]
 楼主| 发表于 2007-12-9 08:33 | 显示全部楼层

你好

我也不知道那错了qdzb是EI-centro波,现在明白了,但是k0=z(k,cn),说z没有定义,怎么定义呀?
发表于 2007-12-10 15:53 | 显示全部楼层
原帖由 冬哲 于 2007-12-7 20:59 发表
qdzb
k=[1,1,1,1,1]*4.2e+7;
m0=[1,1,1,1,1]*2.0e+5;
g1=1.5e+7;
g2=2.01e+7;
nn=2;
cn=5;
as=3e-2;
t=1.3e-2;
O=zeros(cn);
I=eye(cn);
k0=z(k cn);
m=diag(mo);
[x,d]=eig(k0,m);
d=sqrt(d);
...

楼主把我们当神仙了,请用相反的角度来思考问题
发表于 2010-10-15 10:23 | 显示全部楼层
回复 冬哲 的帖子

想问问以下什么意思
kd=matrixju(kd1*nd,cn);
cd=matrixju(cd1*nd,cn);
matrixju这个啥意思?
发表于 2010-10-18 07:06 | 显示全部楼层
冬哲 发表于 2007-12-9 08:33
我也不知道那错了qdzb是EI-centro波,现在明白了,但是k0=z(k,cn),说z没有定义,怎么定义呀?

从程序逻辑上来看,K0应该是程序的结构的等效刚度,所以z应该是等效刚度计算函数
发表于 2010-10-18 07:10 | 显示全部楼层
土木年华 发表于 2010-10-15 10:23
回复 冬哲 的帖子

想问问以下什么意思

matrixju好像类似于一个刚度矩阵整合方面的代码
没有任何注释,且不完成的程序还是别读了,费神而且大多没有意义

评分

1

查看全部评分

发表于 2012-5-16 15:35 | 显示全部楼层
回复 7 # happy 的帖子

想问一下,qdzb的输入问题,simulink输入地震波怎么输入啊,我的simulink模型怎么也不通,说是dzb的输入有问题
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-5 18:40 , Processed in 0.076047 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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