ContContactCLR Struct Reference

Inheritance diagram for ContContactCLR:

Collaboration diagram for ContContactCLR:

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 239 of file module-cont-contact.cc.

Member Function Documentation

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

inlinevirtual

Implements ConstitutiveLawRead< doublereal, doublereal >.

Definition at line 241 of file module-cont-contact.cc.

References ContContactCL::CC_FLORES_ET_AL, ContContactCL::CC_HUNT_CROSSLEY, ContContactCL::CC_LANKARANI_NIKRAVESH, 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("ContContactCLR:\n"
                                "        continuous contact ,\n"
                                "                [ , sign , { negative | positive | <sign> } , ]\n"
                                "                [ , formulation, { flores | lankarani nikravesh | hunt crossley } , ]\n"
                                "                restitution , <restitution_coefficient> , (0 -> 1)\n"
                                "                kappa , <stiffness> , (> 0)\n"
                                "                exp , <exponent> , (>= 1)\n"
                                "                [ , EpsPrimeTol , <EpsPrimeTol> ]\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("ContContactCLR: invalid sign " << d
                                                << " at line " << HP.GetLineData() << std::endl);
                                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                }
                                dSign = copysign(1., d);
                        }
                }
                
                ContContactCL::Type type(ContContactCL::CC_FLORES_ET_AL);
                if (HP.IsKeyWord("formulation")) {
                        if (HP.IsKeyWord("flores")) {
                                type = ContContactCL::CC_FLORES_ET_AL;

                        } else if (HP.IsKeyWord("hunt" "crossley")) {
                                type = ContContactCL::CC_HUNT_CROSSLEY;

                        } else if (HP.IsKeyWord("lankarani" "nikravesh")) {
                                type = ContContactCL::CC_LANKARANI_NIKRAVESH;

                        } else {
                                silent_cerr("ContContactCLR: invalid formulation"
                                        << " at line " << HP.GetLineData() << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }
                }       

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

                } else if ((type == ContContactCL::CC_FLORES_ET_AL) && (dRest <= 0.)) {
                        silent_cerr("ContContactCLR: null \"restitution\" incompatible with \"Flores\" at line "
                                << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }

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

                if (!HP.IsKeyWord("exp")) {
                        silent_cerr("ContContactCLR: \"exp\" expected at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                doublereal dExp = HP.GetReal();
                if (dExp < 1.) {
                        silent_cerr("ContContactCLR: invalid \"exp\" at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                
                doublereal dInitialEpsPrimeTol = 1.e-6;
                if (HP.IsKeyWord("EpsPrimeTol")) {
                        dInitialEpsPrimeTol = HP.GetReal();
                        if (dInitialEpsPrimeTol < 0.) {
                                silent_cerr("ContContactCLR: invalid \"EpsPrimeTol\" at line " << HP.GetLineData() << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }
                }

                /* Prestrain */
                TplDriveCaller<doublereal> *pTplDC = GetPreStrain<doublereal>(pDM, HP);

                SAFENEWWITHCONSTRUCTOR(pCL, ContContactCL,
                        ContContactCL(pTplDC, dSign, type, dRest, dK, dExp, dInitialEpsPrimeTol));

                return pCL;
        };

Here is the call graph for this function:

---

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

• module-cont-contact.cc