Вспомогательный класс Beam для описания упругой хорды деформируемой балки More...

#include <Mechanics2DDeformable.h>

Collaboration diagram for VM2D::Beam:

[legend]

Public Member Functions
double	shape (int n, double x) const

	Beam (const World2D &W_, bool fsi_, double x0_, double L_, int R_)

double	phi (int n, double t) const

void	solveDU (int n, double dt)

void	solveDU_RK (int n, double dt)

double	getTotalDisp (double x, double t) const

double	getGivenLaw (double x, double t, const std::vector< double > &deformParam) const

Public Attributes
const bool	fsi

const World2D &	W

double	x0

double	L

double	rho

double	F

double	EJ

int	R

const std::vector< double >	unitLambda = { 1.87510407, 4.69409113, 7.854757, 10.99554073, 14.1371683 }

std::vector< double >	qCoeff

const double	intSqUnitShape = 0.25

std::vector< double >	currentPhi

std::vector< double >	currentDPhi

const size_t	nLastSteps = 1

std::vector< std::vector< double > >	presLastSteps

Detailed Description

Вспомогательный класс Beam для описания упругой хорды деформируемой балки

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

Version: 1.14

Date: 6 марта 2026 г.

Definition at line 61 of file Mechanics2DDeformable.h.

Constructor & Destructor Documentation

◆ Beam()

Beam::Beam	(	const World2D &	W_,
		bool	fsi_,
		double	x0_,
		double	L_,
		int	R_
	)

Definition at line 56 of file Mechanics2DDeformable.cpp.

                                                                      :
    W(W_),
    fsi(fsi_),
    //filament
    /*
    rho(10000.0),
    F(0.04),
    EJ(1.0 * 1.0),
    */
 
    //Turek
    rho(10000.0),
    F(0.02),
    EJ(1.17032),
 
    R(R_),
    x0(x0_),
    L(L_)
{
    qCoeff.resize(R);
    currentPhi.resize(R);
    currentDPhi.resize(R);
    presLastSteps.reserve(nLastSteps);
};

Member Function Documentation

◆ getGivenLaw()

double Beam::getGivenLaw	(	double	x,
		double	t,
		const std::vector< double > &	deformParam
	)		const

Definition at line 152 of file Mechanics2DDeformable.cpp.

{
    const double c1 = deformParam[3];
    const double c2 = deformParam[4];
    double alpha;
    
    if (fsi) //Turek
        alpha = 0.0; //Turek
    else
        alpha = deformParam[1]; //Fish
    
 
    const double lambda = deformParam[2];
    const double length = 1.0;
    const double f = deformParam[0];
    if (!fsi && (f == 0))
        alpha = 0.0;
 
    auto A = [c1, c2](double xi) {return 1.0 + (xi - 1.0) * c1 + (xi * xi - 1.0) * c2;};
 
    //double result = alpha * A(x + 0.5) * sin(DPI * ((x + 0.5) / (length * lambda) - f * t));
    double result = W.getPassport().physicalProperties.accelCft(t) * alpha * A(x + 0.5) * sin(DPI * ((x + 0.5) / (length * lambda) - f * t));
 
    return result;
}

◆ getTotalDisp()

double Beam::getTotalDisp	(	double	x,
		double	t
	)		const

Definition at line 141 of file Mechanics2DDeformable.cpp.

{
    double result = 0.0;
    for (int i = 0; i < R; ++i)
        result += phi(i, t) * shape(i, x);
    return result;
}

Here is the call graph for this function:

◆ phi()

double Beam::phi	(	int	n,
		double	t
	)		const

Definition at line 81 of file Mechanics2DDeformable.cpp.

{
    return currentPhi[n];
}

Here is the caller graph for this function:

◆ shape()

double VM2D::Beam::shape	(	int	n,
		double	x
	)		const

inline

Definition at line 83 of file Mechanics2DDeformable.h.

        {
            double lam = unitLambda[n] / L;
            double C4 = (sin(unitLambda[n]) - sinh(unitLambda[n])) / (cos(unitLambda[n]) + cosh(unitLambda[n]));
 
            return 0.5 * ( -cos(lam * (x - x0)) + cosh(lam * (x - x0)) - C4 * sin(lam * (x - x0)) + C4 * sinh(lam * (x - x0)) );
        }

Here is the caller graph for this function:

◆ solveDU()

void Beam::solveDU	(	int	n,
		double	dt
	)

Definition at line 86 of file Mechanics2DDeformable.cpp.

{
    double cm = rho * F;
    double ck = EJ * sqr(sqr(unitLambda[n] / L));
    double cq = qCoeff[n];
 
    double omega = sqrt(ck / cm);
    double cDamp = 0.0;
 
    double phiAst = currentPhi[n] + 0.5 * dt * currentDPhi[n];
    double psiAst = currentDPhi[n] + 0.5 * dt * (-ck * currentPhi[n] - cDamp * omega * currentDPhi[n] - cq) / cm;
 
    currentPhi[n] += dt * psiAst;
    currentDPhi[n] += dt * (-ck * phiAst - cDamp * omega * psiAst - cq) / cm;
 
    //std::cout << "n = " << n << std::endl;
    //std::cout << "ck / cm = " << ck / cm << std::endl;
    //std::cout << "cq / cm = " << cq / cm << std::endl;
    //std::cout << "lam = " << unitLambda[n] / L << std::endl;
}

◆ solveDU_RK()

void Beam::solveDU_RK	(	int	n,
		double	dt
	)

Definition at line 107 of file Mechanics2DDeformable.cpp.

                                      {
    double cm = rho * F;
    double ck = EJ * sqr(sqr(unitLambda[n] / L));
    double cq = qCoeff[n];
 
    double omega = sqrt(ck / cm);
    double cDamp = 0.0;
 
    auto acceleration = [&](double phi, double dphi) {
        return (-ck * phi - cDamp * omega * dphi - cq) / cm;
        };
 
    double phi = currentPhi[n];
    double dphi = currentDPhi[n];
 
    double k1_phi = dphi;
    double k1_dphi = acceleration(phi, dphi);
 
    double k2_phi = dphi + 0.5 * dt * k1_dphi;
    double k2_dphi = acceleration(phi + 0.5 * dt * k1_phi, dphi + 0.5 * dt * k1_dphi);
 
    double k3_phi = dphi + 0.5 * dt * k2_dphi;
    double k3_dphi = acceleration(phi + 0.5 * dt * k2_phi, dphi + 0.5 * dt * k2_dphi);
 
    double k4_phi = dphi + dt * k3_dphi;
    double k4_dphi = acceleration(phi + dt * k3_phi, dphi + dt * k3_dphi);
 
    currentPhi[n] = phi + (dt / 6.0) * (k1_phi + 2 * k2_phi + 2 * k3_phi + k4_phi);
    currentDPhi[n] = dphi + (dt / 6.0) * (k1_dphi + 2 * k2_dphi + 2 * k3_dphi + k4_dphi);
    
    std::cout << "phi = " << currentPhi[n] << std::endl;
}

Here is the call graph for this function:

Member Data Documentation

◆ currentDPhi

std::vector<double> VM2D::Beam::currentDPhi

Definition at line 91 of file Mechanics2DDeformable.h.

◆ currentPhi

std::vector<double> VM2D::Beam::currentPhi

Definition at line 91 of file Mechanics2DDeformable.h.

◆ EJ

double VM2D::Beam::EJ

Definition at line 71 of file Mechanics2DDeformable.h.

◆ F

double VM2D::Beam::F

Definition at line 71 of file Mechanics2DDeformable.h.

◆ fsi

const bool VM2D::Beam::fsi

Definition at line 64 of file Mechanics2DDeformable.h.

◆ intSqUnitShape

const double VM2D::Beam::intSqUnitShape = 0.25

Definition at line 80 of file Mechanics2DDeformable.h.

◆ L

double VM2D::Beam::L

Definition at line 69 of file Mechanics2DDeformable.h.

◆ nLastSteps

const size_t VM2D::Beam::nLastSteps = 1

Definition at line 93 of file Mechanics2DDeformable.h.

◆ presLastSteps

std::vector<std::vector<double> > VM2D::Beam::presLastSteps

Definition at line 94 of file Mechanics2DDeformable.h.

◆ qCoeff

std::vector<double> VM2D::Beam::qCoeff

Definition at line 77 of file Mechanics2DDeformable.h.

◆ R

int VM2D::Beam::R

Definition at line 73 of file Mechanics2DDeformable.h.

◆ rho

double VM2D::Beam::rho

Definition at line 71 of file Mechanics2DDeformable.h.

◆ unitLambda

const std::vector<double> VM2D::Beam::unitLambda = { 1.87510407, 4.69409113, 7.854757, 10.99554073, 14.1371683 }

Definition at line 76 of file Mechanics2DDeformable.h.

76{ 1.87510407, 4.69409113, 7.854757, 10.99554073, 14.1371683 };

◆ W

const World2D& VM2D::Beam::W

Definition at line 66 of file Mechanics2DDeformable.h.

◆ x0

double VM2D::Beam::x0

Definition at line 68 of file Mechanics2DDeformable.h.

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

VM2D/Mechanics2D/Mechanics2DDeformable.h
VM2D/Mechanics2D/Mechanics2DDeformable.cpp

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Beam()

Member Function Documentation

◆ getGivenLaw()

◆ getTotalDisp()

◆ phi()

◆ shape()

◆ solveDU()

◆ solveDU_RK()

Member Data Documentation

◆ currentDPhi

◆ currentPhi

◆ EJ

◆ F

◆ fsi

◆ intSqUnitShape

◆ L

◆ nLastSteps

◆ presLastSteps

◆ qCoeff

◆ R

◆ rho

◆ unitLambda

◆ W

◆ x0