#include "myassert.h"
#include "except.h"

Include dependency graph for shapefnc.h:

This graph shows which files directly or indirectly include this file:

Enumerations
enum	Order { ORD_ALG = 0, ORD_D1 = 1, ORD_D2 = 2 }

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

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

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

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

Variables
const doublereal	dN3 [2][3]

const doublereal	dN3P [2][3]

const doublereal	dN3PP [2][3]

const doublereal	dN2 [2]

const doublereal	dN2P [2]

const doublereal	dN2PP [2]

Enumeration Type Documentation

enum Order

Enumerator
ORD_ALG
ORD_D1
ORD_D2

Definition at line 42 of file shapefnc.h.

42 { ORD_ALG = 0, ORD_D1 = 1, ORD_D2 = 2 };

ORD_D1

Definition: shapefnc.h:42

ORD_D2

Definition: shapefnc.h:42

ORD_ALG

Definition: shapefnc.h:42

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 = `ORD_ALG`
	)

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 = `ORD_ALG`
	)

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.;
 }