#include <shelleasans.h>

Inheritance diagram for Shell4EASANS:

Collaboration diagram for Shell4EASANS:

Public Types
enum	IntegrationPoint { IP_1_1 = 0, IP_1_2 = 1, IP_1_3 = 2, IP_2_1 = 3, IP_2_2 = 4, IP_2_3 = 5, IP_3_1 = 6, IP_3_2 = 7, IP_3_3 = 8, NUMIP = 4 }

enum	ShearStrainEvaluationPoint { SSEP_1 = 0, SSEP_2 = 1, SSEP_3 = 2, SSEP_4 = 3, NUMSSEP = 4 }

enum	NodeName { NODE1 = 0, NODE2 = 1, NODE3 = 2, NODE4 = 3, NUMNODES = 4 }

enum	Deformations { STRAIN = 0, CURVAT = 1, NUMDEFORM = 2 }

Public Types inherited from Elem
enum	Type { UNKNOWN = -1, AIRPROPERTIES = 0, INDUCEDVELOCITY, AUTOMATICSTRUCTURAL, GRAVITY, BODY, JOINT, JOINT_REGULARIZATION, BEAM, PLATE, FORCE, INERTIA, ELECTRICBULK, ELECTRIC, THERMAL, HYDRAULIC, BULK, LOADABLE, DRIVEN, EXTERNAL, AEROMODAL, AERODYNAMIC, GENEL, SOCKETSTREAM_OUTPUT, RTAI_OUTPUT = SOCKETSTREAM_OUTPUT, LASTELEMTYPE }

Public Types inherited from SimulationEntity
typedef std::vector< Hint * >	Hints

Public Types inherited from ToBeOutput
enum	{ OUTPUT = 0x1U, OUTPUT_MASK = 0xFU, OUTPUT_PRIVATE = 0x10U, OUTPUT_PRIVATE_MASK = ~OUTPUT_MASK }

Public Types inherited from Shell
enum	Type { UNKNOWN = -1, ELASTIC = 0, VISCOELASTIC, LASTSHELLTYPE }

typedef MyVectorHandler	vh

typedef std::vector< vh >	vvh

typedef FullMatrixHandler	fmh

typedef std::vector< fmh >	vfmh

typedef ConstitutiveLawOwner < vh, fmh >	ConstitutiveLawOwnerType

Public Member Functions
	Shell4EASANS (unsigned int uL, const DofOwner pDO, const StructNode pN1, const StructNode pN2, const StructNode pN3, const StructNode *pN4, const Mat3x3 &R1, const Mat3x3 &R2, const Mat3x3 &R3, const Mat3x3 &R4, const fmh &pDTmp, const vh &PreStress, flag fOut)

virtual	~Shell4EASANS (void)

virtual Shell::Type	GetShellType (void) const

virtual Elem::Type	GetElemType (void) const

virtual unsigned int	iGetNumDof (void) const

virtual DofOrder::Order	GetDofType (unsigned int i) const

virtual DofOrder::Order	GetEqType (unsigned int i) const

virtual std::ostream &	Restart (std::ostream &out) const

virtual void	WorkSpaceDim (integer piNumRows, integer piNumCols) const

void	SetValue (DataManager pDM, VectorHandler &, VectorHandler &, SimulationEntity::Hints ph=0)

virtual SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)

virtual VariableSubMatrixHandler &	AssJac (VariableSubMatrixHandler &WorkMat, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)

virtual void	Output (OutputHandler &OH) const

virtual unsigned int	iGetInitialNumDof (void) const

virtual void	InitialWorkSpaceDim (integer piNumRows, integer piNumCols) const

virtual void	SetInitialValue (VectorHandler &)

virtual VariableSubMatrixHandler &	InitialAssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)

virtual SubVectorHandler &	InitialAssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr)

virtual const StructNode *	pGetNode (unsigned int i) const

virtual void	GetConnectedNodes (std::vector< const Node * > &connectedNodes) const

Public Member Functions inherited from Elem
	Elem (unsigned int uL, flag fOut)

virtual	~Elem (void)

virtual std::ostream &	DescribeDof (std::ostream &out, const char *prefix="", bool bInitial=false) const

virtual void	DescribeDof (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const

virtual std::ostream &	DescribeEq (std::ostream &out, const char *prefix="", bool bInitial=false) const

virtual void	DescribeEq (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const

virtual void	AssMats (VariableSubMatrixHandler &WorkMatA, VariableSubMatrixHandler &WorkMatB, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)

virtual bool	bInverseDynamics (void) const

void	SetInverseDynamicsFlags (unsigned uIDF)

unsigned	GetInverseDynamicsFlags (void) const

bool	bIsErgonomy (void) const

bool	bIsRightHandSide (void) const

virtual VariableSubMatrixHandler &	AssJac (VariableSubMatrixHandler &WorkMat, const VectorHandler &XCurr)

virtual SubVectorHandler &	AssRes (SubVectorHandler &WorkVec, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr, const VectorHandler &XPrimePrimeCurr, InverseDynamics::Order iOrder=InverseDynamics::INVERSE_DYNAMICS)

virtual int	GetNumConnectedNodes (void) const

Public Member Functions inherited from WithLabel
	WithLabel (unsigned int uL=0, const std::string &sN="")

virtual	~WithLabel (void)

void	PutLabel (unsigned int uL)

void	PutName (const std::string &sN)

unsigned int	GetLabel (void) const

const std::string &	GetName (void) const

Public Member Functions inherited from SimulationEntity
	SimulationEntity (void)

virtual	~SimulationEntity (void)

virtual bool	bIsValidIndex (unsigned int i) const

virtual Hint *	ParseHint (DataManager pDM, const char s) const

virtual void	BeforePredict (VectorHandler &, VectorHandler &, VectorHandler &, VectorHandler &) const

virtual void	AfterPredict (VectorHandler &X, VectorHandler &XP)

virtual void	Update (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)

virtual void	DerivativesUpdate (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)

virtual void	Update (const VectorHandler &XCurr, InverseDynamics::Order iOrder)

virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP)

virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP, const VectorHandler &XPP)

virtual unsigned int	iGetNumPrivData (void) const

virtual unsigned int	iGetPrivDataIdx (const char *s) const

virtual doublereal	dGetPrivData (unsigned int i) const

virtual std::ostream &	OutputAppend (std::ostream &out) const

virtual void	ReadInitialState (MBDynParser &HP)

Public Member Functions inherited from ToBeOutput
	ToBeOutput (flag fOut=fDefaultOut)

virtual	~ToBeOutput (void)

virtual void	OutputPrepare (OutputHandler &OH)

virtual void	Output (OutputHandler &OH, const VectorHandler &X, const VectorHandler &XP) const

virtual flag	fToBeOutput (void) const

virtual bool	bToBeOutput (void) const

virtual void	SetOutputFlag (flag f=flag(1))

Public Member Functions inherited from Shell
	Shell (unsigned uLabel, const DofOwner *pDO, flag fOut)

virtual	~Shell (void)

Public Member Functions inherited from ElemGravityOwner
	ElemGravityOwner (unsigned int uL, flag fOut)

virtual	~ElemGravityOwner (void)

virtual doublereal	dGetM (void) const

Vec3	GetS (void) const

Mat3x3	GetJ (void) const

Vec3	GetB (void) const

Vec3	GetG (void) const

Public Member Functions inherited from GravityOwner
	GravityOwner (void)

virtual	~GravityOwner (void)

void	PutGravity (const Gravity *pG)

virtual bool	bGetGravity (const Vec3 &X, Vec3 &Acc) const

Public Member Functions inherited from ElemWithDofs
	ElemWithDofs (unsigned int uL, const DofOwner *pDO, flag fOut)

virtual	~ElemWithDofs (void)

Public Member Functions inherited from DofOwnerOwner
	DofOwnerOwner (const DofOwner *pDO)

virtual	~DofOwnerOwner ()

virtual const DofOwner *	pGetDofOwner (void) const

virtual integer	iGetFirstIndex (void) const

Public Member Functions inherited from InitialAssemblyElem
	InitialAssemblyElem (unsigned int uL, flag fOut)

virtual	~InitialAssemblyElem (void)

Public Member Functions inherited from SubjectToInitialAssembly
	SubjectToInitialAssembly (void)

virtual	~SubjectToInitialAssembly (void)

Static Public Attributes
static doublereal	xi_i [NUMIP][2]

static doublereal	w_i [NUMIP]

static doublereal	xi_A [NUMSSEP][2]

static doublereal	xi_n [NUMNODES][2]

static doublereal	xi_0 [2] = {0., 0.}

Protected Member Functions
virtual void	AddInternalForces (Vec6 &, unsigned int)

Protected Member Functions inherited from ElemGravityOwner
virtual Vec3	GetS_int (void) const

virtual Mat3x3	GetJ_int (void) const

virtual Vec3	GetB_int (void) const

virtual Vec3	GetG_int (void) const

Protected Attributes
const StructNode *	pNode [NUMNODES]

Vec3	xa_0 [NUMNODES]

Vec3	xa [NUMNODES]

Mat3x3	iTa [NUMNODES]

Mat3x3	iTa_i [NUMIP]

Mat3x3	iTa_A [NUMSSEP]

Vec3	phi_tilde_n [NUMNODES]

Vec3	phi_tilde_i [NUMIP]

Vec3	phi_tilde_A [NUMSSEP]

Vec3	phi_tilde_0

Mat3x3	T0_overline

Mat3x3	T_overline

Mat3x3	T_0_0

Mat3x3	T_0_i [NUMIP]

Mat3x3	T_0_A [NUMSSEP]

Mat3x3	T_0

Mat3x3	T_i [NUMIP]

Mat3x3	T_A [NUMSSEP]

Mat3x3	Phi_Delta_i [NUMIP][NUMNODES]

Mat3x3	Phi_Delta_A [NUMIP][NUMNODES]

Mat3x3	Kappa_delta_i_1 [NUMIP][NUMNODES]

Mat3x3	Kappa_delta_i_2 [NUMIP][NUMNODES]

Mat3x3	Q_i [NUMIP]

Mat3x3	Q_A [NUMSSEP]

Vec3	k_1_i [NUMIP]

Vec3	k_2_i [NUMIP]

Vec3	eps_tilde_1_0_i [NUMIP]

Vec3	eps_tilde_2_0_i [NUMIP]

Vec3	eps_tilde_1_0_A [NUMSSEP]

Vec3	eps_tilde_2_0_A [NUMSSEP]

Vec3	eps_tilde_1_i [NUMIP]

Vec3	eps_tilde_2_i [NUMIP]

Vec3	eps_tilde_1_A [NUMSSEP]

Vec3	eps_tilde_2_A [NUMSSEP]

Vec3	k_tilde_1_0_i [NUMIP]

Vec3	k_tilde_2_0_i [NUMIP]

Vec3	k_tilde_1_i [NUMIP]

Vec3	k_tilde_2_i [NUMIP]

fmh	S_alpha_beta_0

vfmh	S_alpha_beta_i

vfmh	S_alpha_beta_A

doublereal	alpha_0

doublereal	alpha_i [NUMIP]

vfmh	L_alpha_beta_i

vfmh	L_alpha_beta_A

vfmh	B_overline_i

vfmh	D_overline_i

vfmh	G_i

vfmh	P_i

vfmh	K_beta_beta_i

Vec3	y_i_1 [NUMIP]

Vec3	y_i_2 [NUMIP]

vh	beta

vh	epsilon_hat

vh	epsilon

vvh	FRef

bool	bPreStress

vh	PreStress

vfmh	DRef

vvh	stress_i

bool	bFirstRes

Protected Attributes inherited from WithLabel
unsigned int	uLabel

std::string	sName

Protected Attributes inherited from ToBeOutput
flag	fOutput

Protected Attributes inherited from GravityOwner
Gravity *	pGravity

Private Member Functions
void	UpdateNodalAndAveragePosAndOrientation ()

void	ComputeInitialNodeOrientation ()

void	InterpolateOrientation ()

void	ComputeIPCurvature ()

Detailed Description

Definition at line 70 of file shelleasans.h.

Member Enumeration Documentation

enum Shell4EASANS::Deformations

Enumerator
STRAIN
CURVAT
NUMDEFORM

Definition at line 144 of file shelleasans.h.

                           {
                 STRAIN = 0,
                 CURVAT = 1,
 
                 NUMDEFORM = 2
         };

enum Shell4EASANS::IntegrationPoint

Enumerator
IP_1_1
IP_1_2
IP_1_3
IP_2_1
IP_2_2
IP_2_3
IP_3_1
IP_3_2
IP_3_3
NUMIP

Definition at line 102 of file shelleasans.h.

                               {
                 IP_1_1 = 0,
                 IP_1_2 = 1,
                 IP_1_3 = 2,
                 IP_2_1 = 3,
                 IP_2_2 = 4,
                 IP_2_3 = 5,
                 IP_3_1 = 6,
                 IP_3_2 = 7,
                 IP_3_3 = 8,
 
                 NUMIP = 4
         };

enum Shell4EASANS::NodeName

Enumerator
NODE1
NODE2
NODE3
NODE4
NUMNODES

Definition at line 131 of file shelleasans.h.

                       {
                 NODE1 = 0,
                 NODE2 = 1,
                 NODE3 = 2,
                 NODE4 = 3,
 
                 NUMNODES = 4
         };

enum Shell4EASANS::ShearStrainEvaluationPoint

Enumerator
SSEP_1
SSEP_2
SSEP_3
SSEP_4
NUMSSEP

Definition at line 120 of file shelleasans.h.

                                         {
                 SSEP_1 = 0,
                 SSEP_2 = 1,
                 SSEP_3 = 2,
                 SSEP_4 = 3,
 
                 NUMSSEP = 4
         };

Constructor & Destructor Documentation

Shell4EASANS::Shell4EASANS	(	unsigned int	uL,
		const DofOwner *	pDO,
		const StructNode *	pN1,
		const StructNode *	pN2,
		const StructNode *	pN3,
		const StructNode *	pN4,
		const Mat3x3 &	R1,
		const Mat3x3 &	R2,
		const Mat3x3 &	R3,
		const Mat3x3 &	R4,
		const fmh &	pDTmp,
		const vh &	PreStress,
		flag	fOut
	)

Definition at line 253 of file shelleasans.cc.

References alpha_0, alpha_i, ComputeInitialNodeOrientation(), ComputeIPCurvature(), DRef, eps_tilde_1_0_A, eps_tilde_1_0_i, eps_tilde_2_0_A, eps_tilde_2_0_i, Mat3x3::GetCol(), iGetNumDof(), k_1_i, k_2_i, k_tilde_1_0_i, k_tilde_2_0_i, L_alpha_beta_A, L_alpha_beta_i, MatrixHandler::MatMatMul(), MatrixHandler::MatTMatMul(), Mat3x3::MulTV(), NODE1, NODE2, NODE3, NODE4, NUMIP, NUMNODES, NUMSSEP, P_i, pNode, S_alpha_beta_0, S_alpha_beta_A, S_alpha_beta_i, SAFENEWWITHCONSTRUCTOR, T0_overline, T_0, T_0_0, T_0_A, T_0_i, T_A, T_i, T_overline, xa, xa_0, xi_0, xi_A, xi_i, y_i_1, and y_i_2.

 : 
 Elem(uL, fOut), 
 Shell(uL, pDO, fOut),
 
 S_alpha_beta_0(2, 2),
 S_alpha_beta_i(NUMIP, fmh(2, 2) ),
 S_alpha_beta_A(NUMSSEP, fmh(2, 2) ),
 
 L_alpha_beta_i(NUMIP, fmh(4, 2) ),
 L_alpha_beta_A(NUMSSEP, fmh(4, 2) ),
 
 B_overline_i(NUMIP, fmh(12, 24) ),
 D_overline_i(NUMIP, fmh(15, 24) ),
 G_i(NUMIP, fmh(15, 15) ),
 
 P_i(NUMIP, fmh(12, iGetNumDof()) ),
 
 K_beta_beta_i(NUMIP, fmh(iGetNumDof(), iGetNumDof()) ),
 
 beta(iGetNumDof()),
 epsilon_hat(12),
 epsilon(12),
 
 #ifndef USE_CL_IN_SHELL
 bPreStress(PreStress.Norm() > 0.),
 PreStress(PreStress),
 #endif // ! USE_CL_IN_SHELL
 
 DRef(NUMIP, fmh(12, 12)),
 stress_i(NUMIP, vh(12))
 {
 #ifdef USE_CL_IN_SHELL
         for (integer i = 0; i < NUMIP; i++) {
                 pD[i] = 0;
                 SAFENEWWITHCONSTRUCTOR(pD[i],
                         ConstitutiveLawOwnerType,
                         ConstitutiveLawOwnerType(pDTmp[i]));    
         }
 #else // ! USE_CL_IN_SHELL
         for (unsigned i = 0; i < NUMIP; i++) {
                 DRef[i] = pDTmp;
         }
 #endif // ! USE_CL_IN_SHELL
 
         pNode[NODE1] = pN1;
         pNode[NODE2] = pN2;
         pNode[NODE3] = pN3;
         pNode[NODE4] = pN4;
 //      f[NODE1] = f1;
 //      f[NODE2] = f2;
 //      f[NODE3] = f3;
 //      f[NODE4] = f4;
         ComputeInitialNodeOrientation();
 //      UpdateNodalAndAveragePosAndOrientation();
 //      InterpolateOrientation();
         // copy ref values
         T0_overline = T_overline;
         T_0_0 = T_0;
         for (integer i = 0; i < NUMNODES; i++) {
                 xa_0[i] = xa[i];
         }
         for (integer i = 0; i < NUMIP; i++) {
                 T_0_i[i] = T_i[i];
         }
         for (integer i = 0; i < NUMSSEP; i++) {
                 T_0_A[i] = T_A[i];
         }
 
 
         fmh M_0(4, 4);
         fmh M_0_Inv(4, 4);
         {
                 Vec3 x_1 = InterpDeriv_xi1(xa, xi_0);
                 Vec3 x_2 = InterpDeriv_xi2(xa, xi_0);
                 S_alpha_beta_0(1, 1) = T_0_0.GetCol(1) * x_1;
                 S_alpha_beta_0(2, 1) = T_0_0.GetCol(2) * x_1;
                 S_alpha_beta_0(1, 2) = T_0_0.GetCol(1) * x_2;
                 S_alpha_beta_0(2, 2) = T_0_0.GetCol(2) * x_2;
                 alpha_0 = S_alpha_beta_0(1, 1) * S_alpha_beta_0(2, 2) -
                         S_alpha_beta_0(1, 2) * S_alpha_beta_0(2, 1);
 
                 M_0(1, 1) = S_alpha_beta_0(1, 1) * S_alpha_beta_0(1, 1);
                 M_0(1, 2) = S_alpha_beta_0(1, 2) * S_alpha_beta_0(1, 2);
                 M_0(1, 3) = S_alpha_beta_0(1, 2) * S_alpha_beta_0(1, 1);
                 M_0(1, 4) = S_alpha_beta_0(1, 1) * S_alpha_beta_0(1, 2);        
 
                 M_0(2, 1) = S_alpha_beta_0(2, 1) * S_alpha_beta_0(2, 1);
                 M_0(2, 2) = S_alpha_beta_0(2, 2) * S_alpha_beta_0(2, 2);
                 M_0(2, 3) = S_alpha_beta_0(2, 2) * S_alpha_beta_0(2, 1);
                 M_0(2, 4) = S_alpha_beta_0(2, 1) * S_alpha_beta_0(2, 2);
         
                 M_0(3, 1) = S_alpha_beta_0(1, 1) * S_alpha_beta_0(2, 1);
                 M_0(3, 2) = S_alpha_beta_0(1, 2) * S_alpha_beta_0(2, 2);
                 M_0(3, 3) = S_alpha_beta_0(1, 1) * S_alpha_beta_0(2, 2);
                 M_0(3, 4) = S_alpha_beta_0(1, 2) * S_alpha_beta_0(2, 1);
 
                 M_0(4, 1) = S_alpha_beta_0(1, 1) * S_alpha_beta_0(2, 1);
                 M_0(4, 2) = S_alpha_beta_0(1, 2) * S_alpha_beta_0(2, 2);
                 M_0(4, 3) = S_alpha_beta_0(1, 2) * S_alpha_beta_0(2, 1);
                 M_0(4, 4) = S_alpha_beta_0(1, 1) * S_alpha_beta_0(2, 2);
                 
                 Inv4x4(M_0, M_0_Inv);
         }
 
         for (integer i = 0; i < NUMIP; i++) {
                 fmh L_alpha_B_i(4, 2);
                 Vec3 x_1 = InterpDeriv_xi1(xa, xi_i[i]);
                 Vec3 x_2 = InterpDeriv_xi2(xa, xi_i[i]);
                 S_alpha_beta_i[i](1, 1) = T_0_i[i].GetCol(1) * x_1;
                 S_alpha_beta_i[i](2, 1) = T_0_i[i].GetCol(2) * x_1;
                 S_alpha_beta_i[i](1, 2) = T_0_i[i].GetCol(1) * x_2;
                 S_alpha_beta_i[i](2, 2) = T_0_i[i].GetCol(2) * x_2;
                 // alpha_i = det(S_alpha_beta_i)
                 alpha_i[i] = S_alpha_beta_i[i](1, 1) * S_alpha_beta_i[i](2, 2) -
                         S_alpha_beta_i[i](1, 2) * S_alpha_beta_i[i](2, 1);
 
                 // xi_i_i = S_alpha_beta_i^{-1}
                 FullMatrixHandler xi_i_i(2, 2);
                 Inv2x2(S_alpha_beta_i[i], xi_i_i);
                 
                 for (integer n = 0; n < NUMNODES; n++) {
                         for (integer ii = 0; ii < 2; ii++) {
                                 L_alpha_B_i(n + 1, ii + 1) = LI_J[n][ii](xi_i[i]);
                         }
                 }
                 
                 L_alpha_B_i.MatMatMul(L_alpha_beta_i[i], xi_i_i); 
 
                 fmh H(4, iGetNumDof());
                 doublereal t = xi_i[i][0] * xi_i[i][1];
                 H(1, 1) = xi_i[i][0];
                 H(1, 2) = t;
 
                 H(2, 3) = xi_i[i][1];
                 H(2, 4) = t;
                 
                 H(3, 5) = xi_i[i][0];
                 H(3, 6) = t;
                 
                 H(4, 7) = xi_i[i][1];
                 H(4, 6) = t;
                 
 //              M_0_Inv.MatTMatMul(P_i[i], H);
 //              P_i[i].ScalarMul(alpha_0 / alpha_i[i]);
 // 
 // //////////////
                 fmh Perm(12, 4), tmpP(4, iGetNumDof());
                         // 1, 5, 4, 2, 3, 6
                         Perm(1, 1) = 1.;
                         Perm(2, 3) = 1.;
                         Perm(4, 4) = 1.;
                         Perm(5, 2) = 1.;
                 
                 M_0_Inv.MatTMatMul(tmpP, H);
                 Perm.MatMatMul(P_i[i], tmpP);
                 P_i[i].ScalarMul(alpha_0 / alpha_i[i]);
 
         }
         // save initial axial values
         ComputeIPCurvature();
         for (integer i = 0; i < NUMIP; i++) {
                 InterpDeriv(xa, L_alpha_beta_i[i], y_i_1[i], y_i_2[i]);
                 eps_tilde_1_0_i[i] = T_i[i].MulTV(y_i_1[i]);
                 eps_tilde_2_0_i[i] = T_i[i].MulTV(y_i_2[i]);
                 k_tilde_1_0_i[i] = T_i[i].MulTV(k_1_i[i]);
                 k_tilde_2_0_i[i] = T_i[i].MulTV(k_2_i[i]);
         }
         for (integer i = 0; i < NUMSSEP; i++) {
                 fmh L_alpha_B_A(4, 2);
                 Vec3 x_1 = InterpDeriv_xi1(xa, xi_A[i]);
                 Vec3 x_2 = InterpDeriv_xi2(xa, xi_A[i]);
                 S_alpha_beta_A[i](1, 1) = T_0_A[i].GetCol(1) * x_1;
                 S_alpha_beta_A[i](2, 1) = T_0_A[i].GetCol(2) * x_1;
                 S_alpha_beta_A[i](1, 2) = T_0_A[i].GetCol(1) * x_2;
                 S_alpha_beta_A[i](2, 2) = T_0_A[i].GetCol(2) * x_2;
 
                 Vec3 y_A_1;
                 Vec3 y_A_2;
                 InterpDeriv(xa, L_alpha_beta_A[i], y_A_1, y_A_2);
                 eps_tilde_1_0_A[i] = T_A[i].MulTV(y_A_1);
                 eps_tilde_2_0_A[i] = T_A[i].MulTV(y_A_2);
 
                 // xi_A_i = S_alpha_beta_A^{-1}
                 FullMatrixHandler xi_A_i(2, 2);
                 Inv2x2(S_alpha_beta_A[i], xi_A_i);
 
                 for (integer n = 0; n < NUMNODES; n++) {
                         for (integer ii = 0; ii < 2; ii++) {
                                 L_alpha_B_A(n + 1, ii + 1) = LI_J[n][ii](xi_A[i]);
                         }
                 }
                 
                 L_alpha_B_A.MatMatMul(L_alpha_beta_A[i], xi_A_i); 
 
         }
         {
                 Vec3 y_A_1;
                 Vec3 y_A_2;
                 for (integer i = 0; i < NUMSSEP; i++) {
                         InterpDeriv(xa, L_alpha_beta_A[i], y_A_1, y_A_2);
                         eps_tilde_1_0_A[i] = T_A[i].MulTV(y_A_1);
                         eps_tilde_2_0_A[i] = T_A[i].MulTV(y_A_2);
                 }
         }
 }

Here is the call graph for this function:

Shell4EASANS::~Shell4EASANS ( void )

virtual

Definition at line 476 of file shelleasans.cc.

References ASSERT, NUMIP, and SAFEDELETE.

 {
 #ifdef USE_CL_IN_SHELL
         for (integer i = 0; i < NUMIP; i++) {
                 ASSERT(pD[i] != 0);
                 if (pD[i] != 0) {
                         SAFEDELETE(pD[i]);
                 }
         }
 #endif // USE_CL_IN_SHELL
 }

Member Function Documentation

virtual void Shell4EASANS::AddInternalForces	(	Vec6 &	,
		unsigned	int
	)

inlineprotectedvirtual

Definition at line 294 of file shelleasans.h.

References NO_OP.

                                                      {
                 NO_OP;
         };

VariableSubMatrixHandler & Shell4EASANS::AssJac	(	VariableSubMatrixHandler &	WorkMat,
		doublereal	dCoef,
		const VectorHandler &	XCurr,
		const VectorHandler &	XPrimeCurr
	)

virtual

Implements Elem.

Definition at line 715 of file shelleasans.cc.

References FullSubMatrixHandler::Add(), FullSubMatrixHandler::AddT(), alpha_i, B_overline_i, D_overline_i, DRef, G_i, DofOwnerOwner::iGetFirstIndex(), StructDispNode::iGetFirstMomentumIndex(), StructDispNode::iGetFirstPositionIndex(), iGetNumDof(), MatrixHandler::MatMatMul(), NUMIP, P_i, pNode, FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), VariableSubMatrixHandler::SetFull(), w_i, and WorkSpaceDim().

 {
         FullSubMatrixHandler& WM = WorkMat.SetFull();
 
         /* Dimensiona e resetta la matrice di lavoro */
         integer iNumRows = 0;
         integer iNumCols = 0;
         WorkSpaceDim(&iNumRows, &iNumCols);
         WM.ResizeReset(iNumRows, iNumCols);
 
         /* Recupera e setta gli indici */
         integer iFirstReactionIndex = iGetFirstIndex();
         for (integer i = 0; i < 4; i++) {
                 integer iNodeFirstMomIndex = pNode[i]->iGetFirstMomentumIndex();
                 integer iNodeFirstPosIndex = pNode[i]->iGetFirstPositionIndex();
                 for (int iCnt = 1; iCnt <= 6; iCnt++) {
                         WM.PutRowIndex(iCnt + 6 * i, iNodeFirstMomIndex + iCnt);
                         WM.PutColIndex(iCnt + 6 * i, iNodeFirstPosIndex + iCnt);
                 }
         }
         for (unsigned int iCnt = 1; iCnt <= iGetNumDof(); iCnt++) {
                 WM.PutRowIndex(24 + iCnt, iFirstReactionIndex + iCnt);
                 WM.PutColIndex(24 + iCnt, iFirstReactionIndex + iCnt);
         }
 
         // tangente
 
         FullMatrixHandler Kg(24, 24);
         FullMatrixHandler Km(24, 24);
         FullMatrixHandler K_beta_q(iGetNumDof(), 24);
         FullMatrixHandler K_beta_beta(iGetNumDof(), iGetNumDof());
 
         FullMatrixHandler Ktg(24, 15);
         FullMatrixHandler Ktm(24, 12);
         
         FullMatrixHandler Ktbetaq(iGetNumDof(), 12);
         
         FullMatrixHandler C(12, 12);
         for (integer i = 0; i < NUMIP; i++) {
                 D_overline_i[i].MatTMatMul(Ktg, G_i[i]);
                 Ktg.MatMatMul(Kg, D_overline_i[i]);
 
 #ifdef USE_CL_IN_SHELL
                 C = pD[i]->GetFDE();
 #else // ! USE_CL_IN_SHELL
                 C = DRef[i];
 #endif // ! USE_CL_IN_SHELL
 
                 B_overline_i[i].MatTMatMul(Ktm, C);
 
                 Ktm.MatMatMul(Km, B_overline_i[i]);
                 
                 // extract Cm matrix from C
                 
                 P_i[i].MatTMatMul(Ktbetaq, C);
                 Ktbetaq.MatMatMul(K_beta_q, B_overline_i[i]);
                 
                 Ktbetaq.MatMatMul(K_beta_beta, P_i[i]);
                 
 //              std::cerr << "Kg:\n" << std::fixed << std::setprecision(12) << Kg << std::endl;
 //              std::cerr << "Km:\n" << std::fixed << std::setprecision(12) << Km << std::endl;
 
                 // AssembleMatrix(WM, 1, 1, Kg, alpha_i[i] * w_i[i] * dCoef);
 
                 // AssembleMatrix(WM, 1, 1, Km, alpha_i[i] * w_i[i] * dCoef);
 
                 // AssembleTransposeMatrix(WM, 1, 25, K_beta_q, alpha_i[i] * w_i[i]);
                 // AssembleMatrix(WM, 25, 1, K_beta_q, alpha_i[i] * w_i[i]);
                 // AssembleMatrix(WM, 25, 25, K_beta_beta, alpha_i[i] * w_i[i] / dCoef);
 
                 WM.Add(1, 1, Kg, alpha_i[i] * w_i[i] * dCoef);
 
                 WM.Add(1, 1, Km, alpha_i[i] * w_i[i] * dCoef);
 
                 WM.AddT(1, 25, K_beta_q, alpha_i[i] * w_i[i]);
                 WM.Add(25, 1, K_beta_q, alpha_i[i] * w_i[i]);
                 WM.Add(25, 25, K_beta_beta, alpha_i[i] * w_i[i] / dCoef);
         }
         return WorkMat;
 
 }

Here is the call graph for this function:

SubVectorHandler & Shell4EASANS::AssRes	(	SubVectorHandler &	WorkVec,
		doublereal	dCoef,
		const VectorHandler &	XCurr,
		const VectorHandler &	XPrimeCurr
	)

virtual

Implements Elem.

Definition at line 488 of file shelleasans.cc.

References alpha_i, B_overline_i, beta, bPreStress, ComputeIPCurvature(), FullMatrixHandler::CopyMatrixRow(), D_overline_i, Vec3::Dot(), DRef, eps_tilde_1_0_A, eps_tilde_1_0_i, eps_tilde_1_A, eps_tilde_1_i, eps_tilde_2_0_A, eps_tilde_2_0_i, eps_tilde_2_A, eps_tilde_2_i, epsilon, epsilon_hat, G_i, DofOwnerOwner::iGetFirstIndex(), StructDispNode::iGetFirstMomentumIndex(), iGetNumDof(), InterpolateOrientation(), k_1_i, k_2_i, k_tilde_1_0_i, k_tilde_1_i, k_tilde_2_0_i, k_tilde_2_i, Kappa_delta_i_1, Kappa_delta_i_2, L_alpha_beta_A, L_alpha_beta_i, MatCross, MatrixHandler::MatMatMul(), mb_deye< Mat3x3 >(), Mat3x3::MulTV(), NUMIP, NUMNODES, NUMSSEP, P_i, Phi_Delta_A, Phi_Delta_i, phi_tilde_n, pNode, PreStress, MyVectorHandler::Put(), FullMatrixHandler::Put(), SubVectorHandler::PutRowIndex(), FullMatrixHandler::PutT(), FullMatrixHandler::Reset(), VectorHandler::ResizeReset(), stress_i, T_A, T_i, Vec3::Tens(), UpdateNodalAndAveragePosAndOrientation(), w_i, WorkSpaceDim(), xa, xi_A, xi_i, y_i_1, y_i_2, and Zero3.

 {
         /* Dimensiona e resetta la matrice di lavoro */
         integer iNumRows = 0;
         integer iNumCols = 0;
         WorkSpaceDim(&iNumRows, &iNumCols);
         WorkVec.ResizeReset(iNumRows);
 
         /* Recupera e setta gli indici */
         for (integer i = 0; i < 4; i++) {
                 integer iNodeFirstMomIndex = pNode[i]->iGetFirstMomentumIndex();
                 for (int iCnt = 1; iCnt <= 6; iCnt++) {
                         WorkVec.PutRowIndex(iCnt + 6 * i, iNodeFirstMomIndex + iCnt);
                 }
         }
 
         integer iFirstReactionIndex = iGetFirstIndex();
         for (unsigned int iCnt = 1; iCnt <= iGetNumDof(); iCnt++) {
                 WorkVec.PutRowIndex(24 + iCnt, iFirstReactionIndex + iCnt);
         }
 
         UpdateNodalAndAveragePosAndOrientation();
         InterpolateOrientation();
 //      ComputeIPSEPRotations();
         for (unsigned int i = 1; i <= iGetNumDof(); i++) {
                 beta(i) = XCurr(iFirstReactionIndex + i);
         }
 
 
         ComputeIPCurvature();
         for (integer i = 0; i < NUMIP; i++) {
                 InterpDeriv(xa, L_alpha_beta_i[i], y_i_1[i], y_i_2[i]);
                 eps_tilde_1_i[i] = T_i[i].MulTV(y_i_1[i]) - eps_tilde_1_0_i[i];
                 eps_tilde_2_i[i] = T_i[i].MulTV(y_i_2[i]) - eps_tilde_2_0_i[i];
                 k_tilde_1_i[i] = T_i[i].MulTV(k_1_i[i]) - k_tilde_1_0_i[i];
                 k_tilde_2_i[i] = T_i[i].MulTV(k_2_i[i]) - k_tilde_2_0_i[i];
                 // parte variabile di B_overline_i
                 for (integer n = 0; n < NUMNODES; n++) {
                         Mat3x3 Phi_Delta_i_n_LI_i = Phi_Delta_i[i][n] * LI[n](xi_i[i]);
 
                         // delta epsilon_tilde_1_i
                         B_overline_i[i].PutT(1, 1 + 6 * n, T_i[i] * L_alpha_beta_i[i](n + 1, 1));
                         B_overline_i[i].Put(1, 4 + 6 * n, T_i[i].MulTM(Mat3x3(MatCross, y_i_1[i])) * Phi_Delta_i_n_LI_i);
 
                         // delta epsilon_tilde_2_i
                         B_overline_i[i].PutT(4, 1 + 6 * n, T_i[i] * L_alpha_beta_i[i](n + 1, 2));
                         B_overline_i[i].Put(4, 4 + 6 * n, T_i[i].MulTM(Mat3x3(MatCross, y_i_2[i])) * Phi_Delta_i_n_LI_i);
 
                         // delta k_tilde_1_i
                         Vec3 phi_tilde_1_i(Zero3);
                         Vec3 phi_tilde_2_i(Zero3);
                         InterpDeriv(phi_tilde_n, L_alpha_beta_i[i], phi_tilde_1_i, phi_tilde_2_i);
                         B_overline_i[i].Put(7, 4 + 6 * n,
                                 T_i[i].MulTM(Mat3x3(MatCross, k_1_i[i])) * Phi_Delta_i_n_LI_i
                                 +
                                 T_i[i].MulTM(Kappa_delta_i_1[i][n])
                         );
 
                         // delta k_tilde_2_i
                         B_overline_i[i].Put(10, 4 + 6 * n, 
                                 T_i[i].MulTM(Mat3x3(MatCross, k_2_i[i])) * Phi_Delta_i_n_LI_i
                                 +
                                 T_i[i].MulTM(Kappa_delta_i_2[i][n])
                         );
 
                         // delta y_alpha_1
                         D_overline_i[i].Put(1, 1 + n * 6, mb_deye<Mat3x3>(L_alpha_beta_i[i](n + 1, 1)));
 
                         // delta y_alpha_2
                         D_overline_i[i].Put(4, 1 + n * 6, mb_deye<Mat3x3>(L_alpha_beta_i[i](n + 1, 2)));
 
                         // delta k_1_i
                         D_overline_i[i].Put(7, 4 + n * 6, Kappa_delta_i_1[i][n]);
 
                         // delta k_1_i
                         D_overline_i[i].Put(10, 4 + n * 6, Kappa_delta_i_2[i][n]);
 
                         // phi_delta
                         D_overline_i[i].Put(13, 4 + n * 6, Phi_Delta_i_n_LI_i);
                 }
         }
         // ANS
         {
                 fmh B_overline_A(6, 24);
                 Vec3 y_A_1;
                 Vec3 y_A_2;
                 fmh B_overline_3_ABCD(4, 24);
                 fmh B_overline_6_ABCD(4, 24);
                 for (integer i = 0; i < NUMSSEP; i++) {
                         B_overline_A.Reset();
                         InterpDeriv(xa, L_alpha_beta_A[i], y_A_1, y_A_2);
                         eps_tilde_1_A[i] = T_A[i].MulTV(y_A_1) - eps_tilde_1_0_A[i];
                         eps_tilde_2_A[i] = T_A[i].MulTV(y_A_2) - eps_tilde_2_0_A[i];
                         for (integer n = 0; n < NUMNODES; n++) {
                                 Mat3x3 Phi_Delta_A_n_LI_i = Phi_Delta_A[i][n] * LI[n](xi_A[i]);
 
                                 // delta epsilon_tilde_1_A
                                 B_overline_A.PutT(1, 1 + 6 * n, T_A[i] * L_alpha_beta_A[i](n + 1, 1));
                                 B_overline_A.Put(1, 4 + 6 * n, 
                                         T_A[i].MulTM(Mat3x3(MatCross, y_A_1)) * Phi_Delta_A_n_LI_i
                                 );
 
                                 // delta epsilon_tilde_2_A
                                 B_overline_A.PutT(4, 1 + 6 * n, T_A[i] * L_alpha_beta_A[i](n + 1, 2));
                                 B_overline_A.Put(4, 4 + 6 * n, 
                                         T_A[i].MulTM(Mat3x3(MatCross, y_A_2)) * Phi_Delta_A_n_LI_i
                                 );
                         }
 #if 0
                         CopyMatrixRow(B_overline_3_ABCD, i + 1, B_overline_A, 3);
                         CopyMatrixRow(B_overline_6_ABCD, i + 1, B_overline_A, 6);
 #endif
                         B_overline_3_ABCD.CopyMatrixRow(i + 1, B_overline_A, 3);
                         B_overline_6_ABCD.CopyMatrixRow(i + 1, B_overline_A, 6);
                 }
                 fmh tmp_B_ANS(1, 24);
                 for (integer i = 0; i < NUMIP; i++) {
                         FullMatrixHandler sh1(1, 4);
                         FullMatrixHandler sh2(1, 4);
                                 sh1(1, 1) = (1. + xi_i[i][1]) * 0.5;
                                 sh1(1, 3) = (1. - xi_i[i][1]) * 0.5;
                                 sh2(1, 4) = (1. + xi_i[i][0]) * 0.5;
                                 sh2(1, 2) = (1. - xi_i[i][0]) * 0.5;
 
                         eps_tilde_1_i[i](3) = 
                                 sh1(1, 1) * eps_tilde_1_A[0](3) + 
                                 sh1(1, 3) * eps_tilde_1_A[2](3)
                         ;
                         eps_tilde_2_i[i](3) = 
                                 sh2(1, 2) * eps_tilde_2_A[1](3) + 
                                 sh2(1, 4) * eps_tilde_2_A[3](3)
                         ;
                         
                         sh1.MatMatMul(tmp_B_ANS, B_overline_3_ABCD);
 #if 0
                         CopyMatrixRow(B_overline_i[i], 3, tmp_B_ANS, 1);
 #endif
                         B_overline_i[i].CopyMatrixRow(3, tmp_B_ANS, 1);
                         
                         sh2.MatMatMul(tmp_B_ANS, B_overline_6_ABCD);
 #if 0
                         CopyMatrixRow(B_overline_i[i], 6, tmp_B_ANS, 1);
 #endif
                         B_overline_i[i].CopyMatrixRow(6, tmp_B_ANS, 1);
                 }
         }
         
         // EAS: B membranali
 //      {
 //              integer tmpidx1[5] = {0, 1, 5, 4, 2};
 //              for (integer i = 0; i < NUMIP; i++) {
 //                      for (integer n = 1; n <= 4; n++) {
 //#if 0
 //                              CopyMatrixRow(B_overline_m_i[i], n, B_overline_i[i], tmpidx1[n]);
 //#endif
 //                              B_overline_m_i[i].CopyMatrixRow(n, B_overline_i[i], tmpidx1[n]);
 //                      }
 //              }
 //      }
 
         /* Calcola le azioni interne */
         for (integer i = 0; i < NUMIP; i++) {
                 epsilon.Put(1, eps_tilde_1_i[i]);
                 epsilon.Put(4, eps_tilde_2_i[i]);
                 epsilon.Put(7, k_tilde_1_i[i]);
                 epsilon.Put(10, k_tilde_2_i[i]);
                 // TODO: recupera epsilon_hat con l'ordine giusto per qua
                 P_i[i].MatVecMul(epsilon_hat, beta);
 
                 epsilon += epsilon_hat;
 #ifdef USE_CL_IN_SHELL
                 pD[i]->Update(epsilon);
                 stress_i[i] = pD[i]->GetF();
 #else // ! USE_CL_IN_SHELL
                 DRef[i].MatVecMul(stress_i[i], epsilon);
                 if (bPreStress) {
                         stress_i[i] += PreStress;
                 }
 #endif // ! USE_CL_IN_SHELL
                 
                 Vec3 n1, n2, m1, m2;
                 ExtractVec3(n1, stress_i[i],  1);
                 ExtractVec3(n2, stress_i[i],  4);
                 ExtractVec3(m1, stress_i[i],  7);
                 ExtractVec3(m2, stress_i[i], 10);
                 
                 Mat3x3 Hh;
                 Vec3 Tn1 = T_i[i] * n1;
                 Vec3 Tn2 = T_i[i] * n2;
                 Vec3 Tm1 = T_i[i] * m1;
                 Vec3 Tm2 = T_i[i] * m2;
                 Hh = Tn1.Tens(y_i_1[i]) - mb_deye<Mat3x3>(Tn1.Dot(y_i_1[i]))
                         + Tn2.Tens(y_i_2[i]) - mb_deye<Mat3x3>(Tn2.Dot(y_i_2[i]))
                         + Tm1.Tens(k_1_i[i]) - mb_deye<Mat3x3>(Tm1.Dot(k_1_i[i]))
                         + Tm2.Tens(k_2_i[i]) - mb_deye<Mat3x3>(Tm2.Dot(k_2_i[i]))
                         ;
 
                 // NOTE: use PutCross()?
                 G_i[i].PutT(1, 13, Mat3x3(MatCross, Tn1));
                 G_i[i].PutT(4, 13, Mat3x3(MatCross, Tn2));
                 G_i[i].Put(13, 7, Mat3x3(MatCross, Tm1));
                 G_i[i].Put(13, 10, Mat3x3(MatCross, Tm2));
                 G_i[i].Put(13, 1, Mat3x3(MatCross, Tn1));
                 G_i[i].Put(13, 4, Mat3x3(MatCross, Tn2));
                 G_i[i].Put(13, 13, Hh);
         }
         
         //Residuo
         //forze di volume
         MyVectorHandler rd(24);
         MyVectorHandler rbeta(iGetNumDof());
         for (integer i = 0; i < NUMIP; i++) {
                 B_overline_i[i].MatTVecMul(rd, stress_i[i]);
                 P_i[i].MatTVecMul(rbeta, stress_i[i]);
                 
                 AssembleVector(WorkVec,  1, rd, -alpha_i[i] * w_i[i]);
                 AssembleVector(WorkVec, 25, rbeta, -alpha_i[i] * w_i[i] / dCoef);
         }
 
         return WorkVec;
 }

Here is the call graph for this function:

void Shell4EASANS::ComputeInitialNodeOrientation ( void )

private

Definition at line 138 of file shelleasans.cc.

References Vec3::Cross(), Eye3, StructNode::GetRCurr(), StructDispNode::GetXCurr(), InterpolateOrientation(), iTa, iTa_A, iTa_i, Vec3::Norm(), NUMIP, NUMNODES, NUMSSEP, pNode, T_A, T_i, UpdateNodalAndAveragePosAndOrientation(), xa, xi_A, xi_i, and xi_n.

Referenced by Shell4EASANS().

 {
         for (integer i = 0; i < NUMNODES; i++) {
                 xa[i] = pNode[i]->GetXCurr();
         }
         for (integer i = 0; i < NUMNODES; i++) {
                 Vec3 t1 = InterpDeriv_xi1(xa, xi_n[i]);
                 t1 = t1 / t1.Norm();
                 Vec3 t2 = InterpDeriv_xi2(xa, xi_n[i]);
                 t2 = t2 / t2.Norm();
                 Vec3 t3 = t1.Cross(t2);
                 t3 = t3 / t3.Norm();
                 t2 = t3.Cross(t1);
                 iTa[i] = (pNode[i]->GetRCurr()).MulTM(Mat3x3(t1, t2, t3));
         }
         for (integer i = 0; i < NUMIP; i++) {
                 iTa_i[i] = Eye3;
         }
         for (integer i = 0; i < NUMSSEP; i++) {
                 iTa_A[i] = Eye3;
         }
         UpdateNodalAndAveragePosAndOrientation();
         InterpolateOrientation();
         for (integer i = 0; i < NUMIP; i++) {
                 Vec3 t1 = InterpDeriv_xi1(xa, xi_i[i]);
                 t1 = t1 / t1.Norm();
                 Vec3 t2 = InterpDeriv_xi2(xa, xi_i[i]);
                 t2 = t2 / t2.Norm();
                 Vec3 t3 = t1.Cross(t2);
                 t3 = t3 / t3.Norm();
                 t2 = t3.Cross(t1);
                 iTa_i[i] = (T_i[i]).MulTM(Mat3x3(t1, t2, t3));
         }
         for (integer i = 0; i < NUMSSEP; i++) {
                 Vec3 t1 = InterpDeriv_xi1(xa, xi_A[i]);
                 t1 = t1 / t1.Norm();
                 Vec3 t2 = InterpDeriv_xi2(xa, xi_A[i]);
                 t2 = t2 / t2.Norm();
                 Vec3 t3 = t1.Cross(t2);
                 t3 = t3 / t3.Norm();
                 t2 = t3.Cross(t1);
                 iTa_A[i] = (T_A[i]).MulTM(Mat3x3(t1, t2, t3));
         }
         InterpolateOrientation();
 }

Here is the call graph for this function:

void Shell4EASANS::ComputeIPCurvature ( void )

private

Definition at line 229 of file shelleasans.cc.

References RotManip::DRot(), RotManip::DRot_I(), RotManip::Elle(), k_1_i, k_2_i, Kappa_delta_i_1, Kappa_delta_i_2, L_alpha_beta_i, Mat3x3::MulMT(), NUMIP, NUMNODES, Phi_Delta_i, phi_tilde_i, phi_tilde_n, T_overline, xi_i, and Zero3.

Referenced by AssRes(), and Shell4EASANS().

 {
         Mat3x3 Gamma_I_n_MT_T_overline[NUMNODES];
         for (integer n = 0; n < NUMNODES; n++) {
                 Gamma_I_n_MT_T_overline[n] = RotManip::DRot_I(phi_tilde_n[n]).MulMT(T_overline);
         }
         for (integer i = 0; i < NUMIP; i++) {
                 Vec3 phi_tilde_1_i(Zero3);
                 Vec3 phi_tilde_2_i(Zero3);
                 InterpDeriv(phi_tilde_n, L_alpha_beta_i[i], phi_tilde_1_i, phi_tilde_2_i);
                 Mat3x3 T_overlineGamma_tilde_i(T_overline * RotManip::DRot(phi_tilde_i[i]));
                 k_1_i[i] = T_overlineGamma_tilde_i * phi_tilde_1_i;
                 k_2_i[i] = T_overlineGamma_tilde_i * phi_tilde_2_i;
                 Mat3x3 tmp1 = T_overline * RotManip::Elle(phi_tilde_i[i], phi_tilde_1_i);
                 Mat3x3 tmp2 = T_overline * RotManip::Elle(phi_tilde_i[i], phi_tilde_2_i);
                 for (int n = 0; n < NUMNODES; n++) {
                         Kappa_delta_i_1[i][n] = tmp1 * Gamma_I_n_MT_T_overline[n] * LI[n](xi_i[i]) + 
                                 Phi_Delta_i[i][n] * L_alpha_beta_i[i](n + 1, 1);
                         Kappa_delta_i_2[i][n] = tmp2 * Gamma_I_n_MT_T_overline[n] * LI[n](xi_i[i]) + 
                                 Phi_Delta_i[i][n] * L_alpha_beta_i[i](n + 1, 2);
                 }
         }
 }

Here is the call graph for this function:

virtual void Shell4EASANS::GetConnectedNodes ( std::vector< const Node * > & connectedNodes ) const

inlinevirtual

Reimplemented from Elem.

Definition at line 452 of file shelleasans.h.

References NUMNODES, and pNode.

                                                                          {
                 connectedNodes.resize(NUMNODES);
                 for (int i = 0; i < NUMNODES; i++) {
                         connectedNodes[i] = pNode[i];
                 }
         };

virtual DofOrder::Order Shell4EASANS::GetDofType ( unsigned int i ) const

inlinevirtual

Reimplemented from Elem.

Definition at line 341 of file shelleasans.h.

References DofOrder::ALGEBRAIC, ASSERT, and iGetNumDof().

                                                                {
                 ASSERT(i >= 0 && i < iGetNumDof());
                 return DofOrder::ALGEBRAIC;
         };

Here is the call graph for this function:

virtual Elem::Type Shell4EASANS::GetElemType ( void ) const

inlinevirtual

Implements Elem.

Definition at line 332 of file shelleasans.h.

References Elem::PLATE.

                                                  {
                 return Elem::PLATE;
         };

virtual DofOrder::Order Shell4EASANS::GetEqType ( unsigned int i ) const

inlinevirtual

Reimplemented from SimulationEntity.

Definition at line 346 of file shelleasans.h.

References DofOrder::ALGEBRAIC, ASSERT, and iGetNumDof().

                                                               {
                 ASSERT(i >= 0 && i < iGetNumDof());
                 return DofOrder::ALGEBRAIC;
         };

Here is the call graph for this function:

virtual Shell::Type Shell4EASANS::GetShellType ( void ) const

inlinevirtual

Definition at line 327 of file shelleasans.h.

References Shell::ELASTIC.

                                                    {
                 return Shell::ELASTIC;
         };

virtual unsigned int Shell4EASANS::iGetInitialNumDof ( void ) const

inlinevirtual

Implements SubjectToInitialAssembly.

Definition at line 417 of file shelleasans.h.

                                                            {
                 return 0;
         };

virtual unsigned int Shell4EASANS::iGetNumDof ( void ) const

inlinevirtual

Reimplemented from Elem.

Definition at line 336 of file shelleasans.h.

Referenced by AssJac(), AssRes(), GetDofType(), GetEqType(), Shell4EASANS(), and WorkSpaceDim().

                                                     { 
 //              return 8;
                 return 7;
         };

VariableSubMatrixHandler & Shell4EASANS::InitialAssJac	(	VariableSubMatrixHandler &	WorkMat,
		const VectorHandler &	XCurr
	)

virtual

Implements SubjectToInitialAssembly.

Definition at line 817 of file shelleasans.cc.

References VariableSubMatrixHandler::SetNullMatrix().

 {
         WorkMat.SetNullMatrix();
         return WorkMat;
 }

Here is the call graph for this function:

SubVectorHandler & Shell4EASANS::InitialAssRes	(	SubVectorHandler &	WorkVec,
		const VectorHandler &	XCurr
	)

virtual

Implements SubjectToInitialAssembly.

Definition at line 825 of file shelleasans.cc.

References VectorHandler::Resize().

 {
         WorkVec.Resize(0);
         return WorkVec;
 }

Here is the call graph for this function:

virtual void Shell4EASANS::InitialWorkSpaceDim	(	integer *	piNumRows,
		integer *	piNumCols
	)		const

inlinevirtual

Implements SubjectToInitialAssembly.

Definition at line 422 of file shelleasans.h.

                                                                           {
                 *piNumRows = 6*4;
                 *piNumCols = 6*4;
         };

void Shell4EASANS::InterpolateOrientation ( void )

private

Definition at line 185 of file shelleasans.cc.

References RotManip::DRot_I(), iTa_A, iTa_i, Mat3x3::MulMT(), NUMIP, NUMNODES, NUMSSEP, Phi_Delta_A, Phi_Delta_i, phi_tilde_0, phi_tilde_A, phi_tilde_i, phi_tilde_n, RotManip::Rot(), RotManip::RotAndDRot(), T_0, T_A, T_i, T_overline, xi_0, xi_A, and xi_i.

Referenced by AssRes(), and ComputeInitialNodeOrientation().

 {
         Mat3x3 DRot_I_phi_tilde_n_MT_T_overline[NUMNODES];
         Mat3x3 Ri, Gammai;
         for (integer n = 0; n < NUMNODES; n++) {
                 DRot_I_phi_tilde_n_MT_T_overline[n] = 
                         RotManip::DRot_I(phi_tilde_n[n]).MulMT(T_overline);
         }
         for (integer i = 0; i < NUMIP; i++) {
                 phi_tilde_i[i] = Interp(phi_tilde_n, xi_i[i]);
                 RotManip::RotAndDRot(phi_tilde_i[i], Ri, Gammai);
                 T_i[i] = T_overline * Ri * iTa_i[i];
                 Mat3x3 T_overline_Gamma_tilde_i(T_overline * Gammai);
                 for (int n = 0; n < NUMNODES; n++) {
                         Phi_Delta_i[i][n] = T_overline_Gamma_tilde_i * 
                                 DRot_I_phi_tilde_n_MT_T_overline[n];
                 }
         }
         Vec3 phi_tilde_0 = Interp(phi_tilde_n, xi_0);
         T_0 = T_overline * RotManip::Rot(phi_tilde_0);
         for (integer i = 0; i < NUMSSEP; i++) {
                 phi_tilde_A[i] = Interp(phi_tilde_n, xi_A[i]);
                 RotManip::RotAndDRot(phi_tilde_A[i], Ri, Gammai);
                 T_A[i] = T_overline * Ri * iTa_A[i];
                 Mat3x3 T_overline_Gamma_tilde_A(T_overline * Gammai);
                 for (int n = 0; n < NUMNODES; n++) {
                         Phi_Delta_A[i][n] = T_overline_Gamma_tilde_A * 
                                 DRot_I_phi_tilde_n_MT_T_overline[n];
                 }
         }
 }

Here is the call graph for this function:

void Shell4EASANS::Output ( OutputHandler & OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 848 of file shelleasans.cc.

References ToBeOutput::bToBeOutput(), WithLabel::GetLabel(), NUMIP, OutputHandler::PLATES, OutputHandler::Plates(), stress_i, and OutputHandler::UseText().

 {
         if (bToBeOutput()) {
                 if (OH.UseText(OutputHandler::PLATES)) {
                         std::ostream& out = OH.Plates();
                         out << std::setw(8) << GetLabel();
                                 // TODO: complete
                         for (integer i = 0; i < NUMIP; i++) {
                                 for (integer r = 1; r <= 12; r++) {
                                         out << " " << stress_i[i](r);
                                 }
                         }
                         out << std::endl;
                 }
         }
 }

Here is the call graph for this function:

const StructNode * Shell4EASANS::pGetNode ( unsigned int i ) const

virtual

Definition at line 833 of file shelleasans.cc.

References ASSERT, MBDYN_EXCEPT_ARGS, and pNode.

 {
         ASSERT(i >= 1 && i <= 4);
         switch (i) {
         case 1:
         case 2:
         case 3:
         case 4:
                 return pNode[i - 1];
         default:
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 }

std::ostream & Shell4EASANS::Restart ( std::ostream & out ) const

virtual

Implements Elem.

Definition at line 802 of file shelleasans.cc.

 {
         return out << "# not implemented yet" << std::endl;
 }

virtual void Shell4EASANS::SetInitialValue ( VectorHandler & X )

inlinevirtual

Initialize state vector used in initial assembly. May set internal states of the element. Do not rely on being always called, because initial assembly could be implicitly or explicitly skipped

Reimplemented from DofOwnerOwner.

Definition at line 427 of file shelleasans.h.

References NO_OP.

                                                         {
                 NO_OP;
         };

void Shell4EASANS::SetValue	(	DataManager *	pDM,
		VectorHandler &	,
		VectorHandler &	,
		SimulationEntity::Hints *	ph = `0`
	)

virtual

Reimplemented from SimulationEntity.

Definition at line 809 of file shelleasans.cc.

References NO_OP.

 {
         NO_OP;
 }

void Shell4EASANS::UpdateNodalAndAveragePosAndOrientation ( void )

private

Definition at line 116 of file shelleasans.cc.

References StructNode::GetRCurr(), StructNode::GetRRef(), StructDispNode::GetXCurr(), iTa, Mat3x3::MulTM(), NUMNODES, phi_tilde_n, pNode, RotManip::Rot(), T_overline, RotManip::VecRot(), xa, and Zero3x3.

Referenced by AssRes(), and ComputeInitialNodeOrientation().

 {
         Mat3x3 T_avg(Zero3x3);
         Mat3x3 Tn[NUMNODES];
         Mat3x3 R_tilde_n[NUMNODES];
         for (integer i = 0; i < NUMNODES; i++) {
 //              xa[i] = pNode[i]->GetXCurr()+pNode[i]->GetRCurr()*f[i];
                 xa[i] = pNode[i]->GetXCurr();
                 Tn[i] = pNode[i]->GetRCurr() * iTa[i];
 //              T_a_avg += Ta[i];
                 T_avg += pNode[i]->GetRRef() * iTa[i];
         }
         T_avg /= 4.;
         //FIXME; alternative solution: polar decomposition of T_a_avg = T0*U
         T_overline = RotManip::Rot(RotManip::VecRot(T_avg));
         for (integer i = 0; i < NUMNODES; i++) {
                 R_tilde_n[i] = T_overline.MulTM(Tn[i]);
                 phi_tilde_n[i] = RotManip::VecRot(R_tilde_n[i]);
         }
 }

Here is the call graph for this function:

virtual void Shell4EASANS::WorkSpaceDim	(	integer *	piNumRows,
		integer *	piNumCols
	)		const

inlinevirtual

Implements Elem.

Definition at line 366 of file shelleasans.h.

References iGetNumDof().

Referenced by AssJac(), and AssRes().

                                                                    {
                 *piNumRows = 6*4 + iGetNumDof();
                 *piNumCols = 6*4 + iGetNumDof();
         };

Here is the call graph for this function:

Member Data Documentation

doublereal Shell4EASANS::alpha_0

protected

Definition at line 252 of file shelleasans.h.

Referenced by Shell4EASANS().

doublereal Shell4EASANS::alpha_i[NUMIP]

protected

Definition at line 253 of file shelleasans.h.

Referenced by AssJac(), AssRes(), and Shell4EASANS().

vfmh Shell4EASANS::B_overline_i

protected

Definition at line 257 of file shelleasans.h.

Referenced by AssJac(), and AssRes().

vh Shell4EASANS::beta

protected

Definition at line 269 of file shelleasans.h.

Referenced by AssRes().

bool Shell4EASANS::bFirstRes

protected

Definition at line 289 of file shelleasans.h.

bool Shell4EASANS::bPreStress

protected

Definition at line 278 of file shelleasans.h.

Referenced by AssRes().

vfmh Shell4EASANS::D_overline_i

protected

Definition at line 259 of file shelleasans.h.

Referenced by AssJac(), and AssRes().

vfmh Shell4EASANS::DRef

protected

Definition at line 283 of file shelleasans.h.

Referenced by AssJac(), AssRes(), and Shell4EASANS().

Vec3 Shell4EASANS::eps_tilde_1_0_A[NUMSSEP]

protected

Definition at line 215 of file shelleasans.h.

Referenced by AssRes(), and Shell4EASANS().

Vec3 Shell4EASANS::eps_tilde_1_0_i[NUMIP]

protected

Definition at line 213 of file shelleasans.h.

Referenced by AssRes(), and Shell4EASANS().

Vec3 Shell4EASANS::eps_tilde_1_A[NUMSSEP]

protected

Definition at line 220 of file shelleasans.h.

Referenced by AssRes().

Vec3 Shell4EASANS::eps_tilde_1_i[NUMIP]

protected

Definition at line 218 of file shelleasans.h.

Referenced by AssRes().

Vec3 Shell4EASANS::eps_tilde_2_0_A[NUMSSEP]

protected

Definition at line 216 of file shelleasans.h.

Referenced by AssRes(), and Shell4EASANS().

Vec3 Shell4EASANS::eps_tilde_2_0_i[NUMIP]

protected

Definition at line 214 of file shelleasans.h.

Referenced by AssRes(), and Shell4EASANS().

Vec3 Shell4EASANS::eps_tilde_2_A[NUMSSEP]

protected

Definition at line 221 of file shelleasans.h.

Referenced by AssRes().

Vec3 Shell4EASANS::eps_tilde_2_i[NUMIP]

protected

Definition at line 219 of file shelleasans.h.

Referenced by AssRes().

vh Shell4EASANS::epsilon

protected

Definition at line 271 of file shelleasans.h.

Referenced by AssRes().

vh Shell4EASANS::epsilon_hat

protected

Definition at line 270 of file shelleasans.h.

Referenced by AssRes().

vvh Shell4EASANS::FRef

protected

Definition at line 277 of file shelleasans.h.

vfmh Shell4EASANS::G_i

protected

Definition at line 260 of file shelleasans.h.

Referenced by AssJac(), and AssRes().

Mat3x3 Shell4EASANS::iTa[NUMNODES]

protected

Definition at line 165 of file shelleasans.h.

Referenced by ComputeInitialNodeOrientation(), and UpdateNodalAndAveragePosAndOrientation().

Mat3x3 Shell4EASANS::iTa_A[NUMSSEP]

protected

Definition at line 167 of file shelleasans.h.

Referenced by ComputeInitialNodeOrientation(), and InterpolateOrientation().

Mat3x3 Shell4EASANS::iTa_i[NUMIP]

protected

Definition at line 166 of file shelleasans.h.

Referenced by ComputeInitialNodeOrientation(), and InterpolateOrientation().

Vec3 Shell4EASANS::k_1_i[NUMIP]

protected

Definition at line 206 of file shelleasans.h.

Referenced by AssRes(), ComputeIPCurvature(), and Shell4EASANS().

Vec3 Shell4EASANS::k_2_i[NUMIP]

protected

Definition at line 207 of file shelleasans.h.

Referenced by AssRes(), ComputeIPCurvature(), and Shell4EASANS().

vfmh Shell4EASANS::K_beta_beta_i

protected

Definition at line 264 of file shelleasans.h.

Vec3 Shell4EASANS::k_tilde_1_0_i[NUMIP]

protected

Definition at line 225 of file shelleasans.h.

Referenced by AssRes(), and Shell4EASANS().

Vec3 Shell4EASANS::k_tilde_1_i[NUMIP]

protected

Definition at line 228 of file shelleasans.h.

Referenced by AssRes().

Vec3 Shell4EASANS::k_tilde_2_0_i[NUMIP]

protected

Definition at line 226 of file shelleasans.h.

Referenced by AssRes(), and Shell4EASANS().

Vec3 Shell4EASANS::k_tilde_2_i[NUMIP]

protected

Definition at line 229 of file shelleasans.h.

Referenced by AssRes().

Mat3x3 Shell4EASANS::Kappa_delta_i_1[NUMIP][NUMNODES]

protected

Definition at line 193 of file shelleasans.h.

Referenced by AssRes(), and ComputeIPCurvature().

Mat3x3 Shell4EASANS::Kappa_delta_i_2[NUMIP][NUMNODES]

protected

Definition at line 194 of file shelleasans.h.

Referenced by AssRes(), and ComputeIPCurvature().

vfmh Shell4EASANS::L_alpha_beta_A

protected

Definition at line 255 of file shelleasans.h.

Referenced by AssRes(), and Shell4EASANS().

vfmh Shell4EASANS::L_alpha_beta_i

protected

Definition at line 254 of file shelleasans.h.

Referenced by AssRes(), ComputeIPCurvature(), and Shell4EASANS().

vfmh Shell4EASANS::P_i

protected

Definition at line 262 of file shelleasans.h.

Referenced by AssJac(), AssRes(), and Shell4EASANS().

Mat3x3 Shell4EASANS::Phi_Delta_A[NUMIP][NUMNODES]

protected

Definition at line 192 of file shelleasans.h.

Referenced by AssRes(), and InterpolateOrientation().

Mat3x3 Shell4EASANS::Phi_Delta_i[NUMIP][NUMNODES]

protected

Definition at line 191 of file shelleasans.h.

Referenced by AssRes(), ComputeIPCurvature(), and InterpolateOrientation().

Vec3 Shell4EASANS::phi_tilde_0

protected

Definition at line 174 of file shelleasans.h.

Referenced by InterpolateOrientation().

Vec3 Shell4EASANS::phi_tilde_A[NUMSSEP]

protected

Definition at line 173 of file shelleasans.h.

Referenced by InterpolateOrientation().

Vec3 Shell4EASANS::phi_tilde_i[NUMIP]

protected

Definition at line 172 of file shelleasans.h.

Referenced by ComputeIPCurvature(), and InterpolateOrientation().

Vec3 Shell4EASANS::phi_tilde_n[NUMNODES]

protected

Definition at line 169 of file shelleasans.h.

Referenced by AssRes(), ComputeIPCurvature(), InterpolateOrientation(), and UpdateNodalAndAveragePosAndOrientation().

const StructNode* Shell4EASANS::pNode[NUMNODES]

protected

Definition at line 153 of file shelleasans.h.

Referenced by AssJac(), AssRes(), ComputeInitialNodeOrientation(), GetConnectedNodes(), pGetNode(), Shell4EASANS(), and UpdateNodalAndAveragePosAndOrientation().

vh Shell4EASANS::PreStress

protected

Definition at line 279 of file shelleasans.h.

Referenced by AssRes().

Mat3x3 Shell4EASANS::Q_A[NUMSSEP]

protected

Definition at line 203 of file shelleasans.h.

Mat3x3 Shell4EASANS::Q_i[NUMIP]

protected

Definition at line 202 of file shelleasans.h.

fmh Shell4EASANS::S_alpha_beta_0

protected

Definition at line 249 of file shelleasans.h.

Referenced by Shell4EASANS().

vfmh Shell4EASANS::S_alpha_beta_A

protected

Definition at line 251 of file shelleasans.h.

Referenced by Shell4EASANS().

vfmh Shell4EASANS::S_alpha_beta_i

protected

Definition at line 250 of file shelleasans.h.

Referenced by Shell4EASANS().

vvh Shell4EASANS::stress_i

protected

Definition at line 286 of file shelleasans.h.

Referenced by AssRes(), and Output().

Mat3x3 Shell4EASANS::T0_overline

protected

Definition at line 177 of file shelleasans.h.

Referenced by Shell4EASANS().

Mat3x3 Shell4EASANS::T_0

protected

Definition at line 187 of file shelleasans.h.

Referenced by InterpolateOrientation(), and Shell4EASANS().

Mat3x3 Shell4EASANS::T_0_0

protected

Definition at line 183 of file shelleasans.h.

Referenced by Shell4EASANS().

Mat3x3 Shell4EASANS::T_0_A[NUMSSEP]

protected

Definition at line 185 of file shelleasans.h.

Referenced by Shell4EASANS().

Mat3x3 Shell4EASANS::T_0_i[NUMIP]

protected

Definition at line 184 of file shelleasans.h.

Referenced by Shell4EASANS().

Mat3x3 Shell4EASANS::T_A[NUMSSEP]

protected

Definition at line 189 of file shelleasans.h.

Referenced by AssRes(), ComputeInitialNodeOrientation(), InterpolateOrientation(), and Shell4EASANS().

Mat3x3 Shell4EASANS::T_i[NUMIP]

protected

Definition at line 188 of file shelleasans.h.

Referenced by AssRes(), ComputeInitialNodeOrientation(), InterpolateOrientation(), and Shell4EASANS().

Mat3x3 Shell4EASANS::T_overline

protected

Definition at line 179 of file shelleasans.h.

Referenced by ComputeIPCurvature(), InterpolateOrientation(), Shell4EASANS(), and UpdateNodalAndAveragePosAndOrientation().

doublereal Shell4EASANS::w_i

static

Initial value:

Definition at line 117 of file shelleasans.h.

Referenced by AssJac(), and AssRes().

Vec3 Shell4EASANS::xa[NUMNODES]

protected

Definition at line 163 of file shelleasans.h.

Referenced by AssRes(), ComputeInitialNodeOrientation(), Shell4EASANS(), and UpdateNodalAndAveragePosAndOrientation().

Vec3 Shell4EASANS::xa_0[NUMNODES]

protected

Definition at line 162 of file shelleasans.h.

Referenced by Shell4EASANS().

doublereal Shell4EASANS::xi_0 = {0., 0.}

static

Definition at line 142 of file shelleasans.h.

Referenced by InterpolateOrientation(), and Shell4EASANS().

doublereal Shell4EASANS::xi_A

static

Initial value:

=  {
        { 0.,  1.},
        {-1.,  0.},
        { 0., -1.},
        { 1.,  0.}
}

Definition at line 129 of file shelleasans.h.

Referenced by AssRes(), ComputeInitialNodeOrientation(), InterpolateOrientation(), and Shell4EASANS().

doublereal Shell4EASANS::xi_i

static

Initial value:

= {
        {-1. / std::sqrt(3.), -1. / std::sqrt(3.)},
        { 1. / std::sqrt(3.), -1. / std::sqrt(3.)},
        { 1. / std::sqrt(3.),  1. / std::sqrt(3.)},
        {-1. / std::sqrt(3.),  1. / std::sqrt(3.)}
}

Definition at line 116 of file shelleasans.h.

Referenced by AssRes(), ComputeInitialNodeOrientation(), ComputeIPCurvature(), InterpolateOrientation(), and Shell4EASANS().

doublereal Shell4EASANS::xi_n

static

Initial value:

= {
        { 1.,  1.},
        {-1.,  1.},
        {-1., -1.},
        { 1., -1.}
}

Definition at line 140 of file shelleasans.h.

Referenced by ComputeInitialNodeOrientation().

Vec3 Shell4EASANS::y_i_1[NUMIP]

protected

Definition at line 266 of file shelleasans.h.

Referenced by AssRes(), and Shell4EASANS().

Vec3 Shell4EASANS::y_i_2[NUMIP]

protected

Definition at line 267 of file shelleasans.h.

Referenced by AssRes(), and Shell4EASANS().

The documentation for this class was generated from the following files:

Public Types

Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation