clear;clc;formatshorteE1=2*10^5;E2=2*10^4;u=0.3;G12=1*10^4;S=[1/E1,-u/E1,0;-u/E1,1/E2,0;0,0,1/G12];Q=S^-1;k=linspace(0,pi/2,100);o=k;m=cos(o);n=sin(o);T=[m.^2,n.^2,2.*m.*...
clear;clc;
format short e
E1=2*10^5;
E2=2*10^4;
u=0.3;
G12=1*10^4;
S=[1/E1,-u/E1,0;
-u/E1,1/E2,0;
0,0,1/G12];
Q=S^-1;
k=linspace(0,pi/2,100);
o=k;
m=cos(o);
n=sin(o);
T=[m.^2,n.^2,2.*m.*n;
n.^2,m.^2,-2.*m.*n;
-m.*n,m.*n,m.^2-n.^2];
gS=T'*S*T;
gQ=gS.^-1;
A=gQ*0.001;
B=0.5*gQ*((-0.0015)^2-(-0.0005)^2);
D=1/3*gQ*((-0.0015)^3-(-0.0005)^3);
o=-k;
m=cos(o);
n=sin(o);
T=[m.^2,n.^2,2.*m.*n;
n.^2,m.^2,-2.*m.*n;
-m.*n