19 lines
653 B
Matlab
19 lines
653 B
Matlab
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;%¾ØÐνØÃ棬ԲÐÎ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];
|
|
|
|
|
|
|
|
|
|
|