rotorcraft dynamics
TRANSCRIPT
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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