5 #if !defined(KRATOS_RUNGE_KUTTA_SCHEME_H_INCLUDED )
6 #define KRATOS_RUNGE_KUTTA_SCHEME_H_INCLUDED
47 void SetTranslationalIntegrationSchemeInProperties(Properties::Pointer pProp,
bool verbose =
true)
const override;
48 void SetRotationalIntegrationSchemeInProperties(Properties::Pointer pProp,
bool verbose =
true)
const override;
50 void UpdateTranslationalVariables(
59 const double force_reduction_factor,
62 const bool Fix_vel[3])
override;
64 void CalculateNewRotationalVariablesOfSpheres(
67 const double moment_of_inertia,
70 const double moment_reduction_factor,
74 const bool Fix_Ang_vel[3])
override;
76 void CalculateNewRotationalVariablesOfRigidBodyElements(
82 const double moment_reduction_factor,
87 const bool Fix_Ang_vel[3])
override;
89 virtual void UpdateRotationalVariables(
92 const double& moment_of_inertia,
99 const bool Fix_Ang_vel[3])
override;
101 void UpdateRotationalVariables(
111 const bool Fix_Ang_vel[3])
override;
113 void UpdateAngularVelocity(
115 const double LocalTensorInv[3][3],
119 void CalculateLocalAngularAcceleration(
120 const double moment_of_inertia,
122 const double moment_reduction_factor,
125 void CalculateLocalAngularAccelerationByEulerEquations(
129 const double moment_reduction_factor,
132 void CalculateAngularVelocityRK(
134 const double& moment_of_inertia,
138 const bool Fix_Ang_vel[3])
override;
140 void CalculateAngularVelocityRK(
146 const bool Fix_Ang_vel[3])
override;
148 void QuaternionCalculateMidAngularVelocities(
150 const double LocalTensorInv[3][3],
158 virtual std::string
Info()
const override{
159 std::stringstream buffer;
160 buffer <<
"RungeKuttaScheme";
166 virtual void PrintInfo(std::ostream& rOStream)
const override{
167 rOStream <<
"RungeKuttaScheme";
172 virtual void PrintData(std::ostream& rOStream)
const override{
189 RungeKuttaScheme(RungeKuttaScheme
const& rOther) {
207 rOStream << std::endl;
PeriodicInterfaceProcess & operator=(const PeriodicInterfaceProcess &)=delete
Definition: dem_integration_scheme.h:24
MeshType::NodesContainerType NodesContainerType
Nodes container. Which is a vector set of nodes with their Id's as key.
Definition: model_part.h:128
This class defines the node.
Definition: node.h:65
Definition: runge_kutta_scheme.h:23
DEMIntegrationScheme * CloneRaw() const override
Definition: runge_kutta_scheme.h:37
RungeKuttaScheme()
Default constructor.
Definition: runge_kutta_scheme.h:32
virtual void PrintInfo(std::ostream &rOStream) const override
Print information about this object.
Definition: runge_kutta_scheme.h:166
virtual std::string Info() const override
Turn back information as a string.
Definition: runge_kutta_scheme.h:158
KRATOS_CLASS_POINTER_DEFINITION(RungeKuttaScheme)
Pointer definition of RungeKuttaScheme.
virtual ~RungeKuttaScheme()
Destructor.
Definition: runge_kutta_scheme.h:35
virtual void PrintData(std::ostream &rOStream) const override
Print object's data.
Definition: runge_kutta_scheme.h:172
ModelPart::NodesContainerType NodesArrayType
Definition: runge_kutta_scheme.h:26
DEMIntegrationScheme::Pointer CloneShared() const override
Definition: runge_kutta_scheme.h:42
dt
Definition: DEM_benchmarks.py:173
REF: G. R. Cowper, GAUSSIAN QUADRATURE FORMULAS FOR TRIANGLES.
Definition: mesh_condition.cpp:21
std::istream & operator>>(std::istream &rIStream, LinearMasterSlaveConstraint &rThis)
input stream function
std::ostream & operator<<(std::ostream &rOStream, const LinearMasterSlaveConstraint &rThis)
output stream function
Definition: linear_master_slave_constraint.h:432
vel
Definition: pure_conduction.py:76
float delta_t
Definition: rotatingcone_PureConvectionBenchmarking.py:129
integer i
Definition: TensorModule.f:17