GRAALLDamperConstitutiveLaw Class Reference

Inheritance diagram for GRAALLDamperConstitutiveLaw:

Collaboration diagram for GRAALLDamperConstitutiveLaw:

Public Member Functions
	GRAALLDamperConstitutiveLaw (const DriveCaller *pDC, const doublereal &rla, const char *const filename)
virtual	~GRAALLDamperConstitutiveLaw (void)
ConstLawType::Type	GetConstLawType (void) const
virtual ConstitutiveLaw < doublereal, doublereal > *	pCopy (void) const
virtual std::ostream &	Restart (std::ostream &out) const
virtual void	Update (const doublereal &Eps, const doublereal &EpsPrime=0.)
Public Member Functions inherited from ConstitutiveLaw< doublereal, doublereal >
	ConstitutiveLaw (void)
virtual	~ConstitutiveLaw (void)
virtual void	AfterConvergence (const doublereal &Eps, const doublereal &EpsPrime=mb_zero< doublereal >())
virtual const doublereal &	GetEpsilon (void) const
virtual const doublereal &	GetEpsilonPrime (void) const
virtual const doublereal &	GetF (void) const
virtual const doublereal &	GetFDE (void) const
virtual const doublereal &	GetFDEPrime (void) const
virtual unsigned int	iGetNumDof (void) const
virtual std::ostream &	DescribeDof (std::ostream &out, const char *prefix="", bool bInitial=false) const
virtual void	DescribeDof (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const
virtual std::ostream &	DescribeEq (std::ostream &out, const char *prefix="", bool bInitial=false) const
virtual void	DescribeEq (std::vector< std::string > &desc, bool bInitial=false, int i=-1) const
virtual DofOrder::Order	GetDofType (unsigned int i) const
Public Member Functions inherited from WithLabel
	WithLabel (unsigned int uL=0, const std::string &sN="")
virtual	~WithLabel (void)
void	PutLabel (unsigned int uL)
void	PutName (const std::string &sN)
unsigned int	GetLabel (void) const
const std::string &	GetName (void) const
Public Member Functions inherited from SimulationEntity
	SimulationEntity (void)
virtual	~SimulationEntity (void)
virtual bool	bIsValidIndex (unsigned int i) const
virtual DofOrder::Order	GetEqType (unsigned int i) const
virtual void	SetValue (DataManager *pDM, VectorHandler &X, VectorHandler &XP, SimulationEntity::Hints *h=0)
virtual Hint *	ParseHint (DataManager *pDM, const char *s) const
virtual void	BeforePredict (VectorHandler &, VectorHandler &, VectorHandler &, VectorHandler &) const
virtual void	AfterPredict (VectorHandler &X, VectorHandler &XP)
virtual void	Update (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	DerivativesUpdate (const VectorHandler &XCurr, const VectorHandler &XPrimeCurr)
virtual void	Update (const VectorHandler &XCurr, InverseDynamics::Order iOrder)
virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP)
virtual void	AfterConvergence (const VectorHandler &X, const VectorHandler &XP, const VectorHandler &XPP)
virtual unsigned int	iGetNumPrivData (void) const
virtual unsigned int	iGetPrivDataIdx (const char *s) const
virtual doublereal	dGetPrivData (unsigned int i) const
virtual std::ostream &	OutputAppend (std::ostream &out) const
virtual void	ReadInitialState (MBDynParser &HP)

Private Types
enum	{   RLA = 1, PNTR = 2, P0 = 3, A0 = 4,   RL0 = 5, PLTREX = 6, TCOST = 7, TLIN = 8,   TEMPORIF = 9, VELRIF = 10, CORONA = 11, TRIFVUOTO = 12,   T = 13, TRIFTEMP = 14, RNURIF = 15, ASTNT = 16,   R0 = 17, TIME = 18 }

Private Attributes
doublereal	v [TIME]
integer	NPDMR
integer	NPDMA
integer	NPCVR
integer	NPCVA
integer	NPVT
doublereal *	TBDMR
doublereal *	TBDMA
doublereal *	TBCVR
doublereal *	TBCVA
doublereal *	VETVIS
doublereal	dT0
DriveOwner	Time
flag	fWorking

Additional Inherited Members
Public Types inherited from ConstitutiveLaw< doublereal, doublereal >
typedef ConstitutiveLaw < doublereal, doublereal > ::ErrNotAvailable	Err
Public Types inherited from SimulationEntity
typedef std::vector< Hint * >	Hints
Protected Attributes inherited from ConstitutiveLaw< doublereal, doublereal >
doublereal	Epsilon
doublereal	EpsilonPrime
doublereal	F
doublereal	FDE
doublereal	FDEPrime
Protected Attributes inherited from WithLabel
unsigned int	uLabel
std::string	sName

Detailed Description

Definition at line 63 of file module-damper-graall.cc.

Member Enumeration Documentation

anonymous enum

private

Enumerator
RLA
PNTR
P0
A0
RL0
PLTREX
TCOST
TLIN
TEMPORIF
VELRIF
CORONA
TRIFVUOTO
T
TRIFTEMP
RNURIF
ASTNT
R0
TIME

Definition at line 66 of file module-damper-graall.cc.

             {
                RLA = 1,
                PNTR = 2,
                P0 = 3,
                A0 = 4,
                RL0 = 5,
                PLTREX = 6,
                TCOST = 7,
                TLIN = 8,
                TEMPORIF = 9,
                VELRIF = 10,
                CORONA = 11,
                TRIFVUOTO = 12,
                T = 13,
                TRIFTEMP = 14,
                RNURIF = 15,
                ASTNT = 16,
                R0 = 17,
                TIME = 18
        };

Constructor & Destructor Documentation

GRAALLDamperConstitutiveLaw::GRAALLDamperConstitutiveLaw ( const DriveCaller *  pDC, const doublereal &  rla, const char *const  filename  )

inline

Definition at line 130 of file module-damper-graall.cc.

References A0, ASTNT, CORONA, DEBUGLCOUT, DriveOwner::dGet(), MYDEBUG_INPUT, NPCVA, NPCVR, NPDMA, NPDMR, NPVT, P0, PLTREX, PNTR, R0, RL0, RLA, RNURIF, SAFENEWARR, T, TBCVA, TBCVR, TBDMA, TBDMR, TCOST, TEMPORIF, TIME, Time, TLIN, TRIFTEMP, TRIFVUOTO, v, VELRIF, and VETVIS.

        : ConstitutiveLaw<doublereal, doublereal>(),
        NPDMR(0), NPDMA(0), NPCVR(0), NPCVA(0), NPVT(0),
        TBDMR(NULL), TBDMA(NULL), TBCVR(NULL), TBCVA(NULL), VETVIS(NULL),
        dT0(0.), Time(pDC), fWorking(0)
        {
                char s[1024];

                /* Si assume che il drive che viene passato sia in realta'
                 * un drive di tipo Time, che verra' usato per conoscere
                 * il tempo della simulazione */

                /* apre il file di ingresso in formato graal */
                std::ifstream in(filename);

                /* legge alcuni parametri */
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Skipping ... (remark: \"" << s
                        << "\")" << std::endl);

                in >> v[P0-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading P0 (=" << v[P0-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[A0-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading A0 (=" << v[A0-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[PNTR-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading PNTR (=" << v[PNTR-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[RL0-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading RL0 (=" << v[RL0-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[PLTREX-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading PLTREX (=" << v[PLTREX-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Skipping ... (remark: \"" << s
                        << "\")" << std::endl);

                in >> v[ASTNT-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading ASTNT (=" << v[ASTNT-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[R0-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading R0 (=" << v[R0-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[TCOST-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading TCOST (=" << v[TCOST-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[TLIN-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading TLIN (=" << v[TLIN-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[TEMPORIF-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading A0 (=" << v[TEMPORIF-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[VELRIF-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading VELRIF (=" << v[VELRIF-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[CORONA-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading CORONA (=" << v[CORONA-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[TRIFVUOTO-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading TRIFVUOTO (=" << v[TRIFVUOTO-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[RNURIF-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading RNURIF (=" << v[RNURIF-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                in >> v[TRIFTEMP-1];
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading TRIFTEMP (=" << v[TRIFTEMP-1]
                        << ", remark: \"" << s << "\")" << std::endl);

                /* legge le tabelle */

                /* legge TBDMR */
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading TBDMR (remark: \"" << s
                        << "\")" << std::endl);

                in >> NPDMR;
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "NPDMR=" << NPDMR
                        << " (remark: \"" << s << "\")" << std::endl);

                SAFENEWARR(TBDMR, doublereal, 2*NPDMR);

                for (integer i = 0; i < NPDMR; i++) {
                        in >> TBDMR[i];
                        in >> TBDMR[NPDMR+i];
                        in.getline(s, sizeof(s));
                        DEBUGLCOUT(MYDEBUG_INPUT, "Reading TBDMR[" << i+1
                                << "]=(" << TBDMR[i]
                                << "," << TBDMR[NPDMR+i]
                                << ") (remark: \"" << s << "\")" << std::endl);
                }

                /* legge TBDMA */
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading TBDMA (remark: \"" << s
                        << "\")" << std::endl);

                in >> NPDMA;
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "NPDMA=" << NPDMA
                        << " (remark: \"" << s << "\")" << std::endl);

                SAFENEWARR(TBDMA, doublereal, 2*NPDMA);

                for (integer i = 0; i < NPDMA; i++) {
                        in >> TBDMA[i];
                        in >> TBDMA[NPDMA+i];
                        in.getline(s, sizeof(s));
                        DEBUGLCOUT(MYDEBUG_INPUT, "Reading TBDMA[" << i+1
                                << "]=(" << TBDMA[i]
                                << "," << TBDMA[NPDMA+i]
                                << ") (remark: \"" << s << "\")" << std::endl);
                }

                /* legge TBCVR */
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading TBCVR (remark: \"" << s
                        << "\")" << std::endl);

                in >> NPCVR;
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "NPCVR=" << NPCVR
                        << " (remark: \"" << s << "\")" << std::endl);

                SAFENEWARR(TBCVR, doublereal, 2*NPCVR);

                for (integer i = 0; i < NPCVR; i++) {
                        in >> TBCVR[i];
                        in >> TBCVR[NPCVR+i];
                        in.getline(s, sizeof(s));
                        DEBUGLCOUT(MYDEBUG_INPUT, "Reading TBCVR[" << i+1
                                << "]=(" << TBCVR[i]
                                << "," << TBCVR[NPCVR+i]
                                << ") (remark: \"" << s << "\")" << std::endl);
                }

                /* legge TBCVA */
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading TBCVA (remark: \"" << s
                        << "\")" << std::endl);

                in >> NPCVA;
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "NPCVA=" << NPCVA
                        << " (remark: \"" << s << "\")" << std::endl);

                SAFENEWARR(TBCVA, doublereal, 2*NPCVA);

                for (integer i = 0; i < NPCVA; i++) {
                        in >> TBCVA[i];
                        in >> TBCVA[NPCVA+i];
                        in.getline(s, sizeof(s));
                        DEBUGLCOUT(MYDEBUG_INPUT, "Reading TBCVA[" << i+1
                                << "]=(" << TBCVA[i]
                                << "," << TBCVA[NPCVA+i]
                                << ") (remark: \"" << s << "\")" << std::endl);
                }

                /* legge VETVIS */
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "Reading VETVIS (remark: \"" << s
                        << "\")" << std::endl);

                in >> NPVT;
                in.getline(s, sizeof(s));
                DEBUGLCOUT(MYDEBUG_INPUT, "NPVT=" << NPVT
                        << " (remark: \"" << s << "\")" << std::endl);

                SAFENEWARR(VETVIS, doublereal, 2*NPVT);

                for (integer i = 0; i < NPVT; i++) {
                        in >> VETVIS[i];
                        in >> VETVIS[NPVT+i];
                        in.getline(s, sizeof(s));
                        DEBUGLCOUT(MYDEBUG_INPUT, "Reading VETVIS[" << i+1
                                << "]=(" << VETVIS[i]
                                << "," << VETVIS[NPVT+i]
                                << ") (remark: \"" << s << "\")" << std::endl);
                }

                /* fine */
                in.close();

                /* Inizializza altri dati */
                v[RLA-1] = rla;   // passato come argomento, per ora
                v[T-1] = 273.16;  // T assoluta ?!?
                v[TIME-1] = Time.dGet();
        };

Here is the call graph for this function:

virtual GRAALLDamperConstitutiveLaw::~GRAALLDamperConstitutiveLaw ( void  )

inlinevirtual

Definition at line 352 of file module-damper-graall.cc.

References SAFEDELETEARR, TBCVA, TBCVR, TBDMA, TBDMR, and VETVIS.

                                                   {
                SAFEDELETEARR(VETVIS);
                SAFEDELETEARR(TBCVA);
                SAFEDELETEARR(TBCVR);
                SAFEDELETEARR(TBDMA);
                SAFEDELETEARR(TBDMR);
        };

Member Function Documentation

ConstLawType::Type GRAALLDamperConstitutiveLaw::GetConstLawType ( void  ) const

inlinevirtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 360 of file module-damper-graall.cc.

References ConstLawType::VISCOELASTIC.

                                                     {
                return ConstLawType::VISCOELASTIC;
        };

virtual ConstitutiveLaw<doublereal, doublereal>* GRAALLDamperConstitutiveLaw::pCopy ( void  ) const

inlinevirtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 364 of file module-damper-graall.cc.

                                                                           {
                return 0;
        };

virtual std::ostream& GRAALLDamperConstitutiveLaw::Restart ( std::ostream &  out ) const

inlinevirtual

Reimplemented from ConstitutiveLaw< doublereal, doublereal >.

Definition at line 368 of file module-damper-graall.cc.

                                                           {
                return out << "GRAALLDamperConstitutiveLaw not implemented yet!" << std::endl;
        };

virtual void GRAALLDamperConstitutiveLaw::Update ( const doublereal &  Eps, const doublereal &  EpsPrime = 0.  )

inlinevirtual

Implements ConstitutiveLaw< doublereal, doublereal >.

Definition at line 372 of file module-damper-graall.cc.

References DriveOwner::dGet(), dmpfr(), dT0, ConstitutiveLaw< doublereal, doublereal >::Epsilon, ConstitutiveLaw< doublereal, doublereal >::EpsilonPrime, ConstitutiveLaw< doublereal, doublereal >::F, grad::fabs(), ConstitutiveLaw< doublereal, doublereal >::FDE, ConstitutiveLaw< doublereal, doublereal >::FDEPrime, fWorking, NPCVA, NPCVR, NPDMA, NPDMR, NPVT, RLA, TBCVA, TBCVR, TBDMA, TBDMR, TIME, Time, v, and VETVIS.

                                                                                    {
                Epsilon = Eps;
                EpsilonPrime = EpsPrime;

                /* se non e' attivo, esegue un test sulla velocita' relativa;
                 * se si sta muovendo, allora setta il flag di stato su attivo,
                 * e calcola il tempo di inizio */
                if (fWorking == 0) {
                        if (fabs(EpsPrime) < 1.e-6) {
                                F = 0.;
                                FDE = 0.;
                                FDEPrime = 0.;
                                return;
                        } else {
                                fWorking = 1;
                                dT0 = Time.dGet();
                        }
                }

                v[TIME-1] = Time.dGet() - dT0;

                doublereal E = Epsilon*v[RLA-1];
                doublereal EP = EpsilonPrime*v[RLA-1];

                __FC_DECL__(dmpfr) (
                        &E,
                        &EP,
                        v,
                        TBDMR,
                        TBDMA,
                        TBCVR,
                        TBCVA,
                        VETVIS,
                        &NPDMR,
                        &NPDMA,
                        &NPCVR,
                        &NPCVA,
                        &NPVT,
                        &F,
                        &FDE,
                        &FDEPrime);
        };

Here is the call graph for this function:

Member Data Documentation

doublereal GRAALLDamperConstitutiveLaw::dT0

private

Definition at line 125 of file module-damper-graall.cc.

Referenced by Update().

flag GRAALLDamperConstitutiveLaw::fWorking

private

Definition at line 127 of file module-damper-graall.cc.

Referenced by Update().

integer GRAALLDamperConstitutiveLaw::NPCVA

private

Definition at line 118 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), and Update().

integer GRAALLDamperConstitutiveLaw::NPCVR

private

Definition at line 118 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), and Update().

integer GRAALLDamperConstitutiveLaw::NPDMA

private

Definition at line 118 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), and Update().

integer GRAALLDamperConstitutiveLaw::NPDMR

private

Definition at line 118 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), and Update().

integer GRAALLDamperConstitutiveLaw::NPVT

private

Definition at line 118 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), and Update().

doublereal* GRAALLDamperConstitutiveLaw::TBCVA

private

Definition at line 122 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), Update(), and ~GRAALLDamperConstitutiveLaw().

doublereal* GRAALLDamperConstitutiveLaw::TBCVR

private

Definition at line 121 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), Update(), and ~GRAALLDamperConstitutiveLaw().

doublereal* GRAALLDamperConstitutiveLaw::TBDMA

private

Definition at line 120 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), Update(), and ~GRAALLDamperConstitutiveLaw().

doublereal* GRAALLDamperConstitutiveLaw::TBDMR

private

Definition at line 119 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), Update(), and ~GRAALLDamperConstitutiveLaw().

DriveOwner GRAALLDamperConstitutiveLaw::Time

private

Definition at line 126 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), and Update().

doublereal GRAALLDamperConstitutiveLaw::v[TIME]

private

Definition at line 87 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), and Update().

doublereal* GRAALLDamperConstitutiveLaw::VETVIS

private

Definition at line 123 of file module-damper-graall.cc.

Referenced by GRAALLDamperConstitutiveLaw(), Update(), and ~GRAALLDamperConstitutiveLaw().

---

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

• module-damper-graall.cc