shelleas.cc File Reference

#include "mbconfig.h"
#include "constltp_impl.h"
#include "tpldrive_impl.h"
#include "shelleas.h"
#include "mynewmem.h"
#include "shell.hc"
#include <iostream>

Include dependency graph for shelleas.cc:

Go to the source code of this file.

Typedefs
typedef LinearElasticGenericConstitutiveLaw < Shell::vh, Shell::fmh >	LEGCLShell

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

Typedef Documentation

typedef LinearElasticGenericConstitutiveLaw<Shell::vh, Shell::fmh> LEGCLShell

Definition at line 1269 of file shelleas.cc.

Function Documentation

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

Definition at line 1272 of file shelleas.cc.

References Eye3, DataManager::fReadOutput(), WithLabel::GetLabel(), DataManager::GetLogFile(), HighParser::GetReal(), MBDynParser::GetRotRel(), HighParser::IsKeyWord(), MBDYN_EXCEPT_ARGS, Elem::PLATE, R, DataManager::ReadNode(), ReadShellConstLaw(), MyVectorHandler::Reset(), SAFENEWWITHCONSTRUCTOR, and Node::STRUCTURAL.

Referenced by DataManager::ReadOneElem().

{
        silent_cout("Shell4EAS(" << uLabel << "): warning, \"shell4eas\" is deprecated; use \"shell4easans\" instead" << std::endl);

        const StructNode* pN[4];
        Mat3x3 R[4];
        for (unsigned i = 0; i < 4; i++) {
                pN[i] = pDM->ReadNode<const StructNode, Node::STRUCTURAL>(HP);
                ReferenceFrame RF(pN[i]);
                if (HP.IsKeyWord("orientation")) {
                        R[i] = HP.GetRotRel(RF);
                } else {
                        R[i] = Eye3;
                }
        }

        /* Prestress and prestrain */
        Shell::vh PreStress(12);
        PreStress.Reset();

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

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

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

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

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

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

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

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

        std::ostream& out = pDM->GetLogFile();
        out << "shell4: " << 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: