FEM-Course-Matlab/1.matlab平面纯弯梁有限元编程/悬臂梁/BeamElementStiffness.m

15 lines
483 B
Mathematica
Raw Permalink Normal View History

2024-01-28 16:46:36 +00:00
function y = BeamElementStiffness(EI,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 = EI/(L*L*L)*[12 6*L -12 6*L ; 6*L 4*L*L -6*L 2*L*L ;
-12 -6*L 12 -6*L ; 6*L 2*L*L -6*L 4*L*L];