TurbulentViscoElasticCLR< T, Tder > Struct Template Reference

Inheritance diagram for TurbulentViscoElasticCLR< T, Tder >:

Collaboration diagram for TurbulentViscoElasticCLR< 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 TurbulentViscoElasticCLR< T, Tder >

Definition at line 1140 of file constltp_impl.cc.

Member Function Documentation

template<class T , class Tder >

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

inlinevirtual

Implements ConstitutiveLawRead< T, Tder >.

Definition at line 1142 of file constltp_impl.cc.

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

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

                CLType = ConstLawType::VISCOELASTIC;

                doublereal dS = HP.GetReal();
                DEBUGCOUT("Visco-Elastic Turbulent Rod Joint, stiffness = "
                                << dS << std::endl);

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

                doublereal dParabStiff = HP.GetReal();
                DEBUGCOUT("stiffness prime = " << dParabStiff << std::endl);

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

                doublereal dTreshold = 0.;
                if (HP.IsArg()) {
                        dTreshold = HP.GetReal(dTreshold);

                        /*
                         * Il legame costitutivo ha la forma seguente:
                         *      F = Kp*e + Kd*(de/dt)
                         * con Kp costante e Kd dato dalla seguente legge:
                         *      Kd = cost2                per fabs(de/dt) < Treshold
                         *      Kd = 2*cost1*fabs(de/dt)  per fabs(de/dt) > Treshold
                         * se non viene inserito il valore di treshold, lo si
                         * assume = 0. e quindi il legame e' sempre del secondo
                         * tipo. Altrimenti, se non viene inserita la seconda
                         * costante cost2, si assume che vi sia raccordo tra
                         * i due tipi di legge, ovvero cost2 = cost1*Treshold
                         * altrimenti e' possibile avere un comportamento,
                         * che in prima approssimazione e' valido
                         * per numerosi fluidi, in cui vi e' un salto tra i due
                         * tipi di legge costitutiva. */
                }

                doublereal dSP = dTreshold*dParabStiff;
                if (HP.IsArg()) {
                        dSP = HP.GetReal(dSP);
                }

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

                typedef TurbulentViscoElasticConstitutiveLaw<T, Tder> L;
                SAFENEWWITHCONSTRUCTOR(pCL,
                                L,
                                L(pTplDC, PreStress,
                                        dS, dSP, dTreshold, dParabStiff));

                return pCL;
        };

Here is the call graph for this function:

---

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

• constltp_impl.cc