OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sfor_visn4.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!|| sfor_visn4 ../engine/source/elements/solid/solide4/sfor_visn4.F
25!||--- called by ------------------------------------------------------
26!|| s4for_distor ../engine/source/elements/solid/solide4/s4for_distor.F
27!||====================================================================
28 SUBROUTINE sfor_visn4(VC , FLD, TOL_V, IFC1,
29 . VX1, VX2, VX3, VX4,
30 . VY1, VY2, VY3, VY4,
31 . VZ1, VZ2, VZ3, VZ4,
32 . FOR_T1, FOR_T2, FOR_T3, FOR_T4,
33 . IFCTL, STIF , MU , NEL ,
34 . E_DISTOR,DT1)
35C-----------------------------------------------
36C I m p l i c i t T y p e s
37C-----------------------------------------------
38#include "implicit_f.inc"
39C-----------------------------------------------
40C G l o b a l P a r a m e t e r s
41C-----------------------------------------------
42#include "mvsiz_p.inc"
43C-----------------------------------------------
44C D u m m y A r g u m e n t s
45C-----------------------------------------------
46 INTEGER, INTENT (IN) :: NEL
47 INTEGER, INTENT (OUT) :: IFCTL
48 INTEGER, DIMENSION(MVSIZ),INTENT (INOUT) :: IFC1
49 my_real, DIMENSION(MVSIZ), INTENT (INOUT) :: STIF
50 my_real, DIMENSION(MVSIZ), INTENT (IN) :: FLD ,
51 4 VX1, VX2, VX3, VX4,
52 6 vy1, vy2, vy3, vy4,
53 8 vz1, vz2, vz3, vz4
54 my_real, DIMENSION(MVSIZ,3), INTENT (IN) :: vc
55 my_real, DIMENSION(MVSIZ,3), INTENT (INOUT) ::
56 . for_t1, for_t2, for_t3, for_t4
57 my_real, INTENT (IN) :: tol_v,mu,dt1
58 my_real, DIMENSION(NEL), INTENT (INOUT) :: e_distor
59C-----------------------------------------------
60C C o m m o n B l o c k s
61C-----------------------------------------------
62C-----------------------------------------------
63C L o c a l V a r i a b l e s
64C-----------------------------------------------
65 INTEGER I,J
66C 12
67 my_real
68 . VXC,VYC,VZC,FX,FY,FZ,FAC,VNJ(4),VL,TOL_V2,
69 . VXMAX,VYMAX,VZMAX,V2MAX,VC2,LAM_MIN,LAM_MIN2,LAM_MAX
70C----------------------------
71 tol_v2 = tol_v*tol_v
72 ifctl = 0
73 DO i=1,nel
74 vc2 = vc(i,1)*vc(i,1)+vc(i,2)*vc(i,2)+vc(i,3)*vc(i,3)
75 IF (vc2 <em20.OR.stif(i)==zero) cycle
76 vl = tol_v2*vc2
77 vnj(1) = vx1(i)*vx1(i) + vy1(i)*vy1(i) + vz1(i)*vz1(i)
78 vnj(2) = vx2(i)*vx2(i) + vy2(i)*vy2(i) + vz2(i)*vz2(i)
79 vnj(3) = vx3(i)*vx3(i) + vy3(i)*vy3(i) + vz3(i)*vz3(i)
80 vnj(4) = vx4(i)*vx4(i) + vy4(i)*vy4(i) + vz4(i)*vz4(i)
81 v2max = max(vnj(1),vnj(2),vnj(3),vnj(4))
82 IF (v2max > vl) ifc1(i) = 1
83 IF (ifc1(i) > 0) ifctl=1
84 END DO
85C
86 IF (ifctl==1) THEN
87 fac = one + two*mu
88 DO i=1,nel
89 IF (ifc1(i)==0) cycle
90 for_t1(i,1) = for_t1(i,1) - fld(i)*(vx1(i)-vc(i,1))
91 for_t1(i,2) = for_t1(i,2) - fld(i)*(vy1(i)-vc(i,2))
92 for_t1(i,3) = for_t1(i,3) - fld(i)*(vz1(i)-vc(i,3))
93 for_t2(i,1) = for_t2(i,1) - fld(i)*(vx2(i)-vc(i,1))
94 for_t2(i,2) = for_t2(i,2) - fld(i)*(vy2(i)-vc(i,2))
95 for_t2(i,3) = for_t2(i,3) - fld(i)*(vz2(i)-vc(i,3))
96 for_t3(i,1) = for_t3(i,1) - fld(i)*(vx3(i)-vc(i,1))
97 for_t3(i,2) = for_t3(i,2) - fld(i)*(vy3(i)-vc(i,2))
98 for_t3(i,3) = for_t3(i,3) - fld(i)*(vz3(i)-vc(i,3))
99 for_t4(i,1) = for_t4(i,1) - fld(i)*(vx4(i)-vc(i,1))
100 for_t4(i,2) = for_t4(i,2) - fld(i)*(vy4(i)-vc(i,2))
101 for_t4(i,3) = for_t4(i,3) - fld(i)*(vz4(i)-vc(i,3))
102 stif(i) = fac*stif(i)
103 e_distor(i)=e_distor(i)- dt1*(for_t1(i,1)*(vx1(i)-vc(i,1))+
104 . for_t1(i,2)*(vy1(i)-vc(i,2))+
105 . for_t1(i,3)*(vz1(i)-vc(i,3))+
106 . for_t2(i,1)*(vx2(i)-vc(i,1))+
107 . for_t2(i,2)*(vy2(i)-vc(i,2))+
108 . for_t2(i,3)*(vz2(i)-vc(i,3))+
109 . for_t3(i,1)*(vx3(i)-vc(i,1))+
110 . for_t3(i,2)*(vy3(i)-vc(i,2))+
111 . for_t3(i,3)*(vz3(i)-vc(i,3))+
112 . for_t4(i,1)*(vx4(i)-vc(i,1))+
113 . for_t4(i,2)*(vy4(i)-vc(i,2))+
114 . for_t4(i,3)*(vz4(i)-vc(i,3)))
115 ENDDO
116 END IF
117C
118 RETURN
119 END
max
#define max(a, b)
Definition macros.h:21
sfor_visn4
subroutine sfor_visn4(vc, fld, tol_v, ifc1, vx1, vx2, vx3, vx4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, for_t1, for_t2, for_t3, for_t4, ifctl, stif, mu, nel, e_distor, dt1)
Definition sfor_visn4.F:35