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

kings2056

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

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

appleseed05

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

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

wendy袁

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

wanyeqing2003

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

yejet

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

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

jlde_639

我现在也是在做Wilson-θ法，程序编出来也有问题！十层的结构，就是当质量由五次方改为4次方的时候，结构响应图的幅值就变得很大！好几百的次方！~~查看了程序也没觉得错:
clear;
clc;
m=25000*speye(10,10);
k=1.00E+08*sparse([1,1,2,2,2,3,3,3,4,4,4,5,5,5,6,6,6,7,7,7,8,8,8,9,9,9,10,10],...
    [1,2,1,2,3,2,3,4,3,4,5,4,5,6,5,6,7,6,7,8,7,8,9,8,9,10,9,10],...
    [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);
[x,d]=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

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

请把dzbo给出来，否则无法运行程序

jlde_639

您好：
我已在“标题: 【结构振动】 Wilson-θ法中质量次方的改变对结果影响很大！求教为什么？”给出我的dzho了！
能否帮忙调试一下？