OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
viscmain.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!|| viscmain ../engine/source/materials/visc/viscmain.F
25!||--- called by ------------------------------------------------------
26!|| mmain ../engine/source/materials/mat_share/mmain.F90
27!|| mulaw ../engine/source/materials/mat_share/mulaw.F90
28!|| usermat_solid ../engine/source/materials/mat_share/usermat_solid.F
29!||--- calls -----------------------------------------------------
30!|| damping_range_solid ../engine/source/general_controls/damping/damping_range_solid.F90
31!|| visc_plas ../engine/source/materials/visc/visc_plas.f90
32!|| visc_prony ../engine/source/materials/visc/visc_prony.f
33!|| visc_prony_lstrain ../engine/source/materials/visc/visc_prony_lstrain.f
34!||--- uses -----------------------------------------------------
35!|| damping_range_solid_mod ../engine/source/general_controls/damping/damping_range_solid.F90
36!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
37!|| visc_param_mod ../common_source/modules/mat_elem/visc_param_mod.F90
38!|| visc_plas_mod ../engine/source/materials/visc/visc_plas.F90
39!||====================================================================
40 SUBROUTINE viscmain(VISC ,NEL ,
41 . NVARVIS ,UVARVIS ,RHO ,VISCMAX ,SOUNDSP ,TIMESTEP,
42 . EPSPXX ,EPSPYY ,EPSPZZ ,EPSPXY ,EPSPYZ ,EPSPZX ,
43 . DEPSXX ,DEPSYY ,DEPSZZ ,DEPSXY ,DEPSYZ ,DEPSZX ,
44 . SVO1 ,SVO2 ,SVO3 ,SVO4 ,SVO5 ,SVO6 ,
45 . SV1 ,SV2 ,SV3 ,SV4 ,SV5 ,SV6 ,
46 . MFXX ,MFXY ,MFXZ ,MFYX ,MFYY ,MFYZ ,
47 . MFZX ,MFZY ,MFZZ ,
48 . S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
49 . DAMP_BUF,IDAMP_FREQ_RANGE ,MVSIZ ,ET ,YOUNG ,
50 . SHEAR_MODULUS,IVISC )
51C-----------------------------------------------
52C M o d u l e s
53C-----------------------------------------------
54 USE visc_param_mod
55 USE elbufdef_mod, only: buf_damp_range_
56 USE damping_range_solid_mod
57 USE visc_plas_mod
58C-----------------------------------------------
59C I m p l i c i t T y p e s
60C-----------------------------------------------
61#include "implicit_f.inc"
62C-----------------------------------------------
63C G l o b a l P a r a m e t e r s
64C-----------------------------------------------
65#include "param_c.inc"
66C-----------------------------------------------
67C D u m m y A r g u m e n t s
68C-----------------------------------------------
69 INTEGER ,INTENT(IN) :: NEL,NVARVIS
70 INTEGER ,INTENT(IN) :: IDAMP_FREQ_RANGE
71 INTEGER ,INTENT(IN) :: MVSIZ
72 INTEGER ,INTENT(IN) :: IVISC
73 my_real ,INTENT(IN) :: TIMESTEP
74 my_real ,INTENT(IN) :: YOUNG
75 my_real, INTENT(IN) :: SHEAR_MODULUS
76 my_real, DIMENSION(NEL), INTENT(INOUT):: S1,S2,S3,S4,S5,S6
77 .
78 my_real, DIMENSION(NEL) ,INTENT(INOUT) :: SOUNDSP, VISCMAX
79 my_real, DIMENSION(NEL) ,INTENT(IN) :: epspxx,epspyy,epspzz,epspxy,epspyz,epspzx,
80 . mfxx,mfxy,mfxz,mfyx,mfyy,mfyz,mfzx,mfzy,mfzz,
81 . depsxx,depsyy,depszz,depsxy,depsyz,depszx,
82 . svo1,svo2,svo3,svo4,svo5,svo6,rho
83 my_real, DIMENSION(NEL) ,INTENT(OUT) :: sv1,sv2,sv3,sv4,sv5,sv6
84
85 my_real, DIMENSION(NEL*NVARVIS) ,INTENT(INOUT) :: uvarvis
86 my_real, DIMENSION(MVSIZ) ,INTENT(IN) :: et
87 TYPE(visc_param_) ,INTENT(IN) :: VISC
88 TYPE(buf_damp_range_) ,INTENT(IN) :: DAMP_BUF
89C-----------------------------------------------
90C L o c a l V a r i a b l e s
91C-----------------------------------------------
92 INTEGER :: NPRONY,J,NVAR_DAMP
93 my_real, DIMENSION(NEL) :: sdamp1,sdamp2,sdamp3,sdamp4,sdamp5,sdamp6
94C=======================================================================
95C
96 IF (idamp_freq_range > 0) THEN
97 nvar_damp = 21
98 sv1(1:nel) = zero
99 sv2(1:nel) = zero
100 sv3(1:nel) = zero
101 sv4(1:nel) = zero
102 sv5(1:nel) = zero
103 sv6(1:nel) = zero
104 ELSE
105 nvar_damp = 0
106 ENDIF
107C
108 IF (ivisc > 0) THEN
109C-- /VISC PRONY
110 IF (visc%ILAW == 1) THEN
111c
112 nprony = visc%IPARAM(1)
113 CALL visc_prony(visc ,nprony ,nel ,nvarvis ,uvarvis ,
114 . epspxx ,epspyy ,epspzz ,epspxy ,epspyz ,epspzx ,
115 . sv1 ,sv2 ,sv3 ,sv4 ,sv5 ,sv6 ,
116 . timestep,rho ,viscmax ,soundsp ,nvar_damp)
117c
118 ELSEIF (visc%ILAW == 2) THEN
119c
120 nprony = visc%IPARAM(1)
121 CALL visc_prony_lstrain(
122 . visc ,nprony ,nel ,nvarvis ,uvarvis ,
123 . timestep,mfxx ,mfxy ,mfxz ,mfyx ,mfyy ,
124 . mfyz ,mfzx ,mfzy ,mfzz ,
125 . s1 ,s2 ,s3 ,s4 ,s5 ,s6)
126 ELSEIF (visc%ILAW == 3) THEN
127 call visc_plas(
128 . visc ,nel , rho ,soundsp ,
129 . depsxx ,depsyy ,depszz ,depsxy ,depsyz ,depszx,
130 . svo1 ,svo2 ,svo3 ,svo4 ,svo5 , svo6,
131 . sv1 ,sv2 ,sv3 ,sv4 ,sv5 ,sv6)
132 ENDIF
133 ENDIF
134C
135C-- Damping frequency range - same formulation as visc/prony with 3 components
136 IF (idamp_freq_range > 0) THEN
137 call damping_range_solid(damp_buf,nel ,nvarvis ,uvarvis ,nvar_damp ,et ,
138 . epspxx ,epspyy ,epspzz ,epspxy ,epspyz ,epspzx ,
139 . sdamp1 ,sdamp2 ,sdamp3 ,sdamp4 ,sdamp5 ,sdamp6 ,
140 . timestep,rho ,soundsp ,young ,shear_modulus)
141 sv1(1:nel) = sv1(1:nel) + sdamp1(1:nel)
142 sv2(1:nel) = sv2(1:nel) + sdamp2(1:nel)
143 sv3(1:nel) = sv3(1:nel) + sdamp3(1:nel)
144 sv4(1:nel) = sv4(1:nel) + sdamp4(1:nel)
145 sv5(1:nel) = sv5(1:nel) + sdamp5(1:nel)
146 sv6(1:nel) = sv6(1:nel) + sdamp6(1:nel)
147 ENDIF
148
149c-----------
150 RETURN
151 END
visc_prony
subroutine visc_prony(visc, nprony, nel, nvarvis, uvarvis, epspxx, epspyy, epspzz, epspxy, epspyz, epspzx, sv1, sv2, sv3, sv4, sv5, sv6, timestep, rho, viscmax, soundsp, nvar_damp)
Definition visc_prony.F:34
visc_prony_lstrain
subroutine visc_prony_lstrain(visc, nprony, nel, nvarvis, uvarvis, timestep, mfxx, mfxy, mfxz, mfyx, mfyy, mfyz, mfzx, mfzy, mfzz, s1, s2, s3, s4, s5, s6)
Definition visc_prony_lstrain.F:35
viscmain
subroutine viscmain(visc, nel, nvarvis, uvarvis, rho, viscmax, soundsp, timestep, epspxx, epspyy, epspzz, epspxy, epspyz, epspzx, depsxx, depsyy, depszz, depsxy, depsyz, depszx, svo1, svo2, svo3, svo4, svo5, svo6, sv1, sv2, sv3, sv4, sv5, sv6, mfxx, mfxy, mfxz, mfyx, mfyy, mfyz, mfzx, mfzy, mfzz, s1, s2, s3, s4, s5, s6, damp_buf, idamp_freq_range, mvsiz, et, young, shear_modulus, ivisc)
Definition viscmain.F:51