kings2056 发表于 2007-9-19 11:29

急!!Wilson' θ法求结构振动响应时θ越大是不是积分的精度越高

用Wilson' θ法求结构振动响应时θ越大是不是积分的精度越高啊?急求指导,或是相关内容查哪些书啊?

[ 本帖最后由 mjhzhjg 于 2007-9-19 14:00 编辑 ]

appleseed05 发表于 2007-9-19 18:01

这个应该不是这样的,我记得邹经湘的结构动力学书上讲的是theta〉1.4无条件稳定,一般优化解是1.42左右。

身边没有结构动力学,你可以查一下

wendy袁 发表于 2007-9-21 10:26

李国豪的《工程结构抗震动力学》有提到,wilson-θ在θ>=1.37时,是无条件稳定的,在大多数情况下,取θ=1.4左右可以给出很好的结果。wilson-θ法的精度取决于荷载情况、结构系统的频率、阻尼以及时间步长的大小。

wanyeqing2003 发表于 2007-9-21 10:47

θ的取值与计算稳定有关,好像与计算精度关系不大。

yejet 发表于 2007-9-25 10:05

θ对计算精度影响还是挺大的
随着θ增加,相对误差和绝对误差都会增加

曾经看过一篇文章,作者的分析结果表明θ=0.79时,误差最小

jlde_639 发表于 2007-9-25 10:18

我现在也是在做Wilson-θ法,程序编出来也有问题!十层的结构,就是当质量由五次方改为4次方的时候,结构响应图的幅值就变得很大!好几百的次方!~~查看了程序也没觉得错:
clear;
clc;
m=25000*speye(10,10);
k=1.00E+08*sparse(,...
    ,...
    [2.8426,-2.0833,-2.0833,4.1667,-2.0833,-2.0833,3.75,-1.6667,-1.6667,3.3333,-1.6667,-1.6667,3.3333,-1.6667,...
    -1.6667,2.52,-0.8533,-0.8533,1.7067,-0.8533,-0.8533,1.7067,-0.8533,-0.8533,1.7067,-0.8533,-0.8533,0.8533]);
k=full(k);
f1=k/m;
f=full(f1);
=eig(f);
w=sqrt(d);
h1=0.05;                   %阻尼比
h2=0.05;                   %阻尼比
time=1:1:10
plot(time,x(1));
a1=(2*w(1,1)*w(2,2)*(h1*w(2,2)-h2*w(1,1)))/(w(2,2)^2-w(1,1)^2);
a2=2*(h2*w(2,2)-h1*w(1,1))/(w(2,2)^2-w(1,1)^2);
c=a1*m+a2*k;
c=full(c);
dzbo;                     %地震波的调用
ag1=a(:,2);
xs=70/max(abs(ag1));
ag=ag1*0.01*xs;
ag=ag';
seita=1.4;
dt=0.02;
T=10.0;
N1=ones(1,10);
N=N1';
n=T/dt+1;
k1=k+6*m/(seita^2*dt^2)+3*c/(seita*dt);
u=zeros(10,n);
v=zeros(10,n);
ju=zeros(10,n);
for j=2:n
      G1=-m*N*seita*(ag(j)-ag(j-1));
      G2=m*(6*v(:,j-1)/seita/dt+3*ju(:,j-1));
      G3=c*(3*v(:,j-1)+seita*dt*ju(:,j-1)/2);
      dp=G1+G2+G3;
      du1=inv(k1)*dp;
      du=du1/seita;
      dju1=6.0*du1/(seita^2*dt^2)-6.0*v(:,j-1)/(seita*dt)-3.0*ju(:,j-1);
      dju=dju1/seita;
      du=dt*v(:,j-1)+dt^2*ju(:,j-1)/2.0+dt^2*dju;
      dv=dt*ju(:,j-1)+dt*dju/2.0;
      u(:,j)=u(:,j-1)+du;
      v(:,j)=v(:,j-1)+dv;
      ju(:,j)=-N*ag(j)-inv(m)*k*u(:,j)-inv(m)*c*v(:,j);
end
t=0:0.02:10;
plot(t,u(10,:));
不知哪位高手可帮忙解决下!
谢谢

yejet 发表于 2007-9-26 16:20

本帖最后由 VibInfo 于 2016-5-25 16:28 编辑

原帖由 jlde_639 于 2007-9-25 10:18 发表
我现在也是在做Wilson-θ法,程序编出来也有问题!十层的结构,就是当质量由五次方改为4次方的时候,结构响应图的幅值就变得很大!好几百的次方!~~查看了程序也没觉得错:
clear;
clc;
m=25000*speye(10,10) ...
请把dzbo给出来,否则无法运行程序

jlde_639 发表于 2007-9-27 15:02

回复 #7 yejet 的帖子

您好:
我已在“标题: 【结构振动】 Wilson-θ法中质量次方的改变对结果影响很大!求教为什么?”给出我的dzho了!
能否帮忙调试一下?
页: [1]
查看完整版本: 急!!Wilson' θ法求结构振动响应时θ越大是不是积分的精度越高