CalculiX  2.13
A Free Software Three-Dimensional Structural Finite Element Program
utemp_ccxlib.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine utemp_ccxlib (temp, msecpt, kstep, kinc, time, node, coords, vold, mi, iponoel, inoel, ipobody, xbody, ibody)
 

Function/Subroutine Documentation

◆ utemp_ccxlib()

subroutine utemp_ccxlib ( real*8, dimension(msecpt), intent(out)  temp,
integer, intent(in)  msecpt,
integer, intent(in)  kstep,
integer, intent(in)  kinc,
real*8, dimension(2), intent(in)  time,
integer, intent(in)  node,
real*8, dimension(3), intent(in)  coords,
real*8, dimension(0:mi(2),*), intent(in)  vold,
integer, dimension(*), intent(in)  mi,
integer, dimension(*), intent(in)  iponoel,
integer, dimension(2,*), intent(in)  inoel,
integer, dimension(2,*), intent(in)  ipobody,
real*8, dimension(7,*), intent(in)  xbody,
integer, dimension(3,*), intent(in)  ibody 
)
21 !
22 ! user subroutine utemp
23 !
24 !
25 ! INPUT:
26 !
27 ! msecpt number of temperature values (for volume elements:1)
28 ! kstep step number
29 ! kinc increment number
30 ! time(1) current step time
31 ! time(2) current total time
32 ! node node number
33 ! coords(1..3) global coordinates of the node
34 ! vold(0..4,1..nk) solution field in all nodes
35 ! (not available for CFD-calculations)
36 ! 0: temperature
37 ! 1: displacement in global x-direction
38 ! 2: displacement in global y-direction
39 ! 3: displacement in global z-direction
40 ! 4: not used
41 ! mi(1) max # of integration points per element (max
42 ! over all elements)
43 ! mi(2) max degree of freedomm per node (max over all
44 ! nodes) in fields like v(0:mi(2))...
45 ! iponoel(i) the network elements to which node i belongs
46 ! are stored in inoel(1,iponoel(i)),
47 ! inoel(1,inoel(2,iponoel(i)))...... until
48 ! inoel(2,inoel(2,inoel(2......)=0
49 ! inoel(1..2,*) field containing the network elements
50 ! ipobody(1,i) points to an entry in fields ibody and xbody
51 ! containing the body load applied to element i,
52 ! if any, else 0
53 ! ipobody(2,i) index referring to the line in field ipobody
54 ! containing a pointer to the next body load
55 ! applied to element i, else 0
56 ! ibody(1,i) code identifying the kind of body load i:
57 ! 1=centrifugal, 2=gravity, 3=generalized gravity
58 ! ibody(2,i) amplitude number for load i
59 ! ibody(3,i) load case number for load i
60 ! xbody(1,i) size of body load i
61 ! xbody(2..4,i) for centrifugal loading: point on the axis,
62 ! for gravity loading with known gravity vector:
63 ! normalized gravity vector
64 ! xbody(5..7,i) for centrifugal loading: normalized vector on the
65 ! rotation axis
66 !
67 ! OUTPUT:
68 !
69 ! temp(1..msecpt) temperature in the node
70 !
71  implicit none
72 !
73  integer msecpt,kstep,kinc,node,mi(*),iponoel(*),inoel(2,*),
74  & ipobody(2,*),ibody(3,*)
75 !
76  real*8 temp(msecpt),time(2),coords(3),vold(0:mi(2),*),xbody(7,*)
77 !
78  intent(in) msecpt,kstep,kinc,time,node,coords,
79  & vold,mi,iponoel,inoel,ipobody,xbody,ibody
80 !
81  intent(out) temp
82 !
83  temp(1)=293.d0
84 !
85  return
Hosted by OpenAircraft.com, (Michigan UAV, LLC)