OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
alewdx_grid_bcs.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!|| alewdx_grid_bcs ../engine/source/ale/grid/alewdx_grid_bcs.F
25!||--- called by ------------------------------------------------------
26!|| alew5 ../engine/source/ale/grid/alew5.f
27!||====================================================================
28 SUBROUTINE alewdx_grid_bcs( SKEW, ISKEW, ICODT, VEC, NALE, NODE_ID)
29C-----------------------------------------------
30C D e s c r i p t i o n
31C-----------------------------------------------
32C This subroutine is imposing BCS for grid points.
33C done at each iteration since skew may evolve during a smoothing iteration
34C-----------------------------------------------
35C I m p l i c i t T y p e s
36C-----------------------------------------------
37#include "implicit_f.inc"
38C-----------------------------------------------
39C C o m m o n B l o c k s
40C-----------------------------------------------
41#include "com04_c.inc"
42#include "param_c.inc"
43C-----------------------------------------------
44C D u m m y A r g u m e n t s
45C-----------------------------------------------
46 INTEGER,INTENT(IN) :: NODE_ID
47 INTEGER,INTENT(INOUT) :: ISKEW(*),ICODT(*),NALE(NUMNOD)
48 my_real,INTENT(INOUT) :: skew(lskew,*),vec(3)
49C-----------------------------------------------
50C L o c a l V a r i a b l e s
51C-----------------------------------------------
52 INTEGER N, ISK, LCOD1,LCOD2, NINDX
53 my_real aa
54 LOGICAL TAG(3)
55C-----------------------------------------------
56C S o u r c e L i n e s
57C-----------------------------------------------
58 lcod1 = 0
59 lcod2 = 0
60C----------------------------------------
61C BOUNDARY CONDITIONS (MESH)
62C----------------------------------------
63 !RETRIEVING CODE FOR BCS et ALE/BCS
64 ! 000 <=> LCOD=0
65 ! 001 <=> LCOD=1 (Z)
66 ! 010 <=> LCOD=2 (Y)
67 ! 011 <=> LCOD=3 (YZ)
68 ! 100 <=> LCOD=4 (X)
69 ! 101 <=> LCOD=5 (XZ)
70 ! 110 <=> LCOD=6 (XY)
71 ! 111 <=> LCOD=7 (XYZ)
72 tag(1:3)=.false.
73 nindx = 0
74 n = node_id
75 IF(nale(n) /= 0)THEN
76 isk=iskew(n)
77 lcod1=icodt(n+numnod) !BCS
78 IF(lcod1 /= 0) THEN
79 nindx = n
80 ENDIF
81 lcod2=icodt(n+numnod+numnod) !ALE/BCS
82 IF(lcod2 /= 0) THEN
83 nindx = n
84 ENDIF
85 ENDIF
86
87 !TEST ENABLED BINARY (check for X,Y or Z in LCOD1 or LCOD2)
88 !
89 IF(btest(lcod1,0).OR.btest(lcod2,0))tag(3)=.true. ! **1 => condition on Z dir
90 IF(btest(lcod1,1).OR.btest(lcod2,1))tag(2)=.true. ! *1* => condition on Y dir
91 IF(btest(lcod1,2).OR.btest(lcod2,2))tag(1)=.true. ! 1** => condition on X dir
92
93 !UPDATE VECTOR DEPENDING ON BOUNDARY CONDITION AT RELATED NODE
94 !
95 IF (nindx /= 0 )THEN
96 isk = iskew(n)
97 IF(isk == 1) THEN
98 !=====REPERE GLOBAL
99 IF(tag(1))vec(1)=zero
100 IF(tag(2))vec(2)=zero
101 IF(tag(3))vec(3)=zero
102 ELSE
103 !=====REPERE OBLIQUE
104 IF(tag(1))THEN
105 aa =skew(1,isk)*vec(1)+skew(2,isk)*vec(2)+skew(3,isk)*vec(3)
106 vec(1)=vec(1)-skew(1,isk)*aa
107 vec(2)=vec(2)-skew(2,isk)*aa
108 vec(3)=vec(3)-skew(3,isk)*aa
109 ENDIF
110 IF(tag(2))THEN
111 aa =skew(4,isk)*vec(1)+skew(5,isk)*vec(2)+skew(6,isk)*vec(3)
112 vec(1)=vec(1)-skew(4,isk)*aa
113 vec(2)=vec(2)-skew(5,isk)*aa
114 vec(3)=vec(3)-skew(6,isk)*aa
115 ENDIF
116 IF(tag(3))THEN
117 aa =skew(7,isk)*vec(1)+skew(8,isk)*vec(2)+skew(9,isk)*vec(3)
118 vec(1)=vec(1)-skew(7,isk)*aa
119 vec(2)=vec(2)-skew(8,isk)*aa
120 vec(3)=vec(3)-skew(9,isk)*aa
121 ENDIF
122 ENDIF
123 ENDIF
124
125 RETURN
126 END
alew5
subroutine alew5(x, v, w, wa, ale_nn_connect, nale, nodft, nodlt, nbrcvois, nbsdvois, lnrcvois, lnsdvois, skew, iskew, icodt)
Definition alew5.F:40
alewdx_grid_bcs
subroutine alewdx_grid_bcs(skew, iskew, icodt, vec, nale, node_id)
Definition alewdx_grid_bcs.F:29
my_real
#define my_real
Definition cppsort.cpp:32