- 积分
- 0
- 注册时间
- 2010-10-1
- 仿真币
-
- 最后登录
- 1970-1-1
|
syms lambda theta0 theta1 theta2 theta3
syms E1 E2 E3 nu1 nu2 nu3
mu1=E1/(2*(1+nu1))
mu2=E2/(2*(1+nu2))
mu3=E3/(2*(1+nu3))
kapa1=3-4*nu1
kapa2=3-4*nu2
kapa3=3-4*nu3
tao21=mu2/mu1
tao31=mu3/mu1
tao32=mu3/mu2
A=[lambda*exp(i*(lambda-1)*theta0),exp(-i*(lambda-1)*theta0),0,0,-lambda*exp(i*(lambda-1)*theta3),-exp(i*(lambda-1)*theta3),exp(i*(lambda+1)*theta0),0,0,0,-exp(i*(lambda+1)*theta3),0;
exp(i*(lambda-1)*theta0),lambda*exp(-i*(lambda-1)*theta0),0,0,-exp(i*(lambda-1)*theta3),-lambda*exp(-i*(lambda-1)*theta3),0,exp(-i*(lambda+1)*theta0),0,0,0,-exp(-i*(lambda+1)*theta3);
lambda*exp(i*(lambda-1)*theta1),exp(-i*(lambda-1)*theta1),-lambda*exp(i*(lambda-1)*theta1),-exp(-i*(lambda-1)*theta1),0,0,exp(i*(lambda+1)*theta1),0,-exp(i*(lambda+1)*theta1),0,0,0;
exp(i*(lambda-1)*theta1),lambda*exp(-i*(lambda-1)*theta1),-exp(i*(lambda-1)*theta1),-lambda*exp(-i*(lambda-1)*theta1),0,0,0,exp(-i*(lambda+1)*theta1),0,-exp(-i*(lambda+1)*theta1),0,0;
0,0,lambda*exp(i*(lambda-1)*theta2),exp(-i*(lambda-1)*theta2),-lambda*exp(i*(lambda-1)*theta2),-exp(i*(lambda-1)*theta2),0,0,exp(i*(lambda+1)*theta2),0,-exp(i*(lambda+1)*theta2),0;
0,0,exp(i*(lambda-1)*theta2),lambda*exp(-i*(lambda-1)*theta2),-exp(i*(lambda-1)*theta2),-lambda*exp(-i*(lambda-1)*theta2),0,0,0,exp(-i*(lambda+1)*theta2),0,-exp(-i*(lambda+1)*theta2);
tao31*kapa1*exp(i*(lambda-1)*theta0),-tao31*lambda*exp(-i*(lambda-1)*theta0),0,0,-kapa3*exp(i*(lambda-1)*theta3),lambda*exp(-i*(lambda-1)*theta3),0,-tao31*exp(-i*(lambda+1)*theta0),0,0,0,exp(-i*(lambda+1)*theta3);
-tao31*lambda*exp(i*(lambda-1)*theta0),tao31*kapa1*exp(-i*(lambda-1)*theta0),0,0,lambda*exp(i*(lambda-1)*theta3),-kapa3*exp(-i*(lambda-1)*theta3),-tao31*exp(i*(lambda+1)*theta0),0,0,0,exp(i*(lambda+1)*theta3),0;
tao21*kapa1*exp(i*(lambda-1)*theta1),-tao21*lambda*exp(-i*(lambda-1)*theta1),-kapa2*exp(i*(lambda-1)*theta1),lambda*exp(-i*(lambda-1)*theta1),0,0,0,-tao21*exp(-i*(lambda+1)*theta1),0,exp(-i*(lambda+1)*theta1),0,0;
-lambda*tao21*exp(i*(lambda-1)*theta1),tao21*kapa1*exp(-i*(lambda-1)*theta1),lambda*exp(i*(lambda-1)*theta1),-kapa2*exp(-i*(lambda-1)*theta1),0,0,-tao21*exp(i*(lambda+1)*theta1),0,exp(i*(lambda+1)*theta1),0,0,0;
0,0,tao32*kapa2*exp(i*(lambda-1)*theta2),-tao32*lambda*exp(-i*(lambda-1)*theta2),-kapa3*exp(i*(lambda-1)*theta2),lambda*exp(-i*(lambda-1)*theta2),0,0,0,-tao32*exp(-i*(lambda+1)*theta2),0,exp(-i*(lambda+1)*theta2);
0,0,-tao32*lambda*exp(i*(lambda-1)*theta2),tao32*kapa2*exp(-i*(lambda-1)*theta2),lambda*exp(i*(lambda-1)*theta2),-kapa3*exp(-i*(lambda-1)*theta2),0,0,-tao32*exp(i*(lambda+1)*theta2),0,exp(i*(lambda+1)*theta2),0] |
|