OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
q4kega2.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine q4kega2 (py1, py2, pz1, pz2, ay, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, sig, air, off, nel)

Function/Subroutine Documentation

◆ q4kega2()

subroutine q4kega2 ( py1,
py2,
pz1,
pz2,
ay,
k11,
k12,
k13,
k14,
k22,
k23,
k24,
k33,
k34,
k44,
sig,
air,
off,
integer nel )

Definition at line 30 of file q4kega2.F.

36C-----------------------------------------------
37C I m p l i c i t T y p e s
38C-----------------------------------------------
39#include "implicit_f.inc"
40C-----------------------------------------------
41C G l o b a l P a r a m e t e r s
42C-----------------------------------------------
43#include "mvsiz_p.inc"
44C-----------------------------------------------
45C C o m m o n B l o c k s
46C-----------------------------------------------
47#include "com01_c.inc"
48C-----------------------------------------------
49C D u m m y A r g u m e n t s
50C-----------------------------------------------
51 INTEGER NEL
52 my_real
53 . py1(*), py2(*), pz1(*), pz2(*), ay(*),
54 . k11(2,2,*), k12(2,2,*), k13(2,2,*), k14(2,2,*), k22(2,2,*),
55 . k23(2,2,*), k24(2,2,*), k33(2,2,*), k34(2,2,*), k44(2,2,*),
56 . sig(nel,6), air(*), off(*)
57C-----------------------------------------------
58C L o c a l V a r i a b l e s
59C-----------------------------------------------
60 INTEGER I,J,K
61 my_real
62 . fyy(mvsiz),fzz(mvsiz),fxx(mvsiz),fyz(mvsiz),ays(mvsiz),
63 . gij(mvsiz), hij(2,2,mvsiz)
64C-----------------------------------------------
65C S o u r c e L i n e s
66C-----------------------------------------------
67 DO i=1,nel
68 fyy(i)=sig(i,1)*air(i)
69 fzz(i)=sig(i,2)*air(i)
70 fyz(i)=sig(i,4)*air(i)
71 ENDDO
72 IF (n2d==1) THEN
73 DO i=1,nel
74 fxx(i)=sig(i,3)*air(i)
75 ays(i)=ay(i)*ay(i)
76 ENDDO
77 ENDIF
78C-----K11,K13,K33
79 CALL q4kegaij2(1, nel,
80 . py1, pz1, py1, pz1, ays, gij, hij,
81 . fyy, fzz, fxx, fyz)
82 DO i=1,nel
83 DO j=1,2
84 k11(j,j,i) = k11(j,j,i) + gij(i)
85 k13(j,j,i) = k13(j,j,i) - gij(i)
86 k33(j,j,i) = k33(j,j,i) + gij(i)
87 ENDDO
88 DO j=1,2
89 DO k=1,2
90 k11(k,j,i) = k11(k,j,i) - hij(k,j,i)
91 k13(k,j,i) = k13(k,j,i) + hij(k,j,i)
92 k33(k,j,i) = k33(k,j,i) - hij(k,j,i)
93 ENDDO
94 ENDDO
95 ENDDO
96C-----K22,K24,K44
97 CALL q4kegaij2(1, nel,
98 . py2, pz2, py2, pz2, ays, gij, hij,
99 . fyy, fzz, fxx, fyz)
100 DO i=1,nel
101 DO j=1,2
102 k22(j,j,i) = k22(j,j,i) + gij(i)
103 k24(j,j,i) = k24(j,j,i) - gij(i)
104 k44(j,j,i) = k44(j,j,i) + gij(i)
105 ENDDO
106 DO j=1,2
107 DO k=1,2
108 k22(k,j,i) = k22(k,j,i) - hij(k,j,i)
109 k24(k,j,i) = k24(k,j,i) + hij(k,j,i)
110 k44(k,j,i) = k44(k,j,i) - hij(k,j,i)
111 ENDDO
112 ENDDO
113 ENDDO
114C-----K12,K14,K34,K23
115 CALL q4kegaij2(1, nel,
116 . py1, pz1, py2, pz2, ays, gij, hij,
117 . fyy, fzz, fxx, fyz)
118 DO i=1,nel
119 DO j=1,2
120 k12(j,j,i) = k12(j,j,i) + gij(i)
121 k14(j,j,i) = k14(j,j,i) - gij(i)
122 k34(j,j,i) = k34(j,j,i) + gij(i)
123 k23(j,j,i) = k23(j,j,i) - gij(i)
124 ENDDO
125 DO j=1,2
126 DO k=1,2
127 k12(k,j,i) = k12(k,j,i) - hij(k,j,i)
128 k14(k,j,i) = k14(k,j,i) + hij(k,j,i)
129 k34(k,j,i) = k34(k,j,i) - hij(k,j,i)
130 ENDDO
131 ENDDO
132 k23(1,1,i) = k23(1,1,i) + hij(1,1,i)
133 k23(2,2,i) = k23(2,2,i) + hij(2,2,i)
134 k23(1,2,i) = k23(1,2,i) + hij(2,1,i)
135 k23(2,1,i) = k23(2,1,i) + hij(1,2,i)
136 ENDDO
137C
138 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
q4kegaij2
subroutine q4kegaij2(lft, llt, pyi, pzi, pyj, pzj, ays, gij, hij, fyy, fzz, fxx, fyz)
Definition q4kegaij2.F:31