OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
int18_law151_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!|| int18_law151_init ../engine/source/interfaces/int18/int18_law151_init.F
25!||--- called by ------------------------------------------------------
26!|| resol ../engine/source/engine/resol.F
27!||--- uses -----------------------------------------------------
28!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
29!|| groupdef_mod ../common_source/modules/groupdef_mod.f
30!|| message_mod ../engine/share/message_module/message_mod.F
31!|| multi_fvm_mod ../common_source/modules/ale/multi_fvm_mod.F90
32!||====================================================================
33 SUBROUTINE int18_law151_init(MULTI_FVM,IGRBRIC,IPARI,IXS,
34 1 IGROUPS ,IPARG ,ELBUF_TAB ,FORCE_INT ,
35 2 X ,V ,MS ,KINET ,
36 3 X_APPEND,V_APPEND,MASS_APPEND,KINET_APPEND,
37 4 FORCE_INT_PON)
38
39!$COMMENT
40! INT18_LAW151_INIT description
41! initialization of array for /INT18 + LAW 151
42!
43! INT18_LAW151_INIT organization :
44! - x/V/MASS_append :
45! * (1:NUMNOD) : classical x/v/mass
46! * (NUMNOD+1:NUMNOD+NUMELS) : phantom node x/v/mass
47! - kinet_append :
48! * (1:NUMNOD) : classical kinet
49! * (NUMNOD+1:NUMNOD+NUMELS) : =0 (no kinet for phantom node)
50! - force_int array :
51! * for parith/off :
52! dimension = 3,NUMELS*NTHREADS
53! each itask thread works in (1+itask*numels:(1+itask)*numels)
54!$ENDCOMMENT
55C-----------------------------------------------
56C M o d u l e s
57C-----------------------------------------------
58 USE multi_fvm_mod
59 USE groupdef_mod
60 USE elbufdef_mod
61 USE message_mod
62C-----------------------------------------------
63C I m p l i c i t T y p e s
64C-----------------------------------------------
65#include "implicit_f.inc"
66#include "comlock.inc"
67#include "com04_c.inc"
68#include "param_c.inc"
69#include "com01_c.inc"
70C-----------------------------------------------
71C D u m m y A r g u m e n t s
72C-----------------------------------------------
73 INTEGER, DIMENSION(NPARI,*), INTENT(in) :: IPARI
74 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), INTENT(in) :: ELBUF_TAB
75 INTEGER, DIMENSION(NUMELS), INTENT(in) ::IGROUPS
76 INTEGER, DIMENSION(NPARG,*), INTENT(in) ::IPARG
77 my_real, DIMENSION(3,*), INTENT(in) :: x,v
78 my_real, DIMENSION(3,*), INTENT(inout) :: x_append,v_append
79 my_real, DIMENSION(*), INTENT(in) :: ms
80 INTEGER, DIMENSION(*), INTENT(in) :: KINET
81 my_real, DIMENSION(*), INTENT(inout) :: mass_append
82 INTEGER, DIMENSION(*), INTENT(inout) :: KINET_APPEND
83 my_real, DIMENSION(3,*), INTENT(inout) :: force_int
84 REAL(kind=8), dimension(3,6,*), INTENT(inout) :: force_int_pon
85 INTEGER, DIMENSION(NIXS, *), INTENT(in) :: IXS
86 TYPE(multi_fvm_struct) :: MULTI_FVM
87 TYPE (GROUP_) , DIMENSION(NGRBRIC), INTENT(in) :: IGRBRIC
88C-----------------------------------------------
89C L o c a l V a r i a b l e s
90C-----------------------------------------------
91 INTEGER :: N,NN,II,JJ,MY_SIZE
92 INTEGER :: NFT,GROUP_ID,ILOC
93 INTEGER :: ISU1,NSN,NTY,INACTI,NODE_ID,IBRIC
94 LOGICAL, DIMENSION(:), ALLOCATABLE :: TAG
95C-----------------------------------------------
96 ! X/V/MASS extended
97 ! 1:NUMNOD --> classical x/v/mass
98 ! NUMNOD+1:NUMNOD+NUMELS --> x/v/mass of phantom nodes (located to the center of
99 ! the ALE elements)
100 ! INT18_GLOBAL_LIST --> true if /INT18 + LAW151 for the NIN interface
101 multi_fvm%INT18_GLOBAL_LIST(1:ninter) = .false.
102 ALLOCATE(tag(numels))
103 tag(1:numels)=.false.
104 IF( multi_fvm%IS_INT18_LAW151 ) THEN
105 DO n=1,multi_fvm%NUMBER_INT18
106 nn = multi_fvm%INT18_LIST(n)
107 multi_fvm%INT18_GLOBAL_LIST(nn) = .true.
108 ENDDO
109! ! 1:NUMNOD --> classical x/v/mass
110 x_append(1:3,1:numnod) = x(1:3,1:numnod)
111 v_append(1:3,1:numnod) = v(1:3,1:numnod)
112 mass_append(1:numnod) = ms(1:numnod)
113 kinet_append(1:numnod) = kinet(1:numnod)
114! ! initialization of phantom nodes
115 x_append( 1:3,numnod+1:numnod+numels ) = zero
116 v_append( 1:3,numnod+1:numnod+numels ) = zero
117 ! KINET : not used for phantom nodes, always equal to 0
118 mass_append(numnod+1:numnod+numels) = zero
119 kinet_append(numnod+1:numnod+numels) = 0
120 force_int(1:multi_fvm%SIZE_FORCE_INT_1,1:multi_fvm%SIZE_FORCE_INT_2) = zero
121
122 my_size = multi_fvm%SIZE_FORCE_INT_PON
123 force_int_pon(1:3,1:6,1:my_size) = zero
124
125 DO nn=1,multi_fvm%NUMBER_INT18
126 n = multi_fvm%INT18_LIST(nn)
127 isu1 = ipari(45,n)
128 nsn = ipari(5,n) ! number of secondary nodes
129 DO ii = 1,nsn
130 ibric = igrbric(isu1)%ENTITY(ii) ! id of the phantom element
131 IF(.NOT. tag(ibric)) THEN ! otherwise already done with a previous interface
132 group_id = igroups(ibric) ! id of the element group
133 nft = iparg(3,group_id)
134 iloc = ibric - nft
135 ! mass
136 mass_append(numnod + ibric) = elbuf_tab(group_id)%GBUF%RHO(iloc) * elbuf_tab(group_id)%GBUF%VOL(iloc)
137 ! position
138 DO jj = 2, 9
139 node_id = ixs(jj, ibric) ! id of node of the phantom element
140 x_append(1, numnod + ibric) = x_append(1, numnod + ibric) + one_over_8 * x(1, node_id)
141 x_append(2, numnod + ibric) = x_append(2, numnod + ibric) + one_over_8 * x(2, node_id)
142 x_append(3, numnod + ibric) = x_append(3, numnod + ibric) + one_over_8 * x(3, node_id)
143 ENDDO
144 ! velocity
145 v_append(1, numnod + ibric) = multi_fvm%VEL(1, ibric)
146 v_append(2, numnod + ibric) = multi_fvm%VEL(2, ibric)
147 v_append(3, numnod + ibric) = multi_fvm%VEL(3, ibric)
148 tag(ibric)=.true.
149 ENDIF
150 ENDDO
151 ENDDO
152 ENDIF
153 IF(ALLOCATED(tag))DEALLOCATE(tag)
154
155 RETURN
156 END SUBROUTINE int18_law151_init
157C===============================================================================
my_real
#define my_real
Definition cppsort.cpp:32
the
end diagonal values have been computed in the(sparse) matrix id.SOL
int18_law151_init
subroutine int18_law151_init(multi_fvm, igrbric, ipari, ixs, igroups, iparg, elbuf_tab, force_int, x, v, ms, kinet, x_append, v_append, mass_append, kinet_append, force_int_pon)
Definition int18_law151_init.F:38
id
initmumps id
Definition multiplerhs_example.m:4
groupdef_mod
Definition groupdef_mod.F:662
message_mod
Definition message_mod.F:1249