HuntCrossleyCLR Struct Reference

Inheritance diagram for HuntCrossleyCLR:

Collaboration diagram for HuntCrossleyCLR:

Public Member Functions
virtual ConstitutiveLaw < doublereal, doublereal > *	Read (const DataManager *pDM, MBDynParser &HP, ConstLawType::Type &CLType)
Public Member Functions inherited from ConstitutiveLawRead< doublereal, doublereal >
virtual	~ConstitutiveLawRead (void)

Detailed Description

Definition at line 154 of file module-hunt-crossley.cc.

Member Function Documentation

virtual ConstitutiveLaw<doublereal, doublereal>* HuntCrossleyCLR::Read ( const DataManager *  pDM, MBDynParser &  HP, ConstLawType::Type &  CLType  )

inlinevirtual

Implements ConstitutiveLawRead< doublereal, doublereal >.

Definition at line 156 of file module-hunt-crossley.cc.

References copysign(), IncludeParser::GetLineData(), HighParser::GetReal(), HighParser::IsArg(), HighParser::IsKeyWord(), MBDYN_EXCEPT_ARGS, SAFENEWWITHCONSTRUCTOR, and ConstLawType::VISCOELASTIC.

                                                                                {
                ConstitutiveLaw<doublereal, doublereal>* pCL = 0;

                CLType = ConstLawType::VISCOELASTIC;

                if (HP.IsKeyWord("help")) {
                        silent_cerr("HuntCrossleyCL:\n"
                                "        hunt crossley,\n"
                                "                [ , sign, { negative | positive | <sign> } , ]\n"
                                "                alpha, <alpha>,\n"
                                "                kappa, <kappa>,\n"
                                "                exp, <exp>,\n"
                                "                [ , prestress, <prestress> ]\n"
                                "                [ , prestrain, (DriveCaller)<prestrain> ]\n"
                                << std::endl);

                        if (!HP.IsArg()) {
                                throw NoErr(MBDYN_EXCEPT_ARGS);
                        }
                }

                doublereal dSign = -1.;
                if (HP.IsKeyWord("sign")) {
                        if (HP.IsKeyWord("positive")) {
                                dSign = 1.;
                        } else if (HP.IsKeyWord("negative")) {
                                dSign = -1.;
                        } else {
                                doublereal d = HP.GetReal();
                                if (d == 0.) {
                                        silent_cerr("HuntCrossleyCLR: invalid sign " << d
                                                << " at line " << HP.GetLineData() << std::endl);
                                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                }
                                dSign = copysign(1., d);
                        }
                }

                if (!HP.IsKeyWord("alpha")) {
                        silent_cerr("HuntCrossleyCLR: \"alpha\" expected at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                doublereal dAlpha = HP.GetReal();
                if (dAlpha < 0.) {
                        silent_cerr("HuntCrossleyCLR: invalid \"alpha\" at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

                if (!HP.IsKeyWord("kappa")) {
                        silent_cerr("HuntCrossleyCLR: \"kappa\" expected at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                doublereal dK = HP.GetReal();
                if (dK <= 0.) {
                        silent_cerr("HuntCrossleyCLR: invalid \"kappa\" at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

                if (!HP.IsKeyWord("exp")) {
                        silent_cerr("HuntCrossleyCLR: \"exp\" expected at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                doublereal dExp = HP.GetReal();
                if (dExp < 1.) {
                        silent_cerr("HuntCrossleyCLR: invalid \"exp\" at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

                /* Prestress and prestrain */
                // doublereal PreStress(0.);
                // GetPreStress(HP, PreStress);
                TplDriveCaller<doublereal> *pTplDC = GetPreStrain<doublereal>(pDM, HP);

                SAFENEWWITHCONSTRUCTOR(pCL, HuntCrossleyCL,
                        HuntCrossleyCL(pTplDC, dSign, dAlpha, dK, dExp));

                return pCL;
        };

Here is the call graph for this function:

---

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

• module-hunt-crossley.cc