OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
uroto_tens2d_aniso.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!|| uroto_tens2d_aniso ../engine/source/materials/tools/uroto_tens2d_aniso.F
25!||--- called by ------------------------------------------------------
26!|| h3d_shell_tensor ../engine/source/output/h3d/h3d_results/h3d_shell_tensor.F
27!|| tensorc ../engine/source/output/anim/generate/tensorc.F
28!||====================================================================
29 SUBROUTINE uroto_tens2d_aniso(NEL,TENS,DIR_A,DIR_B)
30C-----------------------------------------------
31C I m p l i c i t T y p e s
32C-----------------------------------------------
33#include "implicit_f.inc"
34C-----------------------------------------------
35C D u m m y A r g u m e n t s
36C-----------------------------------------------
37 INTEGER ,INTENT(IN) :: NEL
38 my_real ,DIMENSION(NEL,2) ,INTENT(IN) :: dir_a,dir_b
39 my_real ,DIMENSION(NEL,3) ,INTENT(INOUT) :: tens
40C-----------------------------------------------
41C L o c a l V a r i a b l e s
42C-----------------------------------------------
43 INTEGER :: I
44 my_real :: r1,r2,s1,s2,rs1,rs2,rs3,r12a,r22a,s12b,s22b,r3r3,s3s3,
45 . t1,t2,t3
46c=======================================================================
47 DO i = 1,nel
48 r1 = dir_a(i,1)
49 s1 = dir_a(i,2)
50 r2 = dir_b(i,1)
51 s2 = dir_b(i,2)
52 rs1 = r1*s1
53 rs2 = r2*s2
54 r12a = r1*r1
55 r22a = r2*r2
56 s12b = s1*s1
57 s22b = s2*s2
58 rs3 = s1*s2-r1*r2
59 r3r3 = one+s1*r2+r1*s2
60 r3r3 = half*r3r3
61 s3s3 = one-s1*r2-r1*s2
62 s3s3 = half*s3s3
63 t1 = tens(i,1)
64 t2 = tens(i,2)
65 t3 = tens(i,3)
66 tens(i,1) = r12a*t1 + r22a*t2 - rs3*t3
67 tens(i,2) = s12b*t1 + s22b*t2 + rs3*t3
68 tens(i,3) = rs1*t1 + rs2*t2 + (r3r3 - s3s3)*t3
69 ENDDO
70c-----------
71 RETURN
72 END
my_real
#define my_real
Definition cppsort.cpp:32
uroto_tens2d_aniso
subroutine uroto_tens2d_aniso(nel, tens, dir_a, dir_b)
Definition uroto_tens2d_aniso.F:30