VM2D::Beam Class Reference

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

#include <Mechanics2DDeformable.h>

Collaboration diagram for VM2D::Beam:

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, 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 55 of file Mechanics2DDeformable.cpp.

                                                                      :
    W(W_),
    fsi(fsi_),
    rho(10000.0),
    F(0.02),
    EJ(1.1703239289446188),
    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,
		double	deformParam
	)		const

Definition at line 137 of file Mechanics2DDeformable.cpp.

{
    const double c1 = -0.825;
    const double c2 = 1.625;
    double alpha;
    if (fsi) //Turek
        alpha = 0.0; //Turek
    else
        alpha = 0.1; //Fish
    
 
    const double lambda = 1.0;
    const double length = 1.0;
    const double f = deformParam;
 
    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 129 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 71 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 76 of file Mechanics2DDeformable.cpp.

{
    double cm = rho * F;
    double ck = EJ * sqr(sqr(unitLambda[n] / L));
    double cq = qCoeff[n];
    //std::cout << "ck = " << ck << std::endl;
    //std::cout << "lam = " << unitLambda[n] << std::endl;
    //std::cout << "lam = " << unitLambda[n] / L << std::endl;
 
    //double omega = sqrt(ck / cm);
    //double cDamp = 0.0; //0.005;
 
    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;
}

solveDU_RK()

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

Definition at line 95 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.

{ 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