Класс, определяющий вид механической системы More...

#include <Mechanics2DRigidOscillPart.h>

Inheritance diagram for VM2D::MechanicsRigidOscillPart:

Collaboration diagram for VM2D::MechanicsRigidOscillPart:

[legend]

Public Member Functions
Point2D &	getV ()
	текущая скорость профиля

Point2D &	getVOld ()

Point2D &	getR ()
	текущее отклонение профиля

Point2D &	getROld ()

double &	getW ()
	текущая угловая скорость профиля

double &	getWOld ()

double &	getPhi ()
	текущий угол поворота профиля

double &	getPhiOld ()

bool &	getStrongCoupling ()

	MechanicsRigidOscillPart (const World2D &W_, size_t numberInPassport_)
	Конструктор

	~MechanicsRigidOscillPart ()
	Деструктор

virtual void	GetHydroDynamForce () override
	Вычисление гидродинамической силы, действующей на профиль

virtual Point2D	VeloOfAirfoilRcm (double currTime) override
	Вычисление скорости центра масс профиля

virtual Point2D	PositionOfAirfoilRcm (double currTime) override
	Вычисление положения центра масс профиля

virtual double	AngularVelocityOfAirfoil (double currTime) override
	Вычисление угловой скорости профиля

virtual double	AngleOfAirfoil (double currTime) override
	Вычисление угла поворота профиля

virtual void	VeloOfAirfoilPanels (double currTime) override
	Вычисление скоростей начал панелей

virtual void	ReadSpecificParametersFromDictionary () override
	Чтение параметров конкретной механической системы

virtual void	Move () override
	Перемещение профиля в соответствии с законом

void	MoveKinematic ()

void	MoveOnlyVelo ()

void	Initialize (Point2D Vcm0_, Point2D Rcm0_, double Wcm0_, double Phi0_)
	Задание начального положения и начальной скорости

void	GenerateForcesHeader ()
	Генерация заголовка файла нагрузок

void	GeneratePositionHeader ()
	Генерация заголовка файла положения профиля

void	GenerateForcesString ()
	Сохранение строки со статистикой в файл нагрузок

void	GeneratePositionString ()
	Сохранение строки со статистикой в файл нагрузок

Public Attributes
Point2D	Vcm0
	Начальная скорость центра и угловая скорость

double	Wcm0

Point2D	Rcm0
	Начальное положение профиля

double	Phi0

Point2D	Vcm
	Текущие скорость центра и угловая скорость

double	Wcm

Point2D	Rcm
	Текущие положение профиля

double	Phi

Point2D	VcmOld
	Скорость и отклонение с предыдущего шага

double	WcmOld

Point2D	RcmOld
	Текущие положение профиля

double	PhiOld

double	circulation
	Текущая циркуляция скорости по границе профиля

double	circulationOld
	Циркуляция скорости по границе профиля с предыдущего шага

const bool	isMoves
	Переменная, отвечающая за то, двигается профиль или нет

const bool	isDeform
	Переменная, отвечающая за то, деформируется профиль или нет

Point2D	hydroDynamForce
	Вектор гидродинамической силы и момент, действующие на профиль

double	hydroDynamMoment

Point2D	viscousForce
	Вектор силы и момент вязкого трения, действующие на профиль

double	viscousMoment

Protected Attributes
const World2D &	W
	Константная ссылка на решаемую задачу

const size_t	numberInPassport
	Номер профиля в паспорте

Airfoil &	afl

const Boundary &	boundary

const VortexesParams &	virtVortParams

std::unique_ptr< VMlib::StreamParser >	mechParamsParser
	Умный указатель на парсер параметров механической системы

Private Member Functions
void	ReadParametersFromDictionary ()
	Парсинг списка параметров механической системы

Private Attributes
double	m
	масса профиля

double	J
	момент инерции профиля

Point2D	b
	параметр демпфирования механической системы

double	bw

Point2D	k
	параметр жесткости механической системы

double	kw

Point2D	initDisplacement
	начальное отклонение

double	initAngularDisplacement

Point2D	initVelocity
	начальные скорости

double	initAngularVelocity

bool	strongCoupling
	признак полунеявной схемы связывания

Detailed Description

Класс, определяющий вид механической системы

Упруго закрепленное тело, метод расщепления

Author: Марчевский Илья Константинович; Сокол Ксения Сергеевна; Рятина Евгения Павловна; Колганова Александра Олеговна

\Version 1.14

Date: 6 марта 2026 г.

Definition at line 64 of file Mechanics2DRigidOscillPart.h.

Constructor & Destructor Documentation

◆ MechanicsRigidOscillPart()

MechanicsRigidOscillPart::MechanicsRigidOscillPart	(	const World2D &	W_,
		size_t	numberInPassport_
	)

Конструктор

Parameters

[in]	W_	константная ссылка на решаемую задачу
[in]	numberInPassport_	номер профиля в паспорте задачи

Definition at line 52 of file Mechanics2DRigidOscillPart.cpp.

    :
    Mechanics(W_, numberInPassport_, true, false)
    //, V0({ 0.0, 0.0 })
    //, r0({ W_.getAirfoil(numberInPassport_).rcm[0], W_.getAirfoil(numberInPassport_).rcm[1] })
{
    Vcm0 = { 0.0, 0.0 };
    Rcm0 = { W_.getAirfoil(numberInPassport_).rcm[0], W_.getAirfoil(numberInPassport_).rcm[1] };
    Vcm = Vcm0;
    Rcm = Rcm0;
    VcmOld = Vcm0;
    RcmOld = Rcm;
 
    strongCoupling = false;
        
    MechanicsRigidOscillPart::ReadSpecificParametersFromDictionary();
    Initialize(initVelocity, W_.getAirfoil(numberInPassport_).rcm + initDisplacement, initAngularVelocity, W_.getAirfoil(numberInPassport_).phiAfl + initAngularDisplacement);
};

Here is the call graph for this function:

◆ ~MechanicsRigidOscillPart()

VM2D::MechanicsRigidOscillPart::~MechanicsRigidOscillPart ( )

inline

Деструктор

Definition at line 120 of file Mechanics2DRigidOscillPart.h.

120{};

Member Function Documentation

◆ AngleOfAirfoil()

double MechanicsRigidOscillPart::AngleOfAirfoil ( double currTime )

overridevirtual

Вычисление угла поворота профиля

Parameters

[in] currTime текущее время

Implements VM2D::Mechanics.

Definition at line 158 of file Mechanics2DRigidOscillPart.cpp.

{
    if (afl.phiAfl != Phi)
    {
        std::cout << "afl.phiAfl != Phi" << std::endl;
        exit(100600);
    }
 
    return afl.phiAfl;
}//AngleOfAirfoil(...)

◆ AngularVelocityOfAirfoil()

double MechanicsRigidOscillPart::AngularVelocityOfAirfoil ( double currTime )

overridevirtual

Вычисление угловой скорости профиля

Parameters

[in] currTime текущее время

Implements VM2D::Mechanics.

Definition at line 153 of file Mechanics2DRigidOscillPart.cpp.

{
    return Wcm;
}//AngularVelocityOfAirfoil(...)

◆ GenerateForcesHeader()

void Mechanics::GenerateForcesHeader ( )

inherited

Генерация заголовка файла нагрузок

Definition at line 125 of file Mechanics2D.cpp.

{
    std::stringstream forceFileName, forceFileNameCsv;
    forceFileName << W.getPassport().dir << "forces-airfoil-" << numberInPassport;
    forceFileNameCsv << W.getPassport().dir << "forces-airfoil-" << numberInPassport << ".csv";
 
    std::ofstream newForcesFile(forceFileName.str());
    std::ofstream newForcesFileCsv(forceFileNameCsv.str());
 
    VMlib::PrintLogoToTextFile(newForcesFile, forceFileName.str(), "Hydrodynamic loads for the airfoil " + W.getPassport().airfoilParams[numberInPassport].fileAirfoil);
 
    if (!W.getPassport().calcCoefficients)
        VMlib::PrintHeaderToTextFile(newForcesFile, "currentStep     currentTime     Fx     Fy     Mz     Ftaux     Ftauy     Mtau");
    else
        VMlib::PrintHeaderToTextFile(newForcesFile, "currentStep     currentTime     CFx    CFy    CMz    CFtaux    CFtauy    CMtau");
 
    if (!W.getPassport().calcCoefficients)
        newForcesFileCsv << "step,time,Fx,Fy,Mz,Ftaux,Ftauy,Mtau";
    else
        newForcesFileCsv << "step,time,CFx,CFy,CMz,CFtaux,CFtauy,CMtau";
 
    newForcesFile.close();
    newForcesFile.clear();
 
    newForcesFileCsv.close();
    newForcesFileCsv.clear();
 
}//GenerateForcesHeader()

Here is the call graph for this function:

◆ GenerateForcesString()

void Mechanics::GenerateForcesString ( )

inherited

Сохранение строки со статистикой в файл нагрузок

Definition at line 183 of file Mechanics2D.cpp.

{
    std::stringstream forceFileName, forceFileNameCsv;
    forceFileName << W.getPassport().dir << "forces-airfoil-" << numberInPassport;
    forceFileNameCsv << W.getPassport().dir << "forces-airfoil-" << numberInPassport << ".csv";
 
    //double cShock = (W.getPassport().physicalProperties.getCurrTime() > W.getPassport().physicalProperties.timeAccel + 2.0 * W.getPassport().timeDiscretizationProperties.dt) ? 1.0 : 0.0;
    double cShock = 1.0;
 
    Point2D cartesianHydroForce = cShock * hydroDynamForce;
    Point2D cartesianViscoForce = cShock * viscousForce;
 
    double cartesianHydroMoment = cShock * hydroDynamMoment;
    double cartesianViscoMoment = cShock * viscousMoment;
 
    double rotationAngle = W.getPassport().airfoilParams[numberInPassport].angle;
    
    //Для обдува ветром, когда углы считаются по компасу
    //if (W.getPassport().geographicalAngles)
    //{
    //  cartesianHydroForce = cartesianHydroForce.rotated(-0.5 * PI - W.getPassport().airfoilParams[numberInPassport].angle);
    //  cartesianViscoForce = cartesianViscoForce.rotated(-0.5 * PI - W.getPassport().airfoilParams[numberInPassport].angle);
    //  //rotationAngle += 0.5 * PI;
    //}
 
    //if (W.getPassport().rotateForces)
    //{
    //  cartesianHydroForce = cartesianHydroForce.rotated(rotationAngle);
    //  cartesianViscoForce = cartesianViscoForce.rotated(rotationAngle);
    //}
 
    double flowPress = 0.5 * W.getPassport().physicalProperties.rho * sqr(W.getPassport().physicalProperties.vRef);
    if (W.getPassport().calcCoefficients)
    {
        const double& chord = W.getPassport().airfoilParams[numberInPassport].chord;
        cartesianHydroForce /= (flowPress * chord);
        cartesianViscoForce /= (flowPress * chord);
 
        cartesianHydroMoment /= (flowPress * sqr(chord));
        cartesianViscoMoment /= (flowPress * sqr(chord));
    }
 
 
    std::ofstream forcesFile(forceFileName.str(), std::ios::app);
    forcesFile << std::endl << W.getCurrentStep() << "  " << W.getCurrentTime() << "    " << cartesianHydroForce[0] << "    " << cartesianHydroForce[1] << "    " << cartesianHydroMoment << "  " << cartesianViscoForce[0] << "    " << cartesianViscoForce[1] << "    " << cartesianViscoMoment;
    forcesFile.close();
 
    std::ofstream forcesFileCsv(forceFileNameCsv.str(), std::ios::app);
    forcesFileCsv << std::endl << W.getCurrentStep() << "," << W.getCurrentTime() << "," << cartesianHydroForce[0] << "," << cartesianHydroForce[1] << "," << cartesianHydroMoment << "," << cartesianViscoForce[0] << "," << cartesianViscoForce[1] << "," << cartesianViscoMoment;
    forcesFileCsv.close();
 
}//GenerateForcesString()

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GeneratePositionHeader()

void Mechanics::GeneratePositionHeader ( )

inherited

Генерация заголовка файла положения профиля

Definition at line 156 of file Mechanics2D.cpp.

{
    if (isMoves)
    {
        std::stringstream positionFileName, positionFileNameCsv;
        positionFileName << W.getPassport().dir << "position-airfoil-" << numberInPassport;
        positionFileNameCsv << W.getPassport().dir << "position-airfoil-" << numberInPassport << ".csv";
 
        std::ofstream newPositionFile(positionFileName.str());
        std::ofstream newPositionFileCsv(positionFileNameCsv.str());
 
        VMlib::PrintLogoToTextFile(newPositionFile, positionFileName.str(), "Position of the airfoil " + W.getPassport().airfoilParams[numberInPassport].fileAirfoil);
 
        VMlib::PrintHeaderToTextFile(newPositionFile, "currentStep     currentTime     x     y     phi     Vx     Vy     w");
 
        newPositionFileCsv << "step,time,x,y,phi,Vx,Vy,w";
 
        newPositionFile.close();
        newPositionFile.clear();
 
        newPositionFileCsv.close();
        newPositionFileCsv.clear();
    }
}//GeneratePositionHeader()

Here is the call graph for this function:

◆ GeneratePositionString()

void Mechanics::GeneratePositionString ( )

inherited

Сохранение строки со статистикой в файл нагрузок

Definition at line 238 of file Mechanics2D.cpp.

{
    if (isMoves)
    {
        std::stringstream positionFileName, positionFileNameCsv;
        positionFileName << W.getPassport().dir << "position-airfoil-" << numberInPassport;
        positionFileNameCsv << W.getPassport().dir << "position-airfoil-" << numberInPassport << ".csv";
 
        std::ofstream positionFile(positionFileName.str(), std::ios::app);
        positionFile << std::endl << W.getCurrentStep() << "    " << W.getCurrentTime() << "    " << afl.rcm[0] << "    " << afl.rcm[1] << "    " << Phi << "   " << Vcm[0] << "    " << Vcm[1] << "    " << Wcm;
        positionFile.close();
 
        std::ofstream positionFileCsv(positionFileNameCsv.str(), std::ios::app);
        positionFileCsv << std::endl << W.getCurrentStep() << "," << W.getCurrentTime() << "," << afl.rcm[0] << "," << afl.rcm[1] << "," << Phi << "," << Vcm[0] << "," << Vcm[1] << "," << Wcm;
        positionFileCsv.close();
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetHydroDynamForce()

void MechanicsRigidOscillPart::GetHydroDynamForce ( )

overridevirtual

Вычисление гидродинамической силы, действующей на профиль

Todo:: Учитываем только нулевой момент решения. Надо ли учитывать остальные?

Implements VM2D::Mechanics.

Definition at line 72 of file Mechanics2DRigidOscillPart.cpp.

{
    W.getTimers().start("Force");
 
    const double& dt = W.getPassport().timeDiscretizationProperties.dt;
 
    hydroDynamForce = { 0.0, 0.0 };
    hydroDynamMoment = 0.0;
 
    viscousForce = { 0.0, 0.0 };
    viscousMoment = 0.0;
 
    Point2D hDFGam = { 0.0, 0.0 };      //гидродинамические силы, обусловленные присоед.завихренностью
    Point2D hDFdelta = { 0.0, 0.0 };    //гидродинамические силы, обусловленные приростом завихренности
    Point2D hDFQ = { 0.0, 0.0 };        //гидродинамические силы, обусловленные присоед.источниками
 
    double hDMGam = 0.0;                //гидродинамический момент, обусловленный присоед.завихренностью
    double hDMdelta = 0.0;              //гидродинамический момент, обусловленный приростом завихренности
    double hDMQ = 0.0;                  //гидродинамический момент, обусловленный присоед.источниками
    for (size_t i = 0; i < afl.getNumberOfPanels(); ++i)
    {
        Point2D rK = 0.5 * (afl.getR(i + 1) + afl.getR(i)) - afl.rcm;
 
        Point2D velK = 0.5 * (afl.getV(i) + afl.getV(i + 1));
        double gAtt = (velK & afl.tau[i]);
 
        double gAttOld = 0.0;
        if (W.getCurrentStep() > 0)
        {
            auto oldAfl = W.getOldAirfoil(numberInPassport);
            gAttOld = ((0.5 * (oldAfl.getV(i) + oldAfl.getV(i + 1))) & oldAfl.tau[i]);
        }
 
        double deltaGAtt = gAtt - gAttOld;
 
        double qAtt = (velK & afl.nrm[i]);
 
        double deltaK = boundary.sheets.freeVortexSheet(i, 0) * afl.len[i] - afl.gammaThrough[i] + deltaGAtt * afl.len[i];
 
        /*1*/
        hDFdelta += deltaK * Point2D({ -rK[1], rK[0] });
        hDMdelta += 0.5 * deltaK * rK.length2();
 
        /*2*/
        hDFGam += 0.5 * velK.kcross() * gAtt * afl.len[i];
        hDMGam += 0.5 * (rK ^ velK.kcross()) * gAtt * afl.len[i];
 
        /*3*/
        hDFQ -= 0.5 * velK * qAtt * afl.len[i];
        hDMQ -= 0.5 * (rK ^ velK) * qAtt * afl.len[i];
    }
 
    const double rho = W.getPassport().physicalProperties.rho;
 
    hydroDynamForce = rho * (hDFGam + hDFdelta * (1.0 / dt) + hDFQ);
    hydroDynamMoment = rho * (hDMGam + hDMdelta / dt + hDMQ);
 
    if ((W.getPassport().physicalProperties.nu > 0.0)/* && (W.currentStep > 0)*/)
        for (size_t i = 0; i < afl.getNumberOfPanels(); ++i)
        {
            Point2D rK = 0.5 * (afl.getR(i + 1) + afl.getR(i)) - afl.rcm;
            viscousForce += rho * afl.viscousStress[i] * afl.tau[i];
            viscousMoment += rho * (afl.viscousStress[i] * afl.tau[i]) & rK;
        }
 
    W.getTimers().stop("Force");
}// GetHydroDynamForce()

Here is the call graph for this function:

◆ getPhi()

double & VM2D::MechanicsRigidOscillPart::getPhi ( )

inline

текущий угол поворота профиля

Definition at line 108 of file Mechanics2DRigidOscillPart.h.

108{ return Phi; };

Here is the caller graph for this function:

◆ getPhiOld()

double & VM2D::MechanicsRigidOscillPart::getPhiOld ( )

inline

Definition at line 109 of file Mechanics2DRigidOscillPart.h.

109{ return PhiOld; };

VM2D::Mechanics::PhiOld

double PhiOld

Definition Mechanics2D.h:118

Here is the caller graph for this function:

◆ getR()

Point2D & VM2D::MechanicsRigidOscillPart::getR ( )

inline

текущее отклонение профиля

Definition at line 100 of file Mechanics2DRigidOscillPart.h.

100{ return Rcm; };

Here is the caller graph for this function:

◆ getROld()

Point2D & VM2D::MechanicsRigidOscillPart::getROld ( )

inline

Definition at line 101 of file Mechanics2DRigidOscillPart.h.

101{ return RcmOld; };

Here is the caller graph for this function:

◆ getStrongCoupling()

bool & VM2D::MechanicsRigidOscillPart::getStrongCoupling ( )

inline

Definition at line 111 of file Mechanics2DRigidOscillPart.h.

111{ return strongCoupling; };

Here is the caller graph for this function:

◆ getV()

Point2D & VM2D::MechanicsRigidOscillPart::getV ( )

inline

текущая скорость профиля

Definition at line 96 of file Mechanics2DRigidOscillPart.h.

96{ return Vcm; };

Here is the caller graph for this function:

◆ getVOld()

Point2D & VM2D::MechanicsRigidOscillPart::getVOld ( )

inline

Definition at line 97 of file Mechanics2DRigidOscillPart.h.

97{ return VcmOld; };

Here is the caller graph for this function:

◆ getW()

double & VM2D::MechanicsRigidOscillPart::getW ( )

inline

текущая угловая скорость профиля

Definition at line 104 of file Mechanics2DRigidOscillPart.h.

104{ return Wcm; };

Here is the caller graph for this function:

◆ getWOld()

double & VM2D::MechanicsRigidOscillPart::getWOld ( )

inline

Definition at line 105 of file Mechanics2DRigidOscillPart.h.

105{ return WcmOld; };

VM2D::Mechanics::WcmOld

double WcmOld

Definition Mechanics2D.h:115

Here is the caller graph for this function:

◆ Initialize()

void Mechanics::Initialize	(	Point2D	Vcm0_,
		Point2D	Rcm0_,
		double	Wcm0_,
		double	Phi0_
	)

inherited

Задание начального положения и начальной скорости

Parameters

[in]	Vcm0_	- скорость центра масс
[in]	Rcm0_	- положение центра масс
[in]	Wcm0_	- угловая скорость центра масс
[in]	Phi0_	- угол поворота центра масс

Definition at line 92 of file Mechanics2D.cpp.

{
    Vcm0 = Vcm0_;
    Wcm0 = Wcm0_;
    Rcm0 = Rcm0_;
    Phi0 = Phi0_;
 
    Vcm = Vcm0;
    Wcm = Wcm0;
    Rcm = Rcm0;
    Phi = Phi0;
    VcmOld = Vcm0;
    WcmOld = Wcm0;
    RcmOld = Rcm0;
    PhiOld = Phi0;
 
    circulation = 2.0 * afl.area * Wcm;
    circulationOld = circulation;
 
    afl.Move(Rcm - W.getAirfoil(numberInPassport).rcm);
    afl.Rotate(Phi - W.getAirfoil(numberInPassport).phiAfl);
}//Initialize(...)

Here is the call graph for this function:

Here is the caller graph for this function:

◆ Move()

void MechanicsRigidOscillPart::Move ( )

overridevirtual

Перемещение профиля в соответствии с законом

Implements VM2D::Mechanics.

Definition at line 186 of file Mechanics2DRigidOscillPart.cpp.

{
    Point2D meff;
    //if (W.getPassport().airfoilParams[numberInPassport].addedMass.length2() > 0)
    //  meff = Point2D{ m + W.getPassport().airfoilParams[numberInPassport].addedMass[0], m + W.getPassport().airfoilParams[numberInPassport].addedMass[1] };
    //else
        meff = Point2D{ m, m };
 
        double Jeff = J;
 
    VcmOld = Vcm;
    RcmOld = Rcm;
    PhiOld = Phi;
    WcmOld = Wcm;   
 
    Point2D dr, dV;
    double dphi, dw;
 
    //W.getInfo('t') << "k = " << k << std::endl;
 
    if (k[0] > 0)
    {
        double dt = W.getPassport().timeDiscretizationProperties.dt;
        Point2D kk[4];
 
        kk[0] = { Vcm[0], (hydroDynamForce[0] - 2.0 * b[0] * Vcm[0] - k[0] * Rcm[0]) / meff[0] };
        kk[1] = { Vcm[0] + 0.5 * dt * kk[0][1], (hydroDynamForce[0] - 2.0 * b[0] * (Vcm[0] + 0.5 * dt * kk[0][1]) - k[0] * (Rcm[0] + 0.5 * dt * kk[0][0])) / meff[0] };
        kk[2] = { Vcm[0] + 0.5 * dt * kk[1][1], (hydroDynamForce[0] - 2.0 * b[0] * (Vcm[0] + 0.5 * dt * kk[1][1]) - k[0] * (Rcm[0] + 0.5 * dt * kk[1][0])) / meff[0] };
        kk[3] = { Vcm[0] + dt * kk[2][1], (hydroDynamForce[0] - 2.0 * b[0] * (Vcm[0] + dt * kk[2][1]) - k[0] * (Rcm[0] + dt * kk[2][0])) / meff[0] };
 
        dr[0] = dt * (kk[0][0] + 2. * kk[1][0] + 2. * kk[2][0] + kk[3][0]) / 6.0;
        dV[0] = dt * (kk[0][1] + 2. * kk[1][1] + 2. * kk[2][1] + kk[3][1]) / 6.0;
    }
    else
    {
        dr[0] = 0.0;
        dV[0] = 0.0;
    }
 
 
 
    if (k[1] > 0)
    {
        double dt = W.getPassport().timeDiscretizationProperties.dt;
        Point2D kk[4];
        kk[0] = { Vcm[1], (hydroDynamForce[1] - 2.0 * b[1] * Vcm[1] - k[1] * Rcm[1]) / meff[1]};
        kk[1] = { Vcm[1] + 0.5 * dt * kk[0][1], (hydroDynamForce[1] - 2.0 * b[1] * (Vcm[1] + 0.5 * dt * kk[0][1]) - k[1] * (Rcm[1] + 0.5 * dt * kk[0][0])) / meff[1]};
        kk[2] = { Vcm[1] + 0.5 * dt * kk[1][1], (hydroDynamForce[1] - 2.0 * b[1] * (Vcm[1] + 0.5 * dt * kk[1][1]) - k[1] * (Rcm[1] + 0.5 * dt * kk[1][0])) / meff[1]};
        kk[3] = { Vcm[1] + dt * kk[2][1], (hydroDynamForce[1] - 2.0 * b[1] * (Vcm[1] + dt * kk[2][1]) - k[1] * (Rcm[1] + dt * kk[2][0])) / meff[1]};
                
        dr[1] = dt * (kk[0][0] + 2. * kk[1][0] + 2. * kk[2][0] + kk[3][0]) / 6.0;
        dV[1] = dt * (kk[0][1] + 2. * kk[1][1] + 2. * kk[2][1] + kk[3][1]) / 6.0;
    }
    else
    {
        dr[1] = 0.0;
        dV[1] = 0.0;
    }
 
 
 
    if (kw > 0)
    {
        double dt = W.getPassport().timeDiscretizationProperties.dt;
        Point2D kk[4];
 
        kk[0] = { Wcm, (hydroDynamMoment - 2.0 * bw * Wcm - kw * Phi) / Jeff };
        kk[1] = { Wcm + 0.5 * dt * kk[0][1], (hydroDynamMoment - 2.0 * bw * (Wcm + 0.5 * dt * kk[0][1]) - kw * (Phi + 0.5 * dt * kk[0][0])) / Jeff };
        kk[2] = { Wcm + 0.5 * dt * kk[1][1], (hydroDynamMoment - 2.0 * bw * (Wcm + 0.5 * dt * kk[1][1]) - kw * (Phi + 0.5 * dt * kk[1][0])) / Jeff };
        kk[3] = { Wcm + dt * kk[2][1], (hydroDynamMoment - 2.0 * bw * (Wcm + dt * kk[2][1]) - kw * (Phi + dt * kk[2][0])) / Jeff };
 
        dphi = dt * (kk[0][0] + 2. * kk[1][0] + 2. * kk[2][0] + kk[3][0]) / 6.0;
        dw = dt * (kk[0][1] + 2. * kk[1][1] + 2. * kk[2][1] + kk[3][1]) / 6.0;
    }
    else
    {
        dphi = 0.0;
        dw = 0.0;
    }
 
    afl.Move(dr);
    afl.Rotate(dphi);
 
    Rcm += dr;
    Vcm += dV;
 
    Phi += dphi;
    Wcm += dw;
}//Move()

Here is the call graph for this function:

◆ MoveKinematic()

void MechanicsRigidOscillPart::MoveKinematic ( )

Definition at line 278 of file Mechanics2DRigidOscillPart.cpp.

{
    VcmOld = Vcm;
    RcmOld = Rcm;
 
    Point2D dr, dV;
    double dt = W.getPassport().timeDiscretizationProperties.dt;
 
    if (k[1] > 0)
    {           
        dr[1] = Vcm[1] * dt;
        dV[1] = 0.0;
    }
    else
    {
        dr[1] = 0.0;
        dV[1] = 0.0;
    }
 
 
    if (k[0] > 0)
    {       
        dr[0] = Vcm[0] * dt;
        dV[0] = 0.0;
    }
    else
    {
        dr[0] = 0.0;
        dV[0] = 0.0;
    }
 
    afl.Move(dr);   
            
    Rcm += dr;
    Vcm += dV;
}//MoveKinematic()

Here is the call graph for this function:

Here is the caller graph for this function:

◆ MoveOnlyVelo()

void MechanicsRigidOscillPart::MoveOnlyVelo ( )

Definition at line 316 of file Mechanics2DRigidOscillPart.cpp.

{
    Point2D meff;
    if (W.getPassport().airfoilParams[numberInPassport].addedMass.length2() > 0)
        meff = Point2D{ m + W.getPassport().airfoilParams[numberInPassport].addedMass[0], m + W.getPassport().airfoilParams[numberInPassport].addedMass[1] };
    else
        meff = Point2D{ m, m };
 
    Point2D dr, dV; 
 
    if (k[1] > 0)
    {
        double dt = W.getPassport().timeDiscretizationProperties.dt;
        Point2D kk[4];
    
        kk[0] = { Vcm[1], (hydroDynamForce[1] - 2.0 * b[1] * Vcm[1] - k[1] * Rcm[1]) / meff[1]};
        kk[1] = { Vcm[1] + 0.5 * dt * kk[0][1], (hydroDynamForce[1] - 2.0 * b[1] * (Vcm[1] + 0.5 * dt * kk[0][1]) - k[1] * (Rcm[1] + 0.5 * dt * kk[0][0])) / meff[1]};
        kk[2] = { Vcm[1] + 0.5 * dt * kk[1][1], (hydroDynamForce[1] - 2.0 * b[1] * (Vcm[1] + 0.5 * dt * kk[1][1]) - k[1] * (Rcm[1] + 0.5 * dt * kk[1][0])) / meff[1]};
        kk[3] = { Vcm[1] + dt * kk[2][1], (hydroDynamForce[1] - 2.0 * b[1] * (Vcm[1] + dt * kk[2][1]) - k[1] * (Rcm[1] + dt * kk[2][0])) / meff[1]};
 
        dr[1] = 0.0;// dt* (kk[0][0] + 2.0 * kk[1][0] + 2.0 * kk[2][0] + kk[3][0]) / 6.0;
        dV[1] = dt * (kk[0][1] + 2.0 * kk[1][1] + 2.0 * kk[2][1] + kk[3][1]) / 6.0;
    }
    else
    {
        dr[1] = 0.0;
        dV[1] = 0.0;
    }
 
 
    if (k[0] > 0)
    {
        double dt = W.getPassport().timeDiscretizationProperties.dt;
        Point2D kk[4];
        
        kk[0] = { Vcm[0], (hydroDynamForce[0] - 2.0 * b[0] * Vcm[0] - k[0] * Rcm[0]) / meff[0]};
        kk[1] = { Vcm[0] + 0.5 * dt * kk[0][1], (hydroDynamForce[0] - 2.0 * b[0] * (Vcm[0] + 0.5 * dt * kk[0][1]) - k[0] * (Rcm[0] + 0.5 * dt * kk[0][0])) / meff[0]};
        kk[2] = { Vcm[0] + 0.5 * dt * kk[1][1], (hydroDynamForce[0] - 2.0 * b[0] * (Vcm[0] + 0.5 * dt * kk[1][1]) - k[0] * (Rcm[0] + 0.5 * dt * kk[1][0])) / meff[0]};
        kk[3] = { Vcm[0] + dt * kk[2][1], (hydroDynamForce[0] - 2.0 * b[0] * (Vcm[0] + dt * kk[2][1]) - k[0] * (Rcm[0] + dt * kk[2][0])) / meff[0]};
 
        dr[0] = 0.0; // dt* (kk[0][0] + 2.0 * kk[1][0] + 2.0 * kk[2][0] + kk[3][0]) / 6.0;
        dV[0] = dt * (kk[0][1] + 2.0 * kk[1][1] + 2.0 * kk[2][1] + kk[3][1]) / 6.0;
    }
    else
    {
        dr[0] = 0.0;
        dV[0] = 0.0;
    }
 
    //afl.Move({ dx, dy });
    Vcm += dV;  
    
}//MoveOnlyVelo()

Here is the call graph for this function:

Here is the caller graph for this function:

◆ PositionOfAirfoilRcm()

Point2D MechanicsRigidOscillPart::PositionOfAirfoilRcm ( double currTime )

overridevirtual

Вычисление положения центра масс профиля

Parameters

[in] currTime текущее время

Implements VM2D::Mechanics.

Definition at line 148 of file Mechanics2DRigidOscillPart.cpp.

{
    return Rcm;
}//PositionOfAirfoilRcm(...)

◆ ReadParametersFromDictionary()

void Mechanics::ReadParametersFromDictionary ( )

privateinherited

Парсинг списка параметров механической системы

Definition at line 116 of file Mechanics2D.cpp.

{
    std::stringstream ss(W.getPassport().airfoilParams[afl.numberInPassport].mechanicalSystemParameters);
    
    mechParamsParser.reset(new VMlib::StreamParser(W.getInfo(), "mechanical parser", ss));
}//ReadParametersFromDictionary()

Here is the call graph for this function:

◆ ReadSpecificParametersFromDictionary()

virtual void VM2D::MechanicsRigidOscillPart::ReadSpecificParametersFromDictionary ( )

overridevirtual

Чтение параметров конкретной механической системы

Implements VM2D::Mechanics.

Here is the caller graph for this function:

◆ VeloOfAirfoilPanels()

void MechanicsRigidOscillPart::VeloOfAirfoilPanels ( double currTime )

overridevirtual

Вычисление скоростей начал панелей

Parameters

[in] currTime текущее время

Implements VM2D::Mechanics.

Definition at line 170 of file Mechanics2DRigidOscillPart.cpp.

{
    Point2D veloRcm = VeloOfAirfoilRcm(currTime);
 
    std::vector<Point2D> veloW(afl.getNumberOfPanels());
    for (size_t i = 0; i < afl.getNumberOfPanels(); ++i)
        veloW[i] = veloRcm + Wcm * (afl.getR(i) - Rcm).kcross();
 
    afl.setV(veloW);
 
    circulationOld = circulation;
    circulation = 2.0 * afl.area * Wcm;
 
}//VeloOfAirfoilPanels(...)

Here is the call graph for this function:

◆ VeloOfAirfoilRcm()

Point2D MechanicsRigidOscillPart::VeloOfAirfoilRcm ( double currTime )

overridevirtual

Вычисление скорости центра масс профиля

Parameters

[in] currTime текущее время

Implements VM2D::Mechanics.

Definition at line 142 of file Mechanics2DRigidOscillPart.cpp.

{
    return Vcm;
}//VeloOfAirfoilRcm(...)

Here is the caller graph for this function:

Member Data Documentation

◆ afl

Airfoil& VM2D::Mechanics::afl

protectedinherited

Константная ссылка на профиль
инициализируется автоматом в конструкторе, при помощи const_cast

Warning: использует const_cast для получения неконстантной ссылки

Definition at line 87 of file Mechanics2D.h.

◆ b

Point2D VM2D::MechanicsRigidOscillPart::b

private

параметр демпфирования механической системы

Definition at line 76 of file Mechanics2DRigidOscillPart.h.

◆ boundary

const Boundary& VM2D::Mechanics::boundary

protectedinherited

Константная ссылка на граничное условие
инициализируется автоматом в конструкторе

Definition at line 91 of file Mechanics2D.h.

◆ bw

double VM2D::MechanicsRigidOscillPart::bw

private

Definition at line 77 of file Mechanics2DRigidOscillPart.h.

◆ circulation

double VM2D::Mechanics::circulation

inherited

Текущая циркуляция скорости по границе профиля

Definition at line 121 of file Mechanics2D.h.

◆ circulationOld

double VM2D::Mechanics::circulationOld

inherited

Циркуляция скорости по границе профиля с предыдущего шага

Definition at line 124 of file Mechanics2D.h.

◆ hydroDynamForce

Point2D VM2D::Mechanics::hydroDynamForce

inherited

Вектор гидродинамической силы и момент, действующие на профиль

Definition at line 134 of file Mechanics2D.h.

◆ hydroDynamMoment

double VM2D::Mechanics::hydroDynamMoment

inherited

Definition at line 135 of file Mechanics2D.h.

◆ initAngularDisplacement

double VM2D::MechanicsRigidOscillPart::initAngularDisplacement

private

Definition at line 85 of file Mechanics2DRigidOscillPart.h.

◆ initAngularVelocity

double VM2D::MechanicsRigidOscillPart::initAngularVelocity

private

Definition at line 89 of file Mechanics2DRigidOscillPart.h.

◆ initDisplacement

Point2D VM2D::MechanicsRigidOscillPart::initDisplacement

private

начальное отклонение

Definition at line 84 of file Mechanics2DRigidOscillPart.h.

◆ initVelocity

Point2D VM2D::MechanicsRigidOscillPart::initVelocity

private

начальные скорости

Definition at line 88 of file Mechanics2DRigidOscillPart.h.

◆ isDeform

const bool VM2D::Mechanics::isDeform

inherited

Переменная, отвечающая за то, деформируется профиль или нет

Definition at line 131 of file Mechanics2D.h.

◆ isMoves

const bool VM2D::Mechanics::isMoves

inherited

Переменная, отвечающая за то, двигается профиль или нет

Definition at line 128 of file Mechanics2D.h.

◆ J

double VM2D::MechanicsRigidOscillPart::J

private

момент инерции профиля

Definition at line 73 of file Mechanics2DRigidOscillPart.h.

◆ k

Point2D VM2D::MechanicsRigidOscillPart::k

private

параметр жесткости механической системы

Definition at line 80 of file Mechanics2DRigidOscillPart.h.

◆ kw

double VM2D::MechanicsRigidOscillPart::kw

private

Definition at line 81 of file Mechanics2DRigidOscillPart.h.

◆ m

double VM2D::MechanicsRigidOscillPart::m

private

масса профиля

Definition at line 70 of file Mechanics2DRigidOscillPart.h.

◆ mechParamsParser

std::unique_ptr<VMlib::StreamParser> VM2D::Mechanics::mechParamsParser

protectedinherited

Умный указатель на парсер параметров механической системы

Definition at line 98 of file Mechanics2D.h.

◆ numberInPassport

const size_t VM2D::Mechanics::numberInPassport

protectedinherited

Номер профиля в паспорте

Definition at line 82 of file Mechanics2D.h.

◆ Phi

double VM2D::Mechanics::Phi

inherited

Definition at line 112 of file Mechanics2D.h.

◆ Phi0

double VM2D::Mechanics::Phi0

inherited

Definition at line 106 of file Mechanics2D.h.

◆ PhiOld

double VM2D::Mechanics::PhiOld

inherited

Definition at line 118 of file Mechanics2D.h.

◆ Rcm

Point2D VM2D::Mechanics::Rcm

inherited

Текущие положение профиля

Definition at line 112 of file Mechanics2D.h.

◆ Rcm0

Point2D VM2D::Mechanics::Rcm0

inherited

Начальное положение профиля

Definition at line 106 of file Mechanics2D.h.

◆ RcmOld

Point2D VM2D::Mechanics::RcmOld

inherited

Текущие положение профиля

Definition at line 118 of file Mechanics2D.h.

◆ strongCoupling

bool VM2D::MechanicsRigidOscillPart::strongCoupling

private

признак полунеявной схемы связывания

Definition at line 92 of file Mechanics2DRigidOscillPart.h.

◆ Vcm

Point2D VM2D::Mechanics::Vcm

inherited

Текущие скорость центра и угловая скорость

Definition at line 109 of file Mechanics2D.h.

◆ Vcm0

Point2D VM2D::Mechanics::Vcm0

inherited

Начальная скорость центра и угловая скорость

Definition at line 103 of file Mechanics2D.h.

◆ VcmOld

Point2D VM2D::Mechanics::VcmOld

inherited

Скорость и отклонение с предыдущего шага

Definition at line 115 of file Mechanics2D.h.

◆ virtVortParams

const VortexesParams& VM2D::Mechanics::virtVortParams

protectedinherited

Константная ссылка на структуру с параметрами виртуального вихревого слоя для профиля
инициализируется автоматом в конструкторе

Definition at line 95 of file Mechanics2D.h.

◆ viscousForce

Point2D VM2D::Mechanics::viscousForce

inherited

Вектор силы и момент вязкого трения, действующие на профиль

Definition at line 138 of file Mechanics2D.h.

◆ viscousMoment

double VM2D::Mechanics::viscousMoment

inherited

Definition at line 139 of file Mechanics2D.h.

◆ W

const World2D& VM2D::Mechanics::W

protectedinherited

Константная ссылка на решаемую задачу

Definition at line 79 of file Mechanics2D.h.

◆ Wcm

double VM2D::Mechanics::Wcm

inherited

Definition at line 109 of file Mechanics2D.h.

◆ Wcm0

double VM2D::Mechanics::Wcm0

inherited

Definition at line 103 of file Mechanics2D.h.

◆ WcmOld

double VM2D::Mechanics::WcmOld

inherited

Definition at line 115 of file Mechanics2D.h.

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

VM2D/Mechanics2D/Mechanics2DRigidOscillPart.h
VM2D/Mechanics2D/Mechanics2DRigidOscillPart.cpp

Public Member Functions

Public Attributes

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MechanicsRigidOscillPart()

◆ ~MechanicsRigidOscillPart()

Member Function Documentation

◆ AngleOfAirfoil()

◆ AngularVelocityOfAirfoil()

◆ GenerateForcesHeader()

◆ GenerateForcesString()

◆ GeneratePositionHeader()

◆ GeneratePositionString()

◆ GetHydroDynamForce()

◆ getPhi()

◆ getPhiOld()

◆ getR()

◆ getROld()

◆ getStrongCoupling()

◆ getV()

◆ getVOld()

◆ getW()

◆ getWOld()

◆ Initialize()

◆ Move()

◆ MoveKinematic()

◆ MoveOnlyVelo()

◆ PositionOfAirfoilRcm()

◆ ReadParametersFromDictionary()

◆ ReadSpecificParametersFromDictionary()

◆ VeloOfAirfoilPanels()

◆ VeloOfAirfoilRcm()

Member Data Documentation

◆ afl

◆ b

◆ boundary

◆ bw

◆ circulation

◆ circulationOld

◆ hydroDynamForce

◆ hydroDynamMoment

◆ initAngularDisplacement

◆ initAngularVelocity

◆ initDisplacement

◆ initVelocity

◆ isDeform

◆ isMoves

◆ J

◆ k

◆ kw

◆ m

◆ mechParamsParser

◆ numberInPassport

◆ Phi

◆ Phi0

◆ PhiOld

◆ Rcm

◆ Rcm0

◆ RcmOld

◆ strongCoupling

◆ Vcm

◆ Vcm0

◆ VcmOld

◆ virtVortParams

◆ viscousForce

◆ viscousMoment

◆ W

◆ Wcm

◆ Wcm0

◆ WcmOld