OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
law158_init.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!|| law158_init ../starter/source/materials/mat/mat158/law158_init.F
25!||--- called by ------------------------------------------------------
26!|| c3init3 ../starter/source/elements/sh3n/coque3n/c3init3.F
27!|| cinit3 ../starter/source/elements/shell/coque/cinit3.F
28!|| cmaini3 ../starter/source/elements/sh3n/coquedk/cmaini3.F
29!||====================================================================
30 SUBROUTINE law158_init(DIR1 ,DIR2 ,
31 . UVAR,ALDT,NEL ,NUVAR ,TAN_PHI ,
32 . X1 ,X2 ,X3 ,X4 ,Y1 ,Y2 ,
33 . Y3 ,Y4 ,Z1 ,Z2 ,Z3 ,Z4 ,
34 . E1X ,E2X ,E3X ,E1Y ,E2Y ,E3Y ,E1Z ,E2Z ,E3Z )
35C-----------------------------------------------
36C I m p l i c i t T y p e s
37C-----------------------------------------------
38#include "implicit_f.inc"
39C-----------------------------------------------
40C C o m m o n B l o c k s
41C-----------------------------------------------
42#include "vect01_c.inc"
43#include "param_c.inc"
44C-----------------------------------------------
45C D u m m y A r g u m e n t s
46C-----------------------------------------------
47 INTEGER ,INTENT(IN) :: NEL,NUVAR
48 my_real ,DIMENSION(NEL) ,INTENT(IN) :: ALDT
49 my_real ,DIMENSION(NEL,2) ,INTENT(IN) :: DIR1,DIR2
50 my_real, DIMENSION(NEL) ,INTENT(IN) :: E1X,E2X,E3X,E1Y,E2Y,
51 . e3y,e1z,e2z,e3z,x1,x2,x3,x4,y1,y2,y3,y4,z1,z2,z3,z4
52 my_real ,DIMENSION(NEL,NUVAR) ,INTENT(OUT) :: uvar
53 my_real ,DIMENSION(NEL) ,INTENT(OUT) :: tan_phi
54C-----------------------------------------------
55C L o c a l V a r i a b l e s
56C-----------------------------------------------
57 INTEGER I
58 my_real
59 . R,S,D1,D2,D11,D12,D21,D22,U1X,U1Y,U2X,U2Y,DET,W1X,W2X,W1Y,W2Y,
60 . PHI,CSP,SNP,TANA,R1,R2,S1,S2,AA,BB,SUMA,V1,V2,V3,WR,WS
61 my_real , DIMENSION(NEL) :: e11,e12,e13,e21,e22,e23
62C======================================================================|
63C--- isoparametric (material) axes
64 IF (ity == 3) THEN
65C--- shell 4N
66 DO i=1,nel
67 e11(i)= x2(i)+x3(i)-x1(i)-x4(i)
68 e12(i)= y2(i)+y3(i)-y1(i)-y4(i)
69 e13(i)= z2(i)+z3(i)-z1(i)-z4(i)
70 e21(i)= x3(i)+x4(i)-x1(i)-x2(i)
71 e22(i)= y3(i)+y4(i)-y1(i)-y2(i)
72 e23(i)= z3(i)+z4(i)-z1(i)-z2(i)
73 ENDDO
74 ELSEIF (ity == 7) THEN
75C--- shell 3N
76 DO i=1,nel
77 e11(i)= x2(i)-x1(i)
78 e12(i)= y2(i)-y1(i)
79 e13(i)= z2(i)-z1(i)
80 e21(i)= x3(i)-x1(i)
81 e22(i)= y3(i)-y1(i)
82 e23(i)= z3(i)-z1(i)
83 ENDDO
84 ENDIF
85C------ Initialise UVAR
86 DO i=1,nel
87 aa = dir1(i,1)
88 bb = dir1(i,2)
89 v1 = aa*e11(i) + bb*e21(i)
90 v2 = aa*e12(i) + bb*e22(i)
91 v3 = aa*e13(i) + bb*e23(i)
92 wr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
93 ws = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
94 suma = max( sqrt(wr*wr + ws*ws), em20)
95 r1 = wr/suma
96 s1 = ws/suma
97 aa = dir2(i,1)
98 bb = dir2(i,2)
99 v1 = aa*e11(i) + bb*e21(i)
100 v2 = aa*e12(i) + bb*e22(i)
101 v3 = aa*e13(i) + bb*e23(i)
102 wr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
103 ws = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
104 suma = max( sqrt(wr*wr + ws*ws), em20)
105 r2 = wr/suma
106 s2 = ws/suma
107 tana = (r1*r2 + s1*s2) / (r1*s2 - r2*s1) ! Tan(alpha_totale)
108c
109 uvar(i,1:nuvar) = zero
110 uvar(i,6) = tana
111 tan_phi(i) = tana
112 uvar(i,14) = aldt(i)
113 ENDDO
114C----
115 RETURN
116 END
law158_init
subroutine law158_init(dir1, dir2, uvar, aldt, nel, nuvar, tan_phi, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition law158_init.F:35
max
#define max(a, b)
Definition macros.h:21