19 lines
653 B
Mathematica
19 lines
653 B
Mathematica
|
function y = BeamElementStiffness(EI,GA,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.
|
|||
|
|
k=6/5;%<EFBFBD><EFBFBD><EFBFBD>ν<EFBFBD><EFBFBD>棬Բ<EFBFBD><EFBFBD>10/9
|
|||
|
|
kb=EI/L*[0 0 0 0;0 1 0 -1;0 0 0 0;0 -1 0 1];
|
|||
|
|
ks=GA/(6*k*L)*[6 3*L -6 3*L;3*L 2*L^2 -3*L L^2;-6 -3*L 6 -3*L;3*L L^2 -3*L 2*L^2];
|
|||
|
|
y=kb+ks;
|
|||
|
|
% 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];
|
|||
|
|
|
|||
|
|
|
|||
|
|
|
|||
|
|
|
|||
|
|
|