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

topspeed21 · 发表于 2006-6-1 10:15

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

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

x

１．在计算消能结构时我采用wilson-θ法，在给定增量微分方程中各项系数矩阵和初值的条件下，用来逐步解方程的。此时［C］不是从结构中组装出来的，而是用原结构阻尼叠加一个附加阻尼．但是当叠加后的总阻尼大到一定程度的时候会出现计算上的错误．（是不是很有可能算发散？） 程序我是用徐赵东老师的《MATLAB语言在建筑抗震工程中的应用》一书上的弹性时程分析修改的。(猜想可能是线性方程组病态的原因?) ２．另外，如果不叠加附加阻尼的话，原结构阻尼c=am+bk　当c同样被认为扩大到一定倍数时，也会出现发散的问题。 ３．不是说θ＞１.３７就无条件收敛么？ 请高手赐教上述问题～～～　有什么办法可以解决? 　万分感谢！

[此贴子已经被cdwxg于2006-6-1 11:52:28编辑过]

coocoo · 发表于 2006-6-15 12:09

你能把那本MATLAB语言在建筑抗震工程中的应用》书发给我看看 如果你有电子版的话

xuebx · 发表于 2006-6-15 13:28

大家踊跃发言啊，知道的坛友多多指点

<SCRIPT language=javascript type=text/javascript> document.write (usercolor('9','topspeed21')); </script>

复制代码

topspeed21。

chuxidahai · 发表于 2006-6-15 19:52

没有问题,可能你的程序抄错了,仔细些

chuxidahai · 发表于 2006-6-15 20:01

我试了一下,阻尼比不能太大,你想想,实际的阻尼应有个范围的,取1.5还可收敛一般取0.05~0.07

bendashao · 发表于 2006-8-24 10:34

振动方程的适用范围就是小阻尼的，所以阻尼不可以很大的，这不是算法的问题，而是方程适用性的问题。

cassett · 发表于 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

cassett · 发表于 2006-11-13 22:06

有那本书上相关matlab程序的能贴出来给我看看吗。谢谢

账号		自动登录	找回密码
密码			我要加入

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

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

我也有这个程序,

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