bilinearquadelementrestrictp.m

四节点矩形单元matlab有限元程序

function y=BilinearQuadElementRestrictP(NF,P,X,Y)
%    LinearTriangleElementRestrict的作用是给结点配置载荷；
%    所输入的载荷是节点所受的等效载荷,是已经经过等效和叠加的
NFZ=[];
NLZ=0;
NL=input('  此次配置结点力的结点个数　NL=');
while abs(NL-0)>1.0e-5
    NFP=input('  加结点力的结点编号　NFP[NL]=');
    p=input('  结点力大小　p[NF]=');
    for I=1:NL
        for J=1:NF
            i=(NFP(I)-1)*NF+J;
            P(i)=p(J);
        end
    end
    NFZ=[NFZ,NFP];
    NLZ=NLZ+NL;
    NL=input('  NL=');
end
text(X(NFZ),Y(NFZ),'  NP');
y=P;
hold off;