FEM-Course-Matlab/10.四节点八节点四边形单元悬臂梁的Matlab有限元编程/8_nodes/UniLoad.m

function Pe=UniLoad(ie,N_ID_p1,q0,Nodes,Elements)     % 3<EFBFBD><EFBFBD><EFBFBD>ӳ<EFBFBD><EFBFBD>򣬷ֲ<EFBFBD><EFBFBD><EFBFBD><EFBFBD>ص<EFBFBD>ЧΪ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>

% q0=1000;%1kN/m
k=-0.625e-3;                            % <EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ֵ N/mm
s = [-sqrt(1/3) sqrt(1/3)];                 % 2*2 <EFBFBD><EFBFBD>˹<EFBFBD><EFBFBD><EFBFBD>ֵ<EFBFBD>
ENodes = N_ID_p1(ie,:);                    %<EFBFBD><EFBFBD>ȡ<EFBFBD><EFBFBD>Ԫ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
Pe=zeros(16,1);                          %<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ʱ<EFBFBD><EFBFBD>Ԫ<EFBFBD>ڵ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
x1=Nodes(ENodes(1),1);
x6=Nodes(ENodes(4),1);
L16=abs(x6-x1);                          %<EFBFBD><EFBFBD>Ԫ<EFBFBD><EFBFBD><EFBFBD><EFBFBD>
for i=1:2                                 %<EFBFBD><EFBFBD><EFBFBD>ڸ<EFBFBD>˹<EFBFBD><EFBFBD><EFBFBD>ֵ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ѭ<EFBFBD><EFBFBD>
    N_q=ShapeFun(s(i),1);                   % 4<EFBFBD><EFBFBD><EFBFBD>ӳ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>ShapeFun(s(i),1)
    q_x=q0;
    Pe=Pe+N_q'*q_x*[0;L16/2];               %<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>κ<EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD><EFBFBD>
end
end