efix.m

是个好东西

function FF=efix(infe)
%---计算单元固端力
global node element material K bc bc_number nfe_number nfe  node_number element_number %（定义全局变量）
for j=1:1:6
    FF(j)=0.0;
    ie=nfe(infe,1);
    mm=element(ie,3);
    E=material(mm,1);
    I=material(mm,3);
    xi = node( element( ie, 1 ), 1 ) ;  %(两端座标)
    yi = node( element( ie, 1 ), 2 ) ;  %(两端座标)
    xj = node( element( ie, 2 ), 1 ) ;  %(两端座标)
    yj = node( element( ie, 2 ), 2 ) ;  %(两端座标)
    SL = ( (xj-xi)^2 + (yj-yi)^2)^(1/2) ;  %(单元长度)
    IND=nfe(infe,2);
    A=nfe(infe,3);
    B=nfe(infe,4);
    P1=nfe(infe,5);
    P2=nfe(infe,6);
    B1=SL-(2*B+A)/3;
    B2=B-A;
    B3=(A+B)/2;
    C1=SL-(2*B+A)/3;
    C2=B2;
    C3=(2*B+A)/3;
    D1=B^3-A^3;
    D2=B^2-A^2;
    switch IND
        case 1
            FF(2)=-P1*(SL-A)^2*(1+2*A/SL)/SL^2;
            FF(3)=P1*A*(SL-A)^2/SL^2;
            FF(5)=-P1-FF(2);
            FF(6)=-P1*A^2*(SL-A)/SL^2;
        case 2
            FF(2)=-P1*B2*(12*B1^2*SL-8*B1^3+B2^2*SL-2*B1*B2^2)/(4*SL^3);
            FF(3)=P1*B2*(12*B3*B1^2-3*B1*B2^2+B2^2*SL)/12/SL^2;
            FF(5)=-P1*B2-FF(2);
            FF(6)=-P1*B2*(12*B3^2*B1+3*B1*B2^2-2*B2^2*SL)/12/SL^2;
        case 3
            FF(2)=-P2*C2*(18*C1^2*SL-12*C1^3+C2^2*SL-2*C1*C2^2-4*C2^3/45)/12/SL^3;
            FF(3)=P2*C2*(18*C3*C1^2-3*C1*C2^2+C2^2*SL-2*C2^3/15)/36/SL^2;
            FF(5)=-5*P2*C2-FF(2);
            FF(6)=-P2*C2*(18*C3^2*C1-3*C1*C2^2-2*C2^2*SL+2*C2^3/15)/36/SL^2;
        case 4
            FF(2)=-6*P1*A*(SL-A)/SL^3;
            FF(3)=P1*(SL-A)*(3*A-SL)/SL^2;
            FF(5)=-FF(2);
            FF(6)=P1*A*(2*SL-3*A)/SL^2;
        case 5
            FF(2)=-P1*(3*SL*D2-2*D1)/SL^3;
            FF(3)=P1*(2*D2+(B-A)*SL-D1/SL)/SL;
            FF(5)=-FF(2);
            FF(6)=P1*(D2-D1/SL)/SL;
        case 6
            FF(1)=-P1*(1-A/SL);
            FF(4)=-P1*A/S;
        case 7
            FF(1)=-P1*(B-A)*(1-(B+A)/(2*SL));
            FF(4)=-P1*D2/2/SL;
        case 8
            FF(3)=-A*(P1-P2)*E(ie)*I(ie)/B;
            FF(6)=-FF(3);
        otherwise
            disp('error')
    end
end