CalculiX  2.13
A Free Software Three-Dimensional Structural Finite Element Program
Functions/Subroutines
lintemp_th.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine lintemp_th (t0, vold, konl, nope, jj, t0l, t1l, mi)
 

Function/Subroutine Documentation

◆ lintemp_th()

subroutine lintemp_th ( real*8, dimension(*), intent(in)  t0,
real*8, dimension(0:mi(2),*), intent(in)  vold,
integer, dimension(20), intent(in)  konl,
integer, intent(in)  nope,
integer, intent(in)  jj,
real*8, intent(inout)  t0l,
real*8, intent(inout)  t1l,
integer, dimension(*), intent(in)  mi 
)
20 !
21 ! calculates a trilinear approximation to the quadratic interpolation
22 ! of the temperatures in a C3D20 element (full integration). A
23 ! quadratic interpolation of the temperatures leads to quadratic
24 ! thermal stresses, which cannot be handled by the elements
25 ! displacement functions (which lead to linear stresses). Thus,
26 ! the temperatures are approximated by a trilinear function.
27 !
28  implicit none
29 !
30  integer konl(20),nope,jj,i1,mi(*)
31 !
32  real*8 t0(*),vold(0:mi(2),*),t0l,t1l,a20l(20,27)
33 !
34  intent(in) t0,vold,konl,nope,jj,mi
35 !
36  intent(inout) t0l,t1l
37 !
38  a20l=reshape((
39  &/-0.088729832,-0.240369600,-0.059630393,-0.240369600,-0.240369600,
40  & -0.059630393,-0.011270159,-0.059630393, 0.524865555, 0.066666663,
41  & 0.066666663, 0.524865555, 0.066666663, 0.008467776, 0.008467776,
42  & 0.066666663, 0.524865555, 0.066666663, 0.008467776, 0.066666663,
43  & -0.164549715,-0.164549715,-0.149999995,-0.149999995,-0.149999995,
44  & -0.149999995,-0.035450279,-0.035450279, 0.524865544, 0.295766106,
45  & 0.066666668, 0.295766106, 0.066666668, 0.037567223, 0.008467777,
46  & 0.037567223, 0.295766106, 0.295766106, 0.037567223, 0.037567223,
47  & -0.240369600,-0.088729832,-0.240369600,-0.059630393,-0.059630393,
48  & -0.240369600,-0.059630393,-0.011270159, 0.524865555, 0.524865555,
49  & 0.066666663, 0.066666663, 0.066666663, 0.066666663, 0.008467776,
50  & 0.008467776, 0.066666663, 0.524865555, 0.066666663, 0.008467776,
51  & -0.164549715,-0.149999995,-0.149999995,-0.164549715,-0.149999995,
52  & -0.035450279,-0.035450279,-0.149999995, 0.295766106, 0.066666668,
53  & 0.295766106, 0.524865544, 0.037567223, 0.008467777, 0.037567223,
54  & 0.066666668, 0.295766106, 0.037567223, 0.037567223, 0.295766106,
55  & -0.157274855,-0.157274855,-0.157274855,-0.157274855,-0.092725137,
56  & -0.092725137,-0.092725137,-0.092725137, 0.295766105, 0.295766105,
57  & 0.295766105, 0.295766105, 0.037567224, 0.037567224, 0.037567224,
58  & 0.037567224, 0.166666664, 0.166666664, 0.166666664, 0.166666664,
59  & -0.149999995,-0.164549715,-0.164549715,-0.149999995,-0.035450279,
60  & -0.149999995,-0.149999995,-0.035450279, 0.295766106, 0.524865544,
61  & 0.295766106, 0.066666668, 0.037567223, 0.066666668, 0.037567223,
62  & 0.008467777, 0.037567223, 0.295766106, 0.295766106, 0.037567223,
63  & -0.240369600,-0.059630393,-0.240369600,-0.088729832,-0.059630393,
64  & -0.011270159,-0.059630393,-0.240369600, 0.066666663, 0.066666663,
65  & 0.524865555, 0.524865555, 0.008467776, 0.008467776, 0.066666663,
66  & 0.066666663, 0.066666663, 0.008467776, 0.066666663, 0.524865555,
67  & -0.149999995,-0.149999995,-0.164549715,-0.164549715,-0.035450279,
68  & -0.035450279,-0.149999995,-0.149999995, 0.066666668, 0.295766106,
69  & 0.524865544, 0.295766106, 0.008467777, 0.037567223, 0.066666668,
70  & 0.037567223, 0.037567223, 0.037567223, 0.295766106, 0.295766106,
71  & -0.059630393,-0.240369600,-0.088729832,-0.240369600,-0.011270159,
72  & -0.059630393,-0.240369600,-0.059630393, 0.066666663, 0.524865555,
73  & 0.524865555, 0.066666663, 0.008467776, 0.066666663, 0.066666663,
74  & 0.008467776, 0.008467776, 0.066666663, 0.524865555, 0.066666663,
75  & -0.164549715,-0.149999995,-0.035450279,-0.149999995,-0.164549715,
76  & -0.149999995,-0.035450279,-0.149999995, 0.295766106, 0.037567223,
77  & 0.037567223, 0.295766106, 0.295766106, 0.037567223, 0.037567223,
78  & 0.295766106, 0.524865544, 0.066666668, 0.008467777, 0.066666668,
79  & -0.157274855,-0.157274855,-0.092725137,-0.092725137,-0.157274855,
80  & -0.157274855,-0.092725137,-0.092725137, 0.295766105, 0.166666664,
81  & 0.037567224, 0.166666664, 0.295766105, 0.166666664, 0.037567224,
82  & 0.166666664, 0.295766105, 0.295766105, 0.037567224, 0.037567224,
83  & -0.149999995,-0.164549715,-0.149999995,-0.035450279,-0.149999995,
84  & -0.164549715,-0.149999995,-0.035450279, 0.295766106, 0.295766106,
85  & 0.037567223, 0.037567223, 0.295766106, 0.295766106, 0.037567223,
86  & 0.037567223, 0.066666668, 0.524865544, 0.066666668, 0.008467777,
87  & -0.157274855,-0.092725137,-0.092725137,-0.157274855,-0.157274855,
88  & -0.092725137,-0.092725137,-0.157274855, 0.166666664, 0.037567224,
89  & 0.166666664, 0.295766105, 0.166666664, 0.037567224, 0.166666664,
90  & 0.295766105, 0.295766105, 0.037567224, 0.037567224, 0.295766105,
91  & -0.124999996,-0.124999996,-0.124999996,-0.124999996,-0.124999996,
92  & -0.124999996,-0.124999996,-0.124999996, 0.166666664, 0.166666664,
93  & 0.166666664, 0.166666664, 0.166666664, 0.166666664, 0.166666664,
94  & 0.166666664, 0.166666664, 0.166666664, 0.166666664, 0.166666664,
95  & -0.092725137,-0.157274855,-0.157274855,-0.092725137,-0.092725137,
96  & -0.157274855,-0.157274855,-0.092725137, 0.166666664, 0.295766105,
97  & 0.166666664, 0.037567224, 0.166666664, 0.295766105, 0.166666664,
98  & 0.037567224, 0.037567224, 0.295766105, 0.295766105, 0.037567224,
99  & -0.149999995,-0.035450279,-0.149999995,-0.164549715,-0.149999995,
100  & -0.035450279,-0.149999995,-0.164549715, 0.037567223, 0.037567223,
101  & 0.295766106, 0.295766106, 0.037567223, 0.037567223, 0.295766106,
102  & 0.295766106, 0.066666668, 0.008467777, 0.066666668, 0.524865544,
103  & -0.092725137,-0.092725137,-0.157274855,-0.157274855,-0.092725137,
104  & -0.092725137,-0.157274855,-0.157274855, 0.037567224, 0.166666664,
105  & 0.295766105, 0.166666664, 0.037567224, 0.166666664, 0.295766105,
106  & 0.166666664, 0.037567224, 0.037567224, 0.295766105, 0.295766105,
107  & -0.035450279,-0.149999995,-0.164549715,-0.149999995,-0.035450279,
108  & -0.149999995,-0.164549715,-0.149999995, 0.037567223, 0.295766106,
109  & 0.295766106, 0.037567223, 0.037567223, 0.295766106, 0.295766106,
110  & 0.037567223, 0.008467777, 0.066666668, 0.524865544, 0.066666668,
111  & -0.240369600,-0.059630393,-0.011270159,-0.059630393,-0.088729832,
112  & -0.240369600,-0.059630393,-0.240369600, 0.066666663, 0.008467776,
113  & 0.008467776, 0.066666663, 0.524865555, 0.066666663, 0.066666663,
114  & 0.524865555, 0.524865555, 0.066666663, 0.008467776, 0.066666663,
115  & -0.149999995,-0.149999995,-0.035450279,-0.035450279,-0.164549715,
116  & -0.164549715,-0.149999995,-0.149999995, 0.066666668, 0.037567223,
117  & 0.008467777, 0.037567223, 0.524865544, 0.295766106, 0.066666668,
118  & 0.295766106, 0.295766106, 0.295766106, 0.037567223, 0.037567223,
119  & -0.059630393,-0.240369600,-0.059630393,-0.011270159,-0.240369600,
120  & -0.088729832,-0.240369600,-0.059630393, 0.066666663, 0.066666663,
121  & 0.008467776, 0.008467776, 0.524865555, 0.524865555, 0.066666663,
122  & 0.066666663, 0.066666663, 0.524865555, 0.066666663, 0.008467776,
123  & -0.149999995,-0.035450279,-0.035450279,-0.149999995,-0.164549715,
124  & -0.149999995,-0.149999995,-0.164549715, 0.037567223, 0.008467777,
125  & 0.037567223, 0.066666668, 0.295766106, 0.066666668, 0.295766106,
126  & 0.524865544, 0.295766106, 0.037567223, 0.037567223, 0.295766106,
127  & -0.092725137,-0.092725137,-0.092725137,-0.092725137,-0.157274855,
128  & -0.157274855,-0.157274855,-0.157274855, 0.037567224, 0.037567224,
129  & 0.037567224, 0.037567224, 0.295766105, 0.295766105, 0.295766105,
130  & 0.295766105, 0.166666664, 0.166666664, 0.166666664, 0.166666664,
131  & -0.035450279,-0.149999995,-0.149999995,-0.035450279,-0.149999995,
132  & -0.164549715,-0.164549715,-0.149999995, 0.037567223, 0.066666668,
133  & 0.037567223, 0.008467777, 0.295766106, 0.524865544, 0.295766106,
134  & 0.066666668, 0.037567223, 0.295766106, 0.295766106, 0.037567223,
135  & -0.059630393,-0.011270159,-0.059630393,-0.240369600,-0.240369600,
136  & -0.059630393,-0.240369600,-0.088729832, 0.008467776, 0.008467776,
137  & 0.066666663, 0.066666663, 0.066666663, 0.066666663, 0.524865555,
138  & 0.524865555, 0.066666663, 0.008467776, 0.066666663, 0.524865555,
139  & -0.035450279,-0.035450279,-0.149999995,-0.149999995,-0.149999995,
140  & -0.149999995,-0.164549715,-0.164549715, 0.008467777, 0.037567223,
141  & 0.066666668, 0.037567223, 0.066666668, 0.295766106, 0.524865544,
142  & 0.295766106, 0.037567223, 0.037567223, 0.295766106, 0.295766106,
143  & -0.011270159,-0.059630393,-0.240369600,-0.059630393,-0.059630393,
144  & -0.240369600,-0.088729832,-0.240369600, 0.008467776, 0.066666663,
145  & 0.066666663, 0.008467776, 0.066666663, 0.524865555, 0.524865555,
146  & 0.066666663, 0.008467776, 0.066666663, 0.524865555, 0.066666663/
147  & ),(/20,27/))
148 !
149  do i1=1,nope
150  t0l=t0l+a20l(i1,jj)*t0(konl(i1))
151  t1l=t1l+a20l(i1,jj)*vold(0,konl(i1))
152  enddo
153 !
154  return
Hosted by OpenAircraft.com, (Michigan UAV, LLC)