#include <membraneeas.h>

Inheritance diagram for Membrane4EAS:

Collaboration diagram for Membrane4EAS:

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	NodeName { NODE1 = 0, NODE2 = 1, NODE3 = 2, NODE4 = 3, NUMNODES = 4 }

enum	Deformations { STRAIN = 0, NUMDEFORM = 1 }

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 Membrane
enum	Type { UNKNOWN = -1, ELASTIC = 0, VISCOELASTIC, LASTMEMBRANETYPE }

typedef MyVectorHandler	vh

typedef std::vector< vh >	vvh

typedef FullMatrixHandler	fmh

typedef std::vector< fmh >	vfmh

typedef ConstitutiveLawOwner < vh, fmh >	ConstitutiveLawOwnerType

Public Member Functions
	Membrane4EAS (unsigned int uL, const DofOwner pDO, const StructDispNode pN1, const StructDispNode pN2, const StructDispNode pN3, const StructDispNode *pN4, const fmh &pDTmp, const vh &PreStress, flag fOut)

virtual	~Membrane4EAS (void)

virtual Membrane::Type	GetMembraneType (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 StructDispNode *	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 Membrane
	Membrane (unsigned uLabel, const DofOwner *pDO, flag fOut)

virtual	~Membrane (void)

virtual Elem::Type	GetElemType (void) const

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_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 StructDispNode *	pNode [NUMNODES]

Vec3	xa_0 [NUMNODES]

Vec3	xa [NUMNODES]

Mat3x3	T_0_0

Mat3x3	T_0_i [NUMIP]

Mat3x3	T_0

Mat3x3	T_i [NUMIP]

Vec3	eps_tilde_0_i [NUMIP]

Vec3	eps_tilde_i [NUMIP]

fmh	S_alpha_beta_0

doublereal	alpha_0

doublereal	alpha_i [NUMIP]

vfmh	L_alpha_beta_i

vfmh	B_overline_i

vfmh	P_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	UpdateNodalAndAveragePos ()

void	ComputeInitialIptOrientation ()

Detailed Description

Definition at line 61 of file membraneeas.h.

Member Enumeration Documentation

enum Membrane4EAS::Deformations

Enumerator
STRAIN
NUMDEFORM

Definition at line 105 of file membraneeas.h.

                           {
                 STRAIN = 0,
 
                 NUMDEFORM = 1
         };

enum Membrane4EAS::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 75 of file membraneeas.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 Membrane4EAS::NodeName

Enumerator
NODE1
NODE2
NODE3
NODE4
NUMNODES

Definition at line 92 of file membraneeas.h.

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

Constructor & Destructor Documentation

Membrane4EAS::Membrane4EAS	(	unsigned int	uL,
		const DofOwner *	pDO,
		const StructDispNode *	pN1,
		const StructDispNode *	pN2,
		const StructDispNode *	pN3,
		const StructDispNode *	pN4,
		const fmh &	pDTmp,
		const vh &	PreStress,
		flag	fOut
	)

Definition at line 110 of file membraneeas.cc.

References alpha_0, alpha_i, ComputeInitialIptOrientation(), DRef, eps_tilde_0_i, Mat3x3::GetCol(), iGetNumDof(), L_alpha_beta_i, NODE1, NODE2, NODE3, NODE4, NUMIP, NUMNODES, P_i, pNode, S_alpha_beta_0, SAFENEWWITHCONSTRUCTOR, T_0, T_0_0, T_0_i, T_i, xa, xa_0, xi_0, xi_i, y_i_1, and y_i_2.

 : 
 Elem(uL, fOut), 
 Membrane(uL, pDO, fOut),
 
 S_alpha_beta_0(2, 2),               // 2x2
 
 L_alpha_beta_i(NUMIP, fmh(4, 2) ),  // 4x2
 
 B_overline_i(NUMIP, fmh(3, 12) ),   // 12x24
 
 P_i(NUMIP, fmh(3, iGetNumDof()) ),  // 12x(iGetNumDof)
 
 beta(iGetNumDof()),                 // 1x(iGetNumDof)
 epsilon_hat(3),                     // 1x12
 epsilon(3),                         // 1x12
 
 #ifndef USE_CL_IN_MEMBRANE
 bPreStress(PreStress.Norm() > 0.),
 PreStress(PreStress),
 #endif // ! USE_CL_IN_MEMBRANE
 
 DRef(NUMIP, fmh(3, 3)),            // 12x12
 stress_i(NUMIP, vh(3))              // 1x12
 {
 #ifdef USE_CL_IN_MEMBRANE
         for (integer i = 0; i < NUMIP; i++) {
                 pD[i] = 0;
                 SAFENEWWITHCONSTRUCTOR(pD[i],
                         ConstitutiveLawOwnerType,
                         ConstitutiveLawOwnerType(pDTmp[i]));    
         }
 #else // ! USE_CL_IN_MEMBRANE
         for (unsigned i = 0; i < NUMIP; i++) {
                 DRef[i] = pDTmp;
         }
 #endif // ! USE_CL_IN_MEMBRANE
 
         pNode[NODE1] = pN1;
         pNode[NODE2] = pN2;
         pNode[NODE3] = pN3;
         pNode[NODE4] = pN4;
 //      f[NODE1] = f1;
 //      f[NODE2] = f2;
 //      f[NODE3] = f3;
 //      f[NODE4] = f4;
         ComputeInitialIptOrientation();
 //      UpdateNodalAndAveragePos();
         // 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];
         }
 
 
         fmh M_0(3, 3);
         fmh M_0_Inv(3, 3);
         // doublereal pd_M_0[9], *ppd_M_0[3]; fmh M_0(pd_M_0, ppd_M_0, 9, 3, 3);
         // doublereal pd_M_0_Inv[9], *ppd_M_0_Inv[3]; fmh M_0_Inv(pd_M_0_Inv, ppd_M_0_Inv, 9, 3, 3);
         {
                 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(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(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) +
                             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);
                 
                 // M_0(5, 5) = S_alpha_beta_0(1, 1);
                 // M_0(5, 6) = S_alpha_beta_0(1, 2);
                 // M_0(6, 5) = S_alpha_beta_0(2, 1);
                 // M_0(6, 6) = S_alpha_beta_0(2, 2);
                 
                 Inv3x3(M_0, M_0_Inv);
         }
 
         for (integer i = 0; i < NUMIP; i++) {
                 fmh L_alpha_B_i(4, 2);
                 fmh S_alpha_beta_i(2, 2);
                 // doublereal pd_L_alpha_B_i[8], *ppd_L_alpha_B_i[2]; fmh L_alpha_B_i(pd_L_alpha_B_i, ppd_L_alpha_B_i, 8, 4, 2);
                 // doublereal pd_S_alpha_beta_i[4], *ppd_S_alpha_beta_i[2]; fmh S_alpha_beta_i(pd_S_alpha_beta_i, ppd_S_alpha_beta_i, 4, 2, 2);
                 Vec3 x_1 = InterpDeriv_xi1(xa, xi_i[i]);
                 Vec3 x_2 = InterpDeriv_xi2(xa, xi_i[i]);
                 S_alpha_beta_i(1, 1) = T_0_i[i].GetCol(1) * x_1;
                 S_alpha_beta_i(2, 1) = T_0_i[i].GetCol(2) * x_1;
                 S_alpha_beta_i(1, 2) = T_0_i[i].GetCol(1) * x_2;
                 S_alpha_beta_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(1, 1) * S_alpha_beta_i(2, 2) -
                         S_alpha_beta_i(1, 2) * S_alpha_beta_i(2, 1);
 
                 // xi_i_i = S_alpha_beta_i^{-1}
                 FullMatrixHandler xi_i_i(2, 2);
                 Inv2x2(S_alpha_beta_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(3, iGetNumDof());
                 // integer i_H = iGetNumDof(); doublereal pd_H[3*i_H], *ppd_H[i_H]; fmh H(pd_H, ppd_H, 3*i_H, 3, i_H);
                 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(3, 7) = xi_i[i][1];
 //              H(4, 6) = t;
                 
                 fmh Perm(3, 3), tmpP(3, iGetNumDof());    // Perm (12x6), tmpP(6xiGetNumDof)
                 // doublereal pd_Perm[9], *ppd_Perm[3]; fmh Perm(pd_Perm, ppd_Perm, 9, 3, 3);
                 // integer i_tmpP = iGetNumDof(); doublereal pd_tmpP[3*i_tmpP], *ppd_tmpP[i_tmpP]; fmh tmpP(pd_tmpP, ppd_tmpP, 3*i_tmpP, 3, i_tmpP);
 // FIXME: CHECK ORDER !!!! apparently no need for permutation
                 Perm(1, 1) = 1.;
                 Perm(2, 2) = 1.;
                 Perm(3, 3) = 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
         for (integer i = 0; i < NUMIP; i++) {
                 InterpDeriv(xa, L_alpha_beta_i[i], y_i_1[i], y_i_2[i]);
                 fmh F(3, 2);
                 // doublereal pd_F[6], *ppd_F[2]; fmh(pd_F, ppd_F, 6, 3, 2);
                 F.Put(1, 1, y_i_1[i]);
                 F.Put(1, 2, y_i_2[i]);
                 fmh FTF(2, 2);
                 // doublereal pd_FTF[4], *ppd_FTF[2]; fmh FTF(pd_FTF, ppd_FTF, 4, 2, 2);
                 F.MatTMatMul(FTF, F);
                 eps_tilde_0_i[i](1) = 0.5*( FTF(1, 1) - 1 );
                 eps_tilde_0_i[i](2) = 0.5*( FTF(2, 2) - 1 );
                 eps_tilde_0_i[i](3) = FTF(1, 2);
 /*
                 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]);
 */
         }
 }

Here is the call graph for this function:

Membrane4EAS::~Membrane4EAS ( void )

virtual

Definition at line 297 of file membraneeas.cc.

References ASSERT, NUMIP, and SAFEDELETE.

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

Member Function Documentation

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

inlineprotectedvirtual

Definition at line 182 of file membraneeas.h.

References NO_OP.

                                                      {
                 NO_OP;
         };

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

virtual

Implements Elem.

Definition at line 429 of file membraneeas.cc.

References FullSubMatrixHandler::Add(), FullSubMatrixHandler::AddT(), alpha_i, B_overline_i, DRef, DofOwnerOwner::iGetFirstIndex(), StructDispNode::iGetFirstMomentumIndex(), StructDispNode::iGetFirstPositionIndex(), iGetNumDof(), FullSubMatrixHandler::IncCoef(), L_alpha_beta_i, MatrixHandler::MatMatMul(), NUMIP, NUMNODES, P_i, pNode, FullSubMatrixHandler::PutColIndex(), FullSubMatrixHandler::PutRowIndex(), FullSubMatrixHandler::ResizeReset(), VariableSubMatrixHandler::SetFull(), stress_i, T_i, 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 iNodeFirstPosIndex = pNode[i]->iGetFirstPositionIndex();
                 integer iNodeFirstMomIndex = pNode[i]->iGetFirstMomentumIndex();
                 for (int iCnt = 1; iCnt <= 3; iCnt++) {
                         WM.PutRowIndex(iCnt + 3 * i, iNodeFirstMomIndex + iCnt);
                         WM.PutColIndex(iCnt + 3 * i, iNodeFirstPosIndex + iCnt);
                 }
         }
         for (unsigned int iCnt = 1; iCnt <= iGetNumDof(); iCnt++) {
                 WM.PutRowIndex(12 + iCnt, iFirstReactionIndex + iCnt);
                 WM.PutColIndex(12 + iCnt, iFirstReactionIndex + iCnt);
         }
 
         // tangente
         FullMatrixHandler Km(12, 12);
         FullMatrixHandler K_beta_q(iGetNumDof(), 12);
         FullMatrixHandler K_beta_beta(iGetNumDof(), iGetNumDof());
 
         FullMatrixHandler Ktm(12, 3);
         
         FullMatrixHandler Ktbetaq(iGetNumDof(), 3);
         
         FullMatrixHandler C(3, 3);
         for (integer i = 0; i < NUMIP; 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]);
                 
                 P_i[i].MatTMatMul(Ktbetaq, C);
                 Ktbetaq.MatMatMul(K_beta_q, B_overline_i[i]);
                 
                 Ktbetaq.MatMatMul(K_beta_beta, P_i[i]);
                 
                 {
                         Vec3 n1;
                         ExtractVec3(n1, stress_i[i], 1);
 
                         Vec3 Tn1 = T_i[i] * n1;
                         Vec3 Tn2 = T_i[i] * n1;
                         Vec3 Tm1 = T_i[i] * n1;
                         Vec3 Tm2 = T_i[i] * n1;
                         for (int n = 0; n < NUMNODES; n++) {
                                 for (int m = 0; m < NUMNODES; m++) {
                                 }
                         }
                 }
 
                 WM.Add(1, 1, Km, alpha_i[i] * w_i[i] * dCoef);
 
                 WM.AddT(1, 13, K_beta_q, alpha_i[i] * w_i[i]);
                 WM.Add(13, 1, K_beta_q, alpha_i[i] * w_i[i] * dCoef);
                 WM.Add(13, 13, K_beta_beta, alpha_i[i] * w_i[i]);
                 
                 for (integer l = 0; l < 3; l++) {
                         doublereal Fl_1[4], Fl_2[4];
                         for (integer n = 0; n < 4; n++) {
                                 Fl_1[n] = L_alpha_beta_i[i](n + 1, 1);
                                 Fl_2[n] = L_alpha_beta_i[i](n + 1, 2);
                         }
                         for (integer n = 0; n < 4; n++) {
                                 for (integer m = 0; m < 4; m++) {
                                         WM.IncCoef(3 * n + l + 1, 3 * m + l + 1, Fl_1[n] * stress_i[i](1) * Fl_1[m] * alpha_i[i] * w_i[i] * dCoef);
                                         WM.IncCoef(3 * n + l + 1, 3 * m + l + 1, Fl_2[n] * stress_i[i](2) * Fl_2[m] * alpha_i[i] * w_i[i] * dCoef);
                                         WM.IncCoef(3 * n + l + 1, 3 * m + l + 1, (Fl_1[n] * stress_i[i](3) * Fl_2[m] + Fl_2[n] * stress_i[i](3) * Fl_1[m])* alpha_i[i] * w_i[i] * dCoef);
                                 }
                         }
                 }
 
         }
 
         return WorkMat;
 
 }

Here is the call graph for this function:

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

virtual

Implements Elem.

Definition at line 309 of file membraneeas.cc.

References alpha_i, B_overline_i, beta, bPreStress, DRef, eps_tilde_0_i, eps_tilde_i, epsilon, epsilon_hat, DofOwnerOwner::iGetFirstIndex(), StructDispNode::iGetFirstMomentumIndex(), iGetNumDof(), L_alpha_beta_i, NUMIP, NUMNODES, P_i, pNode, PreStress, MyVectorHandler::Put(), SubVectorHandler::PutRowIndex(), VectorHandler::ResizeReset(), stress_i, UpdateNodalAndAveragePos(), w_i, WorkSpaceDim(), xa, y_i_1, and y_i_2.

 {
 
         /* 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 <= 3; iCnt++) {
                         WorkVec.PutRowIndex(iCnt + 3 * i, iNodeFirstMomIndex + iCnt);
                 }
         }
 
         integer iFirstReactionIndex = iGetFirstIndex();
         for (unsigned int iCnt = 1; iCnt <= iGetNumDof(); iCnt++) {
                 WorkVec.PutRowIndex(12 + iCnt, iFirstReactionIndex + iCnt);
         }
 
         UpdateNodalAndAveragePos();
         for (unsigned int i = 1; i <= iGetNumDof(); i++) {
                 beta(i) = XCurr(iFirstReactionIndex + i);
         }
 
 
         for (integer i = 0; i < NUMIP; i++) {
                 InterpDeriv(xa, L_alpha_beta_i[i], y_i_1[i], y_i_2[i]);
                 fmh F(3, 2);
                 F.Put(1, 1, y_i_1[i]);
                 F.Put(1, 2, y_i_2[i]);
                 fmh FTF(2, 2);
                 F.MatTMatMul(FTF, F);
                 eps_tilde_i[i](1) = 0.5*( FTF(1, 1) - 1 );
                 eps_tilde_i[i](2) = 0.5*( FTF(2, 2) - 1 );
                 eps_tilde_i[i](3) = FTF(1, 2);
                 eps_tilde_i[i] -= eps_tilde_0_i[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];
 */
 
                 // parte variabile di B_overline_i
                 for (integer n = 0; n < NUMNODES; n++) {
 
                         fmh TMP_1(3, 3), TMP_2(3, 3);
                         // TMP_1.PutT(1, 1, T_i[i] * L_alpha_beta_i[i](n + 1, 1));
                         // TMP_2.PutT(1, 1, T_i[i] * L_alpha_beta_i[i](n + 1, 2));
 
                         B_overline_i[i].PutCoef(1, 1 + 3 * n, F(1, 1) * L_alpha_beta_i[i](n + 1, 1) );
                         B_overline_i[i].PutCoef(1, 2 + 3 * n, F(2, 1) * L_alpha_beta_i[i](n + 1, 1) );
                         B_overline_i[i].PutCoef(1, 3 + 3 * n, F(3, 1) * L_alpha_beta_i[i](n + 1, 1) );
                         B_overline_i[i].PutCoef(2, 1 + 3 * n, F(1, 2) * L_alpha_beta_i[i](n + 1, 2) );
                         B_overline_i[i].PutCoef(2, 2 + 3 * n, F(2, 2) * L_alpha_beta_i[i](n + 1, 2) );
                         B_overline_i[i].PutCoef(2, 3 + 3 * n, F(3, 2) * L_alpha_beta_i[i](n + 1, 2) );
                         B_overline_i[i].PutCoef(3, 1 + 3 * n, F(1, 1) * L_alpha_beta_i[i](n + 1, 2) + F(1, 2) * L_alpha_beta_i[i](n + 1, 1) );
                         B_overline_i[i].PutCoef(3, 2 + 3 * n, F(2, 1) * L_alpha_beta_i[i](n + 1, 2) + F(2, 2) * L_alpha_beta_i[i](n + 1, 1) );
                         B_overline_i[i].PutCoef(3, 3 + 3 * n, F(3, 1) * L_alpha_beta_i[i](n + 1, 2) + F(3, 2) * L_alpha_beta_i[i](n + 1, 1) );
 
 /* 
                         // delta epsilon_tilde_1_i
                         B_overline_i[i].PutT(1, 1 + 3 * n, T_i[i] * L_alpha_beta_i[i](n + 1, 1));
                         // delta epsilon_tilde_2_i
                         // B_overline_i[i].PutT(4, 1 + 3 * n, T_i[i] * L_alpha_beta_i[i](n + 1, 2));
 */
 
                 }
 
         }
         
 
         /* Calcola le azioni interne */
         for (integer i = 0; i < NUMIP; i++) {
                 epsilon.Put(1, eps_tilde_i[i]);
                 // TODO: recupera epsilon_hat con l'ordine giusto per qua
                 P_i[i].MatVecMul(epsilon_hat, beta);
 
 
 
 
 
 // FIXME: add epsilon_hat to enhance strains
                 epsilon += epsilon_hat;
 
 
 
 #ifdef USE_CL_IN_MEMBRANE
                 pD[i]->Update(epsilon);
                 stress_i[i] = pD[i]->GetF();
 #else // ! USE_CL_IN_MEMBRANE
                 DRef[i].MatVecMul(stress_i[i], epsilon);
                 if (bPreStress) {
                         stress_i[i] += PreStress;
                 }
 #endif // ! USE_CL_IN_MEMBRANE
         }
                 
         
         //Residuo
         //forze di volume
         MyVectorHandler rd(12);
         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, 13, rbeta, -alpha_i[i] * w_i[i]);
         }
 
         return WorkVec;
 }

Here is the call graph for this function:

void Membrane4EAS::ComputeInitialIptOrientation ( void )

private

Definition at line 87 of file membraneeas.cc.

References Vec3::Cross(), Vec3::Norm(), NUMIP, T_0, T_i, UpdateNodalAndAveragePos(), xa, xi_0, and xi_i.

Referenced by Membrane4EAS().

 {
         UpdateNodalAndAveragePos();
         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);
                 T_i[i] = Mat3x3(t1, t2, t3);
         }
         Vec3 t1 = InterpDeriv_xi1(xa, xi_0);
         t1 = t1 / t1.Norm();
         Vec3 t2 = InterpDeriv_xi2(xa, xi_0);
         t2 = t2 / t2.Norm();
         Vec3 t3 = t1.Cross(t2);
         t3 = t3 / t3.Norm();
         t2 = t3.Cross(t1);
         T_0 = Mat3x3(t1, t2, t3);
 }

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from Elem.

Definition at line 327 of file membraneeas.h.

References NUMNODES, and pNode.

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

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

inlinevirtual

Reimplemented from Elem.

Definition at line 216 of file membraneeas.h.

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

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

Here is the call graph for this function:

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

inlinevirtual

Reimplemented from SimulationEntity.

Definition at line 221 of file membraneeas.h.

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

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

Here is the call graph for this function:

virtual Membrane::Type Membrane4EAS::GetMembraneType ( void ) const

inlinevirtual

Definition at line 206 of file membraneeas.h.

References Membrane::ELASTIC.

                                                          {
                 return Membrane::ELASTIC;
         };

virtual unsigned int Membrane4EAS::iGetInitialNumDof ( void ) const

inlinevirtual

Implements SubjectToInitialAssembly.

Definition at line 292 of file membraneeas.h.

                                                            {
                 return 0;
         };

virtual unsigned int Membrane4EAS::iGetNumDof ( void ) const

inlinevirtual

Reimplemented from Elem.

Definition at line 210 of file membraneeas.h.

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

                                                     { 
                 //return 14;
                 //return 13;
                 return 7;
         };

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

virtual

Implements SubjectToInitialAssembly.

Definition at line 542 of file membraneeas.cc.

References VariableSubMatrixHandler::SetNullMatrix().

 {
         WorkMat.SetNullMatrix();
         return WorkMat;
 }

Here is the call graph for this function:

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

virtual

Implements SubjectToInitialAssembly.

Definition at line 550 of file membraneeas.cc.

References VectorHandler::Resize().

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

Here is the call graph for this function:

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

inlinevirtual

Implements SubjectToInitialAssembly.

Definition at line 297 of file membraneeas.h.

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

void Membrane4EAS::Output ( OutputHandler & OH ) const

virtual

Reimplemented from ToBeOutput.

Definition at line 573 of file membraneeas.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 <= 3; r++) {
                                         out << " " << stress_i[i](r);
                                 }
                         }
                         out << std::endl;
                 }
         }
 }

Here is the call graph for this function:

const StructDispNode * Membrane4EAS::pGetNode ( unsigned int i ) const

virtual

Definition at line 558 of file membraneeas.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 & Membrane4EAS::Restart ( std::ostream & out ) const

virtual

Implements Elem.

Definition at line 527 of file membraneeas.cc.

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

virtual void Membrane4EAS::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 302 of file membraneeas.h.

References NO_OP.

                                                         {
                 NO_OP;
         };

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

virtual

Reimplemented from SimulationEntity.

Definition at line 534 of file membraneeas.cc.

References NO_OP.

 {
         NO_OP;
 }

void Membrane4EAS::UpdateNodalAndAveragePos ( void )

private

Definition at line 79 of file membraneeas.cc.

References StructDispNode::GetXCurr(), NUMNODES, pNode, and xa.

Referenced by AssRes(), and ComputeInitialIptOrientation().

 {
         for (integer i = 0; i < NUMNODES; i++) {
                 xa[i] = pNode[i]->GetXCurr();
         }
 }

Here is the call graph for this function:

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

inlinevirtual

Implements Elem.

Definition at line 241 of file membraneeas.h.

References iGetNumDof().

Referenced by AssJac(), and AssRes().

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

Here is the call graph for this function:

Member Data Documentation

doublereal Membrane4EAS::alpha_0

protected

Definition at line 146 of file membraneeas.h.

Referenced by Membrane4EAS().

doublereal Membrane4EAS::alpha_i[NUMIP]

protected

Definition at line 147 of file membraneeas.h.

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

vfmh Membrane4EAS::B_overline_i

protected

Definition at line 150 of file membraneeas.h.

Referenced by AssJac(), and AssRes().

vh Membrane4EAS::beta

protected

Definition at line 157 of file membraneeas.h.

Referenced by AssRes().

bool Membrane4EAS::bFirstRes

protected

Definition at line 177 of file membraneeas.h.

bool Membrane4EAS::bPreStress

protected

Definition at line 166 of file membraneeas.h.

Referenced by AssRes().

vfmh Membrane4EAS::DRef

protected

Definition at line 171 of file membraneeas.h.

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

Vec3 Membrane4EAS::eps_tilde_0_i[NUMIP]

protected

Definition at line 129 of file membraneeas.h.

Referenced by AssRes(), and Membrane4EAS().

Vec3 Membrane4EAS::eps_tilde_i[NUMIP]

protected

Definition at line 131 of file membraneeas.h.

Referenced by AssRes().

vh Membrane4EAS::epsilon

protected

Definition at line 159 of file membraneeas.h.

Referenced by AssRes().

vh Membrane4EAS::epsilon_hat

protected

Definition at line 158 of file membraneeas.h.

Referenced by AssRes().

vvh Membrane4EAS::FRef

protected

Definition at line 165 of file membraneeas.h.

vfmh Membrane4EAS::L_alpha_beta_i

protected

Definition at line 148 of file membraneeas.h.

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

vfmh Membrane4EAS::P_i

protected

Definition at line 152 of file membraneeas.h.

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

const StructDispNode* Membrane4EAS::pNode[NUMNODES]

protected

Definition at line 113 of file membraneeas.h.

Referenced by AssJac(), AssRes(), GetConnectedNodes(), Membrane4EAS(), pGetNode(), and UpdateNodalAndAveragePos().

vh Membrane4EAS::PreStress

protected

Definition at line 167 of file membraneeas.h.

Referenced by AssRes().

fmh Membrane4EAS::S_alpha_beta_0

protected

Definition at line 145 of file membraneeas.h.

Referenced by Membrane4EAS().

vvh Membrane4EAS::stress_i

protected

Definition at line 174 of file membraneeas.h.

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

Mat3x3 Membrane4EAS::T_0

protected

Definition at line 124 of file membraneeas.h.

Referenced by ComputeInitialIptOrientation(), and Membrane4EAS().

Mat3x3 Membrane4EAS::T_0_0

protected

Definition at line 121 of file membraneeas.h.

Referenced by Membrane4EAS().

Mat3x3 Membrane4EAS::T_0_i[NUMIP]

protected

Definition at line 122 of file membraneeas.h.

Referenced by Membrane4EAS().

Mat3x3 Membrane4EAS::T_i[NUMIP]

protected

Definition at line 125 of file membraneeas.h.

Referenced by AssJac(), ComputeInitialIptOrientation(), and Membrane4EAS().

doublereal Membrane4EAS::w_i

static

Initial value:

Definition at line 90 of file membraneeas.h.

Referenced by AssJac(), and AssRes().

Vec3 Membrane4EAS::xa[NUMNODES]

protected

Definition at line 117 of file membraneeas.h.

Referenced by AssRes(), ComputeInitialIptOrientation(), Membrane4EAS(), and UpdateNodalAndAveragePos().

Vec3 Membrane4EAS::xa_0[NUMNODES]

protected

Definition at line 116 of file membraneeas.h.

Referenced by Membrane4EAS().

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

static

Definition at line 103 of file membraneeas.h.

Referenced by ComputeInitialIptOrientation(), and Membrane4EAS().

doublereal Membrane4EAS::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 89 of file membraneeas.h.

Referenced by ComputeInitialIptOrientation(), and Membrane4EAS().

doublereal Membrane4EAS::xi_n

static

Initial value:

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

Definition at line 101 of file membraneeas.h.

Vec3 Membrane4EAS::y_i_1[NUMIP]

protected

Definition at line 154 of file membraneeas.h.

Referenced by AssRes(), and Membrane4EAS().

Vec3 Membrane4EAS::y_i_2[NUMIP]

protected

Definition at line 155 of file membraneeas.h.

Referenced by AssRes(), and Membrane4EAS().

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