有限元软件通用二次开发思路摘要.pptVIP

C MISES STRESSC SMISES=(STRESS(1)-STRESS(2))*(STRESS(1)-STRESS(2)) + (STRESS(2)-STRESS(3))*(STRESS(2)-STRESS(3)) + (STRESS(3)-STRESS(1))*(STRESS(3)-STRESS(1)) DO 90 K1=NDI+1,NTENS SMISES=SMISES+SIX*STRESS(K1)*STRESS(K1) CONTINUE SMISES=SQRT(SMISES/TWO) IF (SMISES.GT.(1.0+TOLER1)*SYIEL0) THEN C FLOW DIRECTIONC SHYDRO=(STRESS(1)+STRESS(2)+STRESS(3))/THREE ONESY=ONE/SMISES DO 110 K1=1,NDI FLOW(K1)=ONESY*(STRESS(K1)-SHYDRO) CONTINUE DO 120 K1=NDI+