OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
viscmain.F File Reference
#include "implicit_f.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines

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)

Function/Subroutine Documentation

◆ viscmain()

subroutine viscmain ( type(visc_param_), intent(in) visc,
integer, intent(in) nel,
integer, intent(in) nvarvis,
intent(inout) uvarvis,
intent(in) rho,
intent(inout) viscmax,
intent(inout) soundsp,
intent(in) timestep,
intent(in) epspxx,
intent(in) epspyy,
intent(in) epspzz,
intent(in) epspxy,
intent(in) epspyz,
intent(in) epspzx,
intent(in) depsxx,
intent(in) depsyy,
intent(in) depszz,
intent(in) depsxy,
intent(in) depsyz,
intent(in) depszx,
intent(in) svo1,
intent(in) svo2,
intent(in) svo3,
intent(in) svo4,
intent(in) svo5,
intent(in) svo6,
intent(out) sv1,
intent(out) sv2,
intent(out) sv3,
intent(out) sv4,
intent(out) sv5,
intent(out) sv6,
intent(in) mfxx,
intent(in) mfxy,
intent(in) mfxz,
intent(in) mfyx,
intent(in) mfyy,
intent(in) mfyz,
intent(in) mfzx,
intent(in) mfzy,
intent(in) mfzz,
intent(inout) s1,
intent(inout) s2,
intent(inout) s3,
intent(inout) s4,
intent(inout) s5,
intent(inout) s6,
type(buf_damp_range_), intent(in) damp_buf,
integer, intent(in) idamp_freq_range,
integer, intent(in) mvsiz,
intent(in) et,
intent(in) young,
intent(in) shear_modulus,
integer, intent(in) ivisc )

Definition at line 40 of file viscmain.F.

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
my_real
#define my_real
Definition cppsort.cpp:32
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