FNetElementFloat.h

 
#ifndef F_NETELEMENT_FLOAT_H
#define F_NETELEMENT_FLOAT_H
#include <cmath>
#include "sfh/constants.h"
#include "sfh/math/array.h"
#include "sfh/util/diagnostic.h"
#include "sfh/sim/kinematics.h"
 
 
/*
 * /// <li> Calculates cross-sectional area of twine; all twines in a net element are equal.
    m_twineCrossSecArea         = sfh::pi*twineDiameter*twineDiameter/4;
 
    for (int i = 0; i<2; i++)
    {
        m_nodesABdist_mesh[i] = m_nodeBPos_mesh[i] - m_nodeAPos_mesh[i];
        m_nodesACdist_mesh[i] = m_nodeCPos_mesh[i] - m_nodeAPos_mesh[i];
        m_nodesBCdist_mesh[i] = m_nodeCPos_mesh[i] - m_nodeBPos_mesh[i];
    }
 
    m_meshDet = sfh::math::Abs(m_nodesACdist_mesh[0]*m_nodesABdist_mesh[1] - m_nodesABdist_mesh[0]*m_nodesACdist_mesh[1]);
    if (sfh::math::Abs(m_meshDet) < 1.0)
        CPrintDuringExec::WriteLog("In CNetElement3N::Initialize() for net "
        + m_sNetName
        + " and section with id "
        + sfh::text::Int2Str(m_iID)
        + ". The mesh determinant in  is less than 1.0. This may be caused by two nodes of a net section being too close. Check your input!", "TempOut.txt");
 
    m_numMeshes     = sfh::math::Abs(m_meshDet/2);  // Number of meshes in the triangular element.
    m_numBars       = sfh::math::Abs(m_meshDet);    // Number of bars equals twice the number of meshes (Priour, 2005, Eq. 7).
    m_numUBars      = sfh::math::Abs(m_meshDet/2);  // Number of u-bars is half the total number of mesh bars.
    m_numVBars      = sfh::math::Abs(m_meshDet/2);  // Number of v-bars is half the total number of mrsh bars.
    m_numKnots      = sfh::math::Abs(m_meshDet/2);  // Number of knots equals the number of meshes;
 
    m_weight = (m_rho_twines - m_rho_water)*m_twineCrossSecArea*m_meshBarLength_unstretched*m_numBars;
    m_NodeWeight = m_weight / 3;
 
    m_Ct = m_Ct_nominal*sfh::pi; // Approximate tangential drag force coefficient, nominal coefficient multiplied by PI due to the typical expression for the drag force.
 
    float CnKnots = m_CnKnots_nominal * sfh::pi * m_knotDiameter * m_knotDiameter / 4;
    m_kKnot = 0.5 * m_rho_water * m_knotDiameter * m_numKnots * CnKnots;
 
    m_kTwineDragNormal[0] = 0.5 * m_rho_water * m_twineDiameter * (m_meshBarLength_unstretched - m_knotDiameter) * m_numUBars ;
    m_kTwineDragNormal[1] = 0.5 * m_rho_water * m_twineDiameter * (m_meshBarLength_unstretched - m_knotDiameter) * m_numVBars ;
 
    m_EATwine = m_EModulus * m_twineCrossSecArea;
    m_dDampingRatio = 10.0;
    float avgMeshesBetweenNodes = (sfh::math::Norm(m_nodesABdist_mesh, 3) + sfh::math::Norm(m_nodesACdist_mesh, 3) + sfh::math::Norm(m_nodesBCdist_mesh, 3)) / 3;
    float avgNodeDist = avgMeshesBetweenNodes * m_meshBarLength_unstretched;
    m_dDampingCoeff  = m_dDampingRatio * sqrt(m_NodeWeight * m_EATwine * m_numBars / avgMeshesBetweenNodes / avgNodeDist) ;
 *
 */
 
// Sets some constants common for all net panels.
namespace netconstantsfloat
{
    const float dHydroForcesRatio       = 1;
    const float dAddedLinearDrag        = 0;
    const float knotDiameterRatio       = 2.0;
    const float m_ny_water              = 1.19e-6;
    const float rho_water               = 1025;
    const float m_knotMomStiff          = 0;// 0.001;
    const float CnLinearLimitAngle      = 4.0 * sfh::pi / 180;
    const float CnLinear                = 0.068 ;
    const float CnTurb                  = 0.8;
    const float CnLam                   = 1.15;
    const float TurbLimit               = 3.4e5;
    const float LamLimit                = 2.0e5;
    const float Ct_nominal              = 0.01;
    const float CnKnots_nominal     = 1.15;
};
 
void AddNetPanelNodeForces(
        int             iElement,
        const float*    adPosAndVel,
        const int*      aiPosIndex,
        const int*      aiVelIndex,
        const float*    adWaterVel_ned,
        float*      adNodeForces,
        int*            aiNetConn,
        float*      adNodesDistMesh,
        float*      adMesdDet,
        float*      adEATwine,
        float*      adMeshBarLength_unstretched,
        float*      adKnotDiameter,
        float*          adKnotMomStiff_contact,
        float*      adNodeWeight,
        float*          adDampingCoeff,
        float*      adKKnot,
        float*      adTwineDragNormal,
        float*      adTwineDiameter,
        const float m_CnLinearLimitAngle,
        const float m_CnLinear,
        const float m_CnTurb,
        const float m_CnLam,
        const float m_TurbLimit,
        const float m_LamLimit,
        const float m_Ct_nominal,
        const float m_CnKnots_nominal);
 
#endif