#include "mbconfig.h"
#include <limits>
#include <cfloat>
#include "matvec3.h"
#include "shapefnc.h"

Include dependency graph for shapefnc.cc:

Functions
doublereal	ShapeFunc2N (doublereal d, integer iNode, enum Order Ord)

doublereal	DxDcsi2N (doublereal d, const Vec3 &X1, const Vec3 &X2)

doublereal	ShapeFunc3N (doublereal d, integer iNode, enum Order Ord)

doublereal	DxDcsi3N (doublereal d, const Vec3 &X1, const Vec3 &X2, const Vec3 &X3)

Variables
const doublereal	dN2_1 = .5

const doublereal	dN2_2 = .5

const doublereal	dN2P_1 = -.5

const doublereal	dN2P_2 = .5

const doublereal	dN2PP_1 = 0.

const doublereal	dN2PP_2 = 0.

const doublereal	dN2 [2]

const doublereal	dN2P [2]

const doublereal	dN2PP [2]

const doublereal	dS = 1./sqrt(3.)

const doublereal	dN1_I = (1.+sqrt(3.))/6.

const doublereal	dN2_I = 2./3.

const doublereal	dN3_I = (1.-sqrt(3.))/6.

const doublereal	dN1P_I = -(2.*sqrt(3.)+3.)/6.

const doublereal	dN2P_I = 2./sqrt(3.)

const doublereal	dN3P_I = -(2.*sqrt(3.)-3.)/6.

const doublereal	dN1PP_I = 1.

const doublereal	dN2PP_I = -2.

const doublereal	dN3PP_I = 1.

const doublereal	dN1II = dN3_I

const doublereal	dN2II = dN2_I

const doublereal	dN3II = dN1_I

const doublereal	dN1PII = -dN3P_I

const doublereal	dN2PII = -dN2P_I

const doublereal	dN3PII = -dN1P_I

const doublereal	dN1PPII = dN3PP_I

const doublereal	dN2PPII = dN2PP_I

const doublereal	dN3PPII = dN1PP_I

const doublereal	dN3 [2][3]

const doublereal	dN3P [2][3]

const doublereal	dN3PP [2][3]

Function Documentation

doublereal DxDcsi2N	(	doublereal	d,
		const Vec3 &	X1,
		const Vec3 &	X2
	)

Definition at line 113 of file shapefnc.cc.

References Vec3::Dot(), ORD_D1, ShapeFunc2N(), and grad::sqrt().

Referenced by AerodynamicBeam2::AssJac(), and AerodynamicBeam2::AssVec().

 {
         doublereal dN1p = ShapeFunc2N(d, 1, ORD_D1);
         doublereal dN2p = ShapeFunc2N(d, 2, ORD_D1);
         Vec3 DXDcsi(X1*dN1p+X2*dN2p);
         doublereal dd = DXDcsi.Dot();
 
         if (dd > std::numeric_limits<doublereal>::epsilon()) {
                 return std::sqrt(dd);
         }
 
         return 0.;
 }

Here is the call graph for this function:

doublereal DxDcsi3N	(	doublereal	d,
		const Vec3 &	X1,
		const Vec3 &	X2,
		const Vec3 &	X3
	)

Definition at line 228 of file shapefnc.cc.

References Vec3::Dot(), ORD_D1, ShapeFunc3N(), and grad::sqrt().

Referenced by AerodynamicBeam::AssJac(), and AerodynamicBeam::AssVec().

 {
         doublereal dN1p = ShapeFunc3N(d, 1, ORD_D1);
         doublereal dN2p = ShapeFunc3N(d, 2, ORD_D1);
         doublereal dN3p = ShapeFunc3N(d, 3, ORD_D1);
         Vec3 DXDcsi(X1*dN1p+X2*dN2p+X3*dN3p);
         doublereal dd = DXDcsi.Dot();
 
         if (dd > std::numeric_limits<doublereal>::epsilon()) {
                 return std::sqrt(dd);
         }
 
         return 0.;
 }

Here is the call graph for this function:

doublereal ShapeFunc2N	(	doublereal	d,
		integer	iNode,
		enum Order	Ord
	)

Definition at line 66 of file shapefnc.cc.

References ASSERT, MBDYN_EXCEPT_ARGS, ORD_ALG, ORD_D1, and ORD_D2.

Referenced by AerodynamicBeam2::AssJac(), AerodynamicBeam2::AssVec(), and DxDcsi2N().

 {
         ASSERT(iNode == 1 || iNode == 2);
         
         switch (Ord) {
         case ORD_ALG:
                 switch (iNode) {
                 case 1:
                         return .5*(1.-d);
                         
                 case 2:
                         return .5*(1.+d);
                         
                 default:
                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                 }
                 
         case ORD_D1:
                 switch (iNode) {
                 case 1:         
                         return -.5;
                         
                 case 2:
                         return .5;
                         
                 default:
                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                 }
                 
         case ORD_D2:
                 switch (iNode) {
                 case 1:         
                         return 0.;
                         
                 case 2:
                         return 0.;
                         
                 default:
                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                 }
         }
         
         /* Per evitare warnings */
         return 0.;
 }

doublereal ShapeFunc3N	(	doublereal	d,
		integer	iNode,
		enum Order	Ord
	)

Definition at line 173 of file shapefnc.cc.

References ASSERT, MBDYN_EXCEPT_ARGS, ORD_ALG, ORD_D1, and ORD_D2.

Referenced by AerodynamicBeam::AssJac(), BeamSliderJoint::AssRes(), AerodynamicBeam::AssVec(), and DxDcsi3N().

 {
         ASSERT(iNode == 1 || iNode == 2 || iNode == 3);
         
         switch (Ord) {
         case ORD_ALG:
                 switch (iNode) {
                 case 1:         
                         return .5*d*(d-1.);
                         
                 case 2:         
                         return 1.-d*d;
                         
                 case 3:         
                         return .5*d*(d+1.);
                         
                 default:
                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                 }
                 
         case ORD_D1:
                 switch (iNode) {
                 case 1:         
                         return d-.5;
                         
                 case 2:
                         return -2.*d;
                         
                 case 3:
                         return d+.5;
                         
                 default:
                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                 }
         case ORD_D2:
                 switch (iNode) {
                 case 1:         
                         return 1.;
                         
                 case 2:
                         return -2.;
                         
                 case 3:
                         return 1.;
                         
                 default:
                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                 }
         }
         
         /* Per evitare warnings */
         return 0.;
 }

Variable Documentation

const doublereal dN1_I = (1.+sqrt(3.))/6.

Definition at line 131 of file shapefnc.cc.

const doublereal dN1II = dN3_I

Definition at line 144 of file shapefnc.cc.

const doublereal dN1P_I = -(2.*sqrt(3.)+3.)/6.

Definition at line 135 of file shapefnc.cc.

const doublereal dN1PII = -dN3P_I

Definition at line 148 of file shapefnc.cc.

const doublereal dN1PP_I = 1.

Definition at line 139 of file shapefnc.cc.

const doublereal dN1PPII = dN3PP_I

Definition at line 152 of file shapefnc.cc.

const doublereal dN2[2]

Initial value:

= {
        dN2_1, dN2_2
}

Definition at line 52 of file shapefnc.cc.

Referenced by AerodynamicBeam::AssJac(), AerodynamicBeam2::AssJac(), Beam2::AssStiffnessMat(), ViscoElasticBeam2::AssStiffnessMat(), AerodynamicBeam::AssVec(), AerodynamicBeam2::AssVec(), Beam2::InterpState(), HBeam::Omega0(), and Beam2::Omega0().

const doublereal dN2_1 = .5

Definition at line 42 of file shapefnc.cc.

const doublereal dN2_2 = .5

Definition at line 43 of file shapefnc.cc.

const doublereal dN2_I = 2./3.

Definition at line 132 of file shapefnc.cc.

const doublereal dN2II = dN2_I

Definition at line 145 of file shapefnc.cc.

const doublereal dN2P[2]

Initial value:

= {
        dN2P_1, dN2P_2
}

Definition at line 56 of file shapefnc.cc.

Referenced by Beam2::AssStiffnessMat(), ViscoElasticBeam2::AssStiffnessMat(), and Beam2::InterpDeriv().

const doublereal dN2P_1 = -.5

Definition at line 45 of file shapefnc.cc.

const doublereal dN2P_2 = .5

Definition at line 46 of file shapefnc.cc.

const doublereal dN2P_I = 2./sqrt(3.)

Definition at line 136 of file shapefnc.cc.

const doublereal dN2PII = -dN2P_I

Definition at line 149 of file shapefnc.cc.

const doublereal dN2PP[2]

Initial value:

= {
        dN2PP_1, dN2PP_2
}

Definition at line 60 of file shapefnc.cc.

const doublereal dN2PP_1 = 0.

Definition at line 48 of file shapefnc.cc.

const doublereal dN2PP_2 = 0.

Definition at line 49 of file shapefnc.cc.

const doublereal dN2PP_I = -2.

Definition at line 140 of file shapefnc.cc.

const doublereal dN2PPII = dN2PP_I

Definition at line 153 of file shapefnc.cc.

const doublereal dN3[2][3]

Initial value:

= {
        { dN1_I, dN2_I, dN3_I },
        { dN1II, dN2II, dN3II }
}

Definition at line 157 of file shapefnc.cc.

Referenced by AerodynamicBeam::AssJac(), Beam::AssStiffnessMat(), ViscoElasticBeam::AssStiffnessMat(), AerodynamicBeam::AssVec(), Beam::InterpState(), and Beam::Omega0().

const doublereal dN3_I = (1.-sqrt(3.))/6.

Definition at line 133 of file shapefnc.cc.

const doublereal dN3II = dN1_I

Definition at line 146 of file shapefnc.cc.

const doublereal dN3P[2][3]

Initial value:

= {
        { dN1P_I, dN2P_I, dN3P_I },
        { dN1PII, dN2PII, dN3PII }
}

Definition at line 162 of file shapefnc.cc.

Referenced by Beam::AssStiffnessMat(), ViscoElasticBeam::AssStiffnessMat(), and Beam::InterpDeriv().

const doublereal dN3P_I = -(2.*sqrt(3.)-3.)/6.

Definition at line 137 of file shapefnc.cc.

const doublereal dN3PII = -dN1P_I

Definition at line 150 of file shapefnc.cc.

const doublereal dN3PP[2][3]

Initial value:

= {
        { dN1PP_I, dN2PP_I, dN3PP_I },
        { dN1PPII, dN2PPII, dN3PPII }
}

Definition at line 167 of file shapefnc.cc.

const doublereal dN3PP_I = 1.

Definition at line 141 of file shapefnc.cc.

const doublereal dN3PPII = dN1PP_I

Definition at line 154 of file shapefnc.cc.

const doublereal dS = 1./sqrt(3.)

Definition at line 128 of file shapefnc.cc.

Functions

Variables

Function Documentation

Variable Documentation