Base class for micromechanical models. More...

#include <micromechanical.hpp>

Inheritance diagram for micromechanical::core::MicromechanicalBase< T >:

[legend]

Collaboration diagram for micromechanical::core::MicromechanicalBase< T >:

[legend]

Public Member Functions
	MicromechanicalBase (const std::variant< ContactLawType, std::string, ContactLawBase< T > * > &contact, const std::unordered_map< std::string, T > &props={}, T npv=1e9, T radius=0.65e-3)
	constructor
	MicromechanicalBase ()=default
	default constructor
void	setContactLaw (const std::variant< ContactLawType, std::string, ContactLawBase< T > * > &contact, const std::unordered_map< std::string, T > &props={})
	Set the contact law if it is not set in the constructor.
virtual Eigen::Matrix< T, 6, 6 >	stiffness1 (const StateVariable< T > &sv0, StateVariable< T > &sv, const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &kn=std::nullopt, const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &ks=std::nullopt, const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &kt=std::nullopt)=0
	Compute the stiffness matrix.
virtual Eigen::Matrix< T, Eigen::Dynamic, 3 >	localizeStress (const Eigen::Matrix< T, 6, 1 > &sig, const StateVariable< T > &sv0, StateVariable< T > &sv)
	localize the macro stress to micro contact force
virtual Eigen::Matrix< T, Eigen::Dynamic, 3 >	localizeStress2 (const Eigen::Matrix< T, 6, 1 > &sig, const StateVariable< T > &sv0, StateVariable< T > &sv, const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &kn=std::nullopt, const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &ks=std::nullopt, const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &kt=std::nullopt)
	localize the macro stress to micro contact force with the kinematic hypothesis, estimated from the strain and contact displacement.
virtual Eigen::Matrix< T, Eigen::Dynamic, 3 >	localizeStrain (const Eigen::Matrix< T, 6, 1 > &eps, const StateVariable< T > &sv0, StateVariable< T > &sv)
	localize the macro strain to micro contact displacement
virtual Eigen::Matrix< T, Eigen::Dynamic, 3 >	localizeStrain2 (const Eigen::Matrix< T, 6, 1 > &eps, const StateVariable< T > &sv0, StateVariable< T > &sv, const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &kn=std::nullopt, const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &ks=std::nullopt, const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &kt=std::nullopt)
	localize the macro strain to micro contact displacement with the static hypothesis, estimated from the stress and contact force.
virtual Eigen::Matrix< T, 6, 1 >	averageStress (const Eigen::Matrix< T, Eigen::Dynamic, 3 > force, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Average the micro contact force to macro stress.
virtual Eigen::Matrix< T, 6, 1 >	averageStrain (const Eigen::Matrix< T, Eigen::Dynamic, 3 > disp, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Average the micro contact displacement to macro strain.
Public Member Functions inherited from constitutive::ElastoplasticModelBase< T >
void	setPlasticMethod (const std::string &method)
	set the plastic method
Eigen::Matrix< T, 6, 6 >	continuumElastoplasticStiffness (const Eigen::Matrix< T, 6, 6 > &K, const StateVariable< T > &sv0, StateVariable< T > &sv, bool check_elastic_state=true, std::optional< bool > check_positive_definite=std::nullopt, std::optional< bool > check_non_singularity=std::nullopt) override
Eigen::Matrix< T, 6, 6 >	symmetricContinuumElastoplasticStiffness (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv, bool check_elastic_state=true, std::optional< bool > check_positive_definite=std::nullopt, std::optional< bool > check_non_singularity=std::nullopt) override
Eigen::Matrix< T, 6, 6 >	consistentElastoplasticStiffness (const Eigen::Matrix< T, 6, 6 > &K, const StateVariable< T > &sv0, StateVariable< T > &sv, bool check_elastic_state=true, std::optional< bool > check_positive_definite=std::nullopt, std::optional< bool > check_non_singularity=std::nullopt) override
virtual T	yieldFunction (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Yield function for elastoplastic model.
Eigen::Matrix< T, 6, 1 >	dfdsig0 (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the yield function with respect to the stress.
virtual Eigen::Matrix< T, 6, 1 >	dfdsig (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the yield function with respect to the stress.
Eigen::Matrix< T, 6, 1 >	dfdsig1 (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the yield function with respect to the stress.
Eigen::Matrix< T, 6, 1 >	dfdeps0 (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the yield function with respect to the strain.
virtual Eigen::Matrix< T, 6, 1 >	dfdeps (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the yield function with respect to the strain.
Eigen::Matrix< T, 6, 1 >	dfdeps1 (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the yield function with respect to the strain.
Eigen::Matrix< T, 6, 1 >	dfdepsp0 (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the yield function with respect to the plastic strain.
virtual Eigen::Matrix< T, 6, 1 >	dfdepsp (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the yield function with respect to the plastic strain.
Eigen::Matrix< T, 6, 1 >	dfdepsp1 (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the yield function with respect to the plastic strain.
Public Member Functions inherited from constitutive::ConstitutiveModelBase< T >
	ConstitutiveModelBase ()=default
	default constructor
virtual	~ConstitutiveModelBase ()=default
	destructor
virtual std::vector< StateVariable< T > >	run ()
	Run the model.
virtual Eigen::Matrix< T, 6, 1 >	increment (const Eigen::Matrix< T, 6, 1 > &dx, Eigen::Matrix< T, 6, 1 > x, const Eigen::Matrix< T, 6, 6 > &S, const Eigen::Matrix< T, 6, 6 > &E, const Eigen::Matrix< T, 6, 6 > &G, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Compute one increment.
void	setElementType (const std::string &element_type)
	Set ABAQUS element type.
void	setStressStrainComponents (const int ndi, const int nshr)
	Set number of stress/strain components.
std::tuple< StateVariable< T >, StateVariable< T > >	umat (const T dload, T stress, T statev, T ddsdde=nullptr, T sse=nullptr, T spd=nullptr, T *scd=nullptr)
	ABAQUS UMAT interface.
void	vumat (const int nBlock, const T loadInc, const T stressOld, T stressNew, const T stateOld, T *stateNew, const int nStatev)
	ABAQUS VUMAT interface.
virtual Eigen::Matrix< T, 6, 6 >	stiffness (const StateVariable< T > &sv0, StateVariable< T > &sv)=0
	Compute the stiffness matrix.
Eigen::Matrix< T, 6, 6 >	elastoplasticStiffness (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > deps, const StateVariable< T > &sv0, StateVariable< T > &sv, bool check_elastic_state=true, std::optional< bool > check_positive_definite=std::nullopt, std::optional< bool > check_non_singularity=std::nullopt)
	Compute the elastoplastic stiffness matrix.

Public Attributes
std::shared_ptr< ContactLawBase< T > >	contact
	contact law
std::shared_ptr< Integration< T > >	integration
	gauss integration points
int	_integration = 0
	number of integration steps
Public Attributes inherited from constitutive::ElastoplasticModelBase< T >
std::string	plasticMethod = "CPA"
	Plastic method.
Public Attributes inherited from constitutive::ConstitutiveModelBase< T >
std::string	name = "Undefined"
	name of the model
std::shared_ptr< Load< T > >	load = std::make_shared<Load<T>>()
	load
std::shared_ptr< Options< T > >	options = std::make_shared<Options<T>>()
	options for the constitutive model
std::shared_ptr< DebugVariable< T > >	debugVariable = std::make_shared<DebugVariable<T>>()
	debug variable
std::vector< StateVariable< T > >	stateVars = {}
	state variables of the current run
Eigen::Matrix< T, 3, 3 >	custom_loading_direction = Eigen::Matrix<T, 3, 3>::Zero()
	custom loading direction
bool	abaqus_initialize = false
	Options for abaqus UMAT/VUMAT interfaces.
bool	abaqus_normalize_ddsdde = true
	whether to normalize stiffness matrix
bool	abaqus_strain_control = true
	the load is controlled by strain (true) or stress (false)
bool	abaqus_disable_unused_components = false
	whether to disable unused components in the flow rule
int	abaqus_ndi = 3
	number of direct stress/strain components
int	abaqus_nshr = 3
	number of shear stress/strain components
std::string	abaqus_stiffness_type = "elastic"
T	abaqus_elastoplastic_stiffness_strain_increment_multiplier = 1.0
	multiplier for strain increment when computing elastoplastic stiffness matrix
bool	abaqus_check_stiffness_matrix_positive_definite = true
	whether to check the positive definiteness of the stiffness matrix
bool	abaqus_check_stiffness_matrix_non_singularity = true
	whether to check the non-singularity of the stiffness matrix
int	_phase = 0
	index of the phase
int	_step = 0
	current step
int	_mixedLoad = 0
	number of mixed load steps

Protected Member Functions
int	nstatev () const
void	fromAbaqus (StateVariable< T > &sv0, T *statev) override
	Sync the state variables sv0 from abaqus statev.
void	toAbaqus (const StateVariable< T > &sv0, StateVariable< T > &sv, T *statev) override
	Sync the state variables sv to abaqus statev.
StateVariable< T >	createStateVariable () override
	Create a new state variable object.
void	initialize (StateVariable< T > &sv) override
	Initialize the state variables.
void	postProcess (StateVariable< T > &sv) override
	Post-process the state variables.
void	preStressStrainIntegrate (const StateVariable< T > &sv0, StateVariable< T > &sv) override
	Pre-processing before the stress-strain integration.
void	postIncrement (const StateVariable< T > &sv0, StateVariable< T > &sv) override
	Post-processing after the increment.
void	contactIntegrates (const Eigen::Matrix< T, Eigen::Dynamic, 3 > &dforce, const Eigen::Matrix< T, Eigen::Dynamic, 3 > &ddisp, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Integrate all the contact laws with cache.
virtual Eigen::Matrix< T, 3, 3 >	contactIntegrate (const int idx, const Eigen::Matrix< T, 3, 1 > &dforce, const Eigen::Matrix< T, 3, 1 > &ddisp, const StateVariable< T > &sv0, StateVariable< T > &sv) const =0
	Integrate the contact law.
Protected Member Functions inherited from constitutive::ElastoplasticModelBase< T >
bool	isExplicit () const
	check if the plastic method is explicit
bool	isCPA () const
	check if the plastic method is CPA
bool	isCPPM () const
	check if the plastic method is CPPM
bool	isRungeKutta () const
	check if the plastic method is RungeKutta
void	stressStrainIntegrate (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &dsig, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv) override
	Stress-strain integration.
void	returnMapping (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &dsig, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv)
	return mapping algorithm to solve plastic problem (plastic strain and stress)
T	elasticTrial (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Perform the elastic trial.
void	revertElasticTrial (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Revert the elastic trial.
virtual void	beginConstantStressStateScope (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Begin the scope of a constant stress state which can be used to cache intermediate variables for the plastic correction, such as derivatives of the yield and plastic potential functions, and hardening law.
void	beginConstantStressStateScope1 (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Begin the scope of a constant stress state which can be used to cache intermediate variables for the plastic correction, such as derivatives of the yield and plastic potential functions, and hardening law.
virtual void	endConstantStressStateScope (const StateVariable< T > &sv0, StateVariable< T > &sv)
	End the scope of a constant stress state.
void	endConstantStressStateScope1 (const StateVariable< T > &sv0, StateVariable< T > &sv)
	End the scope of a constant stress state.
virtual bool	converged (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Whether the current stress state can be accepted as converged, used for elastoplastic substepping.
void	ExplicitCPA (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Explicit/CPA methods for elastoplastic model.
void	CPPM (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv)
	CPPM method for elastoplastic model.
std::tuple< bool, T >	integrateToYieldSurface (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Integrate stress state to the yield surface elastically, solve f(sig0 + alpha * K * deps) = 0.
std::tuple< Eigen::Matrix< T, 6, 6 >, Eigen::Matrix< T, 6, 1 >, Eigen::Matrix< T, 6, 1 >, Eigen::Matrix< T, Eigen::Dynamic, 1 >, Eigen::Matrix< T, Eigen::Dynamic, 6 > >	RKn (int n, const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv, const bool revert=false)
	Perform a Runge-Kutta increment.
void	RungeKutta (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv, std::optional< int > order=std::nullopt)
	Runge-Kutta method for elastoplastic model.
virtual T	CSL (const StateVariable< T > &sv0, StateVariable< T > &sv) const
	calculate the critical state line
T	yieldFunction1 (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Yield function for elastoplastic model.
Eigen::Matrix< T, 6, 6 >	dsigadsig (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the modified stress with respect to the stress.
virtual T	potentialFunction (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Potential function for elastoplastic model.
Eigen::Matrix< T, 6, 1 >	dgdsig0 (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the potential function with respect to the stress.
virtual Eigen::Matrix< T, 6, 1 >	dgdsig (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Derivative of the potential function with respect to the stress.
Eigen::Matrix< T, 6, 1 >	dgdsig1 (const Eigen::Matrix< T, 6, 1 > &dfdsig, const StateVariable< T > &sv0, StateVariable< T > &sv, std::optional< bool > associated=std::nullopt)
	Derivative of the potential function with respect to the stress.
Eigen::Matrix< T, 6, 6 >	d2gdsig20 (const Eigen::Matrix< T, 6, 1 > &dfdsig, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Second-order derivative of the potential function with respect to the stress.
virtual Eigen::Matrix< T, 6, 6 >	d2gdsig2 (const Eigen::Matrix< T, 6, 1 > &dfdsig, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Second-order derivative of the potential function with respect to the stress.
Eigen::Matrix< T, 6, 6 >	d2gdsig21 (const Eigen::Matrix< T, 6, 1 > &dfdsig, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Second-order derivative of the potential function with respect to the stress.
virtual T	plasticModulus (const Eigen::Matrix< T, 6, 1 > &dgdsig, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Plastic modulus for elastoplastic model, including dgdsig.
T	plasticModulus1 (const Eigen::Matrix< T, 6, 1 > &dgdsig, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Plastic modulus for elastoplastic model, including dgdsig.
virtual void	updateHardeningVariables (const T ddlambda, const Eigen::Matrix< T, 6, 1 > &ddepsp, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Update the hardening variables.
virtual void	yieldFunctionCorrection (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &deps, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Yield surface correction to keep the force state on the yield surface.
Eigen::Matrix< T, 3, 3 >	orthogonalDirection (const Eigen::Matrix< T, 3, 3 > &F, const Eigen::Matrix< T, 3, 3 > &n) const
	Orthogonal direction.
Eigen::Matrix< T, 3, 3 >	nonCoaxialDirection (const Eigen::Matrix< T, 3, 3 > &F, const Eigen::Matrix< T, 3, 3 > &n, const Eigen::Matrix< T, 3, 3 > &m) const
	Non-coaxial direction.
void	interpolateFabricTensor (const StateVariable< T > &sv0, StateVariable< T > &sv, const bool clear_caches=true, std::optional< std::string > combined_stress_tensor_type=std::nullopt)
	Interpolate fabric tensor towards critical state.
Protected Member Functions inherited from constitutive::ConstitutiveModelBase< T >
bool	checkStiffnessMatrix (const Eigen::Matrix< T, 6, 6 > &K, std::optional< bool > check_positive_definite=std::nullopt, std::optional< bool > check_non_singularity=std::nullopt)
	Check whether the stiffness matrix is valid.
T	updateVoidRatio (T e0, StateVariable< T > &sv)
	Update the void ratio after the isotropic loading to the void ratio before the isotropic loading.
virtual void	postInitialize (StateVariable< T > &sv)
	Post-initialize the state variables.
void	saveElasticStiffnessMatrix (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Save elastic stiffness matrix.
virtual void	postStressStrainIntegrate (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Post-processing after the stress-strain integration.
virtual void	preIncrement (const StateVariable< T > &sv0, StateVariable< T > &sv)
	Pre-processing before the increment.
std::tuple< Eigen::Matrix< T, 6, 1 >, Eigen::Matrix< T, 6, 1 > >	linearize (const Eigen::Matrix< T, 6, 6 > &K, const Eigen::Matrix< T, 6, 1 > &dx, const Eigen::Matrix< T, 6, 6 > &S, const Eigen::Matrix< T, 6, 6 > &E, const Eigen::Matrix< T, 6, 6 > &G, const StateVariable< T > &sv0, StateVariable< T > &sv)
	Linearize the stress and strain increment.
void	resetCustomStateVariables (StateVariable< T > &sv0, StateVariable< T > &sv)
	Reset the custom state variables to zeros existing in the current state variables but not in the initial state variables.
Eigen::Matrix< T, 3, 3 >	loadingDirection (const StateVariable< T > &sv0, StateVariable< T > &sv) const
	loading direction

Additional Inherited Members
Protected Attributes inherited from constitutive::ElastoplasticModelBase< T >
std::string	state = "none"
bool	cached_ = false
	whether the state variables have been cached
Eigen::Matrix< T, 1, 6 >	dpadsiga = {-1.0 / 3.0, -1.0 / 3.0, -1.0 / 3.0, 0.0, 0.0, 0.0}
Eigen::Matrix< T, 1, 6 >	dqadtaua = Eigen::Matrix<T, 1, 6>::Zero()
Eigen::Matrix< T, 6, 6 >	dtauadsiga
Protected Attributes inherited from constitutive::ConstitutiveModelBase< T >
T	_e0 = 0.0
	original void ratio

Detailed Description

template<typename T>
struct micromechanical::core::MicromechanicalBase< T >

Base class for micromechanical models.

Constructor & Destructor Documentation

◆ MicromechanicalBase() [1/2]

template<typename T>

micromechanical::core::MicromechanicalBase< T >::MicromechanicalBase	(	const std::variant< ContactLawType, std::string, ContactLawBase< T > * > &	contact,
		const std::unordered_map< std::string, T > &	props = {},
		T	npv = 1e9,
		T	radius = 0.65e-3 )

inlineexplicit

constructor

Parameters

contact	contact law
props	material properties
npv	number of integration points
radius	radius of the particle

◆ MicromechanicalBase() [2/2]

template<typename T>

micromechanical::core::MicromechanicalBase< T >::MicromechanicalBase ( )

default

default constructor

Member Function Documentation

◆ averageStrain()

template<typename T>

virtual Eigen::Matrix< T, 6, 1 > micromechanical::core::MicromechanicalBase< T >::averageStrain	(	const Eigen::Matrix< T, Eigen::Dynamic, 3 >	disp,
		const StateVariable< T > &	sv0,
		StateVariable< T > &	sv )

inlinevirtual

Average the micro contact displacement to macro strain.

Parameters

disp	micro contact displacement
sv0	initial state variable
sv	state variable

Returns: macro strain

Reimplemented in micromechanical::core::BestFitAveraging< T >.

◆ averageStress()

template<typename T>

virtual Eigen::Matrix< T, 6, 1 > micromechanical::core::MicromechanicalBase< T >::averageStress	(	const Eigen::Matrix< T, Eigen::Dynamic, 3 >	force,
		const StateVariable< T > &	sv0,
		StateVariable< T > &	sv )

inlinevirtual

Average the micro contact force to macro stress.

Parameters

force	micro contact force
sv0	initial state variable
sv	state variable

Returns: macro stress

Reimplemented in micromechanical::core::LoveWeberAveraging< T >.

◆ contactIntegrate()

template<typename T>

virtual Eigen::Matrix< T, 3, 3 > micromechanical::core::MicromechanicalBase< T >::contactIntegrate	(	const int	idx,
		const Eigen::Matrix< T, 3, 1 > &	dforce,
		const Eigen::Matrix< T, 3, 1 > &	ddisp,
		const StateVariable< T > &	sv0,
		StateVariable< T > &	sv ) const

protectedpure virtual

Integrate the contact law.

Parameters

idx	index of the integration point
dforce	force increment
ddisp	displacement increment
sv0	initial state variable
sv	state variable

Returns: the elastoplastic stiffness matrix

Implemented in micromechanical::ChangHicher< T >.

◆ contactIntegrates()

template<typename T>

void micromechanical::core::MicromechanicalBase< T >::contactIntegrates	(	const Eigen::Matrix< T, Eigen::Dynamic, 3 > &	dforce,
		const Eigen::Matrix< T, Eigen::Dynamic, 3 > &	ddisp,
		const StateVariable< T > &	sv0,
		StateVariable< T > &	sv )

inlineprotected

Integrate all the contact laws with cache.

Parameters

dforce	force increment
ddisp	displacement increment
sv0	initial state variable
sv	state variable

Referenced by micromechanical::ChangHicher< T >::stressStrainIntegrate1().

Here is the caller graph for this function:

◆ createStateVariable()

template<typename T>

StateVariable< T > micromechanical::core::MicromechanicalBase< T >::createStateVariable ( )

inlineoverrideprotectedvirtual

Create a new state variable object.

Implements constitutive::ConstitutiveModelBase< T >.

◆ fromAbaqus()

template<typename T>

void micromechanical::core::MicromechanicalBase< T >::fromAbaqus	(	StateVariable< T > &	sv0,
		T *	statev )

inlineoverrideprotectedvirtual

Sync the state variables sv0 from abaqus statev.

Parameters

sv0	state variables
statev	abaqus state variables

Implements constitutive::ConstitutiveModelBase< T >.

◆ initialize()

template<typename T>

void micromechanical::core::MicromechanicalBase< T >::initialize ( StateVariable< T > & sv )

inlineoverrideprotectedvirtual

Initialize the state variables.

Parameters

sv	A StateVariable object to be initialized

Returns: A StateVariable object with initialized state variables

Implements constitutive::ConstitutiveModelBase< T >.

◆ localizeStrain()

template<typename T>

virtual Eigen::Matrix< T, Eigen::Dynamic, 3 > micromechanical::core::MicromechanicalBase< T >::localizeStrain	(	const Eigen::Matrix< T, 6, 1 > &	eps,
		const StateVariable< T > &	sv0,
		StateVariable< T > &	sv )

inlinevirtual

localize the macro strain to micro contact displacement

Parameters

eps	macro strain
sv0	initial state variable
sv	state variable

Returns: micro contact displacement

Reimplemented in micromechanical::core::KinematicHypothesisLocalization< T >.

◆ localizeStrain2()

template<typename T>

virtual Eigen::Matrix< T, Eigen::Dynamic, 3 > micromechanical::core::MicromechanicalBase< T >::localizeStrain2	(	const Eigen::Matrix< T, 6, 1 > &	eps,
		const StateVariable< T > &	sv0,
		StateVariable< T > &	sv,
		const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &	kn = std::nullopt,
		const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &	ks = std::nullopt,
		const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &	kt = std::nullopt )

inlinevirtual

localize the macro strain to micro contact displacement with the static hypothesis, estimated from the stress and contact force.

Parameters

eps	macro strain
sv0	initial state variable
sv	state variable
kn	contact normal stiffness
ks	contact shear stiffness at the first direction
kt	contact shear stiffness at the second direction

Returns: micro contact displacement

Reimplemented in micromechanical::core::StaticHypothesisLocalization< T >.

◆ localizeStress()

template<typename T>

virtual Eigen::Matrix< T, Eigen::Dynamic, 3 > micromechanical::core::MicromechanicalBase< T >::localizeStress	(	const Eigen::Matrix< T, 6, 1 > &	sig,
		const StateVariable< T > &	sv0,
		StateVariable< T > &	sv )

inlinevirtual

localize the macro stress to micro contact force

Parameters

sig	macro stress
sv0	initial state variable
sv	state variable

Returns: micro contact force

Reimplemented in micromechanical::core::StaticHypothesisLocalization< T >.

◆ localizeStress2()

template<typename T>

virtual Eigen::Matrix< T, Eigen::Dynamic, 3 > micromechanical::core::MicromechanicalBase< T >::localizeStress2	(	const Eigen::Matrix< T, 6, 1 > &	sig,
		const StateVariable< T > &	sv0,
		StateVariable< T > &	sv,
		const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &	kn = std::nullopt,
		const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &	ks = std::nullopt,
		const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &	kt = std::nullopt )

inlinevirtual

localize the macro stress to micro contact force with the kinematic hypothesis, estimated from the strain and contact displacement.

Parameters

sig	macro stress
sv0	initial state variable
sv	state variable
kn	contact normal stiffness
ks	contact shear stiffness at the first direction
kt	contact shear stiffness at the second direction

Returns: micro contact force

Reimplemented in micromechanical::core::KinematicHypothesisLocalization< T >.

◆ nstatev()

template<typename T>

int micromechanical::core::MicromechanicalBase< T >::nstatev ( ) const

inlineprotected

◆ postIncrement()

template<typename T>

void micromechanical::core::MicromechanicalBase< T >::postIncrement	(	const StateVariable< T > &	sv0,
		StateVariable< T > &	sv )

inlineoverrideprotectedvirtual

Post-processing after the increment.

Reimplemented from constitutive::ConstitutiveModelBase< T >.

◆ postProcess()

template<typename T>

void micromechanical::core::MicromechanicalBase< T >::postProcess ( StateVariable< T > & sv )

inlineoverrideprotectedvirtual

Post-process the state variables.

Parameters

sv	state variables

Reimplemented from constitutive::ConstitutiveModelBase< T >.

◆ preStressStrainIntegrate()

template<typename T>

void micromechanical::core::MicromechanicalBase< T >::preStressStrainIntegrate	(	const StateVariable< T > &	sv0,
		StateVariable< T > &	sv )

inlineoverrideprotectedvirtual

Pre-processing before the stress-strain integration.

Reimplemented from constitutive::ConstitutiveModelBase< T >.

◆ setContactLaw()

template<typename T>

void micromechanical::core::MicromechanicalBase< T >::setContactLaw	(	const std::variant< ContactLawType, std::string, ContactLawBase< T > * > &	contact,
		const std::unordered_map< std::string, T > &	props = {} )

inline

Set the contact law if it is not set in the constructor.

Parameters

contact	contact law type
props	material properties

◆ stiffness1()

template<typename T>

virtual Eigen::Matrix< T, 6, 6 > micromechanical::core::MicromechanicalBase< T >::stiffness1	(	const StateVariable< T > &	sv0,
		StateVariable< T > &	sv,
		const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &	kn = std::nullopt,
		const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &	ks = std::nullopt,
		const std::optional< Eigen::Matrix< T, Eigen::Dynamic, 1 > > &	kt = std::nullopt )

pure virtual

Compute the stiffness matrix.

Parameters

sv0	initial state variable
sv	state variable
kn	contact normal stiffness
ks	contact shear stiffness at the first direction
kt	contact shear stiffness at the second direction

Returns: stiffness matrix

Implemented in micromechanical::ChangHicher< T >.

Referenced by micromechanical::core::StaticHypothesisLocalization< T >::localizeStrain2(), and micromechanical::core::KinematicHypothesisLocalization< T >::localizeStress2().

Here is the caller graph for this function:

◆ toAbaqus()

template<typename T>

void micromechanical::core::MicromechanicalBase< T >::toAbaqus	(	const StateVariable< T > &	sv0,
		StateVariable< T > &	sv,
		T *	statev )

inlineoverrideprotectedvirtual

Sync the state variables sv to abaqus statev.

Parameters

sv0	initial state variables
sv	state variables
statev	abaqus state variables

Implements constitutive::ConstitutiveModelBase< T >.

Member Data Documentation

◆ _integration

template<typename T>

int micromechanical::core::MicromechanicalBase< T >::_integration = 0

number of integration steps

Referenced by micromechanical::ChangHicher< T >::stressStrainIntegrate1().

◆ contact

template<typename T>

std::shared_ptr<ContactLawBase<T> > micromechanical::core::MicromechanicalBase< T >::contact

contact law

◆ integration

template<typename T>

std::shared_ptr<Integration<T> > micromechanical::core::MicromechanicalBase< T >::integration

gauss integration points

Referenced by micromechanical::core::BestFitAveraging< T >::averageStrain(), micromechanical::core::LoveWeberAveraging< T >::averageStress(), micromechanical::core::KinematicHypothesisLocalization< T >::localizeStrain(), micromechanical::core::StaticHypothesisLocalization< T >::localizeStrain2(), micromechanical::core::StaticHypothesisLocalization< T >::localizeStress(), micromechanical::core::KinematicHypothesisLocalization< T >::localizeStress2(), micromechanical::ChangHicher< T >::stiffness1(), and micromechanical::ChangHicher< T >::stressStrainIntegrate1().

The documentation for this struct was generated from the following file:

include/micromechanical/core/micromechanical.hpp

Public Member Functions

Public Attributes

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MicromechanicalBase() [1/2]

◆ MicromechanicalBase() [2/2]

Member Function Documentation

◆ averageStrain()

◆ averageStress()

◆ contactIntegrate()

◆ contactIntegrates()

◆ createStateVariable()

◆ fromAbaqus()

◆ initialize()

◆ localizeStrain()

◆ localizeStrain2()

◆ localizeStress()

◆ localizeStress2()

◆ nstatev()

◆ postIncrement()

◆ postProcess()

◆ preStressStrainIntegrate()

◆ setContactLaw()

◆ stiffness1()

◆ toAbaqus()

Member Data Documentation

◆ _integration

◆ contact

◆ integration