Description

General purpose section of an Euler-Bernoulli beam in 3D, not assuming homogeneous density or homogeneous elasticity, given basic material properties.

This is the case where one uses a FEA preprocessor to compute the rigidity of a complex beam made with multi-layered reinforcements with different elasticity and different density - in such a case you could not use ChBeamSectionEulerAdvanced because you do not have a single E or single density, but you rather have collective values of bending rigidities, and collective mass per unit length. This class allows using these values directly, bypassing any knowledge of area, density, Izz Iyy, E young modulus, etc. To be used with ChElementBeamEuler. This material can be shared between multiple beams.

#include <ChBeamSectionEuler.h>

Inheritance diagram for chrono::fea::ChBeamSectionEulerGeneric:

[legend]

Collaboration diagram for chrono::fea::ChBeamSectionEulerGeneric:

[legend]

Public Member Functions
	ChBeamSectionEulerGeneric (const double mAx, const double mTxx, const double mByy, const double mBzz, const double malpha, const double mCy, const double mCz, const double mSy, const double mSz, const double mmu, const double mJxx, const double mMy=0, const double mMz=0)

virtual void	SetAxialRigidity (const double mv)
	Sets the axial rigidity, usually A*E for uniform elasticity, but for nonuniform elasticity here you can put a value ad-hoc from a preprocessor.

virtual void	SetXtorsionRigidity (const double mv)
	Gets the torsion rigidity, for torsion about X axis, at elastic center, usually J*G for uniform elasticity, but for nonuniform elasticity here you can put a value ad-hoc from a preprocessor

virtual void	SetYbendingRigidity (const double mv)
	Gets the bending rigidity, for bending about Y axis, at elastic center, usually Iyy*E for uniform elasticity, but for nonuniform elasticity here you can put a value ad-hoc from a preprocessor.

virtual void	SetZbendingRigidity (const double mv)
	Gets the bending rigidity, for bending about Z axis, at elastic center, usually Izz*E for uniform elasticity, but for nonuniform elasticity here you can put a value ad-hoc from a preprocessor.

virtual void	SetSectionRotation (const double mv)
	Set the rotation in [rad], abour elastic center, of the Y Z axes for which the YbendingRigidity and ZbendingRigidity values are defined.

virtual void	SetCentroidY (const double mv)
	Gets the Y position of the elastic center respect to centerline.

virtual void	SetCentroidZ (const double mv)
	Gets the Z position of the elastic center respect to centerline.

virtual void	SetShearCenterY (const double mv)
	Gets the Y position of the shear center respect to centerline.

virtual void	SetShearCenterZ (const double mv)
	Gets the Z position of the shear center respect to centerline.

virtual void	SetMassPerUnitLength (const double mv)
	Get mass per unit length, ex.SI units [kg/m] For uniform density it would be A*density, but for nonuniform density here you can put a value ad-hoc from a preprocessor.

virtual void	SetInertiaJxxPerUnitLength (const double mv)
	Get the Jxx component of the inertia per unit length (polar inertia), computed at centerline. More...

virtual void	SetInertiaJxxPerUnitLengthInMassReference (const double mv)
	Set inertia moment per unit length Jxx_massref, as assumed computed in the "mass reference" frame, ie. More...

virtual double	GetInertiaJxxPerUnitLengthInMassReference ()
	Get inertia moment per unit length Jxx_massref, as assumed computed in the "mass reference" frame, ie. More...

void	SetCenterOfMass (double my, double mz)
	"mass reference": set the displacement of the center of mass respect to the section centerline reference.

double	GetCenterOfMassY ()

double	GetCenterOfMassZ ()

virtual double	GetAxialRigidity () const
	Gets the axial rigidity, usually A*E, but might be ad hoc.

virtual double	GetXtorsionRigidity () const
	Gets the torsion rigidity, for torsion about X axis at elastic center, usually J*G, but might be ad hoc.

virtual double	GetYbendingRigidity () const
	Gets the bending rigidity, for bending about Y axis at elastic center, usually Iyy*E, but might be ad hoc.

virtual double	GetZbendingRigidity () const
	Gets the bending rigidity, for bending about Z axis at elastic center, usually Izz*E, but might be ad hoc.

virtual double	GetSectionRotation () const
	Set the rotation of the Y Z section axes for which the YbendingRigidity and ZbendingRigidity are defined.

virtual double	GetCentroidY () const
	Gets the Y position of the elastic center respect to centerline.

virtual double	GetCentroidZ () const
	Gets the Z position of the elastic center respect to centerline.

virtual double	GetShearCenterY () const
	Gets the Y position of the shear center respect to centerline.

virtual double	GetShearCenterZ () const
	Gets the Z position of the shear center respect to centerline.

virtual void	ComputeInertiaMatrix (ChMatrixNM< double, 6, 6 > &M) override
	Compute the 6x6 sectional inertia matrix, as in {x_momentum,w_momentum}=[Mm]{xvel,wvel}.

virtual void	ComputeQuadraticTerms (ChVector<> &mF, ChVector<> &mT, const ChVector<> &mW) override
	Compute the centrifugal term and gyroscopic term. More...

virtual double	GetMassPerUnitLength () const
	Get mass per unit length, ex.SI units [kg/m].

virtual double	GetInertiaJxxPerUnitLength () const
	Get the Jxx component of the inertia per unit length (polar inertia), at centerline.

Public Member Functions inherited from chrono::fea::ChBeamSectionEuler
void	SetBeamRaleyghDamping (double mr)
	Set the Rayleigh damping ratio r (as in: R = r * K ), to do: also mass-proportional term.

double	GetBeamRaleyghDamping ()

Public Member Functions inherited from chrono::fea::ChBeamSection
void	SetDrawShape (std::shared_ptr< ChBeamSectionShape > mshape)
	Set the graphical representation for this section. More...

std::shared_ptr< ChBeamSectionShape >	GetDrawShape () const
	Get the drawing shape of this section (i.e.a 2D profile used for drawing 3D tesselation and visualization) By default a thin square section, use SetDrawShape() to change it.

void	SetDrawThickness (double thickness_y, double thickness_z)
	Shortcut: adds a ChBeamSectionShapeRectangular for visualization as a centered rectangular beam, and sets its width/height. More...

void	SetDrawCircularRadius (double draw_rad)
	Shortcut: adds a ChBeamSectionShapeCircular for visualization as a centered circular beam, and sets its radius. More...

void	SetCircular (bool ic)
	OBSOLETE only for backward compability More...

Constructor & Destructor Documentation

◆ ChBeamSectionEulerGeneric()

chrono::fea::ChBeamSectionEulerGeneric::ChBeamSectionEulerGeneric	(	const double	mAx,
		const double	mTxx,
		const double	mByy,
		const double	mBzz,
		const double	malpha,
		const double	mCy,
		const double	mCz,
		const double	mSy,
		const double	mSz,
		const double	mmu,
		const double	mJxx,
		const double	mMy = `0`,
		const double	mMz = `0`
	)

inline

Parameters

mAx	axial rigidity
mTxx	torsion rigidity
mByy	bending regidity about yy
mBzz	bending rigidity about zz
malpha	section rotation about elastic center [rad]
mCy	elastic center y displacement respect to centerline
mCz	elastic center z displacement respect to centerline
mSy	shear center y displacement respect to centerline
mSz	shear center z displacement respect to centerline
mmu	mass per unit length
mJxx	polar inertia Jxx per unit lenght, measured respect to centerline
mMy	mass center y displacement respect to centerline
mMz	mass center z displacement respect to centerline

Member Function Documentation

◆ ComputeQuadraticTerms()

void chrono::fea::ChBeamSectionEulerGeneric::ComputeQuadraticTerms	(	ChVector<> &	mF,
		ChVector<> &	mT,
		const ChVector<> &	mW
	)

overridevirtual

Compute the centrifugal term and gyroscopic term.

Parameters

mF	centrifugal term (if any) returned here
mT	gyroscopic term returned here
mW	current angular velocity of section, in material frame

Implements chrono::fea::ChBeamSectionEuler.

◆ GetInertiaJxxPerUnitLengthInMassReference()

virtual double chrono::fea::ChBeamSectionEulerGeneric::GetInertiaJxxPerUnitLengthInMassReference ( )

inlinevirtual

Get inertia moment per unit length Jxx_massref, as assumed computed in the "mass reference" frame, ie.

centered at the center of mass

◆ SetInertiaJxxPerUnitLength()

virtual void chrono::fea::ChBeamSectionEulerGeneric::SetInertiaJxxPerUnitLength ( const double mv )

inlinevirtual

Get the Jxx component of the inertia per unit length (polar inertia), computed at centerline.

For uniform density it would be Ixx*density or, by polar theorem, (Izz+Iyy)*density, but for nonuniform density here you can put a value ad-hoc from a preprocessor

◆ SetInertiaJxxPerUnitLengthInMassReference()

virtual void chrono::fea::ChBeamSectionEulerGeneric::SetInertiaJxxPerUnitLengthInMassReference ( const double mv )

inlinevirtual

Set inertia moment per unit length Jxx_massref, as assumed computed in the "mass reference" frame, ie.

centered at the center of mass. Call this after you set SetCenterOfMass() and SetMassPerUnitLength()

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

/builds/uwsbel/chrono/src/chrono/fea/ChBeamSectionEuler.h
/builds/uwsbel/chrono/src/chrono/fea/ChBeamSectionEuler.cpp

Description

Public Member Functions

Constructor & Destructor Documentation

◆ ChBeamSectionEulerGeneric()

Member Function Documentation

◆ ComputeQuadraticTerms()

◆ GetInertiaJxxPerUnitLengthInMassReference()

◆ SetInertiaJxxPerUnitLength()

◆ SetInertiaJxxPerUnitLengthInMassReference()