GLOBAL_PROBLEM_PLANESTRESSHARDMIXVAF - solves the global finite element problem

Comments

This function solves the global finite element problem, i.e., the residuals of the discretized weak formulation are minimized. Assumptions: plane stress, isotropic hardening (Voce), kinematic hardening (Armstrong-Frederick).

Input Arguments

n_save_YieldSurface (double) -

n_NR_global (double) - maximum number of Newton-Raphson iterations used for solving the global problem

tol_NR_global (double) - stopping tolerance for the Newton-Raphson iteration used for solving the global problem

n_NR_local (double) - maximum number of Newton-Raphson iterations used for solving the local problem (return mapping algorithm)

tol_NR_local (double) - stopping tolerance for the Newton-Raphson iteration used for solving the local problem (return mapping algorithm)

CPlaneStress (double) - elastic properties (stiffness matrix)

SPlaneStress (double) - elastic properties (compliance matrix)

theta (double) - material parameters

H_isotropic (double) - isotropic hardening parameters

H_kinematic (double) - kinematic hardening parameters

n_dof (double) - see mesh_info

n_element (double) - see mesh_info

n_dof_per_element (double) - see mesh_info

element_dof (double) - see mesh_info

Gauss_weights (double) - see mesh_info

n_Gauss_per_dim (double) - see mesh_info

detJ_GP (double) - see mesh_info

B_GP (double) - see mesh_info

dof_fix (double) - fixed degrees of freedom

dof_displacement (double) - degrees of freedom under displacement control

dof_reaction (double) - degrees of freedom associated with the reaction forces

max_displacement (double) - maximum applied displacement

n_step (double) - number of load steps

displacement_increment (double) - applied displacement increment

Output Arguments

u (double) - displacement at each node for each load step

reaction (double) - reaction forces for each load step

epsilonVp_GP (double) - plastic component of the infinitesimal strain at each Gauss point for the last load step

gamma_GP (double) - plastic multiplier at each Gauss point for the last load step

sigmaV_back_GP (double) - back stress at each Gauss point for the last load step

YieldSurface (double) - stresses at which yielding has occured