1.m

ansys相关的一些简单介绍以及样例介绍

fai=0;E=2e11;A=0.005;I=0.005;L=5
C12=cos(fai);S12=sin(fai);
T12=[C12,S12,0,0,0,0;-S12,C12,0,0,0,0;0,0,1,0,0,0;0,0,0,C12,S12,0;0,0,0,-S12,C12,0;0,0,0,0,0,1]
ke_12=[E*A/L,0,0,-E*A/L,0,0;0,12*E*I/(L^3),6*E*I/(L^2),0,-12*E*I/(L^3),6*E*I/(L^2);0,6*E*I/(L^2),4*E*I/L,0,-6*E*I/(L^2),2*E*I/L;-E*A/L,0,0,E*A/L,0,0;0,-12*E*I/(L^3),-6*E*I/(L^2),0,12*E*I/(L^3),-6*E*I/(L^2);0,6*E*I/(L^2),2*E*I/L,0,-6*E*I/(L^2),4*E*I/L]
ke12=conj(T12)*ke_12*T12


fai=-pi/3;L=6
C23=cos(fai);S23=sin(fai);
T23=[C23,S23,0,0,0,0;-S23,C23,0,0,0,0;0,0,1,0,0,0;0,0,0,C23,S23,0;0,0,0,-S23,C23,0;0,0,0,0,0,1]
ke_23=[E*A/L,0,0,-E*A/L,0,0;0,12*E*I/(L^3),6*E*I/(L^2),0,-12*E*I/(L^3),6*E*I/(L^2);0,6*E*I/(L^2),4*E*I/L,0,-6*E*I/(L^2),2*E*I/L;-E*A/L,0,0,E*A/L,0,0;0,-12*E*I/(L^3),-6*E*I/(L^2),0,12*E*I/(L^3),-6*E*I/(L^2);0,6*E*I/(L^2),2*E*I/L,0,-6*E*I/(L^2),4*E*I/L]
ke23=conj(T23)*ke_23*T23


R=[ke12 zeros(6,3);zeros(3,9)]+[zeros(3,9);zeros(6,3) ke23]

f=inv(R)*[0;0;500;0;100;0;0;0;1200]