rotorcraft dynamics

127
Daniel Guggenheim School of Aerospace Engineering Introduction Kinematical Preliminaries Analysis of 3-D Beam Deformation Stiffness Modeling Geometrically Exact Beam Equations Rotorcraft Dynamics Dewey H. Hodges 1 1 Professor, Daniel Guggenheim School of Aerospace Engineering Georgia Institute of Technology, Atlanta, Georgia AE 6220, Spring 2013 Hodges Composite Rotor Blade Modeling

Upload: others

Post on 05-Dec-2021

26 views

Category:

Documents


2 download

TRANSCRIPT

Page 1: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Rotorcraft Dynamics

Dewey H. Hodges1

1Professor, Daniel Guggenheim School of Aerospace EngineeringGeorgia Institute of Technology, Atlanta, Georgia

AE 6220, Spring 2013

Hodges Composite Rotor Blade Modeling

Page 2: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

BackgroundOutline

By now most rotorcraft dynamicists are accustomed towriting lengthy, complicated equations for analysis of bladedynamicsMany extant sets of such equations contain a host ofoversimplifications

no initial curvatureuniaxial stress field (or alternatively cross-section rigid in itsown plane)only torsional warpingā€œmoderateā€ deflections (or alternatively some sort ofordering scheme)isotropic or transversely isotropic material constructionno shear deformation (invalidates theory for application tocomposite beams)

Hodges Composite Rotor Blade Modeling

Page 3: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

BackgroundOutline

Approximation concepts such as ā€œordering schemesā€ havebeen deeply ingrained in the thinking of most analystsHowever, some research has pointed out problems withthis concept:

It is virtually impossible to apply an ordering scheme in acompletely consistent manner (Stephens, Hodges, Avila,and Kung 1982)Ordering schemes can lead to more lengthy equations dueto expansion of transcendental functions (Crespo da Silvaand Hodges 1986)An ordering scheme that works for one set of configurationparameters may not be suitable for a different set (Hinnantand Hodges 1989)

Hodges Composite Rotor Blade Modeling

Page 4: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

BackgroundOutline

In the present approach, there are several basic and excitingdepartures from the ā€œold schoolā€:

No ordering scheme is needed or usedexact kinematics for beam reference line displacement andcross-sectional rotationgeometrically-exact equations of motion

The beam constitutive law isbased on a separate finite element analysisvalid for anisotropic beams with inhomogeneouscross-sections

Hodges Composite Rotor Blade Modeling

Page 5: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

BackgroundOutline

. . . departures from the ā€œold schoolā€ (continued):A compact matrix notation is usedThe resulting mixed formulation can be put in the weakestform, so

the requirements for the shape functions are minimal,leading to the possibility of shape functions that are assimple as piecewise constantapproximate element quadrature is not required

Hodges Composite Rotor Blade Modeling

Page 6: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

BackgroundOutline

Kinematical PreliminariesAnalysis of 3-D Beam DeformationConstitutive Equations from 2-D Finite Element Analysis1-D Kinematics1-D Equations of Motion1-D Finite Element SolutionExamplesIt should be noted that this material is found in the authorā€™sbook Nonlinear Composite Beam Theory, published byAIAA in 2006.

Hodges Composite Rotor Blade Modeling

Page 7: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

Consider a rigid body B moving in a frame A (frames andrigid bodies are kinematically equivalent)Introduce a dextral unit triad Ai fixed in A (Romansubscripts vary from 1 to 3 unless otherwise specified, andunit vectors are denoted by a bold italic symbol with a ā€œhatā€)Also, introduce a dextral triad Bi fixed in BNow Bi will vary in A as a function of timeA vector is a first-order tensorA vector can always be expressed as a linear combinationof dextral unit vectorsFor example, for an arbitrary vector v with vBi = v Ā· Bi , it isalways true that v = BivBi (note that summation is impliedover any repeated index)

Hodges Composite Rotor Blade Modeling

Page 8: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

Similarly, a dyadic is a second-order Cartesian tensorDyadics are quadratic forms of dextral unit vectorsFor example, consider the relationship of a dyadic T andthe matrix Tij of its components in a mixed set of bases

T = BiTij Aj

The transpose of T is simply

T T = AjTij Bi = AiTji Bj

For simplicity we will not carry names of associated basevectors in the symbol for a particular matrix of dyadiccomponents

Hodges Composite Rotor Blade Modeling

Page 9: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

An exception is the finite rotation tensor for which it is quitehelpful to maintain basis association in naming matricesOne can characterize the rotational motion of B in A inthese two ways:

Bi = CBA Ā· Ai = CBAij Aj

CBA is read as the finite rotation tensor of B in A, given by

CBA = Bi Ai

Both the tensor and its components depend on timeNote the convention for naming the finite rotation tensorand its corresponding matrix of direction cosines

Hodges Composite Rotor Blade Modeling

Page 10: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

The transpose is indicated by reversing the superscripts(CBA

)T= CAB

so thatAi = CAB Ā· Bi

CBA is an orthonormal tensor so that

CBA Ā· CAB = āˆ†

where āˆ† is the identity tensor

Hodges Composite Rotor Blade Modeling

Page 11: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

The matrix of direction cosines CBA is given by

CBAij = Bi Ā· Aj

(= CAB

ji

)CBA is a matrix of the components of the transpose of thefinite rotation tensor CAB so that

CBAij =Ai Ā· CAB Ā· Aj

=Bi Ā· CAB Ā· Bj

The matrix of direction cosines is also orthonormal

CBACAB = āˆ†

where āˆ† is the 3Ɨ 3 identity matrix

Hodges Composite Rotor Blade Modeling

Page 12: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

When an intermediate frame N is involved,

CBA = CBNCNA CBA = CBN Ā· CNA

Consider an arbitrary vector v ; it is always possible to writevZi = v Ā· Z i where Z is an arbitrary frame in which thedextral unit triad Z i is fixedWith the column matrix notation

vZ =

vZ1vZ2vZ3

it is easily demonstrated that

vB = CBAvA vA = CABvB

Hodges Composite Rotor Blade Modeling

Page 13: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

For a rigid body B moving in frame A, there existanalogous vector-dyadic operations

the ā€œpush-forwardā€ operation on v is defined by

CBA Ā· v

the ā€œpull-backā€ operation on v is defined by

CAB Ā· v

As can be demonstrated, these operations rotate thevector by an amount commensurate with the change inorientation from A to B and from B to A respectively.

Hodges Composite Rotor Blade Modeling

Page 14: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

To visualize the pull-back operationimagine the vector v frozen at some instant in time in aframe N which has a dextral triad N i that is coincident withBirotate the frame N so that N i lines up with Aithe rotated image of v is the result of the pull-backoperation

For the push-forward operationimagine the vector v frozen at some instant in time in aframe N which has a dextral triad N i that is coincident withAirotate the frame N so that N i lines up with Bithe rotated image of v is the result of the push-forwardoperation

These operations are useful for describing deformation

Hodges Composite Rotor Blade Modeling

Page 15: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

Note that if a vector is moving in a frame then its timederivative in that frame will be nonzeroHowever, if a vector is fixed in a frame then its timederivative in that frame will be zeroObviously, the concept of the derivative of a vector is thenframe dependentConsider a vector b fixed in B where B is moving in A.Then,

Adbdt6= 0

ButBdbdt

= 0

Hodges Composite Rotor Blade Modeling

Page 16: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

Thus, when the vector v is resolved along unit vectors thatare fixed in the frame in which the derivative is being taken,the derivative of v is easily expressed as

Z dvdt

= Z i vZi

where Z is an arbitrary frame as beforeHowever, as will be clear from material to follow, it is notalways convenient to differentiate a vector in this way

Hodges Composite Rotor Blade Modeling

Page 17: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

For differentiation of a vector in a different frame

Advdt

=Bdvdt

+ Ļ‰BA Ɨ v

where Ļ‰BA is the angular velocity of B in A given by

Ļ‰BA = B1

AdB2

dtĀ· B3 + B2

AdB3

dtĀ· B1 + B3

AdB1

dtĀ· B2 = Ļ‰BA

Bi Bi

The derivatives of the unit vectors are easily expressed interms of the direction cosines

AdBi

dt= CBA

ij Aj = CBAij CAB

jk Bk = āˆ’eijkĻ‰BABj Bk = Ļ‰BA Ɨ Bi

Hodges Composite Rotor Blade Modeling

Page 18: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

Just as the time derivative of a vector depends on theframe in which the derivative is taken, so does the variationof a vectorAs Kane and others have shown, one can express therelationship between variations in two frames as

AĪ“v = BĪ“v + Ī“ĻˆBA Ɨ v

where Ī“ĻˆBA is the virtual rotation of B in A given by

Ī“ĻˆBA

= B1AĪ“B2 Ā· B3 + B2

AĪ“B3 Ā· B1 + B3AĪ“B1 Ā· B2 = Ī“Ļˆ

BABi Bi

andZ Ī“v = Z iĪ“vZi

Hodges Composite Rotor Blade Modeling

Page 19: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

The variations of the unit vectors are easily expressed interms of the direction cosines

AĪ“Bi = Ī“CBAij Aj = Ī“CBA

ij CABjk Bk = āˆ’eijkĪ“Ļˆ

BABj Bk = Ī“Ļˆ

BAƗ Bi

where Ī“( ) is the usual Lagrangean variation

It is evident that the vectors Ļ‰BA and Ī“ĻˆBA can beregarded as operators which produce the time derivativeand variation, respectively, in A of any vector fixed in B

When an additional frame N is involved, Kaneā€™s additiontheorem applies

Ļ‰BA = Ļ‰BN + Ļ‰NA Ī“ĻˆBA

= Ī“ĻˆBN

+ Ī“ĻˆNA

Hodges Composite Rotor Blade Modeling

Page 20: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

Note that to obtain the virtual rotation vector, one needonly replace the dots in the angular velocity vector with Ī“ā€™s,ignoring any other termsThe virtual work in A of an applied torque T acting on abody B is simply

Ī“W = T Ā· Ī“ĻˆBA

Hodges Composite Rotor Blade Modeling

Page 21: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

The velocity of a point P moving in A can be determined bytime differentiation in A of the position vector pP/O whereO is any point fixed in A

vPA =AdpP/O

dt

Often the calculation of velocity in this way is complicatedIn such a case it is helpful to ā€œstepā€ oneā€™s way from theknown to the unknown using the two chain rules:

2 points fixed on a rigid body (or in a frame)1 point moving on a rigid body (or in a frame)

Hodges Composite Rotor Blade Modeling

Page 22: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

For 2 points P and Q fixed on a rigid body B having anangular velocity Ļ‰BA in A, the velocities of these points in Aare related according to

vPA = vQA + Ļ‰BA Ɨ pP/Q

For a point P moving on a rigid body B while B is moving inA, the velocity of P in A is given by

vPA = vPB + v BA

where v BA is the velocity of the point in B that is coincidentwith P at the instant under consideration (this can often beobtained by use of the other theorem)

Hodges Composite Rotor Blade Modeling

Page 23: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

These relationships for the derivative and variation can benicely expressed in matrix notationWeā€™ve already seen how an arbitrary vector v can for anarbitrary frame Z be expressed in terms of

vZ =

vZ1vZ2vZ3

The dual matrix vZ ij = āˆ’eijkvZk has the same measurenumbers but arranged antisymmetrically

vZ =

0 āˆ’vZ3 vZ2vZ3 0 āˆ’vZ1āˆ’vZ2 vZ1 0

Hodges Composite Rotor Blade Modeling

Page 24: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

With the tilde notation identities given in Hodges (1990) arehelpfulWhen Y and Z are 3Ɨ 1 column matrices, it is easilyshown that

(Z )T =āˆ’ Z

ZZ =0

YZ =āˆ’ ZY

Y T Z =āˆ’ Z T Y

Y Z =ZY T āˆ’āˆ†Y T Z

Y Z =Z Y +ĖœYZ

where āˆ† is the 3Ɨ 3 identity matrix

Hodges Composite Rotor Blade Modeling

Page 25: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

This notation also applies to vectorsThe tensor v has components in the Z basis given by thematrix vZ such that

v =v Ɨāˆ†

=Z i vZ ij Z j

where āˆ† is the identity dyadicNote that for any vector w , v Ɨw = v Ā·wNote that for any vectors v and w with measure numbersexpressed in some common basis Z, vZ wZ contains themeasure numbers of the cross product v Ɨw in the Z basisThe ( ) operator is sometimes called a ā€œcross productoperatorā€ for obvious reasons

Hodges Composite Rotor Blade Modeling

Page 26: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

Now with these definitions in mind we note that

Ļ‰BAB = āˆ’CBACAB

andĖœĪ“Ļˆ

BAB = āˆ’Ī“CBACAB

In tensorial form these are

Ļ‰BA = ACBA Ā· CAB

andĖœĪ“Ļˆ

BA= AĪ“CBA Ā· CAB

Hodges Composite Rotor Blade Modeling

Page 27: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

Since CBA is orthonormal, its elements are notindependentEulerā€™s theorem of rotation stipulates that any change oforientation can be characterized as a ā€œsimple rotationā€A motion is a ā€œsimple rotationā€ of B in A if during the motiona line L maintains its orientation in B and in AEuler proved that at least four parameters are necessaryfor singularity free description of finite rotation:

the three measure numbers of a unit vector e along the lineL (ei = eAi = eBi ) andthe magnitude of the rotation Ī±

With e = be1 e2 e3cT , the matrix of direction cosines is

CBA = āˆ† cosĪ± + eeT (1āˆ’ cosĪ±)āˆ’ e sinĪ±

Hodges Composite Rotor Blade Modeling

Page 28: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

According to Kane, Likins, and Levinson (1983), the fourEuler parameters denoted by Īµ0 and

Īµ =

Īµ1Īµ2Īµ3

are defined as

Īµi = ei sin(Ī±/2) Īµ0 = cos(Ī±/2)

They satisfy a constraint

ĪµT Īµ+ Īµ20 = 1

Hodges Composite Rotor Blade Modeling

Page 29: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

The matrix of direction cosines is given by

CBA =(

1āˆ’ 2ĪµT Īµ)

āˆ† + 2(ĪµĪµT āˆ’ Īµ0Īµ

)while the angular velocity is

Ļ‰BAB = 2 [(Īµ0āˆ†āˆ’ Īµ) Īµāˆ’ Īµ0Īµ]

Hodges Composite Rotor Blade Modeling

Page 30: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

It is tempting to eliminate the fourth parameter via theconstraintThis always introduces a singularity, but where thesingularity is can be controlledIntroduce Rodrigues parameters

Īø =

Īø1Īø2Īø3

where Īøi = 2Īµi/Īµ0 = 2ei tan(Ī±/2)

Here the singularity is at Īµ0 = 0 or Ī± = 180

Hodges Composite Rotor Blade Modeling

Page 31: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

Then the matrix of direction cosines is

CBA =

(1āˆ’ ĪøT Īø

4

)āˆ† + ĪøĪøT

2 āˆ’ Īø

1 + ĪøT Īø4

and the angular velocity is

Ļ‰BAB =

(āˆ†āˆ’ Īø

2

)Īø

1 + ĪøT Īø4

Note that in the limit when Ī± is very small, the parametersĪøi are components of the infinitesimal rotation vector

Hodges Composite Rotor Blade Modeling

Page 32: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

In some cases it is useful to characterize rotation via theso-called ā€œfinite rotation vectorā€ the components of whichare Ī¦i = eiĪ± so that

Ī¦ =

Ī¦1Ī¦2Ī¦3

Then the matrix of direction cosines is

CBA = exp(āˆ’Ī¦)

= āˆ†āˆ’ Ī¦ +Ī¦Ī¦

2āˆ’ . . .

Hodges Composite Rotor Blade Modeling

Page 33: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

The angular velocity is

Ļ‰BAB =

[āˆ†

sinĪ±Ī±āˆ’ Ī¦

(1āˆ’ cosĪ±)

Ī±2 + Ī¦Ī¦T (Ī±āˆ’ sinĪ±)

Ī±3

]Ī¦

=

(āˆ†āˆ’ Ī¦

2

)Ī¦ + . . .

where Ī±2 = Ī¦T Ī¦

Hodges Composite Rotor Blade Modeling

Page 34: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Vectors and DyadicsFinite Rotation, Angular Velocity & Differentiation of VectorsVirtual Rotation & Variation of VectorsVelocity Chain RulesTilde NotationImplications of Eulerā€™s Theorem of RotationEpilogue

The material discussed so far has focused on fundamentalrigid-body kinematicsWe must build on this foundation in the lectures following inorder to understand the kinematical foundation of rigorousbeam theoryThe next step is to develop a suitable method fordetermination of strain-displacement relations

Hodges Composite Rotor Blade Modeling

Page 35: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

For describing the deformation, one needs to introduceframe a which is unaffected by beam deformationThe following development is valid even if a is not aninertial frameNext one needs to specify the position vector from somearbitrary point O fixed in a to an arbitrary material point inthe undeformed beamDenote this with r(x1, x2, x3, t) = where x1 is arclengthalong the beam reference line, x2 and x3 arecross-sectional coordinates, and t is time (henceforth,explicit time dependence is ignored)

Hodges Composite Rotor Blade Modeling

Page 36: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

Now introduce a infinite set of frames b along theundeformed beam with a dextral triad bi(x1) fixed thereinThe unit vector b1(x1) is tangent to the undeformed beamreference line at some arbitrary value of x1

In order to describe the geometry of the undeformedbeam, we need to introduce

covariant basis vectors for the undeformed state

g i (x1, x2, x3) =āˆ‚rāˆ‚xi

contravariant base vectors for the undeformed state

g i (x1, x2, x3) =1

2āˆš

geijk g j Ɨ gk

Hodges Composite Rotor Blade Modeling

Page 37: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

Note that g i Ā· g j = Ī“ij and bi = g i(x1,0,0)

Since g = det(g i Ā· g j) > 0, one can always write

r(x1, x2, x3) = r(x1) + Ī¾

where Ī¾ = x2b2 + x3b3 is the position vector to an arbitraryparticle within the reference cross-sectionConsider the position vector R(x1, x2, x3) from O to thesame particle in the deformed beamReference cross-sections undergo two types of motion

rigid-body translations of the order of the beam length andlarge rigid-body rotationssmall deformation of the reference cross-section

Hodges Composite Rotor Blade Modeling

Page 38: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

B 1

B 2

B 3

Deformed State

Undeformed State

r

R

R

s

ru

x1

b 1

b 2

b 3

R Ė†

r Ė†

Figure: Schematic of Beam Kinematics

Hodges Composite Rotor Blade Modeling

Page 39: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

To capture this behavior mathematically, introduce a set offrames B for the deformed beam analogous to b in theundeformed beamGlobal rotation (from b to B) is described by C = CBb

Based on the above, we can express R in the form

R = R + C Ā· (Ī¾ + w)

whereR = r + uu describes the rigid-body translationthe ā€œwarpingā€ w describes cross-sectional deformation

Hodges Composite Rotor Blade Modeling

Page 40: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

Consider the plane of material points that make up thereference cross-section in the undeformed beam:

neither the planar form nor the section shape are preservedin general if w is nonzerothese points will lie very near a plane in the deformedbeam, the orientation of which is determined by sixconstraints on wthe orientation of B is determined by orientation of thisplane

Consider the set of material points that make up thereference line of the undeformed beam:

the reference line of the deformed beam is not the same setof material pointsB1 is not in general tangent to the deformed beamreference line

Hodges Composite Rotor Blade Modeling

Page 41: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

Only the covariant basis vectors for deformed state areneeded

Gi(x1, x2, x3) =āˆ‚Rāˆ‚xi

The deformation is most concisely described in terms ofthe deformation gradient tensor

A = Gig i

The tensor A has a meaning that is heuristically like āˆ‚Rāˆ‚ r

Hodges Composite Rotor Blade Modeling

Page 42: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

A very useful theorem in continuum mechanics called thepolar decomposition theorem (PDT) states that thedeformation gradient can always be decomposed as

A = C Ā· U

whereC is a finite rotation tensor (corresponding to a purerotation)U is the right stretch tensor (corresponding to a pure strain)Thus, the PDT implies

Gi = A Ā· g i = C Ā· U Ā· g i

Hodges Composite Rotor Blade Modeling

Page 43: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

That is, to arrive at the coordinate lines in the deformedbeam, coordinate lines in the undeformed beam (the g i ā€™s)are

first strained (pre-dot-multiplication by U)then rotated (pre-dot-multiplication by C)

C rotates an infinitesimal material element at (x1, x2, x3)

Some of this rotation is due to global rotation and some isdue to local deformation (i.e., warping)Supposing that the warping is small (in some sense), letā€™sdecompose the total rotation so that

C = C(x1) Ā· Clocal(x1, x2, x3)

where Clocal(x1, x2, x3) is the local rotation tensor

Hodges Composite Rotor Blade Modeling

Page 44: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

Recall that any finite rotation tensor can be written in theform

exp(Ļ†) = āˆ† + Ļ†+Ļ†

2

2+Ļ†

3

6+ . . .

Here we let Clocal = exp(Ļ†)

Thus the total finite rotation tensor at some point in thecross-section is

C = C Ā· exp(Ļ†)

Hodges Composite Rotor Blade Modeling

Page 45: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

Now letā€™s introduce the Cauchy strain tensor (also knownas the engineering strain tensor, the Jaumann straintensor, and the Biot strain tensor)

Ī“ = U āˆ’āˆ†

This strain definition is nice because it does not contain thehost of ā€œstrain squaredā€ terms that are a well known part ofother strain tensors (such as the Green strain)Now, from the PDT and the strain definition, we have

Ī“ = CTĀ· Aāˆ’āˆ† = exp

(āˆ’Ļ†)Ā· CT Ā· Aāˆ’āˆ†

Note that Ī“ is a Lagrangean strain tensor

Hodges Composite Rotor Blade Modeling

Page 46: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

Thus, it is appropriate to resolve it along the undeformedbeam reference triad bi yielding

Ī“ = biĪ“ij bj

It is also appropriate to resolve the tensor Ļ† along theundeformed beam triad bi yielding

Ļ† = bi Ļ†ij bj

It now follows that the deformation gradient tensor isresolved along the mixed bases

A = BiAij bj

orAij = (Bi Ā·Gk )(gk Ā· bj)

Hodges Composite Rotor Blade Modeling

Page 47: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

Things are now simple enough that we can finish up inmatrix formIntroducing

Ļ† =

Ļ†1Ļ†2Ļ†3

the matrix of strain components become

Ī“ = exp(āˆ’Ļ†)

Aāˆ’āˆ†

In general an expression for A can be found rather easily,but Ļ† is unknown

Hodges Composite Rotor Blade Modeling

Page 48: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

For the purpose of simplifying this expression for smalllocal deformation, we let

max |Ī“ij (x1, x2, x3)| = Īµ << 1max |Ļ†ij (x1, x2, x3)| = Ļ• < 1

Then the strain becomes

Ī“ = E āˆ’ Ļ†2

2+

12

(E Ļ†āˆ’ Ļ†E

)+ O

(Ļ•4, Ļ•2Īµ

)where

E =A + AT

2āˆ’āˆ†

Ļ† =Aāˆ’ AT

2

Hodges Composite Rotor Blade Modeling

Page 49: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

Assume that Ļ• = O(Īµr )Since Īµ (the strain) is small compared to unity, two casesare of interest:

Small local rotation: r ā‰„ 1:

Ī“ = E

Moderate local rotation: 12 ā‰¤ r < 1.

Ī“ = E āˆ’ Ļ†2

2+

12

(E Ļ†āˆ’ Ļ†E)

For most engineering beam problems, the small localrotation theory is adequateIn most of the rest, incorporation of warping nonlinearitiesexhibited in moderate local rotation theory should beadequate

Hodges Composite Rotor Blade Modeling

Page 50: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

The purpose of this example is to provide a simpleillustration of the theory developed so far which will give usan explicit form of the strainConsider an initially straight beam that is undergoingplanar deformation without warping (w = 0)For the undeformed state, the position vector r is given by

r = x1b1 + x2b2 + x3b3

The covariant and contravariant base vectors are verysimple

g i =āˆ‚rāˆ‚xi

= bi = g i

Hodges Composite Rotor Blade Modeling

Page 51: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

For the deformed state, the position vector R is given by

R = (x1 + ub1)b1 + x2B2 + ub2b2 + x3b3

Since the warping is zero, the reference cross-sectionalplane in the deformed beam is made up of the samematerial points which make up the referencecross-sectional plane of the undeformed beam

The covariant base vectors are Gi = āˆ‚Rāˆ‚xi

such that

G1 =(1 + uā€²b1)b1 + x2(B2)ā€² + uā€²b2b2

G2 =B2

G3 =b3 = B3

Hodges Composite Rotor Blade Modeling

Page 52: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

The matrix of direction cosines C = CBb can be simplyrepresented with a single angle Ī¶

B1

B2

B3

=

cos Ī¶ sin Ī¶ 0āˆ’ sin Ī¶ cos Ī¶ 0

0 0 1

b1

b2

b3

where Ī¶ is the cross-section rotationWith this, the matrix of deformation gradient components inmixed bases can be written as

A =

(1 + uā€²b1) cos Ī¶ + uā€²b2 sin Ī¶ āˆ’ x2Ī¶ā€² 0 0

āˆ’(1 + uā€²b1) sin Ī¶ + uā€²b2 cos Ī¶ 1 00 0 1

Hodges Composite Rotor Blade Modeling

Page 53: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

The strain for small local rotation applies here since thereis no warping

Ī“ =12

(A + AT )āˆ’āˆ†

so that

Ī“11 =(1 + uā€²b1) cos Ī¶ + uā€²b2 sin Ī¶ āˆ’ x2Ī¶ā€² āˆ’ 1

2Ī“12 =2Ī“21 = uā€²b2 cos Ī¶ āˆ’ (1 + uā€²b1) sin Ī¶

Notice that these strains, when linearized, reduce to thoseof a Timoshenko beam

Hodges Composite Rotor Blade Modeling

Page 54: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

For the undeformed state

r(x1, x2, x3) = r(x1) + xĪ±bĪ±(x1)

where r(x1) is the position vector to points on the referencelineIf the reference line is chosen as the locus ofcross-sectional centroids, then

r = 怈r怉 if and only if 怈xĪ±ć€‰ = 0

where the angle brackets stand for the average value overthe cross-sectionThe covariant base vectors are

g1 = r ā€² + xĪ±bā€²Ī± gĪ± = bĪ±

Hodges Composite Rotor Blade Modeling

Page 55: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

It is helpful to derive some special formulae to express g1in a more recognizable form

r ā€² = b1

(bi)ā€² = k Ɨ bi = k Ā· bi

where k = kbi bi (k ) is the curvature vector (tensor) of theundeformed beam and

k =(

CbA)ā€²Ā· CAb = k Ɨāˆ† = bi kij bj = āˆ’bieijlkbl bj

Normally the components of k are known in the b basis:

kb1 = twist per unit length of the undeformed beamkbĪ± = components of undeformed beam curvature

Hodges Composite Rotor Blade Modeling

Page 56: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

With these definitions, the contravariant base vectorsbecome

g1 =b1āˆš

g

g2 =x3kb1b1āˆš

g+ b2

g3 =āˆ’ x2kb1b1āˆšg

+ b3

where āˆšg = 1āˆ’ x2kb3 + x3kb2 > 0

Normally,āˆš

g is near unity

Hodges Composite Rotor Blade Modeling

Page 57: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

Let us now use the more general displacement fieldreferred to in the development

R(x1, x2, x3) = R + xĪ±BĪ± + wi Bi

Here R = r + u, u = ubi bi is the displacement vector of thebeam, and w = wi bi

The push-forward operation altered the bases on the lasttwo terms

Hodges Composite Rotor Blade Modeling

Page 58: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

We must constrain the warping in order for Bi and R to bewell defined

The 1D position is the average position (i.e. averagewarping over the cross-section is zero)

R =āŸØ

RāŸ©

if and only if 怈wi怉 = 0

Orientation of deformed beam cross-sectional frame

B1 Ā·āŸØ

xĪ±RāŸ©

= 0 if and only if 怈xĪ±w1怉 = 0

The average rotation about B1 is zero

B2 Ā·āŸØ

R,2āŸ©

= B3 Ā·āŸØ

R,3āŸ©

if and only if 怈w3,2 āˆ’ w2,3怉 = 0

Hodges Composite Rotor Blade Modeling

Page 59: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

The strain field can be conveniently expressed in terms of1-D variables with the following definitions:

Rā€² =(1 + Ī³11)B1 + 2Ī³1Ī±BĪ±

Bā€²i =K Ɨ Bi

The force strain components are then

Ī³ = C(e1 + uā€²b + kbub)āˆ’ e1

where

Ī³ =

Ī³11

2Ī³122Ī³13

; e1 =

100

Hodges Composite Rotor Blade Modeling

Page 60: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

The moment strains components are given by Īŗ = KB āˆ’ kbso that

Īŗ = āˆ’Cā€²CT + CkbCT āˆ’ kb

In vector-dyadic form

Ī³ =CzB Ā· Rā€² āˆ’ Czb Ā· r ā€²

Īŗ =CzB Ā· K āˆ’ Czb Ā· k

where z is an arbitrary frameLetting the z frame be b, Ī³ and Īŗ have Ī³ and Īŗ asmeasure numbers in the b basis

Hodges Composite Rotor Blade Modeling

Page 61: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

The vector-dyadic form gives more insight as to what theseexpressions mean; see Hodges (1990) for this discussionSometimes it is better to pull back to the a basis and regardmeasure numbers in the a basis as generalized strainsNote that K is the curvature vector of the deformed beamdefined by

K =(

CBa)ā€²Ā· CaB

where

KB1 = the twist per unit length of the deformed beamKBĪ± = components of deformed beam curvature

Hodges Composite Rotor Blade Modeling

Page 62: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

The strain field for small local rotation is given by matrix Eāˆš

gE11 = Ī³11 + x3Īŗ2 āˆ’ x2Īŗ3 + w ā€²1

+ kb1(x3w1,2 āˆ’ x2w1,3) + kb2w3 āˆ’ kb3w2 + Īŗ2w3 āˆ’ Īŗ3w2

2āˆš

gE12 = 2āˆš

gE21 = 2Ī³12 āˆ’ x3Īŗ1 + w ā€²2 +āˆš

gw1,2

+ kb1(x3w2,2 āˆ’ x2w2,3) + kb3w1 āˆ’ kb1w3 + Īŗ3w1 āˆ’ Īŗ1w3

2āˆš

gE13 = 2āˆš

gE31 = 2Ī³13 + x2Īŗ1 + w ā€²3 +āˆš

gw1,3

+ kb1(x3w3,2 āˆ’ x2w3,3) + kb1w2 āˆ’ kb2w1 + Īŗ1w2 āˆ’ Īŗ2w1

E22 = w2,2 2E23 = 2E32 = w2,3 + w3,2 E33 = w3,3

Notice that if O(Īµ2) terms (underlined) are neglected, thenthe strain field is linear in generalized strains Ī³ and Īŗ

Hodges Composite Rotor Blade Modeling

Page 63: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

Undeformed State GeometryDecomposition of the Rotation TensorStrain-Displacement RelationsSimple ExampleRealistic ExampleEpilogue

An example is worked out for a specified warpingIn general one needs to solve for the warping, which isaffected by cross-sectional geometry and materialproperties as well as initial curvature and twistWarping is typically a linear function of the 1-D strainmeasures, leading to a strain energy density of the formU = U(Ī³, Īŗ)

In other words, warping disappears from the problem,affecting only the elastic constants of the beam (as in theSt. Venant problem)Now we will outline a finite element approach tocross-sectional analysis and then proceed as if sectionconstants are known

Hodges Composite Rotor Blade Modeling

Page 64: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Real rotor blades areinternally complex, built-up structuresmore and more likely to be made of composite materials

These provide certain well-known advantagesHigh strength-to-weight ratioLong fatigue lifeDamage toleranceDirectional nature with potential for tailoring

However, they also introduce well-known complexitiesAnisotropicInhomogeneous

Hodges Composite Rotor Blade Modeling

Page 65: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Figure: V-22 Blade Section ā€“ Courtesy Bell Helicopters

Hodges Composite Rotor Blade Modeling

Page 66: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Rotor blades are 3-D bodies and demand a 3-D approachConsider a 3-D representation of the strain energy

U =12

āˆ« āˆ« āˆ«Ī“T DĪ“ dx2dx3dx1

where Ī“ = Ī“(u) and u = u(x1, x2, x3)

This offers:a complete 3-D description of the probleminhomogeneous, anisotropic, nonlinear ā€“ no problem torepresent, but O(106) degrees of freedom may be required!more than adequate motivation to attempt to use beamtheory

Hodges Composite Rotor Blade Modeling

Page 67: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Figure: Schematic of discretized wing

Hodges Composite Rotor Blade Modeling

Page 68: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

To apply beam theory to composite rotor blades, andexpect an accurate answer, requires us to capture 3-Dbehavior with a 1-D model!Problem: as weā€™ve seen, there is a 3-D quantity (warping)present in the energy functionalTherefore, one must start from a general 3-Drepresentation, and solve the problem including

inhomogeneous, anisotropic materialsall possible deformation in the 3-D representationdetermination of the warping as a function of 1-D variables

Hodges Composite Rotor Blade Modeling

Page 69: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Beams have one dimension much larger than the other twoDimensional reduction takes the 3-D body and representsit as a 1-D bodyThis implies that small parameters must be exploited

maximum magnitude of the strain Īµ << 1a < ` (a is a typical cross-sectional diameter and ` is thecharacteristic length of the deformation along the beam)

a < R (R = 1/āˆš

kTb kb)

The result is the strain energy per unit lengthin terms of 1-D measures of strainwith asymptotically exact cross-sectional elastic constantswith asymptotically exact recovering relations

Hodges Composite Rotor Blade Modeling

Page 70: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

We follow this procedure when modeling a beam:Find 2-D (sectional) elastic constants for use in 1-D (beam)theoryFind 1-D (beam) deformation parameters from loading andsectional constantsFind 3-D displacement, strain, and stress in terms of 1-D(beam) deformation parameters

Analogy from elementary beam theory:Constitutive relation: M2 = EI22u3

ā€²ā€²

Equilibrium equation: M2ā€²ā€² = q(x1)

Recovery relation: Ļƒ11 = āˆ’M2x3I22

Approach based on variational-asymptotic method

Hodges Composite Rotor Blade Modeling

Page 71: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

1-D displacements, generalized strains, & stress resultants

2-D cross-sectional analysis (linear)

2-D cross-sectional elastic and inertia

constants

1-D beam analysis (nonlinear)

2-D warping & strain recovery relations

3-D recovery analysis

3-D stress, strain, & displacement fields

cross-sectional geometry, 3-D elastic constants, & density

initial twist & curvature

loads & boundary conditions

Figure: Process of Beam Analysis

Hodges Composite Rotor Blade Modeling

Page 72: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

B 1

B 2

B 3

Deformed State

Undeformed State

r

R

R

s

ru

x1

b 1

b 2

b 3

R Ė†

r Ė†

Figure: Beam kinematics allows for large displacement and rotationwith small strain and local rotation

Hodges Composite Rotor Blade Modeling

Page 73: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Undeformed state:

r(x1, x2, x3) = r(x1) + xĪ±bĪ± = r(x1) + hĪ¶Ī±bĪ±r ā€² = b1 (bi)

ā€² = k Ɨ bi = k Ā· bi 怈怈xĪ±ć€‰ć€‰ = 0

怈怈ā€¢ć€‰ć€‰ =1|A|

āˆ«Sā€¢dx2dx3 怈ā€¢ć€‰ =

1|A|

āˆ«Sā€¢āˆš

gdx2dx3

āˆšg = 1āˆ’ x2kb3 + x3kb2 > 0 g1 =

b1āˆšg

g2 =x3kb1b1āˆš

g+ b2 g3 = āˆ’x2kb1b1āˆš

g+ b3

Hodges Composite Rotor Blade Modeling

Page 74: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Deformed state

R(x1, x2, x3) = R(x1) + xĪ±BĪ±(x1) + wn(x1, Ī¶2, Ī¶3)Bn(x1)

(Bi)ā€² = K Ɨ Bi = K Ā· Bi

Constraints

Rā€² =(1 + Ī³11)B1

怈怈wn(x1, Ī¶2, Ī¶3)怉怉 =0āŸØāŸØw2,3(x1, Ī¶2, Ī¶3)

āŸ©āŸ©=āŸØāŸØ

w3,2(x1, Ī¶2, Ī¶3)āŸ©āŸ©

Hodges Composite Rotor Blade Modeling

Page 75: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Under the condition of small local rotation,Jaumann-Biot-Cauchy strain measures are

Ī“āˆ— =12

(Ļ‡+ Ļ‡T )āˆ’āˆ†

Ļ‡mn = Bm Ā·āˆ‚Rāˆ‚xk

gk Ā· bn

The matrix Ļ‡ contains components of the deformationgradient tensor in mixed basesIn column matrix form they are arranged as

Ī“ = bĪ“āˆ—11 2Ī“āˆ—12 2Ī“āˆ—13 Ī“āˆ—22 2Ī“āˆ—23 Ī“āˆ—33cT

Hodges Composite Rotor Blade Modeling

Page 76: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

The 3D strain is linear in Ī³11, Īŗ, the warping w , and itsderivatives

Ī“ =1a

Ī“aw + Ī“ĪµĪµ+ Ī“Rw + Ī“`w ā€²

Īµ =

Ī³11Īŗ

Ī³11 =e1

T C(e1 + uā€² + ku)āˆ’ 1

0 =eĪ±T C(e1 + uā€² + ku)

Īŗ =āˆ’ Cā€²CT + CkCT āˆ’ k

Hodges Composite Rotor Blade Modeling

Page 77: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

where

Ī“a =

0 0 0āˆ‚āˆ‚Ī¶2

0 0āˆ‚āˆ‚Ī¶3

0 00 āˆ‚

āˆ‚Ī¶20

0 āˆ‚āˆ‚Ī¶3

āˆ‚āˆ‚Ī¶2

0 0 āˆ‚āˆ‚Ī¶3

Ī“Īµ =

1āˆš

g

1 0 āˆ’Ī¶3 Ī¶20 Ī¶3 0 00 āˆ’Ī¶2 0 00 0 0 00 0 0 00 0 0 0

Ī“R =

1āˆš

g

[k + k1

(Ī¶3

āˆ‚āˆ‚Ī¶2āˆ’ Ī¶2

āˆ‚āˆ‚Ī¶3

)āˆ†

0

]Ī“` =

1āˆš

g

āˆ†0

Hodges Composite Rotor Blade Modeling

Page 78: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

The strain energy density for a beam per unit length

U =12

āŸØĪ“T DĪ“

āŸ©The 3-D Jaumann stress Ļƒ, which is conjugate to theJaumann strain Ī“ is

Ļƒ = DĪ“

The basic 3-D problem can be now represented as thefollowing minimization problemāˆ«

U[Īµ(x1),w(x1, Ī¶2, Ī¶3),

āˆ‚w(x1, Ī¶2, Ī¶3)

āˆ‚Ī¶Ī±,w ā€²(x1, Ī¶2, Ī¶3)

]dx1

+ potential energy of external forcesā†’ min

Hodges Composite Rotor Blade Modeling

Page 79: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

In constructing a 1-D beam theory from 3-D elasticity, oneattempts to represent the strain energy stored in the 3-Dbody by finding the strain energy which would be stored inan imaginary 1-D bodyThe warping displacement components wn(x1, Ī¶2, Ī¶3) mustbe written as functions of the 1-D functions R(x1) andBn(x1)

This is too complicated to do exactly due to nonlocaldependenceOne can and must take advantage of small parameters

Hodges Composite Rotor Blade Modeling

Page 80: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Two small parametersThe ratio a

` of the maximum dimension of the cross-section(a) divided by the characteristic wavelength of thedeformation along the beam (`)The maximum dimension of the cross section times themaximum magnitude of initial curvature or twist a

R

Since both of them have the same numerator, expansion ina` and a

R is equivalent to the expansion in a onlyNote the following concerning the maximum strainmagnitude Īµ is

necessary for determining the strain fieldonly needed in the cross-sectional analysis when analyzingthe trapeze effect

Hodges Composite Rotor Blade Modeling

Page 81: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

The variational-asymptotic methodis essentially the work of Berdichevsky and co-workers(1976, 1979, 1981, 1983, 1985, etc.)considers small parameters applied to energy functionalsrather than to differential equations

If Āµ is a typical material modulus, the strain energy is of theform

ĀµĪµ2[O(1) + O

(a`

)+ O

( aR

)+ O

(a`

)2+ O

( aR

)2+ O

(a2

`R

)+ . . .

]

Hodges Composite Rotor Blade Modeling

Page 82: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Finite element discretization of the warping

w(x1, Ī¶2, Ī¶3) = S(Ī¶2, Ī¶3)W (x1)

Strain energy density

2U =

(1a

)2

W T EW

+

(1a

)2W T (DaĪµĪµ+ DaRW + Da`W ā€²)

+(1)(ĪµT DĪµĪµĪµ+ W T DRRW + W ā€²T D``W ā€²

+2W T DRĪµĪµ+ 2W ā€²T D`ĪµĪµ+ 2W T DR`W ā€²)

Hodges Composite Rotor Blade Modeling

Page 83: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

The following definitions were introduced

E =āŸØāŸØ

[Ī“a S]T D [Ī“a S]āŸ©āŸ©

DaĪµ =āŸØāŸØ

[Ī“a S]T D [Ī“Īµ]āŸ©āŸ©

DaR =āŸØāŸØ

[Ī“a S]T D [Ī“R S]āŸ©āŸ©

Da` =āŸØāŸØ

[Ī“a S]T D [Ī“` S]āŸ©āŸ©

DĪµĪµ =āŸØāŸØ

[Ī“Īµ]T D [Ī“Īµ]

āŸ©āŸ©DRR =

āŸØāŸØ[Ī“R S]T D [Ī“R S]

āŸ©āŸ©D`` =

āŸØāŸØ[Ī“` S]T D [Ī“` S]

āŸ©āŸ©DRĪµ =

āŸØāŸØ[Ī“R S]T D [Ī“Īµ]

āŸ©āŸ©D`Īµ =

āŸØāŸØ[Ī“` S]T D [Ī“Īµ]

āŸ©āŸ©DR` =

āŸØāŸØ[Ī“R S]T D [Ī“` S]

āŸ©āŸ©Note: these matrices carry information about the materialproperties and geometry of a given cross section

Hodges Composite Rotor Blade Modeling

Page 84: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Warping field is expanded

W = W0 + aW1 + a2W2

Constraints are discretized

W T HĪØc` = 0

where

H =āŸØāŸØ

ST SāŸ©āŸ©

EĪØc` = 0 ĪØTc`HĪØc` = āˆ†

Hodges Composite Rotor Blade Modeling

Page 85: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Strain energy up to order a2

2U = (1)[ĪµT DĪµĪµĪµ+ W0T EW0 + 2W0

T DaĪµĪµ]+

2(a)[W0T EW1 + W1

T DaĪµĪµ+ W0T DaRW0 + W0

T DRĪµĪµ]+

(a2)[2W0T EW2 + W1

T EW1 + 2W0T (DaR + DT

aR)W1

+ 2W2T DaĪµĪµ+ W0

T DRRW0 + 2W1T DRĪµĪµ]

Hodges Composite Rotor Blade Modeling

Page 86: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Evaluation of W0According to the variational-asymptotic method, one keepsonly the dominant interaction term between W and Īµ andthe dominant quadratic term in W

2U0 =

(1a

)2

W T EW +

(1a

)2W T DaĪµĪµ

The Euler-Lagrange equation (including constraints) is(1a

)EW + DaĪµĪµ = HĪØc`Āµ

The Lagrange multiplier becomes

Āµ = ĪØTc`DaĪµĪµ

Hodges Composite Rotor Blade Modeling

Page 87: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Evaluation of W0 (continued)The warping becomes(

1a

)EW = āˆ’(āˆ†āˆ’ HĪØc`ĪØ

Tc`)DaĪµĪµ

Note the generalized inverse

EE+c` = āˆ†āˆ’ HĪØc`ĪØ

Tc`

E+c`E = āˆ†āˆ’ĪØc`ĪØ

Tc`H

E+c`EE+

c` = E+c`

Hodges Composite Rotor Blade Modeling

Page 88: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Evaluation of W0 (continued)

W = āˆ’aE+c`DaĪµĪµ = W0

Prismatic beam stiffness matrix

2U = ĪµT AĪµ

A = DĪµĪµ āˆ’ [DaĪµ]T E+

c`[DaĪµ]

Hodges Composite Rotor Blade Modeling

Page 89: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Evaluation of W1Perturbing the warping W , one obtains

W = W0 + aW1

The perturbation of the energy then becomes

2U1 =a2W1T EW1+

+2aW1T (DaĪµĪµ+ EW0 + a2DaRW1)

+2a2W1T (DaR + DT

aR)W0

Hodges Composite Rotor Blade Modeling

Page 90: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Evaluation of W1 (continued)Minimization of the perturbed energy (including constraints)yields the Euler-Lagrange equation

EW1 +

(1a

)HĪØc`ĪØ

Tc`DaĪµĪµ+ (DaR + DT

aR)W0 = HĪØc`Āµ

The Lagrange multiplier becomes

Āµ =

(1a

)ĪØT

c`DaĪµĪµ+ ĪØTc`(DaR + DT

aR)W0

The warping then becomes

EW1 = āˆ’(āˆ†āˆ’ HĪØc`ĪØTc`)(DaR + DT

aR)W0

Hodges Composite Rotor Blade Modeling

Page 91: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Evaluation of W1 (continued)The strain energy corrected to first order in a/R is found (allinfluence of the perturbed warping cancels out)

2U = ĪµT Ar Īµ

where

Ar =DĪµĪµ āˆ’ (DaĪµ)T (ĪØc`DaĪµ)

+ a[(ĪØc`DaĪµ)T (DaR + DT

aR)(ĪØc`DaĪµ)

āˆ’ (ĪØc`DaĪµ)T DRĪµ āˆ’ DT

RĪµ(ĪØc`DaĪµ)]

Higher-order corrections in a are presented in detail in thebook by Hodges (2006)

Hodges Composite Rotor Blade Modeling

Page 92: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

The 1-D constitutive law that follows from the energy is ofthe form

FB1

MB

= [S]

Ī³11Īŗ

For isotropic, prismatic beams when x2 and x3 are principalaxes

FB1

MB1

MB2

MB3

=

EA 0 0 00 GJ 0 00 0 EI2 00 0 0 EI3

Ī³11Īŗ1Īŗ2Īŗ3

Hodges Composite Rotor Blade Modeling

Page 93: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Adding modest pretwist and initial curvature, and lettingD = I2 + I3 āˆ’ J and Ī² = 1 + Ī½, one obtains

FB1

MB1

MB2

MB3

=

EA EDk1 āˆ’Ī²EI2k2 āˆ’Ī²EI3k3

EDk1 GJ 0 0āˆ’Ī²EI2k2 0 EI2 0āˆ’Ī²EI3k3 0 0 EI3

Ī³11Īŗ1Īŗ2Īŗ3

For generally anistropic beams, the matrix S becomes fullypopulated (while remaining symmetric and positivedefinite, of course)

Hodges Composite Rotor Blade Modeling

Page 94: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Taking into account the effects of a` , one finds a model for

prismatic, isotropic beams of the form

FB1

FB2

FB3

MB1

MB2

MB3

=

EA 0 0 0 0 00 GK2 0 0 0 00 0 GK3 0 0 00 0 0 GJ 0 00 0 0 0 EI2 00 0 0 0 0 EI3

Ī³112Ī³122Ī³13Īŗ1Īŗ2Īŗ3

Initial twist and curvature add coupling which shifts theneutral and/or shear centersFor generally anistropic beams, the matrix S becomes fullypopulated (while remaining symmetric and positivedefinite, of course)

Hodges Composite Rotor Blade Modeling

Page 95: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Cross-sectional analysis can be undertaken by thecomputer program VABS (commercially available)VABS produces

the 6Ɨ6 matrices for stiffness and inertia propertiesrecovery relations that allow one to use results from thebeam analysis to find 3D stresses, strains anddisplacements

Hodges Composite Rotor Blade Modeling

Page 96: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Figure: Sample cross-sectional PreVABS mesh

Hodges Composite Rotor Blade Modeling

Page 97: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Figure: Sample stiffness output from cross-sectional analysis VABS

Hodges Composite Rotor Blade Modeling

Page 98: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

Figure: Sample stress output from cross-sectional analysis VABS

Hodges Composite Rotor Blade Modeling

Page 99: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

IntroductionApproachKinematicsStrain FieldDimensional Reduction from 3-DEpilogue

A comprehensive beam modeling scheme is presentedwhich

allows for systematic treatment of all possible types ofdeformation in composite beamsis ideal for 2-D finite element sectional analysisgives asymptotically exact section constantsleads to the geometrically exact beam equations, found inHodges (2006)gives formulae needed for recovering strain and stressdistributionsis the basis for the commercial computer code VABS

Hodges Composite Rotor Blade Modeling

Page 100: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

Present an exact set of equations for the dynamics ofbeams in a moving frame suitable for rotorcraft applicationsPresent the equations in a matrix notation so that theentire formulation can be written most conciselyPresent a unified framework in which other less generaldevelopments can be checked for consistencyShow how differences that normally arise in beamanalyses can be interpreted in light of this unifiedframeworkShow how the present unified framework can be used todevelop an elegant basis for accurate, efficient, and robustcomputational solution techniques

Hodges Composite Rotor Blade Modeling

Page 101: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

Published work reveals differences in the way thedisplacement field is represented

different numbers of kinematical variables to describemotion of the cross sectional framedifferent variables to describe finite rotation of this framedifferent orthogonal base vectors for measurement ofdisplacement

In the present approach the kinematical equations areexact and separate from equilibrium and constitutive lawdevelopments

generalized strains are written in simple matrix notationchange of displacement and orientation variables affectsonly kinematics (a relatively small portion of the analysis)

Hodges Composite Rotor Blade Modeling

Page 102: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

Reissner (1973) derived exact intrinsic equations for beamstatic equilibrium (limited to unrestrained warping)

no displacement or orientation variables (ā€œintrinsicā€)reduce to the Kirchhoff-Clebsch-Love equations whenshear deformation is set equal to zerogeometrically exact ā€“ all correct beam equations can bederived from theseintrinsic generalized strains were derived from virtual work

Asymptotic analysis shows that for slender, closed-sectionbeams, a linear 2-D cross-sectional analysis determineselastic constants for use in nonlinear 1-D beam analysis

Hodges Composite Rotor Blade Modeling

Page 103: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

Thus, in the present work we presuppose that an elasticlaw is given as a 1-D strain energy functionPresent analysis is based on exact kinematics and kineticsbut an approximate constitutive law (no ā€œordering schemeā€)Here the exact equilibrium equations are

extended to account for dynamicsderived from Hamiltonā€™s weak principle (HWP) by Hodges(1990) to facilitate development of a finite element method

A mixed finite element approach can be developed fromHWP in which there are many computational advantagesWidely available codes RCAS and DYMORE are based onthese equations, albeit in displacement form in the latter

Hodges Composite Rotor Blade Modeling

Page 104: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

Recall that the displacement field is represented by

r =r + Ī¾ = r + x2b2 + x3b3

R =R + C Ā· (Ī¾ + w) = r + u + x2B2 + x3B3 + wi Bi

with C as the global rotation tensor

Hodges Composite Rotor Blade Modeling

Page 105: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

Force strains

Ī³ =

Ī³11

2Ī³122Ī³13

= C(

e1 + uā€²b + kbub

)āˆ’ e1

where ub is the column matrix whose elements are themeasure numbers of the displacement along bi

Moment strains

Īŗ =

Īŗ1Īŗ2Īŗ3

= KB āˆ’ kb

where KB = āˆ’Cā€²CT + CkbCT

Hodges Composite Rotor Blade Modeling

Page 106: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

Introduce column matrices of known measure numbersalong bi for inertial

velocity of the undeformed beam reference axis vbangular velocity of the undeformed beam cross sectionalframe Ļ‰b

Generalized speeds in the sense of Kane and Levinson(1985) are elements of column matrices that containmeasure numbers along Bi for inertial

velocity of deformed beam reference axis

VB = C(vb + ub + Ļ‰bub)

angular velocity of deformed beam cross-sectional frame

Ī©B = āˆ’CCT + CĻ‰bCT

Hodges Composite Rotor Blade Modeling

Page 107: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

As developed by Borri et al. (1985), HWP for the presentproblem is āˆ« t2

t1

āˆ« `

0

[Ī“(K āˆ’ U) + Ī“W

]dx1dt = Ī“A

Heret1 and t2 are arbitrary fixed timesK and U are the kinetic and strain energy densities per unitlength, respectivelyĪ“A is the virtual action at the ends of the beam and at theends of the time intervalĪ“W is the virtual work of applied loads per unit length

Hodges Composite Rotor Blade Modeling

Page 108: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

St. Venant warping influences the elastic constants in abeam constitutive law written in terms of Ī³ and ĪŗRegarding the strain energy per unit length as U = U(Ī³, Īŗ),one can obtain the variations required in HWP asāˆ« `

0Ī“Udx1 =

āˆ« `

0

[Ī“Ī³T

(āˆ‚Uāˆ‚Ī³

)T

+ Ī“ĪŗT(āˆ‚Uāˆ‚Īŗ

)T]

dx1

The partial derivatives are section force and momentmeasure numbers along Bi

FB =

(āˆ‚Uāˆ‚Ī³

)T

; MB =

(āˆ‚Uāˆ‚Īŗ

)T

Hodges Composite Rotor Blade Modeling

Page 109: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

Introduce a column matrix of virtual displacements definedas Ī“qB = CĪ“ub

Similarly, let the antisymmetric matrix of virtual rotations beĪ“ĻˆB = āˆ’Ī“CCT

Now, one can show that

Ī“Ī³ =Ī“qā€²B + KBĪ“qB + (e1 + Ī³)Ī“ĻˆB

Ī“Īŗ =Ī“Ļˆā€²B + KBĪ“ĻˆB

so that there are neither displacement nor orientationvariables present in the variations

Hodges Composite Rotor Blade Modeling

Page 110: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

The kinetic energy per unit length is K = K (VB,Ī©B)

Thus, the variation required in HWP isāˆ« `

0Ī“Kdx1 =

āˆ« `

0

[Ī“V T

B

(āˆ‚Kāˆ‚VB

)T

+ Ī“Ī©TB

(āˆ‚Kāˆ‚Ī©B

)T]

dx1

Introduce sectional linear and angular momenta, PB andHB, that are conjugate to the generalized speeds

PB =

(āˆ‚Kāˆ‚VB

)T

= m(VB āˆ’ Ī¾BĪ©B)

HB =

(āˆ‚Kāˆ‚Ī©B

)T

= iBĪ©B + mĪ¾BVB

Hodges Composite Rotor Blade Modeling

Page 111: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

With the above definitions of virtual displacement andvirtual rotation, the variations become

Ī“VB = Ė™Ī“qB + Ī©BĪ“qB + VBĪ“ĻˆB

Ī“Ī©B = Ė™Ī“ĻˆB + Ī©BĪ“ĻˆB

which are, as with generalized strain variations,independent of displacement or orientation variablesThese are needed so that contributions of kinetic energy toequilibrium equations can be obtained withoutdisplacement or orientation variables in them

Hodges Composite Rotor Blade Modeling

Page 112: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

The virtual work of external forces per unit length is

Ī“W =

āˆ« `

0

(Ī“qT

B fB + Ī“ĻˆTBmB

)dx1

The virtual action at the ends of the beam and of the timeinterval is

Ī“A =

āˆ« `

0

(Ī“qT

B PB + Ī“ĻˆTBHB

)āˆ£āˆ£āˆ£t2t1

dx1

āˆ’āˆ« t2

t1

(Ī“qT

B FB + Ī“ĻˆTBMB

)āˆ£āˆ£āˆ£`0

dt

Hodges Composite Rotor Blade Modeling

Page 113: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

āˆ« t2

t1

āˆ« `

0

(Ė™Ī“q

T

B āˆ’ Ī“qTBĪ©B āˆ’ Ī“Ļˆ

TB VB

)PB

+

(Ė™Ī“Ļˆ

T

B āˆ’ Ī“ĻˆTBĪ©B

)HB

āˆ’[(Ī“q

ā€²B

)Tāˆ’ Ī“qT

B KB āˆ’ Ī“ĻˆTB(e1 + Ī³

)]FB

āˆ’[(Ī“Ļˆ

ā€²B

)Tāˆ’ Ī“ĻˆT

B KB

]MB + Ī“q

TB fB + Ī“Ļˆ

TBmB

dx1dt

=

āˆ« `

0

(Ī“q

TB PB + Ī“Ļˆ

TBHB

)āˆ£āˆ£āˆ£t2t1

dx1 āˆ’āˆ« t2

t1

(Ī“q

TB FB + Ī“Ļˆ

TBMB

)āˆ£āˆ£āˆ£`0

dt

Hodges Composite Rotor Blade Modeling

Page 114: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

The Euler-Lagrange equations from HWP are thegeometrically-exact partial differential equations of motion

F ā€²B + KBFB+fB = PB + Ī©BPB

M ā€²B + KBMB+(e1 + Ī³

)FB + mB = HB + Ī©BHB + VBPB

HWP also leads to a consistent set of boundary conditionsin which either force or moment can be specified or foundat the ends of the beamThese equations are

geometrically exact equations for the dynamics of a beamin a frame A whose inertial motion is arbitrary and knownidentical to those of Reissner (1973) when specialized tothe static case

Hodges Composite Rotor Blade Modeling

Page 115: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

In order to finalize the development, we use the columnmatrix of Rodrigues parameters Īø as orientation variablesNow the direction cosine matrix can easily be expressed as

C =(1āˆ’ ĪøT Īø

4 )āˆ†āˆ’ Īø + ĪøĪøT

2

1 + ĪøT Īø4

Similarly, Īŗ and Ī© are

Īŗ =

(āˆ†āˆ’ Īø

2

1 + ĪøT Īø4

)Īøā€² + Ckb āˆ’ kb

Ī©B =

(āˆ†āˆ’ Īø

2

1 + ĪøT Īø4

)Īø + CĻ‰b

Hodges Composite Rotor Blade Modeling

Page 116: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

The last two equations should be inverted

Īøā€² =

(āˆ† +

12Īø +

14ĪøĪøT

)(KB āˆ’ Ckb)

and

Īø =

(āˆ† +

12Īø +

14ĪøĪøT

)(Ī©B āˆ’ CĻ‰b)

Also, force strain and velocity equations should be inverted

uā€²b = CT (e1 + Ī³)āˆ’ e1 āˆ’ kbub

andub = CT VB āˆ’ vb āˆ’ Ļ‰bub

Hodges Composite Rotor Blade Modeling

Page 117: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

The column matrices ub, Īø, Ī³, Īŗ, VB, Ī©B, FB, MB, PB, andHB are regarded as independent quantitiesCentral differencing in the spatial variable x1 has beenshown to be equivalent to the mixed finite elementformulation of Hodges (2006)Coefficient matrices are very sparse

Hodges Composite Rotor Blade Modeling

Page 118: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

Nondimensional force

Nondimensional displacements

0 2 4 6 8 100

0.2

0.4

0.6

0.8vertical

horizontal

rotational

Figure: Large displacements of cantilever beam

Hodges Composite Rotor Blade Modeling

Page 119: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

ā€¢ ā€¢ ā€¢ ā€¢

ā€¢ā€¢

ā€¢ā€¢

ā€¢

ā€¢

ā€¢ā€¢

0

10

20

30

40

50

60

70

0 15 30 45

Freq

uenc

y (H

z)

Sweep Angle (degrees)

750 RPM

500 RPM

0 RPM

ā€¢ Experiment (Epps & Chandra) Theory (Epps & Chandra) Theory (Present)

Figure: Frequency of the third bending mode for swept-tip beam

Hodges Composite Rotor Blade Modeling

Page 120: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

0 15 30 45 60 75 900

1

2

3

4

5

Īø (degrees)

MVFDTH

Exp., P = 1Exp., P = 2Exp., P = 3

Ļ†T

IP (d

egre

es)

Ė†

Figure: Princeton beam experiment, torsion for various P versussetting angle Īø

Hodges Composite Rotor Blade Modeling

Page 121: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

0 15 30 45 60 75 900

2

4

6

8

10

Īø (degrees)

wT

IP (i

nche

s)

MVFDTHLinear

Exp., P = 1Exp., P = 2Exp., P = 3

Figure: Princeton beam experiment, out-of-plane bending for variousP versus setting angle Īø

Hodges Composite Rotor Blade Modeling

Page 122: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

0 15 30 45 60 75 900

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

Īø (degrees)

v TIP

(inc

hes)

MVFDTHLinearExp., P = 1Exp., P = 2Exp., P = 3

Figure: Princeton beam experiment, in-plane bending for various Pversus setting angle Īø

Hodges Composite Rotor Blade Modeling

Page 123: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

0 10 20 30 40 50 60 70 80 901.1

1.14

1.18

1.22

1.26

1.3

1.34

1.38

Īø (degrees)

f (H

z)

MVFDTH

LinearExp., P = 2

Figure: Princeton beam experiment, out-of-plane bending frequencyversus P for various setting angles Īø

Hodges Composite Rotor Blade Modeling

Page 124: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

0 10 20 30 40 50 60 70 80 903

3.4

3.8

4.2

4.6

5

5.4

Īø (degrees)

f (H

z)

MVFDTH

LinearExp., P = 2

Figure: Princeton beam experiment, in-plane bending frequencyversus P for various setting angles Īø

Hodges Composite Rotor Blade Modeling

Page 125: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

0 2 4 6 8 10 120

0.4

0.8

1.2

1.6

2

P (lb)

f (H

z)

MVFDTH

LinearExp., Īø = 0

Figure: Princeton beam experiment, out-of-plane frequency versus Pfor setting angle Īø = 0

Hodges Composite Rotor Blade Modeling

Page 126: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

The Euler-Lagrange equations are the kinematical,constitutive, and equations of motion

The kinematical equations are written exactly utilizingso-called intrinsic strain measuresThe equations of motion are written exactly in their intrinsicformThe constitutive law is presumed given and is left in ageneric formThe choice of displacement and rotational variables islocalized in a relatively small portion of the analysis

Hodges Composite Rotor Blade Modeling

Page 127: Rotorcraft Dynamics

Daniel Guggenheim School of Aerospace Engineering

IntroductionKinematical Preliminaries

Analysis of 3-D Beam DeformationStiffness Modeling

Geometrically Exact Beam Equations

ObjectivesPresent Approach1-D KinematicsHamiltonā€™s Weak PrincipleMixed Variational FormulationSample ResultsEpilogue

When specialized, the equations reduce to less generaltreatments in the literatureAlthough the resulting equilibrium equations are identicalto those derived from a Newtonian method, the formulationis variationally consistentThe present development provides substantial insight intorelationships among variational formulations as well asbetween these and Newtonian onesCoefficient matrices are very sparse and the computationalefficiency can be improved by taking advantage thereof

Hodges Composite Rotor Blade Modeling