是别人的程序,一个求悬臂板响应的有限元程序,坐标系建立如图所示
在组装总矩阵的时候,有这样一段我不太明白
for I=1:4
if (YE>=0)
Ksu(5*(ELE(M1,I)-1)+1:5*ELE(M1,I))=Ksu(5*(ELE(M1,I)-1)+1:5*ELE(M1,I))+Kuv((I-1)*5+1:I*5);
else
Ksu(5*(ELE(M1,I)-1)+1:5*ELE(M1,I))=Ksu(5*(ELE(M1,I)-1)+1:5*ELE(M1,I))-Kuv((I-1)*5+1:I*5);
end
end
其中YE为节点的坐标,Ksu为整体的节点载荷,Kuv为单元节点载荷,我不明白的就是为什么YE<0时组装单元节点载荷时用减号