Inheritance diagram for EshelbianPlasticity::HMHStVenantKirchhoff:

Collaboration diagram for EshelbianPlasticity::HMHStVenantKirchhoff:

Public Member Functions
	HMHStVenantKirchhoff (const double lambda, const double mu, const double sigma_y)

MoFEMErrorCode	getOptions ()

virtual OpJacobian *	returnOpJacobian (const std::string &field_name, const int tag, const bool eval_rhs, const bool eval_lhs, boost::shared_ptr< DataAtIntegrationPts > &data_ptr, boost::shared_ptr< PhysicalEquations > &physics_ptr)

MoFEMErrorCode	recordTape (const int tape, DTensor2Ptr *t_h_ptr)

Public Member Functions inherited from EshelbianPlasticity::PhysicalEquations
	PhysicalEquations ()=delete

	PhysicalEquations (const int size_active, const int size_dependent)

virtual	~PhysicalEquations ()

virtual MoFEMErrorCode	recordTape (const int tag, DTensor2Ptr *t_h)=0

virtual OpJacobian *	returnOpJacobian (const std::string &field_name, const int tag, const bool eval_rhs, const bool eval_lhs, boost::shared_ptr< DataAtIntegrationPts > &data_ptr, boost::shared_ptr< PhysicalEquations > &physics_ptr)=0

DTensor2Ptr	get_P ()

DTensor2Ptr	get_Sigma ()

DTensor0Ptr	get_Phi ()

double &	get_PhiRef ()

DTensor2Ptr	get_Flow ()

DTensor3Ptr	get_P_dh0 ()

DTensor3Ptr	get_P_dH0 ()

DTensor3Ptr	get_P_dh1 ()

DTensor3Ptr	get_P_dH1 ()

DTensor3Ptr	get_P_dh2 ()

DTensor3Ptr	get_P_dH2 ()

DTensor3Ptr	get_Sigma_dh0 ()

DTensor3Ptr	get_Sigma_dH0 ()

DTensor3Ptr	get_Sigma_dh1 ()

DTensor3Ptr	get_Sigma_dH1 ()

DTensor3Ptr	get_Sigma_dh2 ()

DTensor3Ptr	get_Sigma_dH2 ()

DTensor2Ptr	get_Phi_dh ()

DTensor2Ptr	get_Phi_dH ()

DTensor3Ptr	get_Flow_dh0 ()

DTensor3Ptr	get_Flow_dH0 ()

DTensor3Ptr	get_Flow_dh1 ()

DTensor3Ptr	get_Flow_dH1 ()

DTensor3Ptr	get_Flow_dh2 ()

DTensor3Ptr	get_Flow_dH2 ()

DTensor2Ptr	get_h ()

DTensor2Ptr	get_H ()

Public Attributes
double	lambda

double	mu

double	sigmaY

ATensor2	th

ATensor2	tH

ATensor2	tP

ATensor2	tSigma

adouble	phi

adouble	s2

adouble	f

adouble	energy

adouble	detH

adouble	detF

adouble	trE

adouble	meanCauchy

ATensor2	tInvH

ATensor2	tF

ATensor2	tInvF

ATensor2	tC

ATensor2	tE

ATensor2	tS

ATensor2	tCauchy

ATensor2	tDevCauchy

ATensor2	tPhi_dDevCauchy

ATensor2	tPhi_dCauchy

ATensor2	tPhi_dSigma

ATensor2	tPulledP

ATensor2	tPulledSigma

ATensor2	tPulledPhi_dSigma

Public Attributes inherited from EshelbianPlasticity::PhysicalEquations
std::vector< double >	activeVariables

std::vector< double >	dependentVariables

std::vector< double >	dependentVariablesDirevatives

Static Public Attributes
static constexpr int	numberOfActiveVariables = 9 + 9

static constexpr int	numberOfDependentVariables = 9 + 9 + 1 + 9

Additional Inherited Members
Public Types inherited from EshelbianPlasticity::PhysicalEquations
typedef FTensor::Tensor1< adouble, 3 >	ATensor1

typedef FTensor::Tensor2< adouble, 3, 3 >	ATensor2

typedef FTensor::Tensor3< adouble, 3, 3, 3 >	ATensor3

typedef FTensor::Tensor0< FTensor::PackPtr< double *, 1 > >	DTensor0Ptr

typedef FTensor::Tensor2< FTensor::PackPtr< double *, 1 >, 3, 3 >	DTensor2Ptr

typedef FTensor::Tensor3< FTensor::PackPtr< double *, 1 >, 3, 3, 3 >	DTensor3Ptr

Static Public Member Functions inherited from EshelbianPlasticity::PhysicalEquations
template<int S>
static DTensor2Ptr	get_VecTensor2 (std::vector< double > &v)

template<int S>
static DTensor0Ptr	get_VecTensor0 (std::vector< double > &v)

template<int S0, int S1, int S2>
static DTensor3Ptr	get_vecTensor3 (std::vector< double > &v)

Detailed Description

Definition at line 231 of file EshelbianADOL-C.cpp.

Constructor & Destructor Documentation

◆ HMHStVenantKirchhoff()

EshelbianPlasticity::HMHStVenantKirchhoff::HMHStVenantKirchhoff	(	const double	lambda,
		const double	mu,
		const double	sigma_y
	)

inline

Definition at line 236 of file EshelbianADOL-C.cpp.

238 : PhysicalEquations(numberOfActiveVariables, numberOfDependentVariables),

239 lambda(lambda), mu(mu), sigmaY(sigma_y) {}

EshelbianPlasticity::HMHStVenantKirchhoff::mu

double mu

Definition: EshelbianADOL-C.cpp:275

EshelbianPlasticity::HMHStVenantKirchhoff::numberOfDependentVariables

static constexpr int numberOfDependentVariables

Definition: EshelbianADOL-C.cpp:234

EshelbianPlasticity::HMHStVenantKirchhoff::sigmaY

double sigmaY

Definition: EshelbianADOL-C.cpp:276

EshelbianPlasticity::HMHStVenantKirchhoff::lambda

double lambda

Definition: EshelbianADOL-C.cpp:274

EshelbianPlasticity::HMHStVenantKirchhoff::numberOfActiveVariables

static constexpr int numberOfActiveVariables

Definition: EshelbianADOL-C.cpp:233

EshelbianPlasticity::PhysicalEquations::PhysicalEquations

PhysicalEquations()=delete

Member Function Documentation

◆ getOptions()

MoFEMErrorCode EshelbianPlasticity::HMHStVenantKirchhoff::getOptions ( )

inline

Definition at line 241 of file EshelbianADOL-C.cpp.

                              {
    MoFEMFunctionBegin;
 
    double E = 1;
    double nu = 0;
 
    CHKERR PetscOptionsBegin(PETSC_COMM_WORLD, "stvenant_", "", "none");
 
    CHKERR PetscOptionsScalar("-young_modulus", "Young modulus", "", E, &E,
                              PETSC_NULL);
    CHKERR PetscOptionsScalar("-poisson_ratio", "poisson ratio", "", nu, &nu,
                              PETSC_NULL);
    CHKERR PetscOptionsScalar("-sigmaY", "plastic strain", "", sigmaY, &sigmaY,
                              PETSC_NULL);
 
    ierr = PetscOptionsEnd();
    CHKERRG(ierr);
 
    lambda = LAMBDA(E, nu);
    mu = MU(E, nu);
 
    MoFEMFunctionReturn(0);
  }

◆ recordTape()

MoFEMErrorCode EshelbianPlasticity::HMHStVenantKirchhoff::recordTape	(	const int	tape,
		DTensor2Ptr *	t_h_ptr
	)

inlinevirtual

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 308 of file EshelbianADOL-C.cpp.

                                                                  {
    MoFEMFunctionBegin;
    FTensor::Index<'i', 3> i;
    FTensor::Index<'j', 3> j;
    FTensor::Index<'I', 3> I;
    FTensor::Index<'J', 3> J;
    FTensor::Number<0> N0;
    FTensor::Number<1> N1;
    FTensor::Number<2> N2;
 
    CHKERR getOptions();
 
    auto ih = get_h();
    auto iH = get_H();
    if (t_h_ptr)
      ih(i, j) = (*t_h_ptr)(i, j);
    else {
      ih(i, j) = 0;
      for (auto ii : {0, 1, 2})
        ih(ii, ii) = 1;
    }
 
    iH(i, j) = 0;
    for (auto ii : {0, 1, 2})
      iH(ii, ii) = 1;
 
    auto r_P = get_P();
    auto r_Sigma = get_Sigma();
    double &r_phi = get_PhiRef();
    auto r_Flow = get_Flow();
 
    enableMinMaxUsingAbs();
    trace_on(tape);
 
    // Set active variables to ADOL-C
    th(i, j) <<= get_h()(i, j);
    tH(i, j) <<= get_H()(i, j);
 
    // Deformation gradient
    CHKERR determinantTensor3by3(tH, detH);
    CHKERR invertTensor3by3(tH, detH, tInvH);
    tF(i, I) = th(i, J) * tInvH(J, I);
    CHKERR determinantTensor3by3(tF, detF);
    CHKERR invertTensor3by3(tF, detF, tInvF);
 
    // Deformation and strain
    tC(I, J) = tF(i, I) * tF(i, J);
    tE(I, J) = tC(I, J);
    tE(N0, N0) -= 1;
    tE(N1, N1) -= 1;
    tE(N2, N2) -= 1;
    tE(I, J) *= 0.5;
 
    // Energy
    trE = tE(I, I);
    energy = 0.5 * lambda * trE * trE + mu * (tE(I, J) * tE(I, J));
 
    // Stress Piola II
    tS(I, J) = tE(I, J);
    tS(I, J) *= 2 * mu;
    tS(N0, N0) += lambda * trE;
    tS(N1, N1) += lambda * trE;
    tS(N2, N2) += lambda * trE;
    // Stress Piola I
    tP(i, J) = tF(i, I) * tS(I, J);
    // Stress Cauchy
    tCauchy(i, j) = tP(i, J) * tF(j, J);
    tCauchy(i, j) *= (1 / detF);
 
    // Stress Eshelby
    tSigma(I, J) = -tF(i, I) * tP(i, J);
    tSigma(N0, N0) += energy;
    tSigma(N1, N1) += energy;
    tSigma(N2, N2) += energy;
 
    // Deviator Cauchy Stress
    meanCauchy = third * tCauchy(i, i);
    tDevCauchy(i, j) = tCauchy(i, j);
    tDevCauchy(N0, N0) -= meanCauchy;
    tDevCauchy(N1, N1) -= meanCauchy;
    tDevCauchy(N2, N2) -= meanCauchy;
 
    // Plastic surface
 
    s2 = tDevCauchy(i, j) * tDevCauchy(i, j); // s:s
    f = sqrt(1.5 * s2);
    phi = f - sigmaY;
 
    // Flow
    tPhi_dDevCauchy(i, j) = tDevCauchy(i, j);
    tPhi_dDevCauchy(i, j) *= 1.5 / f;
    tPhi_dCauchy(i, j) = tPhi_dDevCauchy(i, j);
    tPhi_dCauchy(N0, N0) -= third * tPhi_dDevCauchy(i, i);
    tPhi_dCauchy(N1, N1) -= third * tPhi_dDevCauchy(i, i);
    tPhi_dCauchy(N2, N2) -= third * tPhi_dDevCauchy(i, i);
    tPhi_dSigma(I, J) = tInvF(I, i) * (tPhi_dCauchy(i, j) * tF(j, J));
    tPhi_dSigma(I, J) *= -(1 / detF);
 
    // Pull back
    tPulledP(i, J) = tP(i, I) * tInvH(J, I);
    tPulledP(i, J) *= detH;
    tPulledSigma(i, J) = tSigma(i, I) * tInvH(J, I);
    tPulledSigma(i, J) *= detH;
    tPulledPhi_dSigma(i, J) = tPhi_dSigma(i, I) * tInvH(J, I);
    tPulledPhi_dSigma(i, J) *= detH;
 
    // Set dependent variables to ADOL-C
    tPulledP(i, j) >>= r_P(i, j);
    tPulledSigma(i, j) >>= r_Sigma(i, j);
    phi >>= r_phi;
    tPulledPhi_dSigma(i, j) >>= r_Flow(i, j);
 
    trace_off();
 
    MoFEMFunctionReturn(0);
  }

◆ returnOpJacobian()

virtual OpJacobian * EshelbianPlasticity::HMHStVenantKirchhoff::returnOpJacobian	(	const std::string &	field_name,
		const int	tag,
		const bool	eval_rhs,
		const bool	eval_lhs,
		boost::shared_ptr< DataAtIntegrationPts > &	data_ptr,
		boost::shared_ptr< PhysicalEquations > &	physics_ptr
	)

inlinevirtual

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 266 of file EshelbianADOL-C.cpp.

                                                                    {
    return (
        new OpHMHH(field_name, tag, eval_rhs, eval_lhs, data_ptr, physics_ptr));
  }

Member Data Documentation

◆ detF

adouble EshelbianPlasticity::HMHStVenantKirchhoff::detF

Definition at line 288 of file EshelbianADOL-C.cpp.

◆ detH

adouble EshelbianPlasticity::HMHStVenantKirchhoff::detH

Definition at line 287 of file EshelbianADOL-C.cpp.

◆ energy

adouble EshelbianPlasticity::HMHStVenantKirchhoff::energy

Definition at line 285 of file EshelbianADOL-C.cpp.

◆ f

adouble EshelbianPlasticity::HMHStVenantKirchhoff::f

Definition at line 284 of file EshelbianADOL-C.cpp.

◆ lambda

double EshelbianPlasticity::HMHStVenantKirchhoff::lambda

Definition at line 274 of file EshelbianADOL-C.cpp.

◆ meanCauchy

adouble EshelbianPlasticity::HMHStVenantKirchhoff::meanCauchy

Definition at line 290 of file EshelbianADOL-C.cpp.

◆ mu

double EshelbianPlasticity::HMHStVenantKirchhoff::mu

Definition at line 275 of file EshelbianADOL-C.cpp.

◆ numberOfActiveVariables

constexpr int EshelbianPlasticity::HMHStVenantKirchhoff::numberOfActiveVariables = 9 + 9

staticconstexpr

Definition at line 233 of file EshelbianADOL-C.cpp.

◆ numberOfDependentVariables

constexpr int EshelbianPlasticity::HMHStVenantKirchhoff::numberOfDependentVariables = 9 + 9 + 1 + 9

staticconstexpr

Definition at line 234 of file EshelbianADOL-C.cpp.

◆ phi

adouble EshelbianPlasticity::HMHStVenantKirchhoff::phi

Definition at line 282 of file EshelbianADOL-C.cpp.

◆ s2

adouble EshelbianPlasticity::HMHStVenantKirchhoff::s2

Definition at line 283 of file EshelbianADOL-C.cpp.

◆ sigmaY

double EshelbianPlasticity::HMHStVenantKirchhoff::sigmaY

Definition at line 276 of file EshelbianADOL-C.cpp.

◆ tC

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tC

Definition at line 294 of file EshelbianADOL-C.cpp.

◆ tCauchy

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tCauchy

Definition at line 297 of file EshelbianADOL-C.cpp.

◆ tDevCauchy

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tDevCauchy

Definition at line 298 of file EshelbianADOL-C.cpp.

◆ tE

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tE

Definition at line 295 of file EshelbianADOL-C.cpp.

◆ tF

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tF

Definition at line 292 of file EshelbianADOL-C.cpp.

◆ th

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::th

Definition at line 278 of file EshelbianADOL-C.cpp.

◆ tH

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tH

Definition at line 279 of file EshelbianADOL-C.cpp.

◆ tInvF

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tInvF

Definition at line 293 of file EshelbianADOL-C.cpp.

◆ tInvH

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tInvH

Definition at line 291 of file EshelbianADOL-C.cpp.

◆ tP

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tP

Definition at line 280 of file EshelbianADOL-C.cpp.

◆ tPhi_dCauchy

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tPhi_dCauchy

Definition at line 301 of file EshelbianADOL-C.cpp.

◆ tPhi_dDevCauchy

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tPhi_dDevCauchy

Definition at line 300 of file EshelbianADOL-C.cpp.

◆ tPhi_dSigma

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tPhi_dSigma

Definition at line 302 of file EshelbianADOL-C.cpp.

◆ tPulledP

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tPulledP

Definition at line 304 of file EshelbianADOL-C.cpp.

◆ tPulledPhi_dSigma

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tPulledPhi_dSigma

Definition at line 306 of file EshelbianADOL-C.cpp.

◆ tPulledSigma

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tPulledSigma

Definition at line 305 of file EshelbianADOL-C.cpp.

◆ trE

adouble EshelbianPlasticity::HMHStVenantKirchhoff::trE

Definition at line 289 of file EshelbianADOL-C.cpp.

◆ tS

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tS

Definition at line 296 of file EshelbianADOL-C.cpp.

◆ tSigma

ATensor2 EshelbianPlasticity::HMHStVenantKirchhoff::tSigma

Definition at line 281 of file EshelbianADOL-C.cpp.

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

users_modules/eshelbian_plasticty/src/impl/EshelbianADOL-C.cpp

Public Member Functions

Public Attributes

Static Public Attributes

Additional Inherited Members