beamelementstiffness3.m

来自「有限元程序的matlab工具包」· M 代码 · 共 16 行

function y = BeamElementStiffness3(E,A,I,L)
%BeamElementStiffness   This function returns the element 
%                       stiffness matrix for a beam   
%                       element with modulus of elasticity E,  
%                       moment of inertia I, and length L.
%                       The size of the element stiffness 
%                       matrix is 4 x 4.
y = E/(L*L*L)*[A*L*L 0 0 -A*L*L 0 0 ; 0 12*I 6*I*L 0 -12*I 6*I*L ;
   0 6*I*L 4*I*L*L 0 -6*I*L 2*I*L*L ; -A*L*L 0 0 A*L*L 0 0 ;
   0 -12*I -6*I*L 0 12*I -6*I*L ; 0 6*I*L 2*I*L*L 0 -6*I*L 4*I*L*L];