ContContact3DCLR Struct Reference

Inheritance diagram for ContContact3DCLR:

Collaboration diagram for ContContact3DCLR:

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

Detailed Description

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

Member Function Documentation

virtual ConstitutiveLaw<Vec3, Mat3x3>* ContContact3DCLR::Read ( const DataManager *  pDM, MBDynParser &  HP, ConstLawType::Type &  CLType  )

inlinevirtual

Implements ConstitutiveLawRead< Vec3, Mat3x3 >.

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

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

                                                                                {
                ConstitutiveLaw<Vec3, Mat3x3>* pCL = 0;

                CLType = ConstLawType::VISCOELASTIC;

                if (HP.IsKeyWord("help")) {
                        silent_cerr("ContContact3DCL:\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("ContContact3DCLR: invalid sign " << d
                                                << " at line " << HP.GetLineData() << std::endl);
                                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                }
                                dSign = copysign(1., d);
                        }
                }
                
                ContContact3DCL::Type type(ContContact3DCL::CC_FLORES_ET_AL);
                if (HP.IsKeyWord("formulation")) {
                        if (HP.IsKeyWord("flores")) {
                                type = ContContact3DCL::CC_FLORES_ET_AL;

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

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

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

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

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

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

                if (!HP.IsKeyWord("exp")) {
                        silent_cerr("ContContact3DCLR: \"exp\" expected at line " << HP.GetLineData() << std::endl);
                        throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                }
                doublereal dExp = HP.GetReal();
                if (dExp < 1.) {
                        silent_cerr("ContContact3DCLR: 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("ContContact3DCLR: invalid \"EpsPrimeTol\" at line " << HP.GetLineData() << std::endl);
                                throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                        }
                }

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

                SAFENEWWITHCONSTRUCTOR(pCL, ContContact3DCL,
                        ContContact3DCL(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