📄 a04b.for
字号:
write(66,3404) ICP,ICP2
write(66,3700)
175 ICP=ICP+2
C
IF (IBUG.EQ.0) GO TO 177
IF (MODEX.EQ.0) RETURN
STOP
C
177 ETOLD=PROP(1)
DO 230 J=6,NCON,2
ET=(PROP(J - 1) - PROP(J - 3))/(PROP(J) - PROP(J - 2))
IF (IDATWR.LE.1) write(66,2113) PROP(J-1),PROP(J),ET
C
IF (ET.LT.0.D0) GO TO 179
IF (ET.GE.PROP(1)) GO TO 179
IF (ET.LE.ETOLD) GO TO 178
179 write(66,3600) NG,N
write(66,3405)
write(66,3700)
IF (MODEX.EQ.0) RETURN
STOP
178 ETOLD=ET
C
230 CONTINUE
write(66,2211)
RETURN
C
C
C
C
10 IBUG=0
NPTS=INT(PROP(105))
XCRP=PROP(107)
XINTP=PROP(108)
XSUBM=PROP(109)
XITE=PROP(110)
XALG=PROP(111)
TOLIL=PROP(112)
TOLPC=PROP(113)
C
IF(NPTS.GT.0) GO TO 95
PROP(105)=1.6D1
NPTS=16
C
95 IF(XSUBM.EQ.0.D0) PROP(109)=1.D1
IF(XALG.EQ.2.D0 .AND. XSUBM.LT.3.D0) PROP(109)=3.D0
IF(XITE.EQ.0.D0) PROP(110)=1.5D1
IF(XALG.EQ.0.D0) PROP(111)=1.D0
IF(TOLIL.EQ.0.D0) PROP(112)=5.D-3
IF(TOLPC.EQ.0.D0) PROP(113)=1.D-1
C
C
IF (IDATWR.GT.1) GO TO 120
write(66,2301)
DO 115 K=1,16
IP1=K + 16
IP2=K + 32
IP3=K + 48
IP4=K + 64
IP5=K + 80
115 write(66,2302) PROP(K),PROP(IP1),PROP(IP2),PROP(IP3),PROP(IP4),
1 PROP(IP5)
write(66,2303) (PROP(M),M=97,104)
write(66,2304) (PROP(M),M=105,112)
IF (PROP(111).EQ.1.D0) write(66,2307)
IF (PROP(111).EQ.2.D0) write(66,2308) PROP(113)
C
120 IF (PROP(110).LT.6.D0) write(66,2305)
IF (XCRP.GE.0.D0 .AND. XCRP.LE.3.D0) GO TO 100
write(66,3600) NG,N
write(66,3000)
write(66,3700)
IBUG=1
C
100 IF(XINTP.GE.0.D0.AND.XINTP.LE.1.D0) GO TO 102
write(66,3600) NG,N
write(66,3001)
write(66,3700)
IBUG=1
C
102 IF(NPTS.GE.2.AND.NPTS.LE.16) GO TO 104
IBUG=1
write(66,3600) NG,N
write(66,3002)
write(66,3700)
GO TO 126
C
104 DO 106 J=2,NPTS
JJ=J-1
IF(PROP(J).GT.PROP(JJ)) GO TO 106
IBUG=1
write(66,3600) NG,N
write(66,3003)
write(66,3700)
GO TO 126
106 CONTINUE
C
126 IF (IBUG.EQ.0 .OR. MODEX.EQ.0) GO TO 127
STOP
C
127 JBUG=0
DO 125 I=1,NPTS
IK1=I+16
IK2=I+32
IK3=I+48
IK4=I+64
IF (PROP(IK1).LE.0.D0) JBUG=1
IF (PROP(IK2).LE.-1.D0 .OR. PROP(IK2).GE.0.5D0) JBUG=1
IF (PROP(IK3).LE.0.D0) JBUG=1
IF (PROP(IK4).LT.0.D0) JBUG=1
125 IF (PROP(IK4).GE.PROP(IK1)) JBUG=1
C
IF (JBUG.EQ.0) RETURN
write(66,3600) NG,N
write(66,3850)
write(66,3401)
write(66,3700)
IF (MODEX.EQ.0) RETURN
STOP
C
C
C.... MODEL = 12 EMPTY
C
12 RETURN
C
C
13 write(66,2109) PROP(1),PROP(2)
RETURN
C
C
C
C
C
14 IF (PROP(123).LT.1.D0) PROP(123)=1.D1
write(66,4000)
DO 300 K=1,16
JCP1 =K + 16
JCP2 =K + 32
JCP3 =K + 48
JCP4 =K + 64
JCP5 =K + 80
JCP6 =K + 96
300 write(66,4010) PROP(K),PROP(JCP1),PROP(JCP2),PROP(JCP3),
1 PROP(JCP4),PROP(JCP5),PROP(JCP6)
write(66,4020) (PROP(JC),JC=113,120)
write(66,4030) (PROP(JD),JD=121,128)
C
IBUG=0
NPTS=INT(PROP(121))
C
IF (NPTS.GE.2 .AND. NPTS.LE.16) GO TO 305
IBUG=1
write(66,3600) NG,N
write(66,3002)
write(66,3700)
GO TO 315
C
305 DO 310 J=2,NPTS
JJ=J-1
IF (PROP(J).GT.PROP(JJ)) GO TO 310
IBUG=1
write(66,3600) NG,N
write(66,3003)
write(66,3700)
GO TO 315
310 CONTINUE
C
315 IF (MODEX.EQ.0 .OR. IBUG.EQ.0) RETURN
STOP
C
C
C
1000 FORMAT (I5,4F10.0)
1100 FORMAT (8F10.0)
2100 FORMAT (30H MATERIAL CONSTANTS SET NUMBER,6H .... ,I5//,
1 1H ,4X,29HDEN ..........( DENSITY ).. =, E14.6/)
2101 FORMAT (1H ,4X,29HE ............( PROP(1) ).. =, E14.6/,
1 1H ,4X,29HVNU ..........( PROP(2) ).. =, E14.6///)
2102 FORMAT (1H ,4X,29HE(A) .........( PROP(1) ).. =, E14.6/,
1 1H ,4X,29HE(B) .........( PROP(2) ).. =, E14.6/,
2 1H ,4X,29HE(C) .........( PROP(3) ).. =, E14.6/,
3 1H ,4X,29HVNU(AB) ......( PROP(4) ).. =, E14.6/,
4 1H ,4X,29HVNU(AC) ......( PROP(5) ).. =, E14.6/,
5 1H ,4X,29HVNU(BC) ......( PROP(6) ).. =, E14.6/,
6 1H ,4X,29HG(AB) ........( PROP(7) ).. =, E14.6///)
2301 FORMAT (1H ,4X,17HTEMP (PROP(1-16)),5X,15HE (PROP(17-32)),5X,
1 17HVNU (PROP(33-48)),3X,19HYIELD (PROP(49-64)),3X,
2 16HET (PROP(65-80)),4X,19HALPHA (PROP(81-96)),/)
2302 FORMAT (1H ,4X,6(E14.6,7X))
2303 FORMAT (1H ,//,5X,33HCREEP LAW COEFFICIENTS ..........,//,
1 1H ,4X,30HA0 ............(PROP(97 )).. =,E14.6,/,
2 1H ,4X,30HA1 ............(PROP(98 )).. =,E14.6,/,
3 1H ,4X,30HA2 ............(PROP(99 )).. =,E14.6,/,
4 1H ,4X,30HA3 ............(PROP(100)).. =,E14.6,/,
5 1H ,4X,30HA4 ............(PROP(101)).. =,E14.6,/,
6 1H ,4X,30HA5 ............(PROP(102)).. =,E14.6,/,
7 1H ,4X,30HA6 ............(PROP(103)).. =,E14.6,/,
8 1H ,4X,30HA7 ............(PROP(104)).. =,E14.6,//)
2304 FORMAT (1H ,4X,66HNUMBER OF TEMPERATURE POINTS ...................
1...(PROP(105)).. =,E14.6,/,
2 1H ,4X,66HREFERENCE TEMPERATURE ..........................
3...(PROP(106)).. =,E14.6,/,
4 1H ,4X,66HCREEP LAW KEY ..................................
5...(PROP(107)).. =,E14.6,/,
6 1H ,4X,66HINTEGRATION PARAMETER ..........................
7...(PROP(108)).. =,E14.6,/,
8 1H ,4X,66HMAXIMUM NUMBER OF SUBDIVISIONS .................
9...(PROP(109)).. =,E14.6,/,
A 1H ,4X,66HMAXIMUM NUMBER OF ITERATIONS PER SUBDIVISION ...
B...(PROP(110)).. =,E14.6,/,
C 1H ,4X,66HALGORITHM INDICATOR ............................
D...(PROP(111)).. =,E14.6,/,
E 1H ,4X,66HCONVERGENCE TOLERANCE ..........................
F...(PROP(112)).. =,E14.6)
2307 FORMAT (1H ,4X,80HINELASTIC STRAIN TOLERANCE .....................
G...(PROP(113)).. = ( NOT USED ) //)
2308 FORMAT (1H ,4X,66HINELASTIC STRAIN TOLERANCE .....................
H...(PROP(113)).. =,E14.6,//)
2305 FORMAT (//32H * * * WARNING * * *,//,
1 45H THE USE OF PROP(110) .LT. 6 CAN RESULT IN,/,
2 45H A HIGHLY INACCURATE OR DIVERGENT SOLUTION,//,
3 32H * * * END OF MESSAGE * * *,////)
2103 FORMAT (1H ,4X,17HTEMP (PROP(1-16)),5X,15HE (PROP(17-32)),5X,
1 17HVNU (PROP(33-48)),4X,19HALPHA (PROP(49-64)),/)
2104 FORMAT (1H ,4X,4(E14.6,7X))
2105 FORMAT ( //,4X,46HNUMBER OF TEMPERATURE POINTS ...(PROP(65)).. =,
1 E14.6,/,
2 1H ,3X,46HREFERENCE TEMPERATURE ..........(PROP(66)).. =,
3 E14.6,//)
2106 FORMAT (1H ,4X,29HE ............( PROP(1) ).. =, E14.6/,
1 1H ,4X,29HVNU ..........( PROP(2) ).. =, E14.6/,
2 1H ,4X,29HYIELD ........( PROP(3) ).. =, E14.6/,
3 1H ,4X,29HE (HARDEN) ...( PROP(4) ).. =, E14.6///)
2109 FORMAT (1H ,4X,29HC1 ...........( PROP(1) ).. =, E14.6/,
1 1H ,4X,29HC2 ...........( PROP(2) ).. =, E14.6///)
2110 FORMAT(1H ,4X,29HE ............( PROP(1) ).. =, E14.6/,
1 1H ,4X,29HVNU ..........( PROP(2) ).. =, E14.6/,
2 1H ,4X,29HALFA .........( PROP(3) ).. =, E14.6/,
O 1H ,4X,29HK ............( PROP(4) ).. =, E14.6/,
4 1H ,4X,29HW ............( PROP(5) ).. =, E14.6/,
5 1H ,4X,29HD ............( PROP(6) ).. =, E14.6/,
6 1H ,4X,29HT ............( PROP(7) ).. =, E14.6/,
7 1H ,4X,29HI1A ..........( PROP(8) ).. =, E14.6///)
2111 FORMAT (1H ,4X,29HE ............( PROP(1) ).. =,E14.6,/,
1 1H ,4X,29HVNU ..........( PROP(2) ).. =,E14.6,/,
2 1H ,4X,29HYIELD ........( PROP(3) ).. =,E14.6,//)
2112 FORMAT (1H ,4X,36HPIECEWISE-LINEAR STRESS-STRAIN CURVE,/,
1 1H ,6X,6HSTRESS,10X,6HSTRAIN,12X,2HET,//,
2 6X,E14.6,2X,E14.6)
2113 FORMAT (6X,3(E14.6,2X))
2114 FORMAT (6X,E14.6,2X,E14.6,2X)
2200 FORMAT (///
1 19X,6HVOLUME,8X,7HLOADING,6X,9HUNLOADING,8X,7HLOADING, /
2 19X,6HSTRAIN,2(3X,12HBULK MODULUS),2X,13HSHEAR MODULUS,
3 / 1X)
2210 FORMAT (7H POINT(,I1,2H) ,4E15.4)
2211 FORMAT (////)
2220 FORMAT (35H CRACKING MODE . . . . . (ICRACK) =,I5,/
1 43H EQ.0, CURVE DESCRIPTION NONLINEAR MODEL /,
2 36H EQ.1, SOIL MODEL WITH NO TENSION /,
3 48H EQ.2, SOIL MODEL, NO TENSION, STRESS RELEASE /,
5 31H WEIGHT DENSITY = (,E13.4,1H),/
1 31H STIFFNESS REDUCTION FACTOR = (,E13.4,1H),/
2 31H SHEAR REDUCTION FACTOR = (,E13.4,1H),1X)
2230 FORMAT (//40H (A) UNIAXIAL PARAMETERS ,
1 //49H INITIAL TANGENT MODULUS . . . . . . . (PROP(1))=,E14.6,
2 /49H POISSONS RATIO. . . . . . . . . . . . (PROP(2))=,E14.6,
3 /49H COEFFICIENT OF THERMAL EXPANSION . . .(PROP(3))=,E14.6,
4 /49H UNIAXIAL CUT-OFF TENSILE STRENGTH . . (SIGMAT)=,E14.6,
1 /49H UNIAXIAL MAXIMUM COMPRESSIVE STRESS . .(SIGMAC)=,E14.6,
2 /49H COMPRESSIVE STRAIN AT SIGMAC . . . . . ( EPSC )=,E14.6,
5 /49H UNIAXIAL ULTIMATE COMPRESSIVE STRESS . (SIGMAU)=,E14.6,
8 /49H UNIAXIAL ULTIMATE COMPRESSIVE STRAIN . ( EPSU )=,E14.6)
2235 FORMAT (//40H (B) TRIAXIAL COMPRESSIVE FAILURE CURVES,
1 //4X,10H PRINCIPAL,5X,30X,12HCURVE NUMBER/1X,
2 16H STRESS RATIOS/,9X,3HI=1,10X,1H2,11X,1H3,11X,1H4,11X,
3 1H5,11X,1H6/1X,90(1H-)//1X,6X,6HSP1(I),6X,6F12.4//1X,
4 5X,8HSP3(I,1),5X,6F12.4,/1X,3X,12H(AT SP2=SP1),/1X,
5 5X,8HSP3(I,2),5X,6F12.4,/2X,17H(AT SP2=BETA*SP3),/1X,
6 5X,8HSP3(I,3),5X,6F12.4,/1X,3X,12H(AT SP2=SP3)//1X,90(1H-))
2240 FORMAT (/40H (C) VARIOUS OTHER CONTROL PARAMETERS //
1 ,49H STRESS RATIO FOR FAILURE SURFACE INPUT .(BETA) =,E14.6/
2 ,49H STRAINS SCALING FACTOR - MULTIAXIALITY .(GAMA) =,E14.6/
3 ,49H CONTROL FOR CHANGING MATERIAL LAW . . . (KAPA) =,E14.6/
4 ,49H CONTROL FOR LOADING/UNLOADING CRITERION (ALFA) =,E14.6/
5 ,49H STIFFNESS REDUCTION FACTOR . . . . . .(STIFAC) =,E14.6/
6 ,49H SHEAR REDUCTION FACTOR . . . . . . . .(SHEFAC) =,E14.6)
2241 FORMAT (49H REFERENCE TEMPERATURE . . . . . . . . . (TREF) =,
1 E14.6)
2242 FORMAT (//32H * * * WARNING * * *,//,
1 61H THE TRIAXIAL COMPRESSIVE FAILURE CURVES SHOULD BE
2 CHECKED,//,
2 32H * * * END OF MESSAGE * * *,////)
2243 FORMAT (//32H * * * WARNING * * *,//,
1 29H DEFAULT VALUES ARE USED ,/,
2 35H BECAUSE NO INPUT CURVE IS GIVEN,//,
3 32H * * * END OF MESSAGE * * *,////)
2244 FORMAT (//37H * * * STOP OF SOLUTION * * *,//,
1 58H THE TRIAXIAL COMPRESSIVE FAILURE CURVES SHOULD HA
*VE : ,/,
2 50H SP1(1) .LE. SP1(2) ............. .LE. SP1(6),/,
3 50H SP3(1,1) .LE. SP3(2,1) ....... .LE. SP3(6,1),/,
4 50H SP3(1,2) .LE. SP3(2,2) ....... .LE. SP3(6,2),/,
5 50H SP3(1,3) .LE. SP3(2,3) ....... .LE. SP3(6,3),//
6 37H * * * END OF ERROR MESSAGE * * *,//)
2500 FORMAT (1H ,4X,5HPROP(,I2,10H) ...... =, E14.6)
3000 FORMAT (/ 30H INCORRECT CREEP LAW NUMBER)
3001 FORMAT (/ 35H INCORRECT INTEGRATION PARAMETER)
3002 FORMAT (/ 42H INCORRECT NUMBER OF TEMPERATURE POINTS)
3003 FORMAT (/,52H TEMPERATURE POINTS MUST BE SUCCESSIVELY LARGER )
3399 FORMAT (/,5X,33HPROP(4) MUST BE ZERO OR EQUAL TO ,
1 16HPROP(3)/PROP(1) )
3401 FORMAT (/ 35H YOUNG*S MODULUS MUST BE .GT. 0.,/,
1 49H POISSON*S RATIO MUST BE .GT. -1. AND .LT. 0.5,/,
* 50H HARDENING MODULUS MUST BE .LT. YOUNG*S MODULUS,/,
3 37H HARDENING MODULUS MUST BE .GE. 0.,/,
3 40H INITIAL YIELD STRESS MUST BE .GT. 0.)
3403 FORMAT (//50H INPUT ERROR IN MATERIAL PROPERTIES //
1 15H *** STOP *** //)
3404 FORMAT (/ 48H FOR THE MULTI-LINEAR ELASTIC-PLASTIC MODEL ,/,
1 9H PROP(,I2,16H) EQUAL TO PROP(,I2,16H) IS NOT ALLOW
2ED,/,
3 35H BECAUSE ET CANNOT BE CALCULATED)
3405 FORMAT (/ 46H FOR THE MULTI-LINEAR ELASTIC-PLASTIC MODEL,/,
1 61H WE MUST HAVE E .GT. ET1 .GE. ET2 .GE. ET3 .......
2 .GE. 0.)
3410 FORMAT (/ 35H YOUNG*S MODULUS MUST BE .GT. 0.,/,
1 48H POISSON*S RATIO MUST BE .GE. 0. AND .LT. 0.5,/,
2 35H ALFA MUST BE .GE. 0.,/,
3 35H K MUST BE .GT. 0.,/,
4 35H W AND D MUST BE .LT. 0.,/,
5 35H T MUST BE .GE. 0.,/,
6 35H INITIAL CAP POS.MUST BE .LE. 0.)
3500 FORMAT (/ 45H INITIAL MODULUS (PROP(1)) MUST BE .GT. 0.,/,
1 49H POISSON*S RATIO MUST BE .GT. -1. AND .LT. 0.5,/,
* 26H SIGMAT MUST BE .GE. 0.,/,
* 49H SIGMAC, EPSC, SIGMAU AND EPSU MUST BE .LT. 0.,/,
2 39H ES=SIGMAC/EPSC MUST BE .LE. PROP(1),/,
3 30H SIGMAU MUST BE .GE. SIGMAC,/,
4 26H EPSU MUST BE .LE. EPSC,/,
5 36H BETA MUST BE .GE. 0. AND .LE. 1.,/,
6 24H GAMA MUST BE .GE. 0.,/,
7 36H KAPA MUST BE .GE. 0. AND .LE. 1.,/,
8 59H STIFFNESS REDUCTION FACTORS MUST BE .GE. 0. AND .
9LE. 1.)
3501 FORMAT (//32H * * * WARNING * * *,//,
1 51H THE FOLLOWING CONTROL PARAMETERS ARE RECOMMENDED/,
2 25H FOR THE CONCRETE MODEL,/,
3 30H GAMA .GE. 0. AND .LE. 6.0 ,/,
4 17H ALFA .LE. 0. ,//,
5 32H * * * END OF MESSAGE * * *,///)
3600 FORMAT (//37H * * * STOP OF SOLUTION * * *,//,
1 35H INPUT ERROR IN MATERIAL PROPERTIES,/,
2 19H ELEMENT GROUP NO =,I5,/,
3 27H MATERIAL PROPERTY SET NO =,I5)
3650 FORMAT (/ 35H YOUNG*S MODULUS MUST BE .GT. 0.,/,
1 49H POISSON*S RATIO MUST BE .GT. -1. AND .LT. 0.5)
3700 FORMAT (/ 40H SEE USER*S MANUAL FOR INPUT DESCRIPTION,//,
1 37H * * * END OF ERROR MESSAGE * * *,//)
3750 FORMAT (/ 36H MODULI (E AND G) MUST BE .GT. 0.,/,
1 63H MODULI AND POISSON RATIOS DO NOT SATISFY REQUIRED
2 RELATIONS)
3800 FORMAT (/ 59H MODULI (K AND G) MUST BE .GT. 0. FOR ALL STRAIN
2STATES,/,
1 59H STIFFNESS REDUCTION FACTORS MUST BE .GE. 0. AND .
2LE. 1.)
3850 FORMAT (/ 44H FOR ALL TEMPERATURES THE MODEL REQUIRES :)
4000 FORMAT (/,42H TEMPERATURE-DEPENDENT MATERIAL CONSTANTS ,
A /,6X,10HPROP(1-16),8X,11HPROP(17-32),8X,11HPROP(33-48),8X,
A 11HPROP(49-64),8X,11HPROP(65-80),8X,11HPROP(81-96),7X,
B 12HPROP(97-112),/,9X,4HTEMP,15X,4HALFA,13X,6HCTD(1),
C 13X,6HCTD(2),13X,6HCTD(3),13X,6HCTD(4),13X,6HCTD(5),/)
4010 FORMAT (2H ,6(E14.6,5X),E14.6)
4020 FORMAT (///,50H TEMPERATURE-INDEPENDENT MATERIAL CONSTANTS ,
A //,4X,27H CTI(1)........PROP(113) =,E14.6,
B /,4X,27H CTI(2)........PROP(114) =,E14.6,
C /,4X,27H CTI(3)........PROP(115) =,E14.6,
D /,4X,27H CTI(4)........PROP(116) =,E14.6,
E /,4X,27H CTI(5)........PROP(117) =,E14.6,
F /,4X,27H CTI(6)........PROP(118) =,E14.6,
G /,4X,27H CTI(7)........PROP(119) =,E14.6,
H /,4X,27H CTI(8)........PROP(120) =,E14.6,///)
4030 FORMAT (5X,29HNUMBER OF TEMPERATURE POINTS ,25(1H.),
X 12HPROP(121) =,E14.6,/,
Y 5X,22HREFERENCE TEMPERATURE ,32(1H.),12HPROP(122) =,E14.6,/,
Z 5X,24HNUMBER OF SUBINCREMENTS ,30(1H.),12HPROP(123) =,E14.6,/,
A 5X,37HSOLUTION CONTROL PARAMETER 1 (SCP(1)),17(1H.),
A 12HPROP(124) =,E14.6,/,
B 5X,37HSOLUTION CONTROL PARAMETER 2 (SCP(2)),17(1H.),
B 12HPROP(125) =,E14.6,/,
C 5X,37HSOLUTION CONTROL PARAMETER 3 (SCP(3)),17(1H.),
C 12HPROP(126) =,E14.6,/,
D 5X,37HSOLUTION CONTROL PARAMETER 4 (SCP(4)),17(1H.),
D 12HPROP(127) =,E14.6,/,
E 5X,37HSOLUTION CONTROL PARAMETER 5 (SCP(5)),17(1H.),
E 12HPROP(128) =,E14.6//)
C
C
END
SUBROUTINE PLST3D (XYZ,XX,RE,EDIS,S,IEL,IFLAG)
C
C
C
C***ADD:DPR***
IMPLICIT DOUBLE PRECISION ( A-H,O-Z )
C***END:DPR***
COMMON /VAR/ NG,MODEX,IUPDT,KSTEP,ITEMAX,IEQREF,ITE,KPRI,IREF,
1 IEQUIT,IPRI,KPLOTN,KPLOTE
COMMON /MATMOD/ STRESS(4),STRAIN(4),D(4,4),IPT,NEL,IPS,ISVE
COMMON /XYZLL/ T(3,3),XLOC
COMMON /GNLSTF/ SGNL(45)
C
DIMENSION XYZ(*),XX(*),RE(*),EDIS(*),S(*),DUM(378),ILSK(9),TMA(9)
C
C
X21=XYZ(4)-XYZ(1)
Y21=XYZ(5)-XYZ(2)
Z21=XYZ(6)-XYZ(3)
X31=XYZ(7)-XYZ(1)
Y31=XYZ(8)-XYZ(2)
Z31=XYZ(9)-XYZ(3)
C
XY=X21*Y31 - Y21*X31
YZ=Y21*Z31 - Z21*Y31
ZX=Z21*X31 - X21*Z31
SX=X21*X21 + Y21*Y21 + Z21*Z21
SZ= YZ*YZ + ZX*ZX + XY*XY
C
CALL DETCH (NG,NEL,SX,1)
CALL DETCH (NG,NEL,SZ,1)
C
SX=SQRT(SX)
SZ=SQRT(SZ)
SY=SX*SZ
C
T(1,1)= X21/SX
T(1,2)= Y21/SX
T(1,3)= Z21/SX
T(2,1)=(Z21*ZX - Y21*XY)/SY
T(2,2)=(X21*XY - Z21*YZ)/SY
T(2,3)=(Y21*YZ - X21*ZX)/SY
T(3,1)= YZ/SZ
T(3,2)= ZX/SZ
T(3,3)= XY/SZ
C
IF (IFLAG .GE. 4) GO TO 70
C
C LOCAL PLANE
C
XLOC=0.D0
IF (IFLAG.GT.0) GO TO 45
DO 40 I=1,3
40 XLOC=XLOC + T(3,I)*XYZ(I)
C
45 DO 50 I=1,IEL
I2=2*I
I3=3*I
XX(I2-1)=T(1,1)*XYZ(I3-2) + T(1,2)*XYZ(I3-1) + T(1,3)*XYZ(I3)
XX(I2 )=T(2,1)*XYZ(I3-2) + T(2,2)*XYZ(I3-1) + T(2,3)*XYZ(I3)
IF (IFLAG.LE.1) GO TO 50
TEMPA=T(1,1)*EDIS(I3-2) + T(1,2)*EDIS(I3-1) + T(1,3)*EDIS(I3)
EDIS(I2 )=T(2,1)*EDIS(I3-2) + T(2,2)*EDIS(I3-1) + T(2,3)*EDIS(I3)
EDIS(I2-1)=TEMPA
IF (IFLAG .EQ. 3) GO TO 50
TEMPA=T(1,1)*RE(I3-2) + T(1,2)*RE(I3-1) + T(1,3)*RE(I3)
RE(I2 )=T(2,1)*RE(I3-2) + T(2,2)*RE(I3-1) + T(2,3)*RE(I3)
RE(I2-1)=TEMPA
50 CONTINUE
RETURN
C
C
70 NDOFL=2*IEL
NDOFG=3*IEL
IF (IFLAG .EQ. 5) GO TO 100
DO 75 I=1,NDOFL
75 DUM(I)=RE(I)
DO 80 I=1,IEL
I2=2*I
I3=3*I
RE(I3-2)=T(1,1)*DUM(I2-1) + T(2,1)*DUM(I2)
RE(I3-1)=T(1,2)*DUM(I2-1) + T(2,2)*DUM(I2)
RE(I3 )=T(1,3)*DUM(I2-1) + T(2,3)*DUM(I2)
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -