This subroutine is used for user-defined concentrated load, characterized by the parameter USER on the *CLOAD card. The header and variable description is as follows:
subroutine cload(xload,kstep,kinc,time,node,idof,coords,vold,
& mi,ntrans,trab,inotr,veold,nmethod,nactdof,bcont,fn)
!
! user subroutine cload
!
!
! INPUT:
!
! kstep step number
! kinc increment number
! time(1) current step time
! time(2) current total time
! node node number
! idof degree of freedom
! coords(1..3) global coordinates of the node
! vold(0..mi(2)
! ,1..nk) solution field in all nodes
! 0: temperature
! 1: displacement in global x-direction
! 2: displacement in global y-direction
! 3: displacement in global z-direction
! 4: static pressure
! mi(1) max # of integration points per element (max
! over all elements)
! mi(2) max degree of freedomm per node (max over all
! nodes) in fields like v(0:mi(2))...
! veold(0..3,1..nk) derivative of the solution field w.r.t.
! time in all nodes
! 0: temperature rate
! 1: velocity in global x-direction
! 2: velocity in global y-direction
! 3: velocity in global z-direction
! ntrans number of transform definitions
! trab(1..6,i) coordinates of two points defining transform i
! trab(7,i) -1: cylindrical transformation
! 1: rectangular transformation
! inotr(1,j) transformation number applied to node j
! inotr(2,j) a SPC in a node j in which a transformation
! applied corresponds to a MPC. inotr(2,j)
! contains the number of a new node generated
! for the inhomogeneous part of the MPC
! nmethod kind of procedure
! -1: visco
! 0: no analysis
! 1: static
! 2: frequency
! 3: buckling
! 4: modal dynamic
! 5: modal steady state dynamics
! 6: matrix storage
! nactdof(i,j) number of the degree of freedom in the global
! system of equations of local degree of freedom
! i (0<=i<=mi(2)) in node j; this field is only
! accessible for nmethod=4, else a segmentation
! fault may result
! bcont(i) contact force in global degree of freedom i:
! this option is only available for modal dynamic
! calculations (nmethod=4). In all other cases use
! of this field may lead to a segmentation fault
! fn(0..mi(2)
! ,1..nk) reaction force in all nodes
! 0: concentrated reaction flux
! 1: reaction force in global x-direction
! 2: reaction force in global y-direction
! 3: reaction force in global z-direction
! this option is only available for modal dynamic
! calculations (nmethod=4). In all other cases use
! of this field may lead to a segmentation fault
!
!
! OUTPUT:
!
! xload concentrated load in direction idof of node
! "node" (global coordinates)
!