急!!Wilson' θ法求结构振动响应时θ越大是不是积分的精度越高
用Wilson' θ法求结构振动响应时θ越大是不是积分的精度越高啊?急求指导,或是相关内容查哪些书啊?[ 本帖最后由 mjhzhjg 于 2007-9-19 14:00 编辑 ] 这个应该不是这样的，我记得邹经湘的结构动力学书上讲的是theta〉1.4无条件稳定，一般优化解是1.42左右。
身边没有结构动力学，你可以查一下 李国豪的《工程结构抗震动力学》有提到，wilson-θ在θ>=1.37时，是无条件稳定的，在大多数情况下，取θ=1.4左右可以给出很好的结果。wilson-θ法的精度取决于荷载情况、结构系统的频率、阻尼以及时间步长的大小。 θ的取值与计算稳定有关，好像与计算精度关系不大。 θ对计算精度影响还是挺大的
随着θ增加，相对误差和绝对误差都会增加
曾经看过一篇文章，作者的分析结果表明θ=0.79时，误差最小 我现在也是在做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,:));
不知哪位高手可帮忙解决下！
谢谢 本帖最后由 VibInfo 于 2016-5-25 16:28 编辑
原帖由 jlde_639 于 2007-9-25 10:18 发表
我现在也是在做Wilson-θ法，程序编出来也有问题！十层的结构，就是当质量由五次方改为4次方的时候，结构响应图的幅值就变得很大！好几百的次方！~~查看了程序也没觉得错:
clear;
clc;
m=25000*speye(10,10) ...
请把dzbo给出来，否则无法运行程序
回复 #7 yejet 的帖子
您好：我已在“标题: 【结构振动】 Wilson-θ法中质量次方的改变对结果影响很大！求教为什么？”给出我的dzho了！
能否帮忙调试一下？
页:
[1]