Next: Distributed flux (dflux.f) Up: User-defined loading Previous: Concentrated flux (cflux.f) Contents

###
Concentrated load (cload.f)

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)
    !