#include <aerodata.h>

Inheritance diagram for AeroData:

Collaboration diagram for AeroData:

Public Types
enum	UnsteadyModel { STEADY = 0, HARRIS = 1, BIELAWA = 2, LAST }

enum	{ VX = 0, VY = 1, VZ = 2, WX = 3, WY = 4, WZ = 5, FX = 0, FY = 1, FZ = 2, MX = 3, MY = 4, MZ = 5 }

Public Member Functions
	AeroData (int i_p, int i_dim, UnsteadyModel u=STEADY, DriveCaller *pt=0)

virtual	~AeroData (void)

virtual std::ostream &	Restart (std::ostream &out) const =0

std::ostream &	RestartUnsteady (std::ostream &out) const

virtual void	SetAirData (const doublereal &rho, const doublereal &c)

virtual void	SetSectionData (const doublereal &abscissa, const doublereal &chord, const doublereal &forcepoint, const doublereal &velocitypoint, const doublereal &twist, const doublereal &omega=0.)

virtual int	GetForces (int i, const doublereal W, doublereal TNG, outa_t &OUTA)

virtual int	GetForcesJac (int i, const doublereal W, doublereal TNG, Mat6x6 &J, outa_t &OUTA)

virtual unsigned int	iGetNumDof (void) const

virtual DofOrder::Order	GetDofType (unsigned int i) const

virtual void	AssRes (SubVectorHandler &WorkVec, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr, integer iFirstIndex, integer iFirstSubIndex, int i, const doublereal W, doublereal TNG, outa_t &OUTA)

virtual void	AssJac (FullSubMatrixHandler &WorkMat, doublereal dCoef, const VectorHandler &XCurr, const VectorHandler &XPrimeCurr, integer iFirstIndex, integer iFirstSubIndex, const Mat3xN &vx, const Mat3xN &wx, Mat3xN &fq, Mat3xN &cq, int i, const doublereal W, doublereal TNG, Mat6x6 &J, outa_t &OUTA)

virtual void	AfterConvergence (int i, const VectorHandler &X, const VectorHandler &XP)

AeroData::UnsteadyModel	Unsteady (void) const

Public Member Functions inherited from AeroMemory
	AeroMemory (DriveCaller *pt)

virtual	~AeroMemory (void)

void	Predict (int i, doublereal alpha, doublereal &alf1, doublereal &alf2)

void	Update (int i)

void	SetNumPoints (int i)

int	GetNumPoints (void) const

Protected Member Functions
int	StorageSize (void) const

int	GetForcesJacForwardDiff_int (int i, const doublereal W, doublereal TNG, Mat6x6 &J, outa_t &OUTA)

int	GetForcesJacCenteredDiff_int (int i, const doublereal W, doublereal TNG, Mat6x6 &J, outa_t &OUTA)

Protected Attributes
UnsteadyModel	unsteadyflag

vam_t	VAM

doublereal	Omega

Protected Attributes inherited from AeroMemory
DriveCaller *	pTime

Detailed Description

Definition at line 87 of file aerodata.h.

Member Enumeration Documentation

anonymous enum

Enumerator
VX
VY
VZ
WX
WY
WZ
FX
FY
FZ
MX
MY
MZ

Definition at line 97 of file aerodata.h.

              {
                 VX      = 0,
                 VY      = 1,
                 VZ      = 2,
 
                 WX      = 3,
                 WY      = 4,
                 WZ      = 5,
 
                 FX      = 0,
                 FY      = 1,
                 FZ      = 2,
 
                 MX      = 3,
                 MY      = 4,
                 MZ      = 5
         };

enum AeroData::UnsteadyModel

Enumerator
STEADY
HARRIS
BIELAWA
LAST

Definition at line 89 of file aerodata.h.

                            {
                 STEADY = 0,
                 HARRIS = 1,
                 BIELAWA = 2,
 
                 LAST
         };

Constructor & Destructor Documentation

AeroData::AeroData	(	int	i_p,
		int	i_dim,
		AeroData::UnsteadyModel	u = `STEADY`,
		DriveCaller *	pt = `0`
	)

Definition at line 181 of file aerodata.cc.

References vam_t::bc_position, vam_t::chord, vam_t::density, vam_t::force_position, MBDYN_EXCEPT_ARGS, AeroMemory::SetNumPoints(), vam_t::sound_celerity, STEADY, vam_t::twist, and VAM.

 : AeroMemory(ptime), unsteadyflag(u), Omega(0.)
 {
         // silence static analyzers
         VAM.density = -1.;
         VAM.sound_celerity = -1.;
         VAM.chord = -1.;
         VAM.force_position = 0.;
         VAM.bc_position = 0.;
         VAM.twist = 0.;
 
         if (u != AeroData::STEADY) {
                 if (ptime == 0) {
                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                 }
         }
 
         SetNumPoints(i_p*i_dim);
 }

Here is the call graph for this function:

AeroData::~AeroData ( void )

virtual

Definition at line 202 of file aerodata.cc.

References NO_OP.

 {
         NO_OP;
 }

Member Function Documentation

void AeroData::AfterConvergence	(	int	i,
		const VectorHandler &	X,
		const VectorHandler &	XP
	)

virtual

Reimplemented in TheodorsenAeroData.

Definition at line 416 of file aerodata.cc.

References NO_OP.

 {
         NO_OP;
 }

void AeroData::AssJac	(	FullSubMatrixHandler &	WorkMat,
		doublereal	dCoef,
		const VectorHandler &	XCurr,
		const VectorHandler &	XPrimeCurr,
		integer	iFirstIndex,
		integer	iFirstSubIndex,
		const Mat3xN &	vx,
		const Mat3xN &	wx,
		Mat3xN &	fq,
		Mat3xN &	cq,
		int	i,
		const doublereal *	W,
		doublereal *	TNG,
		Mat6x6 &	J,
		outa_t &	OUTA
	)

virtual

Reimplemented in TheodorsenAeroData.

Definition at line 403 of file aerodata.cc.

References MBDYN_EXCEPT_ARGS.

Referenced by AerodynamicBody::AssJac(), AerodynamicBeam::AssJac(), and AerodynamicBeam2::AssJac().

 {
         silent_cerr("AeroData: AssJac() is undefined" << std::endl);
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 }

void AeroData::AssRes	(	SubVectorHandler &	WorkVec,
		doublereal	dCoef,
		const VectorHandler &	XCurr,
		const VectorHandler &	XPrimeCurr,
		integer	iFirstIndex,
		integer	iFirstSubIndex,
		int	i,
		const doublereal *	W,
		doublereal *	TNG,
		outa_t &	OUTA
	)

virtual

Reimplemented in TheodorsenAeroData.

Definition at line 391 of file aerodata.cc.

References MBDYN_EXCEPT_ARGS.

Referenced by AerodynamicBody::AssVec(), AerodynamicBeam::AssVec(), and AerodynamicBeam2::AssVec().

 {
         silent_cerr("AeroData: AssRes() is undefined" << std::endl);
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 }

DofOrder::Order AeroData::GetDofType ( unsigned int i ) const

virtual

Reimplemented in TheodorsenAeroData.

Definition at line 368 of file aerodata.cc.

References MBDYN_EXCEPT_ARGS.

 {
         silent_cerr("AeroData: "
                 "GetDofType() is undefined because aerodynamic data "
                 "has no degrees of freedom" << std::endl);
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 }

int AeroData::GetForces	(	int	i,
		const doublereal *	W,
		doublereal *	TNG,
		outa_t &	OUTA
	)

virtual

Reimplemented in C81InterpolatedAeroData, C81MultipleAeroData, C81AeroData, and STAHRAeroData.

Definition at line 377 of file aerodata.cc.

References MBDYN_EXCEPT_ARGS.

Referenced by TheodorsenAeroData::AssRes(), AerodynamicBody::AssVec(), AerodynamicBeam::AssVec(), AerodynamicBeam2::AssVec(), GetForcesJacCenteredDiff_int(), and GetForcesJacForwardDiff_int().

 {
         silent_cerr("AeroData: GetForces() is undefined" << std::endl);
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 }

int AeroData::GetForcesJac	(	int	i,
		const doublereal *	W,
		doublereal *	TNG,
		Mat6x6 &	J,
		outa_t &	OUTA
	)

virtual

Reimplemented in C81InterpolatedAeroData, C81MultipleAeroData, C81AeroData, and STAHRAeroData.

Definition at line 384 of file aerodata.cc.

References MBDYN_EXCEPT_ARGS.

Referenced by TheodorsenAeroData::AssJac(), AerodynamicBody::AssJac(), AerodynamicBeam::AssJac(), and AerodynamicBeam2::AssJac().

 {
         silent_cerr("AeroData: GetForcesJac() is undefined" << std::endl);
         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
 }

int AeroData::GetForcesJacCenteredDiff_int	(	int	i,
		const doublereal *	W,
		doublereal *	TNG,
		Mat6x6 &	J,
		outa_t &	OUTA
	)

protected

Definition at line 313 of file aerodata.cc.

References GetForces(), Vec3::Norm(), Omega, and Mat6x6::Put().

 {
         const doublereal epsilon = 1.e-3;
         const doublereal nu = 1.e-9;
 
         Vec3 V(&W[0]);
         Vec3 Omega(&W[3]);
 
         doublereal dv = V.Norm();
         doublereal dw = Omega.Norm();
 
         doublereal TNGp[6];
         doublereal TNGm[6];
 
         GetForces(i, W, TNG0, OUTA);
 
         doublereal *WW = const_cast<doublereal *>(W);
         for (unsigned int iColm1 = 0; iColm1 < 6; iColm1++)     {
                 doublereal dorig;
                 doublereal delta;
 
                 dorig = WW[iColm1];
                 if (iColm1 < 3) {
                         delta = dv*epsilon + nu;
 
                 } else {
                         delta = dw*epsilon + nu;
                 }
 
                 WW[iColm1] = dorig + delta;
 
                 GetForces(i, WW, TNGp, OUTA);
 
                 WW[iColm1] = dorig - delta;
 
                 GetForces(i, WW, TNGm, OUTA);
 
                 doublereal delta2 = 2.*delta;
                 for (unsigned int iRowm1 = 0; iRowm1 < 6; iRowm1++) {
                         J.Put(iRowm1 + 1, iColm1 + 1, (TNGp[iRowm1] - TNGm[iRowm1])/delta2);
                 }
 
                 WW[iColm1] = dorig;
         }
 
         return 0;
 }

Here is the call graph for this function:

int AeroData::GetForcesJacForwardDiff_int	(	int	i,
		const doublereal *	W,
		doublereal *	TNG,
		Mat6x6 &	J,
		outa_t &	OUTA
	)

protected

Definition at line 271 of file aerodata.cc.

References GetForces(), Vec3::Norm(), Omega, and Mat6x6::Put().

Referenced by STAHRAeroData::GetForcesJac(), C81AeroData::GetForcesJac(), C81MultipleAeroData::GetForcesJac(), and C81InterpolatedAeroData::GetForcesJac().

 {
         const doublereal epsilon = 1.e-3;
         const doublereal nu = 1.e-9;
 
         Vec3 V(&W[0]);
         Vec3 Omega(&W[3]);
 
         doublereal dv = V.Norm();
         doublereal dw = Omega.Norm();
 
         doublereal TNG[6];
 
         GetForces(i, W, TNG0, OUTA);
 
         doublereal *WW = const_cast<doublereal *>(W);
         for (unsigned int iColm1 = 0; iColm1 < 6; iColm1++)     {
                 doublereal dorig;
                 doublereal delta;
 
                 dorig = WW[iColm1];
                 if (iColm1 < 3) {
                         delta = dv*epsilon + nu;
 
                 } else {
                         delta = dw*epsilon + nu;
                 }
                 WW[iColm1] = dorig + delta;
 
                 GetForces(i, WW, TNG, OUTA);
 
                 for (unsigned int iRowm1 = 0; iRowm1 < 6; iRowm1++) {
                         J.Put(iRowm1 + 1, iColm1 + 1, (TNG[iRowm1] - TNG0[iRowm1])/delta);
                 }
 
                 WW[iColm1] = dorig;
         }
 
         return 0;
 }

Here is the call graph for this function:

unsigned int AeroData::iGetNumDof ( void ) const

virtual

Reimplemented in TheodorsenAeroData.

Definition at line 362 of file aerodata.cc.

Referenced by AerodynamicBody::AssJac(), AerodynamicBeam::AssJac(), AerodynamicBeam2::AssJac(), AerodynamicBody::AssVec(), AerodynamicBeam::AssVec(), AerodynamicBeam2::AssVec(), AerodynamicBody::InitialAssRes(), ReadAerodynamicBeam(), ReadAerodynamicBeam2(), ReadAerodynamicBody(), and TheodorsenAeroData::TheodorsenAeroData().

 {
         return 0;
 }

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

pure virtual

Implemented in TheodorsenAeroData, C81InterpolatedAeroData, C81MultipleAeroData, C81AeroData, and STAHRAeroData.

Referenced by TheodorsenAeroData::Restart(), AerodynamicBody::Restart(), AerodynamicBeam::Restart(), and AerodynamicBeam2::Restart().

std::ostream & AeroData::RestartUnsteady ( std::ostream & out ) const

Definition at line 252 of file aerodata.cc.

References BIELAWA, HARRIS, and unsteadyflag.

Referenced by STAHRAeroData::Restart(), C81AeroData::Restart(), C81MultipleAeroData::Restart(), and C81InterpolatedAeroData::Restart().

 {
         switch (unsteadyflag) {
         case AeroData::HARRIS:
                 out << ", unsteady, harris";
                 break;
 
         case AeroData::BIELAWA:
                 out << ", unsteady, bielawa";
                 break;
 
         default:
                 break;
         }
 
         return out;
 }

void AeroData::SetAirData	(	const doublereal &	rho,
		const doublereal &	c
	)

virtual

Reimplemented in TheodorsenAeroData.

Definition at line 214 of file aerodata.cc.

References c, vam_t::density, vam_t::sound_celerity, and VAM.

Referenced by AerodynamicBody::AssJac(), AerodynamicBeam::AssJac(), AerodynamicBeam2::AssJac(), AerodynamicBody::AssVec(), AerodynamicBeam::AssVec(), AerodynamicBeam2::AssVec(), and TheodorsenAeroData::SetAirData().

 {
         VAM.density = rho;
         VAM.sound_celerity = c;
 }

void AeroData::SetSectionData	(	const doublereal &	abscissa,
		const doublereal &	chord,
		const doublereal &	forcepoint,
		const doublereal &	velocitypoint,
		const doublereal &	twist,
		const doublereal &	omega = `0.`
	)

virtual

Reimplemented in TheodorsenAeroData, C81InterpolatedAeroData, and C81MultipleAeroData.

Definition at line 237 of file aerodata.cc.

References vam_t::bc_position, vam_t::chord, vam_t::force_position, Omega, vam_t::twist, and VAM.

Referenced by AerodynamicBody::AssJac(), AerodynamicBeam::AssJac(), AerodynamicBeam2::AssJac(), AerodynamicBody::AssVec(), AerodynamicBeam::AssVec(), AerodynamicBeam2::AssVec(), C81MultipleAeroData::SetSectionData(), C81InterpolatedAeroData::SetSectionData(), and TheodorsenAeroData::SetSectionData().

 {
         VAM.chord = chord;
         VAM.force_position = forcepoint;
         VAM.bc_position = velocitypoint;
         VAM.twist = twist;
         Omega = omega;
 }

int AeroData::StorageSize ( void ) const

protectedvirtual

Implements AeroMemory.

Definition at line 221 of file aerodata.cc.

References BIELAWA, HARRIS, MBDYN_EXCEPT_ARGS, STEADY, and unsteadyflag.

 {
         switch (unsteadyflag) {
         case AeroData::HARRIS:
         case AeroData::BIELAWA:
                 return 3;
 
         case AeroData::STEADY:
                 return 0;
 
         default:
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 }