spacetrusselementstress.m

来自「matlab在有限元中应用的源程序及命令 下载后可直接解压」· M 代码 · 共 21 行

function y = SpaceTrussElementStress(E,L,thetax,thetay,thetaz,u)
%SpaceTrussElementStress   This function returns the element stress
%                          given the modulus of elasticity E, the 
%                          length L, the angles thetax, thetay, 
%                          thetaz (in degrees), and the element 
%                          nodal displacement vector u.
x = thetax * pi/180;
w = thetay * pi/180;
v = thetaz * pi/180;
Cx = cos(x);
Cy = cos(w);
Cz = cos(v);
y = E/L*[-Cx -Cy -Cz Cx Cy Cz]*u;