function Pe=UniLoad(ie,N_ID_p1,q0,Nodes,Elements)     % 3级子程序，分布荷载等效为结点荷载

% q0=1000;%1kN/m
k=-0.625e-3;                            % 均布荷载值 N/mm
s = [-sqrt(1/3) sqrt(1/3)];                 % 2*2 高斯积分点
ENodes = N_ID_p1(ie,:);                    %获取单元结点号
Pe=zeros(16,1);                          %生成临时单元节点力零列向量
x1=Nodes(ENodes(1),1);
x6=Nodes(ENodes(4),1);
L16=abs(x6-x1);                          %单元长度
for i=1:2                                 %用于高斯积分的求和循环
    N_q=ShapeFun(s(i),1);                   % 4级子程序：ShapeFun(s(i),1)
    q_x=q0;
    Pe=Pe+N_q'*q_x*[0;L16/2];               %均布形荷载梯形面积
end
end