求解损伤简支梁振型问题

huangying22

>>l=2;
>> x=1;
>>  ch=0.1;
>> ku=0.1;
>>  h=0.05;
>> E=2*10^11;
>> I=ch*ku^3/12;
>>  m=h/ku;
>> J=1.8624*m^2-3.95*m^3+16.37*m^4-37.226*m^5+76.81*m^6-126.9*m^7+172*m^8-143.97*m^9+66.56*m^10;
>> c=(5.346*ku*J)/(E*I);
>> r=1/c;
>> v=7800;
>> p=(pi)^2*((E*I)/(v*l^4))^0.5;
>> k=(p^2*v*(ch*ku)/E/I)^0.25;
>> a=sym('[cosh(k*x),sinh(k*x),cos(k*x), sin(k*x),-cosh(k*x),-sinh(k*x),-cos(k*x),-sin(k*x);k^2*cosh(x),k^2*sinh(x),-k^2*cos(x),k^2*sin(x),sin(k*x);-k^2*cosh(x),-k^2*sinh(x),k^2*cos(x),k^2*sin(x);k^3*sh(x),k^3*cosh(x),k^3*sin(x), -k^3*cos(x),-k^3*sh(x),-k^3*cosh(x),-k^3*sin(x),k^3*cos(x);k*sinh(x)+E*I*k^2*cosh(x)/r,k*cosh(x)+E*I*k^2*sinh(x)/r,-k*sin(x)-E*I*k^2*cos(x)/r,k*cos(x)-E*I*k^2*cosh(l)/r,k*sinh(k*x),k*cosh(k*x),-k*sin(k*x),k*cos(k*x);1,0,1,0,0,0,0,0;0,1,0,1,0,0,0,0;0,0,0,0,cosh(k*l),sinh(k*l),cos(k*l),sin(k*l);0,0,0,0,k*sinh(l),k*cosh(l),-k*sin(l), k*cos(l); ]');
>> b=[0;0;0;0;0;0;0;0;];
>> jacobi(a,b,[0;0;0;0;0;0;0;0;0;0])
??? Error using ==> sym/maple
Error, invalid types in sum
Error in ==> D:\MATLAB6p5\toolbox\symbolic\@sym\minus.m
On line 22  ==> X = maple(A,'-',B);
Error in ==> D:\MATLAB6p5\toolbox\symbolic\@sym\tril.m
On line 36  ==> Y = X - triu(X,offset+1);
Error in ==> d:\MATLAB6p5\work\jacobi.m
On line 4  ==> L=-tril(a,-1);
帮忙看看这个程序错在哪里？怎么修改？谢谢

sigma665

a=sym('[cosh(k*x),sinh(k*x),cos(k*x), sin(k*x),-cosh(k*x),-sinh(k*x),-cos(k*x),-sin(k*x);k^2*cosh(x),k^2*sinh(x),-k^2*cos(x),k^2*sin(x),sin(k*x);-k^2*cosh(x),-k^2*sinh(x),k^2*cos(x),k^2*sin(x);k^3*sh(x),k^3*cosh(x),k^3*sin(x), -k^3*cos(x),-k^3*sh(x),-k^3*cosh(x),-k^3*sin(x),k^3*cos(x);k*sinh(x)+E*I*k^2*cosh(x)/r,k*cosh(x)+E*I*k^2*sinh(x)/r,-k*sin(x)-E*I*k^2*cos(x)/r,k*cos(x)-E*I*k^2*cosh(l)/r,k*sinh(k*x),k*cosh(k*x),-k*sin(k*x),k*cos(k*x);1,0,1,0,0,0,0,0;0,1,0,1,0,0,0,0;0,0,0,0,cosh(k*l),sinh(k*l),cos(k*l),sin(k*l);0,0,0,0,k*sinh(l),k*cosh(l),-k*sin(l), k*cos(l); ]');

这里，k x 都已经是值了，还定义成符号干什么
要么一开始就定义成符号

另外，建议看下sym的用法

[ 本帖最后由 sigma665 于 2008-5-27 15:46 编辑 ]

huangying22

谢谢你了！
那么应该是先定义成符号，再赋值。