OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
r12sumg3.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!|| r12sumg3 ../engine/source/elements/spring/r12sumg3.F
25!||--- called by ------------------------------------------------------
26!|| r12ke3 ../engine/source/elements/spring/r12ke3.F
27!||====================================================================
28 SUBROUTINE r12sumg3(JFT ,JLT ,KX ,DKX ,R11 ,
29 1 R21 ,R31 ,R12 ,R22 ,R32 ,
30 2 KE11 ,KE12 ,KE13 ,KE22 ,KE23 ,
31 3 KE33 )
32C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
33#include "implicit_f.inc"
34#include "mvsiz_p.inc"
35C-----------------------------------------------
36C D U M M Y A R G U M E N T S
37C-----------------------------------------------
38 INTEGER JFT,JLT
39 my_real
40 . KX(*),DKX(*),R11(*),R21(*),R31(*),R12(*),R22(*),R32(*)
41 my_real
42 . ke11(3,3,*),ke22(3,3,*),ke33(3,3,*),ke12(3,3,*),
43 . ke13(3,3,*),ke23(3,3,*)
44C-----------------------------------------------
45C L O C A L V A R I A B L E S
46C-----------------------------------------------
47 INTEGER I, J, EP
48 my_real
49 . K1(3,3,MVSIZ),K2(3,3,MVSIZ),K3(3,3,MVSIZ)
50C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
51C---------------------------------------
52C TRANS LOCAL-->GLOBAL
53C---------------------------------------
54 DO ep=jft,jlt
55 k1(1,1,ep)=r11(ep)*r11(ep)
56 k1(1,2,ep)=r11(ep)*r21(ep)
57 k1(1,3,ep)=r11(ep)*r31(ep)
58 k1(2,2,ep)=r21(ep)*r21(ep)
59 k1(2,3,ep)=r21(ep)*r31(ep)
60 k1(3,3,ep)=r31(ep)*r31(ep)
61 k1(2,1,ep)=k1(1,2,ep)
62 k1(3,2,ep)=k1(2,3,ep)
63 k1(3,1,ep)=k1(1,3,ep)
64 ENDDO
65 DO ep=jft,jlt
66 k2(1,1,ep)=r12(ep)*r12(ep)
67 k2(1,2,ep)=r12(ep)*r22(ep)
68 k2(1,3,ep)=r12(ep)*r32(ep)
69 k2(2,2,ep)=r22(ep)*r22(ep)
70 k2(2,3,ep)=r22(ep)*r32(ep)
71 k2(3,3,ep)=r32(ep)*r32(ep)
72 k2(2,1,ep)=k2(1,2,ep)
73 k2(3,2,ep)=k2(2,3,ep)
74 k2(3,1,ep)=k2(1,3,ep)
75 ENDDO
76 DO ep=jft,jlt
77 k3(1,1,ep)=r11(ep)*r12(ep)
78 k3(1,2,ep)=r11(ep)*r22(ep)
79 k3(1,3,ep)=r11(ep)*r32(ep)
80 k3(2,1,ep)=r21(ep)*r12(ep)
81 k3(2,2,ep)=r21(ep)*r22(ep)
82 k3(2,3,ep)=r21(ep)*r32(ep)
83 k3(3,1,ep)=r31(ep)*r12(ep)
84 k3(3,2,ep)=r31(ep)*r22(ep)
85 k3(3,3,ep)=r31(ep)*r32(ep)
86 ENDDO
87C
88 DO i=1,3
89 DO j=i,3
90 DO ep=jft,jlt
91 ke11(i,j,ep)=kx(ep)*k1(i,j,ep)
92 ke22(i,j,ep)=kx(ep)*(k1(i,j,ep)+k2(i,j,ep)+
93 . k3(i,j,ep)+k3(j,i,ep))
94 ke33(i,j,ep)=kx(ep)*k2(i,j,ep)
95 ENDDO
96 ENDDO
97 ENDDO
98 DO i=1,3
99 DO j=1,3
100 DO ep=jft,jlt
101 ke12(i,j,ep)=-kx(ep)*(k1(i,j,ep)+k3(i,j,ep))
102 ke13(i,j,ep)=kx(ep)*k3(i,j,ep)
103 ke23(i,j,ep)=-kx(ep)*(k2(i,j,ep)+k3(i,j,ep))
104 ENDDO
105 ENDDO
106 ENDDO
107C-----frotement---
108 DO ep=jft,jlt
109 IF (dkx(ep)>zero) THEN
110 DO i=1,3
111 DO j=i,3
112 ke11(i,j,ep)=ke11(i,j,ep)+dkx(ep)*k1(i,j,ep)
113 ke22(i,j,ep)=ke22(i,j,ep)+dkx(ep)*(k1(i,j,ep)-k3(i,j,ep))
114 ke33(i,j,ep)=ke33(i,j,ep)+dkx(ep)*k2(i,j,ep)
115 ke22(i,j,ep)=ke22(i,j,ep)+dkx(ep)*(k2(i,j,ep)-k3(j,i,ep))
116 ENDDO
117 ENDDO
118 DO i=1,3
119 DO j=1,3
120 ke12(i,j,ep)=ke12(i,j,ep)-dkx(ep)*(k1(i,j,ep)-k3(i,j,ep))
121 ke13(i,j,ep)=ke13(i,j,ep)-dkx(ep)*k3(i,j,ep)
122 ke23(i,j,ep)=ke23(i,j,ep)-dkx(ep)*(k2(i,j,ep)-k3(i,j,ep))
123 ENDDO
124 ENDDO
125 ENDIF
126 ENDDO
127C
128 DO i=1,3
129 DO j=i+1,3
130 DO ep=jft,jlt
131 ke11(j,i,ep)=ke11(i,j,ep)
132 ke22(j,i,ep)=ke22(i,j,ep)
133 ke33(j,i,ep)=ke33(i,j,ep)
134 ENDDO
135 ENDDO
136 ENDDO
137C
138 RETURN
139 END
140
r12sumg3
subroutine r12sumg3(jft, jlt, kx, dkx, r11, r21, r31, r12, r22, r32, ke11, ke12, ke13, ke22, ke23, ke33)
Definition r12sumg3.F:32