planeframememberlocalstiffnessmatrix.m

平面刚架（带弹簧支座）有限元程序

function sm=PlaneFrameMemberLocalStiffnessMatrix(pj,pk,ax,iz,e)
dp=pk-pj;
L=norm(dp);
c=e*ax/L;
b1=12*e*iz/L^3;
b2=b1/2*L;
b3=b2/1.5*L;
b4=b3/2;
sm=[c 0 0 -c 0 0 ;...
    0 b1 b2 0 -b1 b2;...
    0 b2 b3 0 -b2 b4;...
    -c 0 0 c 0 0 ;...
    0 -b1 -b2 0 b1 -b2;...
    0 b2 b4 0 -b2 b3;];