*n+2, 2*i: 2*i+2] += Ke[6:8,0:2]
K[2*n: 2*n+2, 2*j: 2*j+2] += Ke[6:8,2:4]
K[2*n: 2*n+...2, 2*m: 2*m+2] += Ke[6:8,4:6]
K[2*n: 2*n+2, 2*n: 2*n+2] += Ke[6:8,6:8]
#将边界条件(力,...loc 取值0,1,2,3和4,分别代表单元的4个节点和单元中心(loc=4)
elem.calculate_Strain_4N(DISPe)
elem.calculate_Stress..._4N()
#print(f"单元{ie}在loc={loc}处的应力矩阵:")
#print(elem.Stress)
#找出每个节点对应的多(或1)个单元及在这些单元上的位置
def...应变和应力的求解是在单元中进行的
应变和应力 在各节处的取值(平均值 or最大值)又需要在 共享该节点的各单元上 取平均 或者取最大值
结果云图绘制又是按单元进行的
所以数据需要按 节点->单元 ->