OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_prop27.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!|| hm_read_prop27 ../starter/source/properties/spring/hm_read_prop27.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_properties ../starter/source/properties/hm_read_properties.F
27!||--- calls -----------------------------------------------------
28!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
29!|| hm_get_floatv_dim ../starter/source/devtools/hm_reader/hm_get_floatv_dim.F
30!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
31!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
32!||--- uses -----------------------------------------------------
33!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
34!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
35!|| message_mod ../starter/share/message_module/message_mod.f
36!|| submodel_mod ../starter/share/modules1/submodel_mod.F
37!||====================================================================
38 SUBROUTINE hm_read_prop27(GEO ,IGEO ,PROP_TAG ,IGTYP ,IG ,
39 . UNITAB ,LSUBMODEL)
40C-----------------------------------------------
41 USE unitab_mod
42 USE elbuftag_mod
43 USE message_mod
44 USE submodel_mod
45 USE hm_option_read_mod
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50C-----------------------------------------------
51C C o m m o n B l o c k s
52C-----------------------------------------------
53#include "units_c.inc"
54#include "param_c.inc"
55#include "tablen_c.inc"
56C-----------------------------------------------
57C D u m m y A r g u m e n t s
58C-----------------------------------------------
59 TYPE (UNIT_TYPE_),INTENT(IN) :: UNITAB
60 INTEGER IGEO(NPROPGI),IGTYP,IG
61C REAL
62 my_real
63 . geo(npropg)
64 TYPE(prop_tag_) , DIMENSION(0:MAXPROP) :: PROP_TAG
65 TYPE(submodel_data),INTENT(IN) :: LSUBMODEL(*)
66C-----------------------------------------------
67C L o c a l V a r i a b l e s
68C-----------------------------------------------
69 INTEGER IFUNC,IFUNC2,ISENS,IFL,IFAIL,J,ILENG,ITENS,Fsmooth
70C REAL
71 my_real
72 . mass,stiff,damp,gap,min_rup,max_rup,
73 . ascale1,ascale2,fscale1,fscale2,
74 . ascale1_unit,ascale2_unit,fscale1_unit,
75 . fscale2_unit,nexp,fcut
76 LOGICAL IS_AVAILABLE, IS_ENCRYPTED
77C=======================================================================
78C
79 is_encrypted = .false.
80 is_available = .false.
81C
82C--------------------------------------------------
83C EXTRACT DATA (IS OPTION CRYPTED)
84C--------------------------------------------------
85 CALL hm_option_is_encrypted(is_encrypted)
86C
87C--------------------------------------------------
88C EXTRACT DATAS (REAL VALUES)
89C-------------------------------------------------
90C 1st card
91 CALL hm_get_floatv('MASS' ,mass ,is_available,lsubmodel,unitab)
92 CALL hm_get_intv('ISENSOR' ,isens ,is_available,lsubmodel)
93 CALL hm_get_intv('ISFLAG' ,ifl ,is_available,lsubmodel)
94 CALL hm_get_intv('Ileng' ,ileng ,is_available,lsubmodel)
95 CALL hm_get_intv('Itens' ,itens ,is_available,lsubmodel)
96 CALL hm_get_intv('Ifail' ,ifail ,is_available,lsubmodel)
97C 2nd card
98 CALL hm_get_floatv('STIFF' ,stiff ,is_available,lsubmodel,unitab)
99 CALL hm_get_floatv('DAMP' ,damp ,is_available,lsubmodel,unitab)
100 CALL hm_get_floatv('NEXP' ,nexp ,is_available,lsubmodel,unitab)
101 CALL hm_get_floatv('MIN_RUP' ,min_rup ,is_available,lsubmodel,unitab)
102 CALL hm_get_floatv('MAX_RUP' ,max_rup ,is_available,lsubmodel,unitab)
103C 3rd card
104 CALL hm_get_floatv('GAP' ,gap ,is_available,lsubmodel,unitab)
105 CALL hm_get_intv('FSMOOTH' ,fsmooth ,is_available,lsubmodel)
106 CALL hm_get_floatv('FCUT' ,fcut ,is_available,lsubmodel,unitab)
107C 4th card
108 CALL hm_get_intv('FUN1' ,ifunc ,is_available,lsubmodel)
109 CALL hm_get_intv('FUN2' ,ifunc2 ,is_available,lsubmodel)
110 CALL hm_get_floatv('ASCALE1' ,ascale1 ,is_available,lsubmodel,unitab)
111 IF (ascale1 == zero) THEN
112 CALL hm_get_floatv_dim('ASCALE1' ,ascale1_unit ,is_available, lsubmodel, unitab)
113 ascale1 = one * ascale1_unit
114 ENDIF
115 CALL hm_get_floatv('FSCALE1' ,fscale1 ,is_available,lsubmodel,unitab)
116 IF (fscale1 == zero) THEN
117 CALL hm_get_floatv_dim('FSCALE1' ,fscale1_unit ,is_available, lsubmodel, unitab)
118 fscale1 = one * fscale1_unit
119 ENDIF
120 CALL hm_get_floatv('ASCALE2' ,ascale2 ,is_available,lsubmodel,unitab)
121 IF (ascale2 == zero) THEN
122 CALL hm_get_floatv_dim('ASCALE2' ,ascale2_unit ,is_available, lsubmodel, unitab)
123 ascale2 = one * ascale2_unit
124 ENDIF
125 CALL hm_get_floatv('FSCALE2' ,fscale2 ,is_available,lsubmodel,unitab)
126 IF (fscale2 == zero) THEN
127 CALL hm_get_floatv_dim('FSCALE2' ,fscale2_unit ,is_available, lsubmodel, unitab)
128 fscale2 = one * fscale2_unit
129 ENDIF
130C
131C-----------------------------
132C CHECKING VALUES
133C-----------------------------
134C
135 ! Checking negative values of MIN_RUP and Gap
136 min_rup = -abs(min_rup)
137 gap = -abs(gap)
138 IF (gap < zero) itens = 0
139 ! Non-linear stiffness exponent
140 IF (nexp <= zero) nexp = one
141 ! Spring force filtering flag
142 IF (fsmooth /= 0) fsmooth = 1
143 ! If functions are activated
144 IF (ifunc /= 0) stiff = zero
145 IF (ifunc2 /= 0) damp = zero
146 ! Sensor flag
147 IF (ifl == 1) isens = -isens
148 ! Default failure limits
149 IF (min_rup == zero) min_rup = -infinity
150 IF (max_rup == zero) max_rup = infinity
151 ! Filtering frequency
152 IF (fcut /= zero) THEN
153 ! If a filtering frequency is given by the user
154 fsmooth = 1
155 ELSE
156 ! If no filtering frequency is given but the flag is activated
157 IF (fsmooth /= 0) THEN
158 fcut = 100000.0d0*unitab%FAC_T_WORK
159 ! If no filtering frequency and no flag is activated => no filtering
160 ELSE
161 fcut = zero
162 ENDIF
163 ENDIF
164C
165C-----------------------------
166C PROPERTY TABLES
167C-----------------------------
168 ! Filling IGEO table
169 igeo(1) = ig
170 igeo(3) = isens
171 igeo(11) = igtyp
172 igeo(101) = ifunc
173 igeo(102) = ifunc2
174C
175 ! Filling GEO table
176 geo(1) = mass
177 geo(2) = stiff
178 geo(3) = damp
179 geo(5) = ep06
180 geo(7) = zero
181 geo(8) = onep1
182 geo(10) = fscale1
183 geo(12) = igtyp+em01
184 geo(15) = min_rup
185 geo(16) = max_rup
186 geo(18) = ascale2
187 geo(19) = gap
188 geo(25) = zero ! NUVAR - Number of user variables
189 geo(39) = ascale1
190 geo(43) = ifail
191 geo(80) = ifl
192 geo(93) = ileng
193 geo(132) = fscale2
194 geo(133) = itens
195 geo(134) = nexp
196 geo(135) = fsmooth
197 geo(136) = fcut
198C
199C----------------------
200C FROM LECGEO - GENERAL
201C----------------------
202 IF (geo(171) /= zero .AND. igeo(10) == 0) igeo(10) = nint(geo(171))
203C
204C-----------------------------
205C PROPERTY BUFFER
206C-----------------------------
207 prop_tag(igtyp)%G_EINT = 1
208 prop_tag(igtyp)%G_FOR = 1
209 prop_tag(igtyp)%G_LENGTH = 1 ! X0 (AL0) - total length
210 prop_tag(igtyp)%G_TOTDEPL = 1 ! DX - total deformation (translation)
211 prop_tag(igtyp)%G_DEP_IN_TENS = 1 ! DPX (DPY,DPZ) - max displacement in tension
212 prop_tag(igtyp)%G_DEP_IN_COMP = 1 ! DPX2 (DPY2,DPZ2) - max displacement in compression
213 prop_tag(igtyp)%G_POSX = 5
214 prop_tag(igtyp)%G_LENGTH_ERR = 1
215 prop_tag(igtyp)%G_NUVAR = max(prop_tag(igtyp)%G_NUVAR,nint(geo(25))) ! additional internal variables for h=6
216 prop_tag(igtyp)%G_DEFINI = 1
217 prop_tag(igtyp)%G_FORINI = 1
218 prop_tag(igtyp)%G_RUPTCRIT = 1
219C
220C-----------------------------
221C PRINTING OUT DATA
222C-----------------------------
223 WRITE(iout,2000)
224 IF (.NOT. is_encrypted) THEN
225 ! Property and spring mass
226 IF (ileng == 0) THEN
227 WRITE(iout,1400) ig,mass
228 ELSE
229 WRITE(iout,1410) ig,mass
230 ENDIF
231 ! Spring stiffness
232 ! -> Linear
233 IF (ifunc == 0) THEN
234 IF (ileng == 0) THEN
235 WRITE(iout,1450) stiff,nexp
236 ELSE
237 WRITE(iout,1460) stiff,nexp
238 ENDIF
239 ! -> Non-linear tabulated
240 ELSE
241 IF (ileng == 0) THEN
242 WRITE(iout,1500) ifunc,ascale1,fscale1
243 ELSE
244 WRITE(iout,1510) ifunc,ascale1,fscale1
245 ENDIF
246 ENDIF
247 ! Spring damping
248 ! -> Linear
249 IF (ifunc2 == 0) THEN
250 IF (ileng == 0) THEN
251 WRITE(iout,1550) damp
252 ELSE
253 WRITE(iout,1560) damp
254 ENDIF
255 ! -> Non-linear tabulated
256 ELSE
257 IF (ileng == 0) THEN
258 WRITE(iout,1600) ifunc2,ascale2,fscale2
259 ELSE
260 WRITE(iout,1610) ifunc2,ascale2,fscale2
261 ENDIF
262 ENDIF
263 ! Gap for activation
264 IF (gap /= zero) THEN
265 IF (ileng == 0) THEN
266 WRITE(iout,1650) gap
267 ELSE
268 WRITE(iout,1660) gap
269 ENDIF
270 ENDIF
271 ! Unit length flag
272 WRITE(iout,1700) ileng
273 ! Tensile behavior flag
274 WRITE(iout,1800) itens
275 ! Velocity filtering
276 WRITE(iout,1900) fsmooth,fcut
277 ! Failure criterion
278 IF (ifail /= 0) THEN
279 IF (ileng == 0) THEN
280 WRITE(iout,1750) ifail,min_rup,max_rup
281 ELSE
282 WRITE(iout,1760) ifail,min_rup,max_rup
283 ENDIF
284 ENDIF
285 ELSE
286 WRITE(iout,1000) ig
287 ENDIF
288C
289C-------------------------------------------------------------------------
2902000 FORMAT(
291 & 5x,'------------------------------------------------------'/,
292 & 5x,' BOUNDED DAMPER SPRING PROPERTY SET '/,
293 & 5x,'------------------------------------------------------'/)
2941000 FORMAT(
295 & 5x,'PROPERTY SET NUMBER . . . . . . . . . . . . . . . . .=',i10/,
296 & 5x,'CONFIDENTIAL DATA'//)
2971400 FORMAT(
298 & 5x,'PROPERTY SET NUMBER . . . . . . . . . . . . . . . . .=',i10/,
299 & 5x,'SPRING MASS . . . . . . . . . . . . . . . . . . . . .=',1pg20.13/)
3001410 FORMAT(
301 & 5x,'PROPERTY SET NUMBER . . . . . . . . . . . . . . . . .=',i10/,
302 & 5x,'SPRING LINEIC MASS . . . . . . . . . . . . . . . . . .=',1pg20.13/)
3031450 FORMAT(
304 & 5x,'LINEAR SPRING STIFFNESS . . . . . . . . . . . . . . .=',1pg20.13/,
305 & 5x,'NON-LINEAR EXPONENT . . . . . . . . . . . . . . . . .=',1pg20.13/)
3061460 FORMAT(
307 & 5x,'LINEAR SPRING STIFFNESS PER UNIT LENGTH. . . . . . . .=',1pg20.13/,
308 & 5x,'NON-LINEAR EXPONENT . . . . . . . . . . . . . . . . .=',1pg20.13/)
3091500 FORMAT(
310 & 5x,'NON LINEAR STIFFNESS FORCE FUNCTION ID . . . . . . . .=',i10/,
311 & 5x,'ELONGATION SCALE FACTOR FOR STIFFNESS FUNCTION . . . .=',1pg20.13/,
312 & 5x,'FORCE SCALE FACTOR FOR STIFFNESS FUNCTION . . . . . .=',1pg20.13/)
3131510 FORMAT(
314 & 5x,'NON LINEAR STIFFNESS FORCE FUNCTION ID . . . . . . . .=',i10/,
315 & 5x,'STRAIN SCALE FACTOR FOR STIFFNESS FUNCTION . . . . . .=',1pg20.13/,
316 & 5x,'FORCE SCALE FACTOR FOR STIFFNESS FUNCTION . . . . . .=',1pg20.13/)
3171550 FORMAT(
318 & 5x,'LINEAR DAMPING COEFFICIENT . . . . . . . . . . . . . .=',1pg20.13/)
3191560 FORMAT(
320 & 5x,'LINEAR DAMPING COEFFICIENT PER UNIT LENGTH . . . . . .=',1pg20.13/)
3211600 FORMAT(
322 & 5x,'NON LINEAR DAMPING FORCE FUNCTION ID . . . . . . . . .=',i10/,
323 & 5x,'VELOCITY SCALE FACTOR FOR DAMPING FUNCTION . . . . . .=',1pg20.13/,
324 & 5x,'FORCE SCALE FACTOR FOR DAMPING FUNCTION . . . . . . .=',1pg20.13/)
3251610 FORMAT(
326 & 5x,'NON LINEAR DAMPING FORCE FUNCTION ID . . . . . . . . .=',i10/,
327 & 5x,'STRAIN-RATE SCALE FACTOR FOR DAMPING FUNCTION . . . .=',1pg20.13/,
328 & 5x,'FORCE SCALE FACTOR FOR DAMPING FUNCTION . . . . . . .=',1pg20.13/)
3291650 FORMAT(
330 & 5x,'MINIMUM COMPRESSION GAP (LENGTH) FOR ACTIVATION . . .=',1pg20.13/)
3311660 FORMAT(
332 & 5x,'MINIMUM COMPRESSION GAP (STRAIN) FOR ACTIVATION . . .=',1pg20.13/)
3331700 FORMAT(
334 & 5x,'UNIT LENGTH FLAG . . . . . . . . . . . . . . . . . . .=',i10/,
335 & 5x,' ILENG = 0 : INPUT VALUES WITH CLASSICAL UNITS ',/,
336 & 5x,' ILENG = 1 : INPUT VALUES PER UNIT LENGTH ',/)
3371750 FORMAT(
338 & 5x,'FAILURE FLAG IFAIL . . . . . . . . . . . . . . . . . .=',i10/,
339 & 5x,' IFAIL = 1 : DISPLACEMENT FAILURE CRITERION ',/,
340 & 5x,' IFAIL = 2 : FORCE FAILURE CRITERION ',/,
341 & 5x,'NEGATIVE FAILURE LIMIT . . . . . . . . . . . . . . . .=',1pg20.13/,
342 & 5x,'POSITIVE FAILURE LIMIT (IF ITENS = 1). . . . . . . . .=',1pg20.13/)
3431760 FORMAT(
344 & 5x,'FAILURE FLAG IFAIL . . . . . . . . . . . . . . . . . .=',i10/,
345 & 5x,' IFAIL = 1 : STRAIN FAILURE CRITERION ',/,
346 & 5x,' IFAIL = 2 : FORCE FAILURE CRITERION ',/,
347 & 5x,'NEGATIVE FAILURE LIMIT . . . . . . . . . . . . . . . .=',1pg20.13/,
348 & 5x,'POSITIVE FAILURE LIMIT (IF ITENS = 1). . . . . . . . .=',1pg20.13/)
3491800 FORMAT(
350 & 5x,'TENSILE BEHAVIOR FLAG . . . . . . . . . . . . . . . .=',i10/,
351 & 5x,' ITENS = 0 : NO STIFFNESS AND DAMPING IN TENSION ',/,
352 & 5x,' ITENS = 1 : STIFFNESS AND DAMPING FOR TENSION ',/)
3531900 FORMAT(
354 & 5x,'SPRING FORCE FILTERING FLAG . . . . . . . . . . . . .=',i10/,
355 & 5x,' FSMOOTH = 0 : NO FILTERING ',/,
356 & 5x,' FSMOOTH = 1 : FILTERING ACTIVATED ',/,
357 & 5x,'CUTOFF FREQUENCY . . . . . . . . . . . . . . . . . . .=',1pg20.13/)
358c-----------
359 END
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_floatv_dim
subroutine hm_get_floatv_dim(name, dim_fac, is_available, lsubmodel, unitab)
Definition hm_get_floatv_dim.F:159
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:284
hm_read_prop27
subroutine hm_read_prop27(geo, igeo, prop_tag, igtyp, ig, unitab, lsubmodel)
Definition hm_read_prop27.F:40
max
#define max(a, b)
Definition macros.h:21
elbuftag_mod
Definition elbuftag_mod.F:227
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
starter
program starter
Definition starter.F:39
elbuftag_mod::prop_tag_
Definition elbuftag_mod.F:374
submodel_mod::submodel_data
Definition submodel_mod.F:723