[求助]线弹性结构wilson-θ法的小问题!
１．在计算消能结构时我采用wilson-θ法，在给定增量微分方程中各项系数矩阵和初值的条件下，用来逐步解方程的。此时［C］不是从结构中组装出来的，而是用原结构阻尼叠加一个附加阻尼．但是当叠加后的总阻尼大到一定程度的时候会出现计算上的错误．（是不是很有可能算发散？）<br>程序我是用徐赵东老师的《MATLAB语言在建筑抗震工程中的应用 》一书上的弹性时程分析修改的。(猜想可能是线性方程组病态的原因?)<br>２．另外，如果不叠加附加阻尼的话，原结构阻尼c=am+bk 当c同样被认为扩大到一定倍数时，也会出现发散的问题。<br>３．不是说θ＞１.３７就无条件收敛么？<br>请高手赐教上述问题～～～ 有什么办法可以解决?<br> 万分感谢！[此贴子已经被cdwxg于2006-6-1 11:52:28编辑过]
<P>你能把那本MATLAB语言在建筑抗震工程中的应用 》书发给我看看<BR>如果你有电子版的话</P> 大家踊跃发言啊，知道的坛友多多指点
<SCRIPT language=javascript type=text/javascript><br>document.write (usercolor('9','topspeed21'));<br></script>
<FONT face=Verdana color=#61b713><B>topspeed21。</B></FONT>
我也有这个程序,
没有问题,可能你的程序抄错了,仔细些kxi1=0.05;kxi2=0.07; %第一振型和第二振型阻尼比以求reily阻尼
<P>我试了一下,阻尼比不能太大,你想想,实际的阻尼应有个范围的,取1.5还可收敛一般取0.05~0.07</P> 振动方程的适用范围就是小阻尼的，所以阻尼不可以很大的，这不是算法的问题，而是方程适用性的问题。 我有点不明白，我自己编的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 有那本书上相关matlab程序的能贴出来给我看看吗。谢谢
页:
[1]