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;