membraneeas.cc File Reference

#include "mbconfig.h"
#include "constltp_impl.h"
#include "tpldrive_impl.h"
#include "membraneeas.h"
#include "mynewmem.h"
#include "shell.hc"

Include dependency graph for membraneeas.cc:

Go to the source code of this file.

Typedefs
typedef LinearElasticGenericConstitutiveLaw < Membrane::vh, Membrane::fmh >	LEGCLMembrane

Functions
Elem *	ReadMembrane4EAS (DataManager *pDM, MBDynParser &HP, const DofOwner *pDO, unsigned int uLabel)

Typedef Documentation

typedef LinearElasticGenericConstitutiveLaw<Membrane::vh, Membrane::fmh> LEGCLMembrane

Definition at line 591 of file membraneeas.cc.

Function Documentation

Elem* ReadMembrane4EAS	(	DataManager *	pDM,
		MBDynParser &	HP,
		const DofOwner *	pDO,
		unsigned int	uLabel
	)

Definition at line 594 of file membraneeas.cc.

References DataManager::fReadOutput(), WithLabel::GetLabel(), DataManager::GetLogFile(), HighParser::GetReal(), MBDYN_EXCEPT_ARGS, Elem::PLATE, ReadMembraneConstLaw(), DataManager::ReadNode(), MyVectorHandler::Reset(), SAFENEWWITHCONSTRUCTOR, and Node::STRUCTURAL.

Referenced by DataManager::ReadOneElem().

{

        const StructDispNode* pN[4];
        for (unsigned i = 0; i < 4; i++) {
                pN[i] = pDM->ReadNode<const StructDispNode, Node::STRUCTURAL>(HP);
        }

        /* Prestress and prestrain */
        Membrane::vh PreStress(3);
        PreStress.Reset();

#ifdef USE_CL_IN_MEMBRANE
        const ConstitutiveLaw<Membrane::vh, Membrane::fmh> *pD[4];

        TplDriveCaller<Membrane::vh>* pTplDC = new ZeroTplDriveCaller<Membrane::vh>;

        Membrane::fmh S(3, 3);
        for (unsigned ir = 1; ir <= 3; ir++) {
                for (unsigned ic = 1; ic <= 3; ic++) {
                        S(ir, ic) = HP.GetReal();
                }
        }

        pD[0] = 0;
        SAFENEWWITHCONSTRUCTOR(pD[0], LEGCLMembrane, LEGCLMembrane(pTplDC, PreStress, S));

        for (unsigned i = 1; i < 4; i++) {
                pD[i] = pD[0]->Copy();
        }
#else // ! USE_CL_IN_MEMBRANE
        Membrane::fmh pD(3, 3);

        if (ReadMembraneConstLaw(HP, pD, PreStress)) {
                silent_cerr("Membrane(" << uLabel << "): unable to read constitutive law" << std::endl);
                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
        }
#endif // ! USE_CL_IN_MEMBRANE

        flag fOut = pDM->fReadOutput(HP, Elem::PLATE);

        Elem *pEl = 0;
        SAFENEWWITHCONSTRUCTOR(pEl, Membrane4EAS,
                Membrane4EAS(uLabel, pDO,
                        pN[0], pN[1], pN[2], pN[3],
#ifdef USE_CL_IN_MEMBRANE
#else // ! USE_CL_IN_MEMBRANE
                        pD, PreStress,
#endif // ! USE_CL_IN_MEMBRANE
                        fOut));

        std::ostream& out = pDM->GetLogFile();
        out << "membrane4: " << uLabel
                << " " << pN[0]->GetLabel()
                << " " << pN[1]->GetLabel()
                << " " << pN[2]->GetLabel()
                << " " << pN[3]->GetLabel()
                << std::endl;

        return pEl;
}

Here is the call graph for this function: