声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1956|回复: 8

[综合讨论] 反推隔震结构的水平刚度的程序

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

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

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

x
本帖最后由 土木年华 于 2010-10-22 16:41 编辑

各位大侠好,
    小弟请教一个问题, 现在正在做一个有意义的问题,就是在正常的环境激励下隔震层的水平刚度怎么理解。前几天我们做了对一隔震房屋做了大地脉动测试,识别出了频率,想与结构的理论振型频率进行比较,如
│sw1-w1│≤0.001且│sw2-w2│≤0.001而反推出隔震层的刚度Kd,想编个MATLAB小程序,请教各位大侠,给指导一下该如何编程,  谢谢大家!  下面是各层的质量,刚度,实测频率,哪位大侠给改改一改?

sw=[1;1;1;]      
sw1=sw(1,:)
sw2=sw(2,:)          % 实测识别前两阶频率

cn=4;
m0=[2762,2762,2300,2300];
k0=[kd;2.485;1.921;1.511]*1e+5;      % 一二三层及隔震层质量,刚度
M=diag(m0);             %生成结构的质量和刚度矩阵
[K]=matrixju(k0,cn);    %刚度矩阵的聚合

[z,d]=eig(K,M);        %计算结构频率和振型
w=diag(sqrt(d))      %理论振型频率

w1=w(1,:)
w2=w(2,:)
回复
分享到:

使用道具 举报

发表于 2010-10-22 21:01 | 显示全部楼层
实在水平专业有限, 看不清楚LZ要的
建议LZ说清楚或简化些
发表于 2010-10-23 08:30 | 显示全部楼层
具体专业问题不太懂
不过如果你这里的kd和最后输出的固有频率之间的关系连续的话
那这个问题不难,是一个简单的参数识别问题
最笨的办法可以采用二分法
如果要做的好一点可以在优化工具箱中找一个合适的函数进行识别
具体的你找相关算法看看吧
 楼主| 发表于 2010-10-23 09:02 | 显示全部楼层
回复 ChaChing 的帖子

我的意思是用我测出来的两阶频率与我用eig函数算出来的频率比较,使他们两个的差值小于某一很小的数,就是想借此看一看隔震层的刚度kd,请问老哥这个程序该怎么弄啊  ,指导一下,哈 谢谢
 楼主| 发表于 2010-10-23 09:04 | 显示全部楼层
回复 happy 的帖子


我的意思是用我测出来的两阶频率与我用eig函数算出来的频率比较,使他们两个的差值小于某一很小的数,就是想借此看一看隔震层的刚度kd ,哈 你就别想这个问题,就给我指导一下这个程序该如何编啊 ,哈 我的MATLAB编程学的不是很好,谢谢你
发表于 2010-10-23 11:58 | 显示全部楼层
简单介绍一下采用二分法的思路吧,程序还是需要你自己去完成的

假设你的kd取值范围为[kd(0),kd(n)]
1. 将该范围分成n份,这样就形成了[kd(0),kd(1)],[kd(1),kd(2)],......,[kd(n-1),kd(n)]共n个区间
2. 计算下列节点位置的的固有频率
[kd(0),kd(1),kd(2),......,kd(n-1),kd(n)]
得到固有频率
[w1(0),w1(1),w1(2),......,w1(n-1),w1(n)]
[w2(0),w2(1),w2(2),......,w2(n-1),w2(n)]
3. 判断解区间
判断寻找满足(w1(i+1)-sw1)(w1(i)-sw1)<0 且 (w1(i+1)-sw1)(w1(i)-sw1)<0的情况
找出对应的存在解得区间为[kd(i),kd(i+1)]
4. 在[kd(i),kd(i+1)]的区间内使用二分法找出你所需要的解
二分法程序参考http://forum.vibunion.com/thread-96621-1-1.html
所不同的是你这里需要同时满足两个条件

该方法整体思路比较简单清晰,但是问题比较多,比如找不到解、效率低下等
如果取得较好的效果还是用优化工具箱里的函数,你这个问题可以表述为多目标优化问题
比如采用fgoalattain函数等
发表于 2010-10-24 01:36 | 显示全部楼层
回复 土木年华 的帖子

还是不清楚LZ要的!:@L
LZ是要做校准(tuning)吗?:@)
 楼主| 发表于 2010-10-24 09:41 | 显示全部楼层
回复 ChaChing 的帖子

就是要做一个迭代, 每次变换Kd一个值求出振型频率w1,w2,然后再与SW1,SW1,分别比较,当两者之差小于某一直时结束
发表于 2010-10-24 10:34 | 显示全部楼层
回复 土木年华 的帖子

help for/break/while/switch ?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-13 16:44 , Processed in 0.072339 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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