Wilson-θ法中质量次方的改变对结果影响很大！求教为什么？
我现在在做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,:));
不知哪位高手可帮忙解决下！急！
谢谢！！！
[ 本帖最后由 jlde_639 于 2007-9-25 10:39 编辑 ] 各位高手谁有Wilson-θ法的程序，可否不吝赐予？！
我的邮箱是：jlde_639@163.com
在写论文中，要用到它，所以挺急的，谢谢各位了！ 本帖最后由 VibInfo 于 2016-5-25 16:26 编辑
原帖由 jlde_639 于 2007-9-25 11:33 发表
各位高手谁有Wilson-θ法的程序，可否不吝赐予？！
我的邮箱是：jlde_639@163.com
在写论文中，要用到它，所以挺急的，谢谢各位了！
http://www.chinavib.com/forum/viewthread.php?tid=34732 由于没有你程序中dzbo的数据，所以没办法帮你调试，你用上面连接中的程序是试试看吧 本帖最后由 VibInfo 于 2016-5-25 16:26 编辑
原帖由 yejet 于 2007-9-26 16:07 发表
由于没有你程序中dzbo的数据，所以没办法帮你调试，你用上面连接中的程序是试试看吧
不好意思，当时心情比较急没有把dzho的数据上传。这个dzho的数据其实就是埃而森特罗波南北向10秒的数据！我现在上传它的数据，到时把txt扩展名改回.m就行咯。那能不能再帮我调试下呢？在此万分感谢！！！
页:
[1]