OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cstraini4.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| cstraini4 ../starter/source/elements/shell/coqueba/cstraini4.F
25!||--- called by ------------------------------------------------------
26!|| cbainit3 ../starter/source/elements/shell/coqueba/cbainit3.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| cg2leps ../starter/source/elements/shell/coqueba/scigini4.F
30!||--- uses -----------------------------------------------------
31!|| message_mod ../starter/share/message_module/message_mod.F
32!||====================================================================
33 SUBROUTINE cstraini4(
34 1 JFT ,JLT ,NFT ,NEL ,NUMSHEL,
35 2 ISTRAIN,GSTR ,SIGSH ,NSIGSH ,NUMEL ,
36 4 IX ,NIX ,PTSH ,THKE ,GSTRM ,
37 7 E1X ,E2X ,E3X ,E1Y ,E2Y ,
38 8 E3Y ,E1Z ,E2Z ,E3Z )
39C-----------------------------------------------
40C M o d u l e s
41C-----------------------------------------------
42 USE message_mod
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C A n a l y s e M o d u l e
49C-----------------------------------------------
50C-----------------------------------------------
51C G l o b a l P a r a m e t e r s
52C-----------------------------------------------
53#include "mvsiz_p.inc"
54C-----------------------------------------------
55C C o m m o n B l o c k s
56C-----------------------------------------------
57#include "com01_c.inc"
58C-----------------------------------------------
59C D u m m y A r g u m e n t s
60C-----------------------------------------------
61 INTEGER, INTENT(IN) :: JFT,JLT,NUMEL,NIX,NFT,
62 . ISTRAIN,NEL,NSIGSH,NUMSHEL
63 INTEGER, DIMENSION(NIX,NUMEL), INTENT(IN) :: IX
64 INTEGER, DIMENSION(NUMEL), INTENT(IN) :: PTSH
65 my_real, DIMENSION(NSIGSH,NUMSHEL),INTENT(IN) :: SIGSH
66 my_real, DIMENSION(NEL,8,4),INTENT(OUT) :: gstr
67 my_real, DIMENSION(NEL,8), INTENT(OUT) :: gstrm
68 my_real, DIMENSION(NEL),INTENT(IN) :: thke
69 my_real, DIMENSION(MVSIZ),INTENT(IN) ::
70 . e1x,e2x,e3x,e1y,e2y,e3y,e1z,e2z,e3z
71C------------------------------------------------------
72C L o c a l V a r i a b l e s
73C-----------------------------------------------
74 INTEGER I,II,J,JJ,N,NPTI,I1,I2,PT,PID1,IPID1,L_PLA,NLAY,
75 . ILAY,LAYNPT_MAX,LAY_MAX,NPTT,NPTMX,IP,PTS,LENS,IPT_ALL,
76 . IPT,PTN,JDIR,ILAW,NPGI,NPG
77 INTEGER LI(MVSIZ)
78 my_real
79 . e1(6),e2(6),z1,z2,z0,aa,e1g(6,4),e2g(6,4),z1g(4),z2g(4),ung,r1
80C=======================================================================
81!
82 npg=4
83 DO i=jft,jlt
84 IF (abs(isigi) /= 3.AND.abs(isigi) /= 4.AND.abs(isigi) /= 5)THEN
85 ii = i+nft
86 n = nint(sigsh(1,ii))
87 IF (n /= ix(nix,ii)) THEN
88 jj = ii
89 DO j = 1,numel
90 ii= j
91 n = nint(sigsh(1,ii))
92 IF (n == 0) GOTO 100
93 IF (n == ix(nix,jj)) GOTO 60
94 ENDDO
95 60 CONTINUE
96 ENDIF
97 ELSE
98 jj=nft+i
99 n =ix(nix,jj)
100 ii=ptsh(jj)
101 IF (ii == 0) GOTO 100
102 ENDIF
103 li(i) = ii
104 npti=nint(sigsh(2,ii))
105 npgi=nint(sigsh(nvshell,ii))
106C--- in global sys NPGI=4 for QEPH like stress now
107 IF(sigsh(17,ii) == one .AND. npgi==npg )THEN
108 pt = inishvar1
109 IF (npti==1) THEN
110 DO ip = 1,npg
111 e1(1:6) = sigsh(pt:pt+5,ii)
112 z1 = sigsh(pt+6,ii)
113 CALL cg2leps(
114 7 e1x(i) ,e2x(i),e3x(i),e1y(i),e2y(i),e3y(i),
115 8 e1z(i) ,e2z(i),e3z(i),e1 )
116 gstr(i,1:5,ip) = e1(1:5)
117 pt = pt + 7
118 END DO
119 ELSE
120 aa = half*thke(i)
121 DO ip = 1,npg
122 e1g(1:6,ip) = sigsh(pt:pt+5,ii)
123 z1g(ip) = sigsh(pt+6,ii)
124 pt = pt + 7
125 END DO
126 DO ip = 1,npg
127 e2g(1:6,ip) = sigsh(pt:pt+5,ii)
128 z2g(ip) = sigsh(pt+6,ii)
129 pt = pt + 7
130 END DO
131 DO ip = 1,npg
132 CALL cg2leps(
133 7 e1x(i) ,e2x(i),e3x(i),e1y(i),e2y(i),e3y(i),
134 8 e1z(i) ,e2z(i),e3z(i),e1g(1,ip))
135 CALL cg2leps(
136 7 e1x(i) ,e2x(i),e3x(i),e1y(i),e2y(i),e3y(i),
137 8 e1z(i) ,e2z(i),e3z(i),e2g(1,ip))
138 IF (z1g(ip)==z2g(ip)) THEN
139c error out
140 CALL ancmsg(msgid=1904,
141 . anmode=aninfo,
142 . msgtype=msgerror,
143 . i1=n,
144 . r1=z1g(ip))
145 ELSEIF (z1g(ip)==zero) THEN
146 gstr(i,1:5,ip)=e1g(1:5,ip)
147 z0 = aa*z2g(ip)
148 gstr(i,6:8,ip)=(e2g(1:3,ip)-e1g(1:3,ip))/z0
149 ELSEIF (z2g(ip)==zero) THEN
150 gstr(i,1:5,ip)=e2g(1:5,ip)
151 z0 = aa*z1g(ip)
152 gstr(i,6:8,ip)=(e1g(1:3,ip)-e2g(1:3,ip))/z0
153 ELSE
154 z0 = aa*(z2g(ip)-z1g(ip))
155 gstr(i,6:8,ip)=(e2g(1:3,ip)-e1g(1:3,ip))/z0
156 gstr(i,1:3,ip)=e1g(1:3,ip)-aa*z1g(ip)*gstr(i,6:8,ip)
157 gstr(i,4:5,ip)= half*(e2g(4:5,ip)+e1g(4:5,ip))
158 END IF
159 END DO
160 END IF
161 gstrm(i,1:8) = zero
162 ung = one/npg
163 DO ip = 1,npg
164 gstrm(i,1:8) = gstrm(i,1:8) + ung*gstr(i,1:8,ip)
165 END DO
166 ENDIF
167 100 CONTINUE
168 ENDDO
169C-----------
170 RETURN
171 END
cstraini4
subroutine cstraini4(jft, jlt, nft, nel, numshel, istrain, gstr, sigsh, nsigsh, numel, ix, nix, ptsh, thke, gstrm, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition cstraini4.F:39
message_mod
Definition message_mod.F:1249
cg2leps
subroutine cg2leps(e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, eps)
Definition scigini4.F:794
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889