OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
read_sensor_dist_surf.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!|| read_sensor_dist_surf ../starter/source/tools/sensor/read_sensor_dist_surf.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_sensors ../starter/source/tools/sensor/hm_read_sensors.F
27!||--- calls -----------------------------------------------------
28!|| anodset ../starter/source/output/analyse/analyse_node.c
29!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
30!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
31!|| usr2sys ../starter/source/system/sysfus.F
32!||--- uses -----------------------------------------------------
33!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
34!|| message_mod ../starter/share/message_module/message_mod.F
35!|| submodel_mod ../starter/share/modules1/submodel_mod.F
36!||====================================================================
37 SUBROUTINE read_sensor_dist_surf(SENSOR_PTR ,SENS_ID ,
38 . ITABM1 ,UNITAB ,LSUBMODEL )
39C-----------------------------------------------
40C M o d u l e s
41C-----------------------------------------------
42 USE unitab_mod
43 USE message_mod
44 USE submodel_mod
45 USE sensor_mod
46 USE hm_option_read_mod
47C-----------------------------------------------
48C I m p l i c i t T y p e s
49C-----------------------------------------------
50#include "implicit_f.inc"
51C-----------------------------------------------
52C A n a l y s e M o d u l e
53C-----------------------------------------------
54#include "analyse_name.inc"
55C-----------------------------------------------
56C C o m m o n B l o c k s
57C-----------------------------------------------
58#include "com04_c.inc"
59#include "units_c.inc"
60C-----------------------------------------------
61C D u m m y A r g u m e n t s
62C-----------------------------------------------
63 INTEGER ,INTENT(IN) :: SENS_ID
64 INTEGER ,DIMENSION(NUMNOD) ,INTENT(IN) :: ITABM1
65 TYPE (SENSOR_STR_) :: SENSOR_PTR
66 TYPE (SUBMODEL_DATA) ,DIMENSION(NSUBMOD) :: LSUBMODEL
67 TYPE (UNIT_TYPE_) ,INTENT(IN) ::UNITAB
68C-----------------------------------------------
69C L o c a l V a r i a b l e s
70C-----------------------------------------------
71 INTEGER :: NOD_ID,ID_N1,ID_N2,ID_N3,NS1,NS2,NS3,NS4,SENS_TYPE,NPARIS,NPARRS,NVAR
72 my_real :: tdel,dmin,dmax,tmin
73 CHARACTER(len=40) :: MESS
74 LOGICAL :: IS_AVAILABLE
75 DATA mess/'SENSOR DEFINITION '/
76C-----------------------------------------------
77C E x t e r n a l F u n c t i o n s
78C-----------------------------------------------
79 INTEGER USR2SYS
80C--------------------------------
81C SENSOR BASED ON RELATIVE DISPLACEMENT OF TWO NODES
82C=======================================================================
83 is_available = .false.
84 sens_type = 15
85c--------------------------------------------------
86card1
87 CALL hm_get_floatv('Tdelay' ,tdel ,is_available,lsubmodel,unitab)
88card2
89 CALL hm_get_intv ('NODE1' ,nod_id ,is_available,lsubmodel)
90 CALL hm_get_intv ('NODE2' ,id_n1 ,is_available,lsubmodel)
91 CALL hm_get_intv ('NODE3' ,id_n2 ,is_available,lsubmodel)
92 CALL hm_get_intv ('NODE4' ,id_n3 ,is_available,lsubmodel)
93card3
94 CALL hm_get_floatv('Dmin' ,dmin ,is_available,lsubmodel,unitab)
95 CALL hm_get_floatv('Dmax' ,dmax ,is_available,lsubmodel,unitab)
96 CALL hm_get_floatv('Tmin' ,tmin ,is_available,lsubmodel,unitab)
97c
98 IF (dmax == zero) dmax = infinity
99c-------------------
100 nparis = 5
101 nparrs = 2
102 nvar = 0
103 ALLOCATE (sensor_ptr%IPARAM(nparis))
104 ALLOCATE (sensor_ptr%RPARAM(nparrs))
105 ALLOCATE (sensor_ptr%VAR(nvar))
106 sensor_ptr%VAR(:) = zero
107c--------------------------------------------------
108c Transform node IDs to system numbers and save
109c--------------------------------------------------
110 sensor_ptr%IPARAM(1) = usr2sys(nod_id,itabm1,mess,ns1)
111 CALL anodset(sensor_ptr%IPARAM(1), check_used)
112 sensor_ptr%IPARAM(3) = usr2sys(id_n1,itabm1,mess,ns2)
113 CALL anodset(sensor_ptr%IPARAM(3), check_used)
114 sensor_ptr%IPARAM(4) = usr2sys(id_n2,itabm1,mess,ns3)
115 CALL anodset(sensor_ptr%IPARAM(4), check_used)
116 sensor_ptr%IPARAM(5) = usr2sys(id_n3,itabm1,mess,ns4)
117 CALL anodset(sensor_ptr%IPARAM(5), check_used)
118c--------------------------------------------------
119 sensor_ptr%NPARI = nparis
120 sensor_ptr%NPARR = nparrs
121 sensor_ptr%NVAR = nvar
122c
123 sensor_ptr%TYPE = sens_type ! sensor type
124 sensor_ptr%SENS_ID = sens_id
125 sensor_ptr%STATUS = 0 ! status = deactivated
126 sensor_ptr%TSTART = infinity
127 sensor_ptr%TCRIT = infinity
128 sensor_ptr%TMIN = tmin
129 sensor_ptr%TDELAY = tdel ! time delay before activation
130 sensor_ptr%VALUE = infinity
131c
132 sensor_ptr%IPARAM(2) = 0 ! ISURF
133 sensor_ptr%RPARAM(1) = dmin
134 sensor_ptr%RPARAM(2) = dmax
135c------------------------------------------------------------
136 WRITE(iout, 1000) sens_id,tdel
137 WRITE (iout,2000) nod_id,id_n1,id_n2,id_n3,dmin,dmax,tmin
138c------------------------------------------------------------
139 1000 FORMAT(
140 . 5x,' SENSOR TYPE 15: DISTANCE NODE TO SURFACE '/,
141 . 5x,' ---------------------------------------- '/,
142 . 5x,'SENSOR ID. . . . . . . . . . . . . . . . .=',i10/
143 . 5x,'TIME DELAY BEFORE ACTIVATION . . . . . . .=',e12.4)
144 2000 FORMAT(
145 . 5x,'NODE ID . . . . . . . . . . . . . . . . .=',i10/
146 . 5x,'SURFACE NODE 1 ID . . . . . . . . . . . .=',i10/
147 . 5x,'SURFACE NODE 2 ID . . . . . . . . . . . .=',i10/
148 . 5x,'SURFACE NODE 3 ID . . . . . . . . . . . .=',i10/
149 . 5x,'MINIMAL DISTANCE . . . . . . . . . . . . .=',e12.4/
150 . 5x,'MAXIMAL DISTANCE . . . . . . . . . . . . .=',e12.4/
151 . 5x,'MIN DURATION . . . . . . . . . . . . . . .=',e12.4/)
152c-----------
153 RETURN
154 END
anodset
void anodset(int *id, int *type)
Definition analyse_node.c:321
my_real
#define my_real
Definition cppsort.cpp:32
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:415
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_option_read_mod
Definition hm_option_read_mod.F:408
message_mod
Definition message_mod.F:1249
submodel_mod
Definition submodel_mod.F:712
unitab_mod
Definition unitab_mod.F:513
read_sensor_dist_surf
subroutine read_sensor_dist_surf(sensor_ptr, sens_id, itabm1, unitab, lsubmodel)
Definition read_sensor_dist_surf.F:39