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) ! ! 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 (for modal ! dynamics: in all nodes for which output ! was requested or a force was applied) ! 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(for modal ! dynamics: in all nodes for which output ! was requested or a force was applied) ! 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 ! ! OUTPUT: ! ! xload concentrated load in direction idof of node ! "node" (global coordinates) !