ANSYS经典版三个主应力代数和云图显示UPFS子程序

c--- 状态变量数大于零时进行if(nsvar .gt. 0) thenc--- 在当前单元各积分点间循环do k = 1, numIntgc--- 获得应力向量call get_ElmData ('SIG ',elemId, k, ncomp, stress)c--- 给stitemp赋值do i=1,ncompstitemp(i) = stress(i)end doc--- 计算三个主应力代数和call prinst(stitemp(1))svar(1)=stitemp(7)+stitemp(8)+stitemp(9)      c--- 将新的状态变量写入ANSYS数据库call put_ElmData ('SVAR',elemId, k, nsvar, svar)end doend if