planeframememberconcentratedload.m

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

function aml=PlaneFrameMemberConcentratedLoad(pj,pk,as,a)
dp=pk-pj;
L=norm(dp);
cx=dp(1)/L;
cy=dp(2)/L;
b=L-a;
R=[cx cy 0;...
    -cy cx 0;...
    0 0 1;];
TML=[-b/L 0 0;...
    0 -b^2*(3*a+b)/L^3 6*a*b/L^3;...
    0 -a*b^2/L^2 b*(2*a-b)/L^2;...
    -a/L 0 0;...
    0 -a^2*(a+3*b)/L^3 -6*a*b/L^3;...
    0 a^2*b/L^2 a*(2*b-a)/L^2;];
aml=TML*R*as;