OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
rsens_nic.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!|| rsens_nic ../engine/source/tools/sensor/rsens_nic.F
25!||--- called by ------------------------------------------------------
26!|| rforc3 ../engine/source/elements/spring/rforc3.F
27!||--- uses -----------------------------------------------------
28!|| sensor_mod ../common_source/modules/sensor_mod.F90
29!||====================================================================
30 SUBROUTINE rsens_nic(NEL ,IXR ,FOR ,MOM ,SKEW ,
31 . NSENSOR,SENSOR_TAB)
32C-----------------------------------------------
33C M o d u l e s
34C-----------------------------------------------
35 USE sensor_mod
36C-----------------------------------------------
37C I m p l i c i t T y p e s
38C-----------------------------------------------
39#include "implicit_f.inc"
40C-----------------------------------------------
41C C o m m o n B l o c k s
42C-----------------------------------------------
43#include "param_c.inc"
44#include "com04_c.inc"
45C-----------------------------------------------
46C D u m m y A r g u m e n t s
47C-----------------------------------------------
48 INTEGER ,INTENT(IN) :: NEL,NSENSOR
49 INTEGER ,DIMENSION(NIXR,NUMELR) ,INTENT(IN) :: IXR
50 my_real ,DIMENSION(NEL) ,INTENT(IN) :: for,mom
51 TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) ,INTENT(INOUT) :: SENSOR_TAB
52C-----------------------------------------------
53C L o c a l V a r i a b l e s
54C-----------------------------------------------
55 INTEGER :: I,J,SPRING_ID,ISENS,SPR_SKW,SEN_SKW,IDIR_F,IDIR_M
56 my_real :: fz,my,ftens,fcomp,mflex,mext,fgx,fgy,fgz,mgx,mgy,mgz
57 my_real ,DIMENSION(3) :: ff,mm
58 my_real ,DIMENSION(LSKEW,NUMSKW) :: skew
59c-----------------------------------------------------------------------
60c send current normalized forces and moments in the spring to sensor for processing
61c=======================================================================
62 isens = 0
63 DO j = 1,nsensor
64 IF (sensor_tab(j)%TYPE == 19) THEN ! sensor NIC_Nij
65 isens = j
66 sen_skw = sensor_tab(isens)%IPARAM(3)
67 spr_skw = sensor_tab(isens)%IPARAM(4)
68 idir_f = sensor_tab(isens)%IPARAM(5)
69 idir_m = sensor_tab(isens)%IPARAM(6)
70c
71 ftens = sensor_tab(isens)%RPARAM(2)
72 fcomp = sensor_tab(isens)%RPARAM(3)
73 mflex = sensor_tab(isens)%RPARAM(4)
74 mext = sensor_tab(isens)%RPARAM(5)
75 EXIT
76 END IF
77 END DO
78c we assume there's only one NIC-Nij sensor in the model
79 IF (isens > 0) THEN
80 my = zero
81 fz = zero
82 DO i = 1,nel
83 spring_id = ixr(6,i)
84 IF (sensor_tab(isens)%IPARAM(1) == spring_id) THEN
85 ff(1) = for(i)
86 ff(2) = for(i+nel)
87 ff(3) = for(i+nel*2)
88 mm(1) = mom(i)
89 mm(2) = mom(i+nel)
90 mm(3) = mom(i+nel*2)
91c rotate to sensor skew if necessary
92 IF (sen_skw > 1) THEN
93 IF (sen_skw /= spr_skw) THEN
94 fgx = skew(1,spr_skw)*ff(1)+skew(4,spr_skw)*ff(2)+skew(7,spr_skw)*ff(3)
95 fgy = skew(2,spr_skw)*ff(1)+skew(5,spr_skw)*ff(2)+skew(8,spr_skw)*ff(3)
96 fgz = skew(3,spr_skw)*ff(1)+skew(6,spr_skw)*ff(2)+skew(9,spr_skw)*ff(3)
97 ff(1) = skew(1,sen_skw)*fgx+skew(2,sen_skw)*fgy+skew(3,sen_skw)*fgz
98 ff(2) = skew(4,sen_skw)*fgx+skew(5,sen_skw)*fgy+skew(6,sen_skw)*fgz
99 ff(3) = skew(7,sen_skw)*fgx+skew(8,sen_skw)*fgy+skew(9,sen_skw)*fgz
100 mgx = skew(1,spr_skw)*mm(1)+skew(4,spr_skw)*mm(2)+skew(7,spr_skw)*mm(3)
101 mgy = skew(2,spr_skw)*mm(1)+skew(5,spr_skw)*mm(2)+skew(8,spr_skw)*mm(3)
102 mgz = skew(3,spr_skw)*mm(1)+skew(6,spr_skw)*mm(2)+skew(9,spr_skw)*mm(3)
103 mm(1) = skew(1,sen_skw)*mgx+skew(2,sen_skw)*mgy+skew(3,sen_skw)*mgz
104 mm(2) = skew(4,sen_skw)*mgx+skew(5,sen_skw)*mgy+skew(6,sen_skw)*mgz
105 mm(3) = skew(7,sen_skw)*mgx+skew(8,sen_skw)*mgy+skew(9,sen_skw)*mgz
106 END IF
107 fz = ff(idir_f)
108 my = mm(idir_m)
109 END IF
110c save normalized FZ and My values
111 IF (fz >= zero) THEN
112 fz = fz / ftens
113 ELSE
114 fz = abs(fz) / fcomp
115 END IF
116 IF (my >= zero) THEN
117 my = my / mflex
118 ELSE
119 my = abs(my) / mext
120 END IF
121 sensor_tab(isens)%VAR(1) = fz
122 sensor_tab(isens)%VAR(2) = my
123 sensor_tab(isens)%VAR(3) = one ! spmd flag for sensor indicating present domain
124 EXIT
125 END IF
126 END DO
127 END IF
128c-----------
129 RETURN
130 END
my_real
#define my_real
Definition cppsort.cpp:32
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
rsens_nic
subroutine rsens_nic(nel, ixr, for, mom, skew, nsensor, sensor_tab)
Definition rsens_nic.F:32