16 lines
675 B
Matlab
16 lines
675 B
Matlab
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 |