#include <aerodata_impl.h>

Inheritance diagram for C81InterpolatedAeroData:

Collaboration diagram for C81InterpolatedAeroData:

Public Member Functions
	C81InterpolatedAeroData (int i_p, int i_dim, AeroData::UnsteadyModel u, std::vector< unsigned > &p, std::vector< doublereal > &ub, std::vector< const c81_data * > &d, doublereal dcltol, DriveCaller *ptime=0)

	~C81InterpolatedAeroData (void)

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

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

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

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

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

virtual	~AeroData (void)

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

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

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 Attributes
std::vector< unsigned >	profiles

std::vector< doublereal >	upper_bounds

std::vector< const c81_data * >	data

std::vector< c81_data >	i_data

Protected Attributes inherited from AeroData
UnsteadyModel	unsteadyflag

vam_t	VAM

doublereal	Omega

Protected Attributes inherited from AeroMemory
DriveCaller *	pTime

Additional Inherited Members
Public Types inherited from AeroData
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 }

Protected Member Functions inherited from AeroData
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)

Detailed Description

Definition at line 121 of file aerodata_impl.h.

Constructor & Destructor Documentation

C81InterpolatedAeroData::C81InterpolatedAeroData	(	int	i_p,
		int	i_dim,
		AeroData::UnsteadyModel	u,
		std::vector< unsigned > &	p,
		std::vector< doublereal > &	ub,
		std::vector< const c81_data * > &	d,
		doublereal	dcltol,
		DriveCaller *	ptime = `0`
	)

Definition at line 244 of file aerodata_impl.cc.

References ASSERT, c81_data_merge(), data, PntWght::dGetPnt(), GaussDataIterator::fGetNext(), GaussDataIterator::GetFirst(), i_data, MBDYN_EXCEPT_ARGS, profiles, grad::sqrt(), and upper_bounds.

 : AeroData(i_p, i_dim, u, ptime),
 profiles(p), upper_bounds(ub), data(d),
 i_data(i_p*i_dim)
 {
         ASSERT(!profiles.empty());
         ASSERT(!upper_bounds.empty());
         ASSERT(!data.empty());
         ASSERT(profiles.size() == upper_bounds.size());
         ASSERT(profiles.size() == data.size());
 
         ASSERT(i_dim >= 1);
         ASSERT(i_dim <= 3);
         ASSERT(!i_data.empty());
 
         const doublereal bs[3][4] = {
                 { -1., 1. },
                 { -1., 0., 1. },
                 { -1., -1./std::sqrt(3.), 1./std::sqrt(3.), 1. }
         };
 
         int i_point = 0;
         for (int b = 0; b < i_dim; b++) {
                 GaussDataIterator GDI(i_p);
                 PntWght PW = GDI.GetFirst();
                 do {
                         doublereal dCsi = PW.dGetPnt();
 
                         doublereal bl = bs[i_dim - 1][b + 1] - bs[i_dim - 1][b];
                         doublereal bm = bs[i_dim - 1][b + 1] + bs[i_dim - 1][b];
                         dCsi = (bm + dCsi*bl)/2.;
 
                         if (c81_data_merge(data.size(), &data[0], &upper_bounds[0],
                                 dCsi, dcltol, &i_data[i_point]))
                         {
                                 // function logs error
                                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                         }
 
                         i_point++;
                 } while (GDI.fGetNext(PW));
         }
 }

Here is the call graph for this function:

C81InterpolatedAeroData::~C81InterpolatedAeroData ( void )

Definition at line 296 of file aerodata_impl.cc.

References c81_data_destroy(), and i_data.

 {
         for (unsigned i = 0; i < i_data.size(); i++) {
                 c81_data_destroy(&i_data[i]);
         }
 }

Here is the call graph for this function:

Member Function Documentation

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

virtual

Reimplemented from AeroData.

Definition at line 330 of file aerodata_impl.cc.

References outa_t::alf1, outa_t::alf2, grad::atan2(), AeroData::BIELAWA, c81_aerod2_u(), AeroData::HARRIS, i_data, AeroMemory::Predict(), AeroData::unsteadyflag, AeroData::VAM, AeroData::VX, and AeroData::VY.

 {
         
         switch (unsteadyflag) {
         case AeroData::HARRIS:
         case AeroData::BIELAWA: {
                 Predict(i, atan2(-W[VY], W[VX]), OUTA.alf1, OUTA.alf2);
                 break;
         }
 
         default:
                 break;
         }
 
         return c81_aerod2_u(const_cast<doublereal *>(W), &VAM, TNG, &OUTA,
                 const_cast<c81_data *>(&i_data[i]), unsteadyflag);
 }

Here is the call graph for this function:

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

virtual

Reimplemented from AeroData.

Definition at line 349 of file aerodata_impl.cc.

References AeroData::GetForcesJacForwardDiff_int().

 {
         return AeroData::GetForcesJacForwardDiff_int(i, W, TNG, J, OUTA);
 }

Here is the call graph for this function:

std::ostream & C81InterpolatedAeroData::Restart ( std::ostream & out ) const

virtual

Implements AeroData.

Definition at line 304 of file aerodata_impl.cc.

References profiles, AeroData::RestartUnsteady(), and upper_bounds.

 {
         out << "C81, interpolated";
         for (unsigned i = 0; i < profiles.size(); i++) {
                 out << ", " << profiles[i] << ", " << upper_bounds[i];
         }
 
         return RestartUnsteady(out);
 }

Here is the call graph for this function:

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

virtual

Reimplemented from AeroData.

Definition at line 315 of file aerodata_impl.cc.

References ASSERT, and AeroData::SetSectionData().

 {
         ASSERT(abscissa >= -1. && abscissa <= 1.);
 
         AeroData::SetSectionData(abscissa, chord, forcepoint, velocitypoint,
                 twist, omega);
 }