LinearViscoElasticCLR< T, Tder > Struct Template Reference

Inheritance diagram for LinearViscoElasticCLR< T, Tder >:

Collaboration diagram for LinearViscoElasticCLR< T, Tder >:

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

Detailed Description

template<class T, class Tder>
struct LinearViscoElasticCLR< T, Tder >

Definition at line 818 of file constltp_impl.cc.

Member Function Documentation

template<class T , class Tder >

virtual ConstitutiveLaw<T, Tder>* LinearViscoElasticCLR< T, Tder >::Read ( const DataManager *  pDM, MBDynParser &  HP, ConstLawType::Type &  CLType  )

inlinevirtual

Implements ConstitutiveLawRead< T, Tder >.

Definition at line 820 of file constltp_impl.cc.

References DEBUGCOUT, dS, IncludeParser::GetLineData(), GetPreStress(), HighParser::GetReal(), HighParser::IsKeyWord(), SAFENEWWITHCONSTRUCTOR, and ConstLawType::VISCOELASTIC.

                                                                                {
                ConstitutiveLaw<T, Tder>* pCL = 0;

                CLType = ConstLawType::VISCOELASTIC;

                doublereal dS = HP.GetReal();
                DEBUGCOUT("Stiffness = " << dS << std::endl);

                if (dS <= 0.) {
                        silent_cerr("warning, null or negative stiffness at line "
                                << HP.GetLineData() << std::endl);
                }

                doublereal dSP = 0.;
                if (HP.IsKeyWord("proportional")) {
                        doublereal k = HP.GetReal();
                        dSP = k*dS;
                } else {
                        dSP = HP.GetReal();
                }
                DEBUGCOUT("stiffness prime = " << dSP << std::endl);

                if (dSP <= 0.) {
                        silent_cerr("warning, null or negative stiffness prime at line "
                                << HP.GetLineData() << std::endl);
                }

                /* Prestress and prestrain */
                T PreStress(mb_zero<T>());
                GetPreStress(HP, PreStress);
                TplDriveCaller<T>* pTplDC = GetPreStrain<T>(pDM, HP);

#if 0   // TODO: implement a "null" constitutive law that does nothing
                if (dS == 0. && dSP == 0.) {

                } else
#endif
                if (dS == 0.) {
                        silent_cerr("warning, null stiffness, "
                                "using linear viscous constitutive law instead"
                                << std::endl);

                        typedef LinearViscousIsotropicConstitutiveLaw<T, Tder> L;
                        SAFENEWWITHCONSTRUCTOR(pCL, L, L(PreStress, dSP));
                        if (pTplDC) {
                                delete pTplDC;
                        }

                } else if (dSP == 0.) {
                        silent_cerr("warning, null stiffness prime, "
                                "using linear elastic constitutive law instead"
                                << std::endl);

                        typedef LinearElasticIsotropicConstitutiveLaw<T, Tder> L;
                        SAFENEWWITHCONSTRUCTOR(pCL, L, L(pTplDC, PreStress, dS));

                } else {
                        typedef LinearViscoElasticIsotropicConstitutiveLaw<T, Tder> L;
                        SAFENEWWITHCONSTRUCTOR(pCL, L, L(pTplDC, PreStress, dS, dSP));
                }

                return pCL;
        };

Here is the call graph for this function:

---

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

• constltp_impl.cc