FNetElement.h

 
#ifndef F_NETELEMENT_H
#define F_NETELEMENT_H
#include <cmath>
#include "sfh/constants.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.
 
    double 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;
    double avgMeshesBetweenNodes = (sfh::math::Norm(m_nodesABdist_mesh, 3) + sfh::math::Norm(m_nodesACdist_mesh, 3) + sfh::math::Norm(m_nodesBCdist_mesh, 3)) / 3;
    double 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 netconstants
{
    const double dHydroForcesRatio      = 1;
    const double dAddedLinearDrag       = 0;
    const double knotDiameterRatio      = 2.0;
    const double m_ny_water             = 1.19e-6;
    const double rho_water              = 1025;
    const double m_knotMomStiff         = 0;// 0.001;
    const double CnLinearLimitAngle     = 4.0 * sfh::pi / 180;
    const double CnLinear               = 0.068 ;
    const double CnTurb                 = 0.8;
    const double CnLam                  = 1.15;
    const double TurbLimit              = 3.4e5;
    const double LamLimit               = 2.0e5;
    const double Ct_nominal             = 0.01;
    const double CnKnots_nominal        = 1.15;
};
 
void AddNetPanelNodeForces(
        int             iElement,
        const double*   adPosAndVel,
        const int*      aiPosIndex,
        const int*      aiVelIndex,
        const double*   adWaterVel_ned,
        double*         adNodeForces,
        int*            aiNetConn,
        double*         adNodesDistMesh,
        double*         adMesdDet,
        double*         adEATwine,
        double*         adMeshBarLength_unstretched,
        double*         adKnotDiameter,
        double*         adKnotMomStiff_contact,
        double*         adNodeWeight,
        double*         adDampingCoeff,
        double*         adKKnot,
        double*         adTwineDragNormal,
        double*         adTwineDiameter,
        const double    m_CnLinearLimitAngle,
        const double    m_CnLinear,
        const double    m_CnTurb,
        const double    m_CnLam,
        const double    m_TurbLimit,
        const double    m_LamLimit,
        const double    m_Ct_nominal,
        const double    m_CnKnots_nominal);
 
#endif