bilinearquadelementrestictu.m

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

function y=BilinearQuadElementRestictU(NF,U,X,Y)
NUZ=[];
NRZ=0;
NR=input('  此次配置位移的结点个数　NR=');
while abs(NR-0)>1.0e-5
    NU=input('  配置结点位移的结点编号　NU[NR]=');
    u=input('  结点位移大小　U[NF]=');
    for I=1:NR
        for J=1:NF
            i=(NU(I)-1)*NF+J;
            U(i)=u(J);
        end
    end
    NUZ=[NUZ,NU];
    NRZ=NRZ+NR;
    NR=input('  NR=');
end;
text(X(NUZ),Y(NUZ),'  NU');
y=U;