planeframememberglobalstiffnessmatrix.m

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

function smd=PlaneFrameMemberGlobalStiffnessMatrix(pj,pk,ax,iz,e)
dp=pk-pj;
L=norm(dp);
cx=dp(1)/L;
cy=dp(2)/L;
scm1=e*ax/L;
scm2=4.0*e*iz/L;
scm3=1.5*scm2/L;
scm4=2.0*scm3/L;
c1=scm1*cx*cx+scm4*cy*cy;
c2=(scm1-scm4)*cx*cy;
c3=-scm3*cy;
c4=scm1*cy*cy+scm4*cx*cx;
c5=scm3*cx;
smd=[c1 c2 c3 -c1 -c2 c3;...
    c2 c4 c5 -c2 -c4 -c5;...
   c3 c5 scm2 -c3 -c5 scm2/2;...
   -c1 -c2 -c3 c1 c2 -c3;...
   -c2 -c4 -c5 c2 c4 -c5;...
   c3 c5 scm2/2 -c3 -c5 scm2;];