声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2259|回复: 7

[综合讨论] [求助]线弹性结构wilson-θ法的小问题!

[复制链接]
发表于 2006-6-1 10:15 | 显示全部楼层 |阅读模式

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

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

x
1.在计算消能结构时我采用wilson-θ法,在给定增量微分方程中各项系数矩阵和初值的条件下,用来逐步解方程的。此时[C]不是从结构中组装出来的,而是用原结构阻尼叠加一个附加阻尼.但是当叠加后的总阻尼大到一定程度的时候会出现计算上的错误.(是不是很有可能算发散?)<br>程序我是用徐赵东老师的《MATLAB语言在建筑抗震工程中的应用 》一书上的弹性时程分析修改的。(猜想可能是线性方程组病态的原因?)<br>2.另外,如果不叠加附加阻尼的话,原结构阻尼c=am+bk 当c同样被认为扩大到一定倍数时,也会出现发散的问题。<br>3.不是说θ>1.37就无条件收敛么?<br>请高手赐教上述问题~~~ 有什么办法可以解决?<br>                                                                            万分感谢!
[此贴子已经被cdwxg于2006-6-1 11:52:28编辑过]

回复
分享到:

使用道具 举报

发表于 2006-6-15 12:09 | 显示全部楼层
<P>你能把那本MATLAB语言在建筑抗震工程中的应用 》书发给我看看<BR>如果你有电子版的话</P>
发表于 2006-6-15 13:28 | 显示全部楼层
大家踊跃发言啊,知道的坛友多多指点
  1. &lt;SCRIPT language=javascript type=text/javascript&gt;<br>document.write (usercolor('9','topspeed21'));<br>&lt;/script&gt;
复制代码

<FONT face=Verdana color=#61b713><B>topspeed21。</B></FONT>
发表于 2006-6-15 19:52 | 显示全部楼层

我也有这个程序,

没有问题,可能你的程序抄错了,仔细些
发表于 2006-6-15 20:01 | 显示全部楼层

kxi1=0.05;kxi2=0.07; %第一振型和第二振型阻尼比以求reily阻尼

<P>我试了一下,阻尼比不能太大,你想想,实际的阻尼应有个范围的,取1.5还可收敛一般取0.05~0.07</P>
发表于 2006-8-24 10:34 | 显示全部楼层
振动方程的适用范围就是小阻尼的,所以阻尼不可以很大的,这不是算法的问题,而是方程适用性的问题。
发表于 2006-11-13 22:05 | 显示全部楼层
我有点不明白,我自己编的wilson-theta法,因为没有《MATLAB语言在建筑抗震工程中的应用 》。theta=1.4   。算平面框架杆系模型的时候,出现响应无限大。阻尼采用c=am+bk ,阻尼比为0.05和0.07,而算层模型时是能够得到收敛的结果,请高手指点。结构为普通的三层平面框架。用有限元方法求的杆系刚度矩阵和一致质量矩阵。梁柱尺寸为柱A=0.25,I=0.0052,l=3,梁A=0.15,I=0.003125,l=4
发表于 2006-11-13 22:06 | 显示全部楼层
有那本书上相关matlab程序的能贴出来给我看看吗。谢谢
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-12 15:27 , Processed in 0.061059 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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