Solid mechanics conforming fractures solver

Introduction

Theory

Under construction

Governing Equations

Under construction

Parameters

In the preceding XML block, The SolidMechanicsLagrangeContact is specified by the title of the subblock of the Solvers block. The following attributes are supported in the input block for SolidMechanicsLagrangeContact:

Name

Type

Default

Description

cflFactor

real64

0.5

Factor to apply to the CFL condition when calculating the maximum allowable time step. Values should be in the interval (0,1]

discretization

groupNameRef

required

Name of discretization object (defined in the Numerical Methods) to use for this solver. For instance, if this is a Finite Element Solver, the name of a Finite Element Discretization should be specified. If this is a Finite Volume Method, the name of a Finite Volume Discretization discretization should be specified.

initialDt

real64

1e+99

Initial time-step value required by the solver to the event manager.

logLevel

integer

0

Log level

massDamping

real64

0

Value of mass based damping coefficient.

maxNumResolves

integer

10

Value to indicate how many resolves may be executed after some other event is executed. For example, if a SurfaceGenerator is specified, it will be executed after the mechanics solve. However if a new surface is generated, then the mechanics solve must be executed again due to the change in topology.

name

groupName

required

A name is required for any non-unique nodes

newmarkBeta

real64

0.25

Value of \beta in the Newmark Method for Implicit Dynamic time integration option. This should be pow(newmarkGamma+0.5,2.0)/4.0 unless you know what you are doing.

newmarkGamma

real64

0.5

Value of \gamma in the Newmark Method for Implicit Dynamic time integration option

stabilizationName

groupNameRef

required

Name of the stabilization to use in the lagrangian contact solver

stabilizationScalingCoefficient

real64

1

It be used to increase the scale of the stabilization entries. A value < 1.0 results in larger entries in the stabilization matrix.

stiffnessDamping

real64

0

Value of stiffness based damping coefficient.

strainTheory

integer

0

Indicates whether or not to use Infinitesimal Strain Theory, or Finite Strain Theory. Valid Inputs are:
0 - Infinitesimal Strain
1 - Finite Strain

targetRegions

groupNameRef_array

required

Allowable regions that the solver may be applied to. Note that this does not indicate that the solver will be applied to these regions, only that allocation will occur such that the solver may be applied to these regions. The decision about what regions this solver will beapplied to rests in the EventManager.

timeIntegrationOption

geos_SolidMechanicsLagrangianFEM_TimeIntegrationOption

ExplicitDynamic

Time integration method. Options are:
* QuasiStatic
* ImplicitDynamic
* ExplicitDynamic

LinearSolverParameters

node

unique

Element: LinearSolverParameters

NonlinearSolverParameters

node

unique

Element: NonlinearSolverParameters

The following data are allocated and used by the solver:

Name

Type

Description

contactRelationName

groupNameRef

Name of contact relation to enforce constraints on fracture boundary.

maxForce

real64

The maximum force contribution in the problem domain.

maxStableDt

real64

Value of the Maximum Stable Timestep for this solver.

meshTargets

geos_mapBase<std_pair<string, string >, LvArray_Array<string, 1, camp_int_seq<long, 0l>, int, LvArray_ChaiBuffer>, std_integral_constant<bool, true> >

MeshBody/Region combinations that the solver will be applied to.

surfaceGeneratorName

string

Name of the surface generator to use

LinearSolverParameters

node

Datastructure: LinearSolverParameters

NonlinearSolverParameters

node

Datastructure: NonlinearSolverParameters

SolverStatistics

node

Datastructure: SolverStatistics