di_fgm_setup
: Domain Integral FGM SETUP¶
Description¶
call di_fgm_setup ( 1, nonode, out, temperatures )
di_fgm_setup:
allocate data structures for two terms used in the calculation of the derivative of the stress work density. these are: nodal values of stress work density and strain.
Output:
j_data.extrap_counts
INTEGER (nonode) ALLOCATABLE SAVE
j_data.swd_at_nodes
(nonode) ALLOCATABLE SAVE
j_data.strain_at_nodes
DOUBLE PRECISION (6,nonode) ALLOCATABLE SAVE
j_data.fgm_e
LOGICAL
j_data.fgm_nu
LOGICAL
Calling Tree¶
c ***************************************************************
c * *
c * -di_fgm_setup *
c * -di_nod_vals *
c * -di_extrap_to_nodes *
c * -ndpts1.f *
c * -oulg1.f (oulgf) *
c * *
c ***************************************************************
Call di_nod_vals
¶
call di_nod_vals( extrap_counts, swd_at_nodes, strain_at_nodes )
build average nodal values of the strain and stress work density for all nodes in the model.
c build average nodal values of the stress work density
c (swd) and strains in the model. for each element,
c shape functions will be used to extrapolate integration-
c point values to each of the element's nodes where they
c are then averaged.
Call di_extrap_to_nodes
¶
subroutine to extrapolate strain and stress work density values from integration points to nodes for one element.
c loop over element nodes. for 8 and 20-noded
c hex elements using 2x2x2 integration, we
c extrapolate to the nodes using the lagrangian
c polynomials. otherwise we just average the
c integration-point values and use that value
c at every node. warp3d uses the following
c arguments to describe elements and integration
c order:
c
c etype = 1: 20-noded brick
c etype = 2: 8-noded brick
c etype = 3: 12-noded brick
c etype = 4: 15-noded brick
c etype = 5: 9-noded brick
c
c etype = 1, int_order = 1: 27 point rule (not used)
c etype = 1, int_order = 8: 8 point rule
c etype = 1, int_order = 9: 14 point rule
c etype = 2, int_order = 1: 8 point rule
c etype = 2, int_order = 2: 6 point rule (not used)