Shimmy Design
Shimmy Design
336853 AQARD-U-800
T A \ $D
ADVISORY GROUP FOR AEROSPACE RESEARCH & DEVELOPMENT
7 RUE ANCELLE, 92200 NEUILLY-SUR-SEINE, FRANCE
-
Papers presented at the 81s Meeting of the AGARD Structum d Materials Panet
held in BanB Canada Octobw lH5.
..........signed.
by DIMS
Papers presented at the 81st Meeting of the AGARD Structures and Materials Panel,
held in Banff, Canada 4-5 October 1995.
According to its Charter, the mission of AGARD is to bring together the leading personalities of the NATO nations in the
fields of science and technology relating to aerospace for the following purposes:
- Recommending effective ways for the member nations to use their research and development capabilities for the
common benefit of the NATO community;
- Providing scientific and technical advice and assistance to the Military Committee in the field of aerospace research
and development (with particular regard to its military application);
- Continuously stimulating advances in the aerospace sciences relevant to strengthening the common defence posture;
- Improving the co-operation among member nations in aerospace research and development;
- Providing assistance to member nations for the purpose of increasing their scientific and technical potential;
- Rendering scientific and technical assistance, as requested, to other NATO bodies and to member nations in
connection with research and development problems in the aerospace field.
The highest authority within AGARD is the National Delegates Board consisting of officially appointed senior
representatives from each member nation. The mission of AGARD is carried out through the Panels which are composed of
experts appointed by the National Delegates, the Consultant and Exchange Programme and the Aerospace Applications
Studies Programme. The results of AGARD work are reported to the member nations and the NATO Authorities through the
AGARD series of publications of which this is one.
Participation in AGARD activities is by invitation only and is normally limited to citizens of the NATO nations.
ISBN 92-836-1032-6
ii
The Design, Qualification and Maintenance of
Vibration-Free Landing Gear
(AGARD R-800)
Executive Summary
Aircraft landing gears are crucial for safety, comfort (both for passengers and pilots) and for weight
considerations. As the element responsible for safely moving the aircraft on the ground, the landing
gear has to fulfill several, sometimes conflicting, requirements.
I
Landing gears that shimmy (shimmy can be defined as a self-excited instability during take-off, landing
or taxiing, involving up to three vibration motions: angular wheel motions about a vertical axis - yaw -,
and a fore and aft axis - roll -, and lateral displacement of the wheel) are unacceptable. In fact, a severe
occurrence of shimmy can damage the landing gear and its attaching structure, resulting in significant
I repair costs and airplane down time. Some assurance is therefore needed that landing gear designs will
be free from shimmy under all operating conditions including the normal wear and tear experienced in
service.
- One of the difficulties of shimmy analysis is that real landing gear systems exhibit many non-linear
r
characteristics. Tests on life-size aircraft are obviously expensive and risky, and tests on test-rigs
(namely drop-test facilities) allow only for limited information about the landing gear's dynamics; the
"
interaction between aircraft and landing gear is especially difficult to assess. On the other hand,
'- simulation offers a means to examine the behaviour of the landing gear as part of a complex system at a
reasonable cost. Both rigid and elastic body motions can be modelled.
In view of the problems involved in correcting landing gear vibrations emerging late in the design
process or even after delivery of the aircraft to service, the Workshop participants concluded that the
state of the art of analyzing landing gear vibrations is not quite up to other similar subjects and that the
cooperation of experts across aircraft, landing gear and tyre industry institutions is badly needed.
Workshop participants did not expect rapid progress unless there was a well planned and coordinated
approach to the problem. They unanimously identified AGARD as the only institution which had
enough power and authority to promote such a coordinated effort in the interest of NATO member
nations.
iii
Letude, lhomologation et la maintenance des
trains datterrissage B amortissement
(AGARD R-800) I1 1
Synthkse
Le train datterrissage dun aCronef est un ClCment determinant pour la sCcuritC et le confort des
passagers et des Cquipages, en plus des considCrations de coefficient de chargement. Cest le train
datterrissage qui autorise le dkplacement de lavion au sol dans des conditions de sCcuritC acceptables.
I1 doit, par condquent, satisfaire B plusieurs critkres, qui sont parfois contradictoires.
I
I
seraient, bien entendu, trop cofiteux et trop dangereux, tandis que les essais effectuCs sur des
installations fixes (et notamment les Cpreuves de chute) ne donnent que des informations limities sur
les caractCristiques dynamiques du train; linteraction entre 1aCronef et le train datterrissage est
particulikrement dklicate B Cvaluer. En revanche, la simulation permet dCtudier le comportement du
train datterrissage en tant quC1Cment constitutif dun systkme complexe, et ce, pour un coot abordable.
La simulation permet la modClisation des mouvements rigides et Clastiques du fuselage.
iv
Contents
Page
Synthbe iv
Preface vi
Reference
Analysis a n d Control of the Flexible Dynamics of Landing Gear in the Presence of Antiskid 4
Control Systems
by E. Denti and D. Fanteria
Unsteady Tire Dynamics and the Application Thereof to Shimmy and Landing Load
Computations
by K. Koenig
A Nonlinear Model for Landing Gear Shimmy with Applications to the McDonnell Douglas 10
FIA-18A
by J. Baumann
V
Preface
Fully reliable procedures for designing vibration-free landing gear still do not exist. This is in large part due to the absence of
accurate dynamic models for describing the tyres in ground contact, the complexity of the (generally non-linear) dynamic
behaviour of the structural systems involved as well as the dynamic interactions with steering and braking systems which are I
contributing stability factors.
The Workshop focused on the various vibrational and stability problems (e.g. shimmy, antiskid induced vibrations) that must
be considered in the early design phase of landing gear systems, thereby especially addressing problems which are related to
vibrations of the combined structural system formed by the landing gear, its tyres and the flexible aircraft structure. The
intention was to indicate the impact of (combined) landing geadaircraft vibration problems on aircraft design and operations.
A further aim of the Workshop was to bring together specialists from aircraft, landing gear and tyre manufacturers to discuss
the state-of-the-art technology in this area and to define possible future steps of development.
Dr. R. Freymann
Workshop Chairman
vi
Structures and Materials Panel
SUB-COMMITTEE MEMBERS
I - -
PANEL EXECUTIVE
vii
Technical Evaluation Report
A. Kraw
Daimler-Benz Aerospace AG
Military Aircraft LME24
8 1663 Munich
GermiWJJ
Present author has spent almost thirty years of his TYPE AND CAUSE OF VIBRATION
professional life in nircraft and landing gear design
and analysis. It appeared to him that to arrive at a In s u m m a r y , m o s t a t t e n t i o n w a s p a i d to
design satisfying straightforward and clearly written conventional shimmy, which implies a combined
physical requirements was hard enough. It was even lateral I torsional I tilting motion of the wheel (3
harder to defend such design from the particular degrrrs of freedom). In particular, [21, [81. nnd t 101
interests of other design disciplines. However, at concentrate on this type of landing gcar vibration.
least in military landing gear performance and Cases of forelaft oscillations (1 d.0.f.) were
design requirements the situation with regard to presented at 131 and 141. Though not a cax of self-
"vibration-free" landing gear was even more inducal oscillatiom, there was a study on the use of
difficult. There was not only a lack of guidance with semi-active landing ge& technology (control of
regard to acceptable methods of design and stroke damping force coefficient) to damp fuselage
analysis, the requirement f o r a vibration-free vertieal bending oscillations o f a passenger
landing gear itself was compromised by accepting transport aircraft [5]. A view on forced (vertical)
"shimmy" to the extent that pilots could still control potentially r e s o n a n t a i r c r a f t l l a n d i n g g e a r
the aircraft, that shimmy loads did not exceed oscillations on rwgb p u n d is given at [I 1. In one
structural limits, and that the phenomenon was ade- of the two classical shimmy cases presented at t71.
quately taken into account in fatigue analysis (see the analytical model included a stick model of the
for instance British requirements). What wonder wing which the landing gear is attached to. At [Z],
that budgets for shimmy analyses were limited at there is also evidence that a predictedly stable
best and that an analytical prediction of landing landing gear exhibited shimmy because of
gear shimmy was not considered a hard fact. Hence neglecting attachment/fuselage flexibility in the
the landing gear designer did not get much support analysis. Staning from isolated h n t wheel shimmy
in striving for a landing gear free of shimmy. of motorcycles, [6] nlso considers oscillationsof the
compl*6 vehiile.
It became evident bat the classical problem of
shimmy is becoming increasingly complex due to model which in its
coupling with other k r a f t systems. [SI indicates measurement of 36
that a fully active landing gear could cause stability
problems in vertical (stroking)d d o n ; in view of
[l]. [3], [41, f61,and 171 pnscat author expects a lot
of additional problmns if a fully active landing par
were to b e p r o v e n t o be v i b r a t i o n - f r e e . of modcl parameters is
notwithstanding ground coupling of flight control
system. The coupling of the flexible dynamics of
landing gear and antiskid feedback dynamics to
cause fore-aftvibrations ("gemalk")is studied at with 9 elementary tnnsfer
[4]. A r e l a t e d c a s e of s t r o n g braking f o r c e manufacturer dynamic tyre
oscillations (involving virtually rigid wheel motorcycles. [6] not only
suspension) is being &Mibed at [3]. At [7lthen is .o.f. but also reviews so-
indicated a case invalving feedback from the
s p n i n g acwam.
i
represent tyre dynamic pmeertics in solution of one Paper [Z]trea the relative influence of model
existing and one predicted sbimmy problem. 171 parameter vari tions on shimmy stability pmdic-
pmfers the von Schlippc medel due 10 difficlilties tions; it also p esents practical experience w.r.t.
encountered in defining the Moreland tyre time pndictionlnali disnepancics and reasons t h d .
constant. Paper [la] refers to the Moreland point- Paper [31 treat a case of apparent shimmy (which,
accordtng to the explanauon given i the paper. is reinforcement for increascd stiffness or in form of
rather a case of camouflaged gear walk), which mass balance for changing vibration modes and/or
vanishes at increased speed presumably due to frequencies.
nonlinearity of tyre circumferential force charac-
teristics. Also in [3] there is descripted a rig test There was not one serious attempt reported to
case (virmally rigid) w h m tyre and braUantiskid modify geometry of a landing gear affe*ed, for in-
system nonlinurities pmhibit linear treatment. [41 e implement a m m suitable
s t ~ c to idnation
deseribcs a similar phenomenon observed on and/or mil m.It appurs that volume end form of
aim& a& embarks on dcvaloptnsnt of a nonlinear landing gear bays are being defined to the laat
(time domain) model f o r calculation of tyre millimeter a t a much too e a r l y time i n t h e
longitudinal force. [5] discusses feasibility of development process. l e a v i n g landing gear
linenrimtion in context witb control optimization of designers at both airframe and l a n d i n g gear
a semi-active nose landing gear. Both linear and manufacturers a t a loss to re-establish a m l y
d idifferential equations are developd at [61 vibdon-free landing gear. This again empbasiiacs
for tbe wobble (shimmy) of a motorcycle front tbc necessity ofprrpsntioa and as m l y as possibly
wheel: nonlinear effects arise from inclusion of application of qualified and commonly accepted
wheel unbalance and periodic tyre radial force methods of design for vibration-free landinggw.
vpriation. In studying shimmy having occurred in
f l i g h trials, nonlinear efkts were also accounted
for by [7lin time domain analyses. Paper [8] p- ROUND TABLE DISCUSSION
sents a complete set of nonlinear equations of
unsteady tyre dynamics. Paper [9] discusses Thc Workshop was concluded with a Round Table
influence of nonlincarities on shimmy, presenting Discussion kd by Dr.R. Fnymum. The discussion
models f o r typical nonlincarities as stick-slip was characterized by deep concern about present
friction and freeplay. hxldon of aoaliaearitics in a state and f u m development of "TEEDESIGN,
simulation model of a cantilevered landing gear QUALIFICATION AND MAMTENA" OF
g.nrrey is prrscntcd at [IO]. VIBRATION-FREE LANDING GEAR'.
C o m p a r e d w i t h e f f o r t s p e n t on a e r o -
From a scan of the sub@ of d i t i e s , p v n t (servo)elaaticity it is obvious that the field of
author got the impression that there is widespread l d i n g g w r v i b r s t i o a s has bccn negkctadbadlyfor
concern about the reliability of linearized landing a long time. Tbcro was also snse discussionabout
gear models. However, increased modelling e f f m the reasona therefor. Nowcver, in view of the cost
a n d cost o f c o m p u t a t i o n appear to i m p e d e involved in correcting landing gear vibration
application of nonlinear moder8 except in cases problems emerging late in the design p r o m s or
when shimmy occurred on existing hardware, even not before aircraft delivery to service,
notwithstanding cases whcre linearization is imp- diseussion turned to future. It became clear tbat the
PrOpriaD (landing, braking). state of the art of mating leading gear vibrations
was not quite up to olber similar subjects and that
Interesting enough to present author none of the cooperation was ncGded of expem widely diatrib-
papers mentioned shimmy under nonsymmetric utEd across aircraft, landiig gear and tyre industry
b a s i c conditions, e.g. a wheel running a t a (the latter was not represented at the workshop),
geometrically or elastically induced pre-set slip imtitutions, and Sum0ritie.s. Workshop participants
mgk ('be-in") or tilt angk. did not expect the rapid progress needed. unless
there wps a well planned and c- appmtdr
to the problem. Unanimously. AGARD was
c o ~ M R Q s u a E s identified the only instituhn which bad enough
power and wtbority to promote such a c m r d b t c d
In spite of all effort spem on design of supposedly effatintheinurratofNATOmMnbcrn;ltiear.The
vibration-free landing gear, the workshop was format of a Working Group was considered most
pmented with some practical e&lmplcs of landing suitable to acbicve the pmgrcss mmssery within a
gear shimmy ([2]. [31. [71, [IO]). It was interesting reasonrbLtim9r.Fnymnrmwasaskcdtoconvey
to note that in most c m the 1andiig gem affected this proposal to tbe appropriate AOARD M y . It
cmld not really be cured. Rather in most instances was conaidemi useful to bave tentative 'Tmms of
the butden had to be put on maintaining close Rckrena"writt.n for this task. John O l w r of de
tolerances w.r.t. f w p l a y and tyre pressure, on Havilland. William E. Krabacher of Wright
costly early tyrc replacement, on frrsucnt checks of Laboratory, and Amulf Krauss of Daimler-Benz
tyre unbalance and out-of-round. and on pilots Aerospace were n o m i d for an informal ''Point
having to h e spccial openting proccdurcs. In of coruct"group.
somc instances there has also mass b.en added to
tbe l a d i n g par, either in form of landing gear
T-4
B. M. Crenshaw
Susan C. Brown
Lockheed Martin Aeronautical Systems
86 South Cobb Drive, Marietta Ga. 30063 USA
The differences in requirements of land based (LB) and g acceleration of gravity -in./sec. squared
I
carrier based (CB) aircraft landing gear are reviewed L aerodynamic lift (per strut) - Ib.
with respect to landing impact and ground surface N load factor - nondimensional
roughness. Frequently the issue of operational S strut stroke - in.
roughness requirements vs. taxi requirements arises. V sink speed - in./sec.
More often, it seems, attention is being focused on W weight per strut - Ib.
operational runway roughness requirements. The MIL- 8, tire deflection - in.
SPEC roughness amplitudes and wavelengths may not 7 strut efficiency - nondimensional
represent the operational capability of an aircraft when 7 tire efficiency - nondimensional
combined sources of loading are considered. There are
severe loads on both main and nose gears during 1. INTRODUCTION
landing rollout on a rough runway surface if braking is
used; however, levels of braking and roughness An effective general approach to minimizing vibration
combinations are not always clearly defined in effects on landing gear and preventing transmission of
procurement specifications. ground generated high frequency loads into other
aircraft structure is to utilize the lowest practical tire
The more robust landing gears sized for carrier pressures, and to keep the strut airspring constants low.
operations are examined to determine their potential To achieve these goals, thorough preliminary analysis
- operational performance for land based roughness is needed to avoid later fixes which can make the
levels when combined loads from rollout braking are gears marginal energy absorbers for ground operations.
considered. While the weight penalty associated with
carrier qualified landing gears is commonly recognized, All too frequently, landing gears and their
and weight efficiency requirements of fighter aircraft retractionlextension mechanisms become overly
may ultimately outweigh cost considerations for complicated or compromised because the implications
commonality of landing gear components, it is of structural requirements were not studied thoroughly
nevertheless worthwhile to consider methods for enough during the initial aircraft layout to lay claim to
reducing costs through multiple application designs and the space and volume necessary to accommodate all
parts usage. Although without great previous success, landing gear needs. Sometimes inevitable aircraft
multi-service application of designs has long been an growth or mission changes can require a tire size
attractive concept. If severe runway roughness change or stroke adjustment, only to have an
capability is considered as an operational requirement interfering bulkhead preclude a straightforward change.
for land based gears, the weight obstacles to common In the tightly packed wheel wells of fighter aircraft,
landing gears may diminish. any unconservative gear sizing may even prevent such
changes as slight castedcamber adjustments to improve
A hypothetical future fighter aircraft is utilized to tire wear life or trail adjustment to optimize shimmy
compare runway roughness capability of gears sized for stability. It is therefore important to exercise all
carrier landing and arrestment loads with gears sized available computing and simulation techniques early in
exclusively for land based use, to compare relative the aircraft preliminary design stage to ensure sufficient
landing gear weights, and to develop concepts of gear stroke and wheel well clearances to support the
landing gears with the potential of multi-service usage. aircraft design maturing process without undue cast
A simple forward retracting vertical post strut penalties to the landing gear or the necessity for
arrangement has been examined from a conceptual cornpromise solut ions.
standpoint for low cost production in two stroke
lengths. For a long period, landing gear load/stroke
characteristics were determined primarily by landing
impact energy absorption requirements. Taxi loads
Paper presented at the 81st Meeting of the AGARD SMP Panel on The Design, Qualification and Maintenance
of Vibration-Free Landing Gear, held in Ban& Canada from 4-5 October 1995, and published in R-800.
1-2
Once preliminary stroke values are obtained, tires must runway roughness" but does not mention braking. In
be chosen. Reference 2 suggests that tires should have addition, an alternative 2.0 x MLG static load and 3.0
a deflection capability from static to fully compressed x NLG static load criteria is allowed by MIL-A-8862.
equal to the step bump height to be traversed; however, These alternate load levels allow very little energy
this is only practical for steps and short wavelength absorption reserve to be used for braking.
bumps. For longer roughness wavelengths with their
large amplitudes, it is believed to be reasonable to MIL-A-87221 is also considered ambiguous because it
choose a tire with a total deflection capability of specifies taxi analysis without braking in the taxi load
around 40 percent of the strut stroke as used in section, but requires consideration of runway roughness
Equation 4. in the braking section. A review of the structural
design criteria for several modem fighter aircraft did
During preliminary designs, it is usually sufficient to not reveal any allowance for combined braking with
consult tire manufacturers' handbooks and select an taxi loads.
existing tire if possible. For single tire struts, the tire
must have the required maximum static load rating. A I t should not come as a surprise, with the pressure to
little more freedom exists for multiple tires in that minimize landing gear weight, that gears designed to
there may be more sizes available which together will these specifications provide little braking capability
have the necessary load rating. It is desirable to utilize over the defined runways, and "operational" roughness
as much tire deflection capability as practical, yet keep capability today remains largely an undefined quantity.
the volume occupied by the tires small enough to fit in It is recommended that procurement documents in the
a minimal size and crowded wheel well. future provide more specific definitions as to the level
of "operational" runway roughness capability required.
For this preliminary study, an F-111 main tire was the
only listed size found with both reasonable tire Any additional roughness capability above that
diameter and the required static load rating for the achieved by current land based aircraft must come with
main gear. Dual C-130 size nose tires were selected, some weight penalty. The next section of this paper
keeping in mind the heavy nose gear loading resulting attempts to quantify this penalty by comparisons with a
from the carrier loading condition of free-flight gear sized for MIL-A-8863 carrier landing criteria
engagement by the arresting cable. To maximize based on a hypothetical aircraft. This comparison is
roughness amplitude capability of the study gears, the done only at a preliminary design level and is intended
same tire sizes are used for both carrier and land based to help quantify the potential weight penalty. It
designs. remains a customer's decision to accept higher landing
gear weight to achieve increased operational capability.
3. RUNWAY ROUGHNESS
4. AIRCRAFT AND LANDING GEAR
Once strut stroke lengths and tire sizes have been CONFIGURATION
selected, runway roughness response must be
addressed. The procurement authority specifies A hypothetical aircraft has been configured for this
roughness requirements, usually in terms of existing landing gear comparative study. The aircraft is not
military specifications. In the future the requirements representative of any existing aircraft; however, an
will likely be those in MIL-A-87221. Figure 2 effort was made to achieve realistic proportions,
compares portions of these specifications for (1-cosine) weights, and size anticipated for future fighters. In
discrete bumps for surfaces defined as either prepared addition to landing gear characteristics, the required
or semi-prepared at wavelengths up to 400 ft. Further parameters are landing and takeoff weight, pitching
breakdowns in speeds below and above 50 knots are moment of inertia, wheel base, and main gear tread
are made in MIL-A-87221; however, only speeds width. Tread width was chosen to ensure turnover
above 50 knots are considered in this study. Takeoff requirements are mel. An artist's concept of the
or landing dictates responses for higher speeds; but at aircraft is shown in Figure 3.
50 knots and lower "taxi" conditions, speed limitations
are feasible for controlling roughness induced loads. Two single tire, forward retracting main gear struts are
Table 1 summarizes several requirements in addition to used. A dual tire nose strut is used for compatibility
amplitude and wavelength. with carrier nose gear catapults. These configurations
are illustrated in Figure 4. Both struts have simple
Of the three specifications considered (References 3 - 9 , forward retraction mechanisms, and are mounted
only MIL-A-8863 explicitly requires combined braking laterally between spherical bearings to reduce bending
simultaneously with loads from the rough surfaces; loads transmitted to the airframe. On each gear, one
however, it does not specify a level of braking to be bearing is locked. The lateral positioning of the side
used. MIL-A-8862 defines the profiles as "acceptable brace should isolate any spin-up/springback and
1-4
antiskid vibrations from transmission into the aircraft Figure 8 shows drop test simulations for nose and main
structure through installed bearing tolerances. No gears and free flight engagement simulation loads for
details were considered for catapult attachment. the nose gear. Since aerodynamic characteristics of the
Pertinent aircraft data is shown in Table 2. hypothetical aircraft were not defined, a conservative
tail down angle of 16 degrees was used for the free
Figure 5 shows the struts in more detail. Gas separator flight engagement simulation.
pistons are used in each strut to minimize any pressure
loss potential from gas going into solution in the 6. LANDING GEAR GROUND ROUGHNESS
hydraulic oil. There is no significant pressure CAPABILITY COMPARISONS
differential across this separator to cause leakage.
Although metering pins are shown, computations were Evaluations of the roughness capability of these gears
carried out assuming a constant orifice size. Metering during landing rollout have been based on techniques
pins should offer some improvement in efficiency, with developed for evaluating existing aircraft to operate
resultant lower loads, but would not be expected to from bomb damage repaired runways as reported in
affect relative comparisons of gear capability. Reference 7. In this method, a matrix of aircraft
Rebound damping using a flap valve below the strut velocity and roughness wavelength combinations is
upper bearing is used in all gears. Rebound damping simulated and the maximum landing gear vertical loads
generally follows recommendations given in Reference are tabulated. For purposes of this comparison, the
6. runway roughness amplitude is defined by the
equations of MIL-A-8862 that express the roughness
Both main and nose struts are assumed to be of the amplitude as a function of the wavelength. Velocity
same diameter for computational convenience. For a increments of five knots were used from 50 to 150
given load/stroke curve, the strut diameter does not knots, along with wavelengths from 50 ft. to 400 ft. in
determine vertical load, it establishes strut internal 50 ft. increments. MIL-A-8862 specifies use of single
pressures. and two bumps. In this study, two bumps were
assumed. Also only dips were simulated as dips are
5. LANDING LOADS usually found in other studies to produce higher loads
than corresponding bumps of equal amplitude.
Landing impact loads for land based gears are
established for 10 fps sink rates. Sink rates for carrier From the resultant 3-dimensional table of velocity,
landing gears are based on the multivariate analysis wavelength, and maximum gear load, contours of
outlined in MIL-A-8863. For the study aircraft, this constant load level can be drawn, producing a graphical
analysis results in a sink rate of 24.5 fps as shown by description of load trends. Contours of loads equal to
Figure 6, based on an approach speed of approximately or higher than limit load define combinations of
125 knots. Besides the landing impact loading, roughness wavelength and airplane speed that must be
allowance must be made for tire encounters with avoided. For full capability over the entire
carrier deck obstructions. Thus, additional tire speed/wavelength range, no load contour can exceed
deflection capability equal to the deck obstruction the limit vertical load of the landing gear. If limit load
height needs to be provided. exceedances are found, the bump amplitude can be
reduced until the gear loads are within the limit load
The main gear drop weight is calculated for an aft c.g. envelope. To minimize the number of computer
at design landing weight. The nose gear drop weight is simulations required, the amplitudes have been reduced
calculated for the most forward c.g. plus a forward by an equal percentage over the entire wavelength
acting c.g. acceleration of 10 fps squared, to account range studied. Thus when the results are discussed in
for pitching following impact with a tail down aircraft terms of a percentage of the specification roughness
attitude. amplitude, the quoted percentage represents the point at
which the highest load occurred, and more capability
Carrier nose gears are very highly loaded by the "free- exists for other noncritical combinations of speed and
flight engagement" condition where the tail hook wavelength.
catches an arresting cable at the instant of main gear
contact while at low sink speed and a high angle of Initial simulations were made using conventional static
attack. For this design point, the nose gear must react to compressed 3:l compression ratio strut sizing. The
impact loads plus absorb a large amount of aircraft land based gears sized for the alternate 2g/3g criteria
rotational energy introduced by hook loads. could not achieve more than approximately 55% of the
sem i-prepared roughness requirement based on single
Figure 7 shows landing impact load/stroke calculations bumps and the carrier based gear could not reach 100%
for both land based and carrier based nose and main of the semi-prepared surface single bump requirements.
gears along with nose and main tire load deflections. The airspring curves were then recomputed using
increased inflation pressures and lower compression Figure 10. For quantifying the weights, it is believed
ratio values. Because of high preload pressures, these that the conventional land based design would be in the
modifications represent what is believed to be the range of 2.2% of aircraft weight. The 0.9% relative
maximum capability for conventional single chamber increment shown for the carrier based design is
gears with these stroke lengths. believed to approximate the additional landing gear
weight necessary for 25 inch stroke gears to achieve
Figure 9 shows strut capability for these gears where full takeoff/landing/RTO capability on MIL-A-8862
sufficiently high roughness levels have been used to paved runways and limited operational capability on
produce at least limit gear load over some portion of semi-prepared surfaces. Full operational capability
the speed/wavelength range studied. Unrestricted with braking over MIL-A-8862 semi-prepared surfaces
operations would require a lower roughness level with two bumps or over the MIL-A-8863 continuous
where no limit loads would occur. (l-cosine) bumps could not be achieved with these
struts as configured on the hypothetical aircraft.
Figure 9 (a) shows the land based 12-inch stroke gear
results without braking. The roughness level was 45% 8. REDUCING ROUGH FIELD LANDING
of the MIL-A-8862 semi-prepared surface. High load GEAR COSTS
levels were produced for both main and nose gears in
the wavelength range above 300 ft. Escalating costs of aircraft procurement requires that
all possible methods of cost reduction be examined.
Since only limited semi-prepared surface roughness For landing gears, this means design simplicity and
capability exists for the land base gear, braking commonality wherever possible. Cost and weight
simulations were made assuming only prepared surface optimization must be considered jointly.
roughness as shown in Figure 9 (b). A braking
coefficient of 0.4 was assumed. Main gear loads were Forging costs represent as much as 5 - 10% of total
well below limit, but nose gear loads did exceed limit cost of modem landing gears. Applying design
at the highest speeds around a wavelength of 100 ft. initiatives to maximize commonality of forgings
between land and carrier based gears, and possibly
As configured in this study, the carrier sized 25-inch between nose and main gear components can reduce
stroke gear could not achieve 2-bump capability to production costs. To do so, the general arrangements
100% of the MIL-A-8862 roughness level at all of the land based and carrier based gears should be as
speeds and wavelengths. Figure 9 (c) shows that limit similar as is practical. Using the same basic geometry
loads are reached in the 300 ft. and above wavelength and trunnion attach points can facilitate commonality
range at 75% of the semi-prepared surface roughness of components and maximize interchangeability.
level without braking. The carrier sized gear does have
substantial braking capability over roughness as Major structural components such as the shock strut
compared to the land based design. Figure 9 (d) shows cylinder and piston could perhaps be machined from
that limit loads are just reached at 50% of the MIL-A- common forgings. Forging dies cut to produce
8862 roughness level. forgings for a 25 inch stroke gear could have excess
material removed when producing 12 inch stroke gears.
Undoubtedly, additional changes could provide some Common strut diameters could be achieved by
improvement in the roughness tolerance of these gears; machining to different wall thicknesses as necessary for
however, it is not believed that the changes would higher carrier landing loads. Navy requirements now
affect the relative capabilities of the carrier and land indicate Aermet 100 material, and this steel is
based designs. One potential method of obtaining presently a significant cost factor. Aermet 100 material
additional roughness capability is the use of double removed during machining can be recovered to recoup
chamber struts. In recognition of objections to double some raw material cost.
chamber struts by maintenance and gear overhaul
personnel, this survey was limited to single chamber Minor structural components such as drag and side
gears only. braces may also be made from common forgings, again
being sized for the heavier gear design. Allowances
7. WEIGHT COMPARISONS for wider lugs should be made during the design
process. Common trunnion pins may be used with a
It is not possible to estimate accurate landing gear larger inner diameter for the land based gears and a
weights until all load components, including drag, side, smaller diameter for the carrier based gears. Retract
and miscellaneous loads are known and the gear design and steering actuators may also be made from common
has progressed at least through the detail layout task. forgings, allowing for larger piston areas for heavier
Consequently, parametric weight estimating techniques gears.
have been used to compare these gears as shown in
1-6
Many nonstructural components may be designed to be maximum stroke lengths and tire sizes that could be
interchangeable from carrier based to land based anticipated for future, perhaps multi service usage of
designs. By using equal piston diameters, common the airframes.
components such as separator pistons, gland nuts, and
strut bearings may be incorporated. Also, system type Future studies should focus on the inherent capabilities
components such as steering valves may be completely of gears designed for shipboard operations to establish
com mon . benefits for use on rough runways.
To benefit from design commonality, both carrier based Additional work is needed to rapidly quantify landing
landing gear requirements and the maximum gear weight as a function of required capability through
operational runway/taxiway roughness should be improvements in automated loads and stress analysis.
considered at the beginning of the design process. Modeling tools of this type must be available and
History has shown that landing gears initially designed utilized from the earliest design stages.
for land base use are impractical for carrier operations
and entirely new designs are required. If carrier It is recommended that cost reduction methods,
requirements are considered initially, the landing gear including simplified gear retraction geometry, be given
may be more easily adapted to multi role uses in the high priority in the initial gear design. Sizeable c a t
future and at the same time much higher roughness savings might be realized, both in development and
tolerance potential can be achieved. manufacturing with more progress toward standardized
gear arrangements, and by reduced gear complexity
9. CONCLUSIONS which enable the use of common parts and forgings.
Specifications for landing gear design ground loads are 11. REFERENCES
not always explicit in describing operational
requirements on rough runways. The lack of clarity of 1. Currie, Norman S., "Aircraft Landing Gear
operational requirements with respect to braking levels Design: Principles and Practices," AIAA
while on rough surfaces can lead to designs with Education Series, American Institute of
reduced operational capability, i. e. limited to unbraked Aeronautics and Astronautics, Inc. ,
taxi or free rollout where the possibility of roughness Washington, D. C., 1988, p.35.
resonance exists.
2. Williams, W. W., Williams, G. K., and
Navy carrier based aircraft specifications (MIL-A- Garrard, W. C. J., "Soft and Rough Field
8863) require landinghakeoff operations from lower Landing Gears," SAE Paper 650844, October, - -
amplitude rough surfaces than the corresponding MIL- 1965.
A-8862 surfaces; however, additional capability is
inherent over corresponding Air Force gears from other 3. Airplane Strength and Rigidity, Landplane
requirements associated with Navy shipboard Landing and Ground Handling Loads, MIL-A-
operations. When severe roughness combined with 8862A, MarchJ971.
other load sources is considered, weight differences in
Navy and Air Force landing gears may diminish. 4. Airplane Strength and Rigidity, Ground Loads
for Navy Acquired Airplanes, MIL-A-
Airfield roughness capability of two landing gears with 8863B(AS), 6 May 1987.
different stroke lengths using MIL-A-8862 roughness
specifications have been compared for a hypothetical 5. Aircraft Structures, General Specification For,
aircraft. The landing gears were sized by Air Force MIL-A-87221, 28 February, 1985.
and Navy Landing Impact Criteria. The landing gear
weight penalty to provide operational capability with 6. Conway, H. G., "Landing Gear Design,"
braking to Air Force gears has been estimated from Chapman and Hall Ltd., London, 1958. , p.
this comparison at approximately 0.9 percent of aircraft 184
weight.
7. Anon, "Aircraft Operation on Repaired
10. Recommendations Runways," Report of Working Group 22 of
the Structures and Materials Panel,
Landing gear design and user technical personnel need NATO/AGARD, Paris France, AGARD-R-
to improve, clarify, and develop realistic ground 731, August, 1990.
roughness performance requirements definitions to
avoid potential problems in landing gear/vehicle
integration. Wheel well sizing should be based on the
1-7
0 5 10 15 20 25
-
SINK RATE FT.lSEG.
30
f
25
/
4
10
4m
5
0 I I I I I I I
0 50 100 150 200 250 300 350 400 0 50 100 150 200 250 300 350 400
TABLE 1
I
High Catapult and Arresting Loads
I 15Degrees I
I Turnover Angle (Forward C. G.) I 53.819 Degrees I
I Minimum Static Nme Gear Load I 8,121 Lb. I
I Maximum Static Nose Gear Load I 18,727Lb. I
- Nose Gear Load with Braking (10 ft./sec? Deal.) 28,211 Lb.
Front View
UPPER BEARING
SEPARATOR PISTON
SPACER \
LOWER BEARING
7,
7
X
20
b.- 20
4 15
w
$
f? 10
2
f n 5
0
0 5 10 15 20 25 0 5 10 15 20 25
STRUTSTROKE - IN. STRUTSTROKE - IN.
(a) Nose Struts (b) Main Struts
12
10
0 5 10 15 20 0 5 10 15 20
TIRE FORCE - LB. x io4 TIRE FORCE -
LB.X lo4
(c) Nme Tire Curves (d) Main Tire. Curves
1.2e+5 1.6~5 I
Free-Flight
Engagement 1.4et5
1.0~5
1.2~5
8.0e+4 m
J
w 0 1.0e+5
n
2o 6.0~4
4
9 8.0et4
U
5 +
E 4.0~4
2+ 6.0et4
v) ..
____....___. v)
4.0e+4
2.0~4
2.0et4
O.OetO
0 5 10 15 20 25 0 5 10 15 20 25
400
\I -i 400
.!2
-
300
350
250
-
-
I@(
u
hl
84 5' 1' I
&
z
350 -92K
300 -
250 -
4 -2
- v)
200 - Nose Limit Load iz 200 - Main Limit Load
v) 57, Lb. 0 92,000 Lb.
W
$ 150 -
3
$ 150-
? U
g 100-
3
g ioo-
50
I I I I I I I I I 50 . I. I. I. l. i. I. I. i. I.
50 60 70 80 90 100110120130140150 50 60 70 80 90 100110120130140150
AIRCRAFTVELOCITY - KNOTS AIRCRAFTVELOCITY - KNOTS
Nose Gear Main Gear
Figure 9 (a). Land Based Gear Capability over Two Bumps at 45 Percent MILA-8862
Semi-prepared Roughness.
1-14
kI 40
-I
350
300
k 350
400
300
4 4
$
-
Nose Limit Load Main Limit Load
250 92,000 Lb.
57,000 Lb.
U
i? 200
I
:
w
I
w:
$ 150 $ 150
8 100
9g 100
50 IIII(IIIIII 50
50 60 70 80 90 100110120130140150 50 60 70 80 90 100110120130140150
Figure 9 @). Land Based Gear Capability over MILA-8862 PrepaRd Surface Roughness with
0.4 Braking Coefficient.
400
k 350
300
4
$ 250 Second Contour: 150K ..
a Nose LimR Load
v)
200 200
Main Limit Load
3 w 180,000 Lb.
z 150 $ 150
g 100
50 11111111111
50 60 70 80 90 100110120130140150 50 60 70 80 90 100110120130140150
Figure 9 (c). Carrier Based Gear Capability over 75 Percent MIL-A-8862 Semi-prepared Roughness.
1-15
400 400
-
350 350
0
i
P 300 E
4 300 -
250 Second Contour: 90K
250 - Second Contour: 1501
d
$
W
200
Nose LImR Load
;
d
w
200 -
50 I I I I I I I I I
50 60 70 80 90 100110120130140150 50 60 70 80 90 100110120130140150
AIRCRAFTVELOCITY - KNOTS AIRCRAFTMLOCITY KNOTS -
Nose Gear Main Gear
Figure 9 (d). Carrier Based Gear Capability over 50 Percent MILA-8862 Semi-prepared Roughness witb
0.4 Braking Coefficient.
"a 1.0%
5w 1.5%
:
L
1.0%
0.5%
0.0%
LANDBASED CARRIu(-BASu)
WILLIAM E. KRABACHER
WL/FIVMA Building 45
WRIGHT LABORATORY
WRIGHT PATTERSON AFB
2130 Eighth Street, Suite 1
Ohio 45433-7542, USA
Paper presented ut the 8lst Meeting of the AGARD SMP Panel on The Design, Qualification and Maintenance
of Vibration-Free Landing Gear, held in Ban8 Canada from 4-5 October 1995, and published in R-800.
2-2
K,, - Torsional Spring Rate of the Steering Actuator critical to gear stability, to review some aircraft
about the Strut Vertical Centerline landing gear data sets indicating in one case the non-
(in-lb/rad) linear variation of the various input parameters as a
function of strut stroke, and to review mathematical
L - Mechanical Trail Length, defined as the model stability predictions for these data sets.
Normal Distance from the Axle Centroid to
the Strut Vertical Centerline (in) 2. THE MATHEMATICAL MODELS
p, - Tire Moment Coefficient (in-lb/rad) From the landing gear structural standpoint, the
two models to be presented are completely identical.
y - Attitude of the Gear, Positive Wheel Forward The only difference between the two models is that
(rad) one uses the Moreland tire model (2) and the other
uses the Von Schlippe-Dietrich tire model (3)(4).
Other symbols used in the analysis The mathematical equations of motion for the
Moreland tire model are given by
F, - Ground Reaction Normal to the Wheel Plane
Ob) 1. m+d%/dt2- (m - mJ+L+d20/dt2=
F, - Ks+y - K,+a - C,+dy/dt
T - Torque Moment about the Strut Vertical Axis
due to Torsional Spring Rate of the Structure 2. Is*dzO/dt2= T + (V/R)+I*dddt - H+\V
(in-lb) -
- FN*(L+ La) + (m m J + L * d ~ / d t z
-
W+(A - R+a)+sin(y)
Q - Tire Roll Rotation about the Fore-AR Axis gear in severe divergent shimmy. The two computer
through the Trail Arm analyses also output the 0. Q,, a,A or yo. and torque
T variables in similar graphical plots.
v, - Tire Torsional Rotation about the Vertical
Axis through the Axle Centerline
> -=
I
obtained from the two computer analyses based upon
the above two sets of differential equations, reference 30
is made to Figures 2, 3, 4. and 5 below. These four
-
graphs are plots of the y Lo@ variable for the cases
of a stable landing gear. a gear in l i t cycle
oscillation, a gear in divergent shimmy, and a
FIGURE 4
2-5
Increases Stability
- Decreases Stability
- Decreases Stability until I, = 2.25 and then
Increases Stability
- Greatly Increases Stability
- Decreases Stability
- Increases Stability
- Increases Stability
- Greatly Decreases Stability
- Greatly Decreases Stability'
- Greatly Increases Stability
- Decreases Stability until m = ,225 and then
Neutral Stability
- Increases Stability
- Slightly Decreases Stability
- Greatly Decreases Stability
FIGURE S - Greatly Increases Stability
- Slightly Increases Stability
- Slightly Increases Stability
- Decreases Stability
3. SENSITIVITY ANALYSIS OF PARAMETER
VARIATION 4. LANDING GEAR STABILW CRITICAL
PARAMETERS
About f h e n years ago, a sensitivity analysis was
conducted using some tighter nose landing gear data In one particular shimmy investigation involving
and the Moreland Tire Model version of the shimmy a trainer nose landing gear a problem was uncovered
analysis. While the reaults of this study are only that identified stability critical parametersu.
totally accurate for the aircraft analysed, it has been Specifically, the vendor supplied a complete set of
found that the general trends in this study have been
a valuable guide in analysing other nose landing gear. ' The some of the original values used in the
The data d e f d g this nose landing gear from a parametric study were accidentally lost "he
shimmy standpoint is presented in Table 1 under the values listed in Table 1 arc approximately correct.
column for Fighter A.' However, they should not ba used in the analysis in
an attempt to reconstruct the results cited here. The
D e f d g the values in this table as the nominal results of this study should be viewed simply as a
values of the gear. the parametric study was general indicator of how dynamic response may
conducted by inormumtally varying one paramctcr at change as a parameter value is altered
a tima and increasing its value until its value was 1.5
times its nominal value. The computer output from ' Since the nominal value of L for this landing gear
the variation of each parameter was reviewed for is 0.8 inches, the maximum value uscd in this
stability effects and graded according to the scale: scnsitivity study is 1.2 inches. Therefore, greatly
greatly increases stability, increases stability, slightly decreasing stability makes physical sense in this case.
increases stability, no change, slightly decreases Generally, positive trail values below 3 inches tend to
stability, decreases stability, greatly decreases be destabilizing.
stability. Using this scale the following results were
obtained. ' The term stability critical parameter is taken to
mean that this parameter will produce greatly varying
C - Increases Stability analysis prediction results for small changes in the
C, - Increases Stability parameter value.
C, - Slightly Increases Stability ' The actual experimentally measured data for this
C, - Increases Stability
C, - Increases Stability trainer is presented in Table 1.
C, - Slightly Decreases Stability
I - Decreases Stability until 1 = 1.6 and then
2-6
landing gear shimmy analysis data which was entered 5. REVIEW OF LANDING GEAR DATA
into a data file for the shimmy analysis. When the SETS AND STABILITY PREDICTIONS
analysis was executed, the analysis indicated the gear
would shear off in the rust quarter-cycle. (Refer to Table 1 below provides the shimmy data sets for
Figure 5 above for a plot of the actual dynamic two different tighter (Fighter B and Fighter C) and
response predicted by the analysis.) Since this level one trainer aircraft that have been analysed for
of unstable dynamic response had never been shimmy behavior. These three aircraft were selected
encountered before, it was determined that something
must be wrong with the data. After the vendor
supplied four additional data sets. all of which also
predicted the gear would-shear off in the first quarter-
cycle, an investigation was conducted by comparing
the values of the data in these five data sets with the
data sets from other aircraft nose landing gears to
determine the physical realism of the data. Using
the above introduced parametric study as a guide,
focus quickly narrowed onto the value of the
parameter K, since its variation greatly alters the
analysis predictions. The result determined was that
in all five data sets supplied the value of K,, relative
to the values of K, for other landing gears was in
error by being two orders of magnitude too large.
landing gear when the values listed in Table 1 arc This procedure was repeated at increments of 0.05
used in the shimmy analysis. An example of the degrees of torsional Geeplay from 0.5 to 1.2 degrees.
typical dynamic response of this landing gear is
provided in Figure 6 below for the 8 variable. (Note: As a comparison. a static analysis of the gear was
8 in this figure is in degrees.) It is noted that this also performed to evaluate the prediction variations.
gear is very stable as there is little dynamic response The static analysis was conducted by selecting the 11
to the initial 5 degree deflection of 0. parameter values at the static stroke position of the
gear and incrcmenting the taxi speed in 5 knot
increments from IS to 150 knots while holding these
Y I I 11 parameters constant. Again the shimmy speed was
taken as the lowest speed at which the analysis
indicated sustained limit cycle oscillation of the gear.
The comparison between the dynamic and static
analyses predictions are shown in chart 14. Also
included in this charr arc the actual experimentally
determined shimmy boundaries of the gear. The
lower experimental boundary shows when shimmy
onset occurred in the gear and the upper experimental
boundary indicates strong limit cycle oscillation
shimmy.
1::
l-
I -
P *-
I 1-
FIGURE 7 FIGURE 11
:
.
.I
1 Y
.
I
FIGURE 8 FIGURE 12
FIGURE 9 FIGURE 13
FIGURE 10
2-9
In comparing these graphs it is noted that as the Referring to Figure 15, it is seen that the lateral
taxi speed increases from 80 to 90 knots a stiffness values obtained for a landing gear rigidly
progressively longer period of time is required for the mounted in a fixture (in the chart, strut only
vibration to damp out. At 95 knots the gear has the (measured)) vary from one and a half to three times
onset of limit cycle shimmy. At 100 h o t s the limit the values of the lateral stiffness actually measured on
cycle shimmy is definitely established. From these a landing gear mounted on the aircraft. This
time histories it is seen that the shimmy oscillation comparison establishes the fact that backup structural
frequency is about 27.06 HZ at 95 knots. When the flexibilities of the airframe must be taken into
actual aircraft was tested, it was found that the nose consideration when the four stability critical stiffness
gear went into a sudden onset of shimmy at 92 knots values are determined.
and had a frequency of 27.5 HZ. Repeated testing of
this landing gear, yielded reasonably consistent In addition to this primary cause of the difference
repeatability at this shimmy speed and frequency. between analysis predictions, a secondary cause was
also established. The vendor's analysis assumed a
In the course of the shimmy investigation for value of 6% of critical damping as the magnitude of
Fighter C, one very important point came to the structural damping available. Whereas, based upon
light concerning the measurement of landing gear personal experience with the variablity of available
shimmy parameters. The vendor's original analysis structural damping, the author decided to assume that
consistently indicated that this landing gear was stable the maximum amount of available damping was 1%
at all taxi speeds. The analysis provided by the of critical. At the time, the main reason for this
author indicated the 95 knot shimmy event at 27.06 selection was a desire to select a damping value that
HZ. The question arose as to what was different would expose any potential sensitivities the gear may
between the two analyses. After review of the input have toward shimmying. Further support for this
data a primary cause of this difference was position was established after the fact when Denis J.
established. The stability critical parameters KT,K,, Feld published his results on an analytical
K,, and K, used in the vendor's analysis were investigation of landing gear shimmy damping (8).
obtained by rigidly mounting the landing gear in a The essential conclusion of the analytical
fixture and measuring the resulting parameters. investigation was that the amount of available
Whereas, the author's analysis took into consideration damping at any instant during a taxi run could vary
- - the softening effects of the landing gear mounted in a from 1% to 10% of critical viscous damping.
- flexible fuselage and reduced the stiffness values Therefore, at the beginning of any landing gear
supplied by the vendor. The amount of this shimmy investigation an assumption of 1% critical
- -
reduction was largely based on experience with other damping is a reasonable starting point.
landing gear.
Before departing from Figure 15, another subject
The reason for taking fuselage flexibility effects is in need of mention. In the discussion of landing
into consideration was based on some data obtained gear stability critical parameters above in section 4, it
in another landing gear shimmy investigation. was pointed out that caution should be exercised
when obtaining stability critical stiffness values
because the finite element analysis predictions for
U T L R A L S C R I N O RATE these values tend to error on the high side. The proof
S R U T ONLY IYEAWREDI
of this error is established in the above figure for the
strut only (calculated) case where it is seen that the
largest error in stiffness values occurs between the
calculated values and the strut mounted on an actual
airframe.
6. CONCLUSIONS
these models uses the Moreland Tire model and the TIRE MOTION AND WHEEL SHIMMY, Robert F.
other uses the Von Schlippe-Dietrich Tire model. Smiley NACA TR 1299 (1 957)
A classification of the various types of analysis
output was given which consisted of a stable gear, a 5. EMPIRICAL ESTIMATION OF TIRE
gear with limit cycle oscillation shimmy, a gear with PARAMETERS, Raymond J. Black, paper
divergent shimmy, and a gear with severe divergent presented at SAE Committee A-5 Meeting, 28 April
shimmy. The results of a sensitivity study based on 1980, Renton, Washington
Fighter A's landing gear data was given. The
parameters K,, K,, K,, and K, were identified as 6. THE EXPERIMENTAL MEASUREMENT OF
being stability critical in the analysis of gear shimmy. THE T-46 NOSE LANDING GEAR SHIMMY
The Trainer data set described a stable landing gear. PARAMETERS, William E. Krabacher, ASDENFSL,
The Fighter B data set described a gear with limit Wright Patterson AFB, OH (1986) (In-House
cycle oscillation shimmi at a frequency of 27 HZ Publication)
which was verified experimentally. Finally, the
I
Fighter C data set analytically described a gear also 7. A DYNAMIC SHIMMY ANALYSIS, William
with limit cycle oscillation shimmy at 95 knots with a E. Krabacher, ASDENFSL, Wright Patterson AFB
frequency of 27.06 HZ. Experimentally, this gear OH (1981) (In-House Publication).
went into limit cycle oscillation at 92 knots with a
frequency of 27.5 HZ. 8. AN ANALYTICAL INVESTIGATION OF
DAMPING OF LANDING GEAR SHIMMY,
In the course of the presentation a number of Denis J. Feld, SAE Technical Paper Series 902015
important points were made. First, the Fighter B data (1 990).
set indicated the nonlinear variation of 11 input
parameters to the shimmy analysis and suggested' that
dynamic modelling of the gear will significantly
improve the accuracy of the overall analysis C T I R E YAW C O E F F I C I E N T
predictions. The Fighter C data set established the
importance of taking into consideration the backup
fuselage flexibilities when determining the stability
critical parameters KT,K,, K,, and K,. This data set
also suggested that an initial assumption of 1% of
critical viscous damping is a good manner in which
to expose the potential senstivities of a landing gear x o .O1 f ; : : : ,
O ~ ~ : : : :
to shimmy. Lastly, the tendency of finite element 2 3 4 5 6 7 8 91011121314
analyses to error on the high side when calculating N L G V E R T I C A L L O A D X10*.3 LB
the stability critical parameters was established.
BIBLIOGRAPHY
1. AIRCRAFT LANDING GEAR DYNAMICS
PRESENT AND FUTURE, William E. Krabacher, C 1 T I R E TIME C O N S T A N T
S A E Technical Paper Series 931400 (1993).
4,
2. THE STORY OF SHIMMY, William J.
Moreland, Journal of the Aeronautical Sciences,
pp. 793-808 (1954).
- . . , , , . , . . . , , .
0 1 2 3 4 5 6 7 8 9101112
0 2 4 6 8 101214 NLG VERTICAL LOAD
LOAD X 10**3 LB XlO"3 LB
K3 TORSIONAL K4 C R O S S COUPLED
STIFFNESS S P R I N G RATE
c
aa 0 2 4 6 8101214 .- .
- - 1 3 5 7 9 1 1 1 3 1 5
NLG STRUT STROKE
(INCHES FROM FULLY NLG VERTICAL LOAD
EXTENDED) XlO'"3 LB
0 3
c
x 2
o h 3 q : ; : ,
0 2 4 6 8101214 0
0 2 4 6 8 10 12 14 16
NGVHTICAL LOAD NLG STROKE (INCHES FROM
X 10-3 LB FULLY EXTENDED)
2-12
5
a
U
8-
0- 0
0 2 -4 6 8 10 12 14 16
8.5 I . . . . . . . .
NLG S T R U T S T R O K E 2 4 6 8 10 12 14 16 18
(INCHES FROM FULLY NLG VERTICAL LOAD
E X T E N DE D ) XlO"3 LES
1-12 I
SHIMMY ONSET
TAXI SPEED - LOAD CURVE \ 7S T R O N G SHIMM'I
I
--DYNAMICANALYSI
c I \
6
i
2000 STATIC ANALYSIS
I
I
' 0
0 40 80 120 160 0 02 0.4 0.6 0.8 1 12
G R O U N D SPEED (KNOTS) TORSIONAL FREEPLAY (DEG)
3- I
Daimler-Benz Aerospace AG
Military Aircraft LME24
8 1663 Munich
Germany
Paper presented at the 81st Meeting of the AGARD SMP Panel on "The Design, Qualification and Maintenance
of Vibration-Free Landing Gear", held in Ban8 Canada from 4-5 October 1995, and published in R-800.
3-2
2.4 EXPLANATION
Rudder Rotation 3.55 The speed of the aircraft was slightly decreased to 32 kts. -
Fin 1 st Bending 7.02 A more detailed sketch of the pure shimmy, only the x-
and y- behaviour of the nose landing gear is depicted in
- ~
The nose landing gear strut is a single axle forward a) Change of trace line
retracting gear with an oil filled, nitrogen charged shock Comparison of the original nose landing gear on T-46A
strut. The gear was originally developed for the Fairchild aircraft shows that the inclination of training aircraft
Republic T-46A Trainer A/C. It was designed for a sink landing gear leg is the same. To change the inclination of
rate of 13 fps limit and 16.2 fps ultimate cases. To the nose landing gear either the actuator rod end or the
maintain the same trail of the nose landing gear as on T - attachments of the strut must be modified.
46A aircraft the original actuator rod end has been
modified by Daimler Benz Aerospace. The nose landing b) Influence of elastic fuselage modes
gear wheel and tyre is from the Pilatus PC-9 and Due to the shimmy frequency of 25.7 Hz it is assumed
delivered via BFGoodrich, Ohio, because the original I
that the eigen-modes of the fuselage, 1st( 12.6 Hz) and
wheel and tyre was not more available.
3-3
2nd (5 1 .d Hz) vertical bending are not involved in this fastest possible brake force rise at high speed. A fast
oscillation problem. brake force rise was considered to yield. at least in
theory. shortest possible landing ground run distance.
C) Tyre Pressure
Under the circumstances that the original tyre was not The test set-up "dynamometer"consisted of one main
available and the used tyre has a different characteristic, wheel with brake running within a heavy drum. where the
an increase of pressure could solve the problem due to rotational inertia of the drum was to represent respective
shorten the relaxation line. mass per braked wheel of the aircraft. This arrangement
not only gives a representative picture of aircraft
d) Installation of shimmy damper deceleration by the brake but also provides correct
The shimmy damper is a hydraulic unit which reduces kinetic energy to be dissipated by the brake. The brake
the tendency of the wheel to oscillate from one side to was actuated by representative hydraulic components.
side. This dampers are usually constructed in one of two Brake control was performed by one lane of the brake
general designs, piston type and vane type, both of which control computer, pilot's pedal input being replaced by a
might be modified to provide power steering as well as synthetic signal. Anti-skid functions of the brake control
shimmy damper actions. computer were implemented and active.
e) Mounting of massbalance
Tuning the eigenfrequencies of the landing gear 3.2 DESCRIPTION OF THE PHENOMENON
massbalance weights have to be mounted on an arm
forwards or rearwards on the strut. Massbalance is a At brake initiation ("pilot" quasi "jumping" onto brake
useful tool if there is enough space to mount the weights. pedal) there was not the expected crisp yet steady rise of
brake force with eventual subsequent anti-skid system
f) Landing gear torsion stiffness activity. Rather there was a sequence of rapid on/off
Possible effects of torsional backlash is a reduced switching of the brake which lasted for about 1.5 seconds
effective torsional stiffness of the strut. Generally there is at a frequency of about 6 Hz. Thereafter the system
an undue sensitivity by reduction in torsional stiffness. stabilised to a normal behaviour.
Taxi tests after increasing the tyre pressure reveals no Fig. 3-1 shows a principal sketch of those three test
shimmy. Fig. 2-6 shows this sufficient improvement measurement traces which are considered essential for
- -. because the nose landing g e a r is almost free of description and explanation of the phenomenon. From
oscillation. top to bottom, these are Wheel Speed. Brake Pressure
No further changes were investigated. (i.e. pressure on brake piston face), and Brake Torque.
- -
It should be noticed that worn and cold tyres made The very first increase of brake pressure is needed to
undercarriages less stable than it was with new tyres. overcome the piston lifting spring. To arrive at the net
This was reported by pilots during taxi tests and is likely pressure acting on the brake pad, that "spring pressure" is
due to the change i n tyre stiffness as the tread wears to be subtracted from the value measured: thus, at the
away. There is also a possibility that out of balance first three troughs in Brake Pressure trace, brake pads are
f o r c e s in the worn tyre induces high frequency practically unloaded. This becomes apparent in Brake
oscillations which nullifies the effect of friction in the Torque trace, where for all practical purposes torque
landing gear. reduces. to zero at the corresponding non-zero troughs of
Brake Pressure.
3. BRAKE TORQUE OSCILLATIONS AT The initial phase of Brake Pressure deserves a further
BRAKE INITIATION comment: While brake servo valve output pressure (not
shown) follows solenoid current with negligible delay.
The case following demonstrates that not only unsuitable this is not true for the brake pressure proper as long as
combinations of structural stiffness, damping, and brake pistons are moving to close the gaps between rotor
pneumatic tyre characteristics may lead to unexpected and stator disks of the brake package. The flow of
vibration problems on landing gears. Rather, an unlucky hydraulic fluid from servo valve to brake is restricted by
combination of brake system design with the peculiarities a quite narrow Restrictor orifice which limits loss of
of circumferential force development by a tyre can also hydraulic fluid in case of e.g. rupture of a flexible hose.
produce a serious vibration problem. The Restrictor leads to a pressure drop between servo
valve output and brake piston as long as the piston is
moving. This pressure drop is the larger the faster the
3.1 THE SITUATION pilot tries to actuate the brake. However, when fluid flow
stops, pressure drop across the Restrictor vanishes and
In a series of development rig tests for a fighter aircraft full servo valve brake pressure gets through to brake
brake system. some test, were scheduled to demonstrate
3-4
pistons within fractions of a second. Thus indirectly the between corresponding tread and ground points at exit
safety feature "Restrictor" is responsible for the from the ground contact patch. This "exit" distance de is.
extremely steep brake pressure rise following on the under stationary speed conditions, calculated to be
brake filling phase.
3.3 EXPLANATIONS
From measurement traces not shown here i t became where ground speed
apparent that oscillations observed were produced by vg
vc wheel circumferential speed
anti-skid system action. However, given the software and
the parameters for skid identification installed in the length of ground contact patch
IP
control algorithm, the anti-skid system worked correctly.
The real reason for that disagreeable and inappropriate Please note that circumferential force in this regime is
series of anti-skid actions at brake initiation is a dependent on a speed ratio rather than on actual speed.
mismatch of brake system control and tyre physics.
b) At increasing speed difference, adhesive friction
3.3.1 Physics of Tyre Braking Force capability will be exceeded in small portions of the
footprint mainly at the exit end of the ground contact
Explanations following aim for a basic understanding of patch. Within the footprint spots affected, rubber will
tyre mechanics with regard to developing circumferential start to slip relative to ground albeit at a very low slip
force. speed and for very short distance. This means that
circumferential force will still increase with increasing
When a tread element of a straight-rolling tyre enters the speed difference. However, rise rate will decrease with
ground contact patch it will stick to that ground element. friction limited areas covering increasing portions of the
On a tyre rolling exactly at the forward speed of the tyre footprint.
wheel axle, the position in space of tyre tread element
and ground element will be practically identical at entry c ) At speed difference approaching ground speed
to and exit from ground contact patch. Nevertheless practically all of the rubber in the footprint will be -
within the patch there will be slight distortion due to the skidding at nearly ground speed and for extended
tyre tread being forced from originally circular shape to distance. Since skidding friction depends on skidding - -
straight (ground) shape. Contact patch distortion speed and distance covered, circumferential force of a .
changing sign at half patch length leads to approximately non-rotating tyre (full skid) will be less at high ground
zero net circumferential force on the wheel. speed than at low ground speed. Furthermore, at fixed - -
ground speed circumferential force of a skidding tyre
Non zero net circumferential force develops when wheel will also be less than the maximum achievable at a lower
circumferential speed differs from axle ground speed. differential speed, where almost all of the footprint area
Due to speed difference the tread point and the ground may also be skidding yet at a higher coefficient of
point which coincided at entry to the contact patch depart friction.
from each other on their travel through the contact patch.
As long as the relative distance of the two points is not
too large. tread rubber will adhere to the original entry Circumferential coefficient of friction (pc) of a tyre
ground point: elastic deformation of the tread element usually is presented as a function of "Slip Ratio" (SR),
increases on its travel through the contact patch, thereby where
increasing the local circumferential load on the tyre. If
Va - v,
wheel circumferential speed is sufficiently less than SR = -
ground speed, at some point along the contact patch
shear force on the tread element will exceed adhesive
friction capacity of the element. As a consequence, the
This presentation compresses all tyre angular rates from
tread element will start to slide on the ground. Sliding
friction of rubber in quite complex manner depends.
-
rolling at ground speed (SR 0) to full skid (SR = 1 )
into an abscissa from 0 to 1. However. as shown above
amongst other parameters. on sliding distance covered
and sliding speed. In general it may gradually fall well pc with an increase of slip ratio becomes increasingly
below adhesive friction with increasing sliding distance dependent on ground speed. Hence there will be different
and speed. This admittedly coarse view of rubber tyre pc vs. SR curves for different values of ground speed.
mechanics leads to the following conjectures about Fig. 3-2 presents a qualitative picture of this speed effect
circumferential force on a braked wheel: on circumferential coefficient of friction.
Figures concerning tyre circumferential force are in
a) At low speed differences, circumferential force on the essence modifications of figures presented at Ref. 3-2.
tyre should be approximately proportional to the distance
3-5
Ref. 3-2 treats this subject in much more detail than it is During perfectly steady braking wheel angular
done in present paper. d e c e l e r a t i o n is proportional to a i r c r a f t l i n e a r
deceleration. Hence, on first view one could assume that
3.3.2 Conclusions with Regard to Sensible Brake any wheel deceleration beyond this value could be
Control interpreted as an incipient skid. In theory this assumption
holds only if the brake is operating at the maximum
3.3.2.1 Brake Control during Braking wheel torque transmittable to ground, because any loss of
transmittable wheel torque will drive the working point
Fig. 3-3 is an somewhat exaggerated variation of Fig.3-2. to the right (unstable) side of the curve (Fig.3-3 and 3-4).
However, the ordinate has been renamed "Wheel
Torque" after, in ideas, having multiplied circumferential If the brake operates on the stable side somewhere below
friction coefficient with wheel load and ground-to-axle maximum transmittable wheel torque a reduction of
distance. transmitted wheel torque will be followed by a "useful"
increment of wheel deceleration which brings transmitted
Best possible deceleration of the aircraft by wheel wheel torque back to brake torque. Therefore, anti-skid
braking could certainly be achieved if retarding moment action should only be triggered on exceedance of this
produced by the brake (!'brake torque") was just equal to "useful" wheel deceleration increment. As stated above,
the maximum driving moment achievable by the tyre admittable trigger level is zero if brake torque equals
("wheel torque"). However, this maximum point will transmittable wheel t o r q u e ; t r i g g e r l e v e l r i s e s
move around very quickly due to wheel load fluctuation (progressively) with the ratio of transmittable wheel
on uneven ground, due to fluctuation of ground to tyre torque to brake torque.
friction on varying ground roughness. due to variation of
ground-to-axle distance. and due to other effects (e.g. Assumed that a brake control system in its anti-skid
side load on a braked wheel). branch contains a fixed wheel deceleration trigger
criterion, the control system should also provide for the
Although brake torque is measurable directly, wheel appropriate torque reserve between brake torque and
torque is not. Hence it is virtually impossible to exploit maximum transmittable wheel torque. Since maximum
100 percent of achievable wheel torque throughout the transmittable wheel torque is not measurable directly,
- braked ground run. Therefore, brake control systems in anti-skid systems of the type considered here contain
general are aiming for a working point below the algorithms which reduce brake pressure output from the
- . maximum wheel torque on the left (stable) branch of the servo valve according to frequency and intensity of
. wheel torque vs. slip ratio curve. However, if by chance
maximum wheel torque falls below brake torque or if
previous anti-skid actions. Brake pressure will be
cautiously re-increased (eventually to the level
- - brake torque is increased beyond maximum achievable corresponding to pilot's command) if no more anti-skid
wheel torque, this statically unbalanced torque difference action was triggered in a sufficiently long time interval.
will decelerate the wheel and thereby increase slip ratio.
Provided torque difference is large enough and is acting 3.3.2.2 Brake Initiation
for sufficiently long time, slip ratio will be increased to
the unstable part of the wheel torque curve. Up to the point of brake initiation, the wheel rolls freely
at negligible wheel torque just balancing rolling drag.
In order to avoid this critical situation it is necessary in Therefore at brake initiation, rise of brake torque cannot
the first place to timely recognise an incipient skid and to immediately be counteracted by wheel torque. Rather,
reduce brake torque fast enough such that slip ratio is momentary difference between wheel torque and brake
kept on the stable side. torque leads to angular deceleration of the wheel. This
effects an increase of slip ratio and wheel torque. If brake
I
However, a backup procedure is needed in case that slip t o r q u e i s l i m i t e d to a v a l u e b e l o w m a x i m u m
ratio has eventually reached the unstable side of the transmittable wheel torque, slip ratio will eventually tune
wheel torque vs. slip ratio curve. If this happened the in to stationary balance of brake torque by wheel torque.
wheel would continue to decelerate to an eventual
standstill as long as brake torque is larger than actual Wheel deceleration occurs during brake initiation as well
wheel torque, even if maximum achievable wheel torque as during a skid. Hence, if during brake initiation wheel
had meanwhile recovered beyond actual brake torque. To deceleration exceeds the threshold set for skid detection,
recover from this "deep skid" the safest way is to lift anti-skid control will unnecessarily lift the brake. The
brake completely until wheel speed has returned to the wheel will recuperate to or nearly to freely rolling
stable part of the wheel torque vs. slip ratio curve. conditions. When brake pressure is re-applied by anti-
skid control another undue skid prevention cycle may
For better understanding of the case presented i t is occur.
necessary to explain the principal method applied here to
recognise an incipient skid: Wheel angular deceleration For sake of a vibration-free brake onset it appears useful
is the most important parameter used for skid detection. to gain insight into the factors influencing wheel angular
3-6
deceleration during brake initiation. To this end quite Fig. 3-5 shows simulation results of a braked wheel
basic considerations may be helpful: during brake initiation. The model used is quite simple. It
contains just one degree of freedom representing rotation
a ) Step increase of brake torque: If brake torque is of a wheel moving at constant ground speed. Brake
assumed to follow a step function then extremum wheel torque is assumed to form a ramp type function of time
deceleration will occur at t = 0 and will equal brake while K is assumed constant. i.e. brake torque is assumed
torque divided by wheel moment of inertia. to not exceed the linear part of wheel torque vs slip ratio
function. Eventual dynamic deviations of wheel torque vs
b) Creeping increase of brake torque: If brake torque is slip ratio from quasi-stationary behaviour were not taken
applied very slowly, this will lead to low wheel into account.
deceleration values as well. because the wheel is being
given ample time to adapt to the slowly increasing Numerical results from simulation confirmed principal
demand on slip ratio. considerations presented above. Hence a systematic
evaluation of simulations was performed and summarised
c) Influence of "brakes on" aircraft speed: In a diagram in Fig. 3-6. This figure presents kind of a design chart.
showing wheel torque versus slip ratio (see e.g. Fig. 3-3) With maximum design ''brake on" speed given it shows
the initial linear part is virtually independent of ground which brake torque rise rate is at best admissible if skid
speed. However. analysis of the equation defining slip detection threshold shall not be exceeded during brake
ratio. viz. initiation. This type of diagram can easily be set up for a
specific aircraft using its tyre and wheel characteristics
(mainly moment of inertia and wheel torque curve) in
c o m b i n a t i o n with b r a k e a n d a n t i - s k i d c o n t r o l
characteristics (primarily skid detection threshold).
indicates that wheel torque rise rate over time will reduce In developing Fig. 3-6 it was tacitly assumed that brake
with an increase of ground speed. torque shall not exceed the straight part of the wheel
torque vs slip ratio curve. However, in reality this is
Assumed wheel torque is linearly dependent on slip ratio, rather improbable, because maximum brake torque
i.e. achievable is mostly close to or even beyond maximum ~
It is also concluded that wheel deceleration may vary Assumptions taken for Fig. 3-6 are therefore on the
with tyre type. e.g. if different tyre construction and/or optimistic side. It is advisable to apply more realistic
tread material changes the slope K of the wheel torque vs assumptions (especially on wheel torque vs slip ratio) in
slip ratio, K may also change with operating conditions. developing brake and anti-skid control system algorithms
e. g. dry or wet ground (Ref. 3-2) and parameters.
Furtheron. wheel deceleration increases on transition Nevertheless, the simplified approach which was taken
from the quasilinear part of the wheel torque slope into here yields sufficient insight into the oscillation problem
the degressive part while approaching maximum observed during dynamometer tests of rapid brake
transmittable wheel torque. application at high aircraft speed.
3-7
From all insight gained into the problem there is no Sustained (e.g. self-induced) landing gear vibrations
solution to be expected from any kind of bolt-on often are considered to form an exotic and rather local
- structural damping device. class of dynamic stability problems. However, the
ramifications of a vibration problem on landing gear may
well extend across the entire aircraft.
- 3.5 POINTS OF IMPROVEMENT IN THE First of all. that vibrations can present severe problems
BRAKE SYSTEM with regard to dynamic loads on landing gear and its
- - attachments.
The brake torque oscillations observed in dynamometer Even if oscillatory landing gear loads do not exceed
tests at extreme "brakes on'' speed and at extreme step loads determined by quasi-static load cases. they may
input to the system can be avoided through improvement well become a problem with regard to structural fatigue
, by design as well as through improvement by operation. (e.g. Ref. 4.1).
Furtheron, landing gear vibrations may impair pilot's
Improvement by operation means that pilots are ability to control the aircraft.
instructed to initiate brakes such that brake torque rise Excessive aircraft on landing gear vibrations deteriorate
rate is kept below critical values. Improvement by passenger comfort.
operation is viable as long as an average pilot is not Early and comprehensive consideration of landing gear
overcharged by the operation required. At best. vibration is therefore an indispensable part in the
oscillations should not occur if the pilot just avoids development of a successful aircraft.
"jumping" onto brake pedals.
"Shimmy" is not the only but one of the most frequent
Nevertheless there are aspects of the problem rendering and most dangerous sources of landing gear vibration.
improvement by design a desirable alternative. Shimmy is an oscillatory wheel mode induced by a
mechanism similar to that of a wing in flutter.
Experience shows that aircraft get heavier in course The shimmy mode is important because it can lead to a
of time. As a consequence maximum "brakes on" dynamic stability problem.
speed shows tendencies to increase, too. If safety Primary source of the phenomenon is the flexibility
margins against oscillations shall be kept constant between tyre footprint and the wheel plane. Other factors
then brake torque rise rate and/or skid detection such as wheel speed, tyre geometry and torsional
threshold must be changed. If brake torque rise constraints have a contributing effect. From flutter
rate is a hardware constant then the necessary calculations we know that a component investigation
lifting of skid detection threshold will lead to show general existence of instability. Total aircraft
reduction of skid detection quality. In turn, analysis reveal always the same flutter mode a s
reduced skid detection quality reduces anti-skid component calculation but very often at lower speed.
3-8
Through interference of aerodynamic and interaction of Third. it turned out that the elastic fuselage modes are
coupled components new flutter modes are possible. not important , if the leg mode frequencies are well
Similar to this flutter behaviour it is believed that the separated from the aircraft mode frequencies.
eigenmode of the flexible aircraft can influence the
stability of landing gear.
The aim of the landing g e a r dynamic behaviour REFERENCES
programme is to obtain dynamic data to support the
design specification of an aircraft. To do this sufficiently Smiley R.F.
the c o d e must be capable to reconcile different Correlation, Evaluation and Extension of
subsystems (structure. flight control system, hydraulic Linearized Theories for Tyre Motion and
etc.) and improve data with test results. Wheel Shimmy
At the beginning of a project the data available for sub NACA TR 1299; 1957
systems are very limited. The analysis can be used to
obtain preliminary results by assumed or old data and Pacejka H.B.
can be updated with more detailed information during the The Wheel Shimmy Phenomenon
development. The analysis at Dasa is now organised as a A Theoretical and Experimental Investigation
series of individual modules for dynamic investigation with Particular Reference to the Non Linear
for components and subsystems. This system allow easy Problem
but time intensive upgrading of modules and trend Technical University Delft, Netherlands
I
I calculation by using generalised modes.
Schweiger J., Weiss F., Dobbs S.
Each full analytical model of the leg has about 50 Aeroelastic Design and Flight Test Evaluation
eigenmodes, depending on the configuration selected. of the RANGER 2000 Training Aircraft
The high order of the model is a direct consequence of International Forum on Aeroelasticity and
including all the features that are necessary to represent Structural Dynamic 1995 Manchester U.K.
the interactions with other sub-systems during landing
and ground roll simulation. Krauss, A. J.
Considerations on Optimality of Landing Gear
Lower order models are required to permit efficient time Arrangement and Design.
response evaluation. Reduction in modal order is Paper #8 from AGARD Conference
achieved by removing the dynamics with the higher Proceedings 484: Landing Gear Design Loads
frequencies modes. This can be achieved by dynamic
condensation of the mass and stiffness matrix. For this 3.2 Brewer, H. K.
application, with relatively small finite element models, Parameters Affecting Aircraft Tire Control
the method is preferable to that of R. Guyan, because the Forces.
dynamic condensation or also the residual flexibility AIAA Paper No. 74-966
methods preserves the mode shapes very well and
permits a simple procedure for the selection of modes. 4.1 McBrearty, J. F.
Up to the lowest five modes of the leg model. selected in A Review of Landing Gear and Ground Loads
the manner described above are retained by transforming Problems.
the equations into modal coordinate system. It is also AGARD Report 1 18
possible to determine the lowest eigenmodes on a free
Ground Vibration Test (GVT) with retracted [ 4.21 Gu yan,R.
undercarriage. With those derived modes a stiffness Reduction of Stiffness and Mass Matrices
matrix can be calculated, which corresponds to the mass AIAA Journal February 1965
matrix but is not useable for static calculations.
A damping matrix corresponding to the generalised co-
ordinates can be easily established. because each mode is
selectable. This generalised matrices are used as input to
the stability criteria of R. Smiley in connection with the
linearised tyre characteristic and the applicable ground
force.
Generally speaking it turned out that the most important
parameter in the shimmy investigation is the relaxation
length of the tyre, because changes of this length results
in large changes of the landing gear stability.
Second, the torsion (friction) damping of the piston
against the cylinder, which is derived from measurement
of the landing gear producer has also an imminent
influence to shimmy.
3-9
81 . b 3
66.60.
C G R O U N D LINE
WL 59.70
STA
144.40
-
Fig. 2 2 Nose Landing Gear
3-11
<
E
. A
-
Fig. 2 3 Time History Nose Undercarriage
1 second
., .. ..
I
1 ....; ...............;.............. j............... ~ ;_.__..._..
_._...........
.,. I , .
-
Fig. 2 4 Shimmy
3-12
X 1
-
Fig. 2 5 Power Spectral Density
.....
...............
....... ...........
f
....
-
Fig. 2 6 Time History after Modification
3-13 I
BRAKE PRESS&
BRAKE TORQUE
I second
U
.. Fig. 3-1 Measurement Traces from Brake Fig. 3-3 Brake Working Points in Relation to
Dynamometer Test Wheel Torque versus Slip Ratio
z I
I I I I I
0 0.1 0.h 0.6 0.8 1.0
U
Fig. 3-2 Circumferential Friction Coefficient Fig. 3-4 Wheel Torque Gradient K
Dependence on Slip Ratio and Speed versus Slip Ratio
W.
200.
100.
0
-
0.5 1.0 1.5 2.0
TIME
Rr
W TORWE RISK RATE
Fig. 3-6 Design Chart for Brake Torque Rise Rate
3-15
Fig. A-1
ANALYSIS AND CONTROL OF THE FLEXIBLE DYNAMICS OF LANDING GEAR
IN THE PRESENCE OF ANTISKID CONTROL SYSTEMS
E. Denti, D. Fanteria
Dipartimento di Ingegneria Aerospaziale, Universits di Pisa
Via Diotisalvi 2. 56126 PISA. Italy
Tel. +39.50.550200 - Fax t39.50.553654
SUMMARY
This work is part of research carried out to investigate the vibration problems due to coupling between rile flrsihle
dynamics of aircraji landing gear and the dynamics of antiskid control systems. The last objective is the dei~elopnrerrt
of new antiskiddevices for the suppression of landing gear vibrations during braking.
In the paper the above mentioned problems are described together with the models and methods available for their
investigation. Preliminary results of the research are also presented. Namelv, an analytical and nwirerical model ofthe
tyres in runway contact dynamics, that is one of the most critical models involved, is developed. A Iineuri,-ed nrodel
of the system composed of the wheel and the landing gear is developed and the system poles are imestigoted.
Paper presented at the 81st Meefrng of the AGARD SMP Panel on The Design, Qualification and Maintenance
of Vibrmon-FreeLnndrng Gear. held in Ban8 Canaak from 4-5 October 1995, and published in R-800.
0-2
3 . SYSTEM MODELLING
The dynamic system to be studied is composed of:
aircraft, landing gear, wheels, the braking system and the
antiskid control system, connected as shown in Fig. 3-1.
Starting from the simulation code described in [DENTI],
a computer program has been developed for the
numerical simulation of the system. This program is
based on the following preliminary choices. The
dynamics of the aircraft is represented by a 3 degrees of
freedom (DOF) rigid model. The elastic dynamics of the
landing gear is represented by a 1 DOF second order
linear model. The vertical and horizontal ground forces
are computed as functions respectively o f tyre
0 IO TIME (s) 20 deformation and slip coefficient; these functions are
Fig. 2-1 An example of the gear walk phenomena supplied to the program by tables as input data. The
dynamic model of the wheel spin takes into account the
horizontal ground force, the rolling friction and the
braking moment. The braking system model takes into
account only the servo-valve dynamics, represented by a
second order linear model, while the braking moment is
assumed to be proportional to the output pressure of the
servo-valve. The antiskid control box implements a
basic control strategy and it is possible to build other
strategies of the users choice. The time history of the
pilots braking command is an input datum.
All these models are discussed in the following.
govern the dynamic generation of the braking moments. wheel hub velocity ( V,) and spin velocity (a), that are
Friction is a non linear phenomenon depending upon both measurable quantities; R , depends upon V C .
temperature, with histeresis and random components. So vertical load on the hub ( F N ) and tyre features. In
a lot of work would seem necessary to improve the generic braking conditions, circumferential velocity ( V , )
simple model chosen. In any case, the way of modelling is defined as the product R R e , while slip velocity is
this system is not essential for the present work because given by the relationship: Vs = VC - V,. (Fig. 4-1 ).
it does not affect the dynamics of the wheel-landing gear
subsystem. However it becomes fundamental when the
antiskid control loop is eventually closed.
Antiskid control box. The antiskid system is a
dynamic box that allows simulation of both basic and
newly designed antiskid control laws.
Recent research on antiskid control systems mainly
concerns ground vehicles (cars, bikes and trucks). For
example new antiskid control laws are studied in
[LINDO, MIYAS, WATAN], where experimental and
simulation results are given, too, and in [TANCH,
YEH-I, YEH-21 from a more theoretical point of view.
In [HOHNO, M A W ] systems based on Fuzzy Logic and Fig. 4-1 Classical model
Neural Networks are presented. Practical data for car
antiskid device simulation are also available in Local sliding velocity and longitudinal global force are
linked to wheel dynamic variables such as V c , R , FN ,
[BOWMA, FANCHI.
In the aeronautical field the development of antiskid, braking torque (T) and time; so, in steady state
from the early antilock on-off systems (1950) to the conditions, we have p = p ( F N , T , V C . R ) w i t h
modem digital antiskid devices, is described in [CURRE, F , = ~ F N .Since F, , T and R are linked by the wheel
spin equation of motion, the previous relationship may
HIRZE] and NASA research on antiskid is reviewed in
[TANNE]. But the strategic importance of the matter b e written again 8s p = p ( F ~ . V c , f 2 ) o r
leads to a virtually total lack of recent research and p=p(FN,VC.Vs). In literature [CLARK, PACEZ.
detailed information and data about commercial antiskid SAKAI] global friction coefficient p is given as a
devices. So, an accurate simulation of these devices may function of FN , V c and slip coefficient s , ~ =Vs/Vc (or
of FN , Vc and theoretical slip coefficient U= V,y/V,.).
~
itself able to predict it; B elastic link between rigid belt & rim
3) belt and tread are considered to be flat in the C elastic link between rigid belt & rim + residual
contact zone and tread is undeformahle stiffness
orthogonally to the ground; A) An elastic ring tyre model was developed by
4) the belt velocity, in the contact zone, is assumed [SGONG, ZEGEL] to solve the problem of correctly
to be uniform and equal to slip velocity Vs; representing the high frequency vibrational phenomena
5 ) at free rolling ( Vs=O) the tread is longitudinally of the tyre belt. Particularly, the need for a model
undeformed and no shear stress exists; capable of simulating the observed Standing Waves
6) in generic braking conditions ( V @ ) , for the phenomena was felt [CLARK, HUANI, HUANZ]. The
adhering part of the contact patch, the horizontal model consists of an elastic ring, treated as a circular
shear stress is assumed to be linearly dependent elastic beam, a mass, representing the wheel rim, and
on the local longitudinal tread deformation, while radial and torsional springs, connecting the ring to the
for the sliding part, shear stress depends on the rim, whose aim is to reproduce tyre side-wall and
local sliding velocity ( V J L )on the basis of the pressurised air behaviour (foundation stiffness). The
dynamic friction law IpAcU]. ring is assumed to have the radius of the belt and ground
To discretize the problem the tread in contact may be loads are concentrated at the intersection between the
regarded as a brush (Fig. 4-2) whose bristles are ring itself and a straight line perpendicular to the ground
attached to the belt and move through the contact zone and passing by the undefomed ring centre. The model
with the relative velocity V,. At least one bristle is gives good results [SGONG, ZECEL] and makes it
always adhering at the beginning of the contact patch possible to take into account the fundamental elastic in-
and the adhering zone extension is determined by the plane dynamics of the wheel. On the other hand. if one
static fnction threshold. is not interested in high frequency modes, only
The brush model, whose fust formulation can be found translational and rotational rigid displacements of the
in [DUGOF](1970), gives good results and is widely belt with respect to the rim are really important and the
used, both in its simplest form [DUGOF. PACEZ] and model becomes too sophisticated.
with experimental or empincal improvements [SAKAI, B) In this case the ring modelling the tyre belt is
SAKAZ, SAKA3, SAKA41. linked to the rim by radial and torsional springs as in the
previous model but it is assumed to be rigid [CLARK,
SAKAI, ZANTI, ZANTZ]. The ring radius is assumed 10
be equal to the effective rolling radius [ZANTI, ZANTZ]
to fit the kinematic relationships at free rolling [PACE21
and the previous hypotheses about load introduction are
retained. The rigid ring model provides good results in
wheel dynamics simulation as demonstrated in [ZANTI]
and [ZANIZ].
C) Recently, an intermediate model has been proposed
in [ZEGEL].The first test results show a good agreement
with experimental data for medium range frequency
phenomena. The ring is rigid and it has the radius of the
undeformed belt, but, between the load introduction
Fig. 4-2 Brush model point and the ring. more elasticity has been introduced.
- Unsready Models of contact force generation are rare Such elasticity (referred to as residual stiffness)
in literature. Examples can be found in [ZANTl] and consists of vertical, longitudinal and torsional springs
[ZANTZ] where the model adopted for computer whose stiffnesses are calibrated on the basis of
experimental wheel elastic properties [ZEGEL].
simulation is an unsteady brush model, implemented
from an Eulerian point of view, which gives very . I ~ ~ ~ ~ ~ ~ ~
, ~
~ ,A, ~
c ::l ~
~ ~ , , ~ : ~
In Fig. 5.1-3 qualitative responses to a step in slip 5 . 2 Brush models and software structure
velocity ( V,) are shown for the three models presented
above. The K Model (Fig. 5.1-3 a) behaves in a rather The software represents a brush built with a set of
simple way even if it is quite "unphysical": as bristle bristles behaving according to one of the three
sliding happens, the sliding velocity V, suddenly previously described models. The basic hypotheses are
reaches its steady state value (V,), while the relative the same as those mentioned in sec. 4 for brush models.
displacement of the bristle shows a discontinuity due to with the assumption of a parabolic distribution of
the fact that the spring force (proportional to the relative vertical forces (see Fig. 4-2). The simulation is
displacement of the bristle) must balance ground friction performed by means of a Lagrangian approach: each
which decreases as VsL increases. The CK Model (Fig. bristle is followed as an independent dynamic system
5.1-3 b) avoids the discontinuity of the kinematic during its motion along the contact patch. The roots of
quantities, by introducing a V, dynamic transient the bristles are attached to the belt and equally spaced:
before the steady state value V, is gained. The MCK they move all together with the slip velocity ( V.7 ) rhat
Model (Fig. 5.1-3 c) shows very complex dynamics is one of the two inputs of the system, the other heing
whose effects in a dynamic brush model simulation are wheel centre velocity ( Vc). When a bristle exits from
not valuable a priori. the contact patch, a new bristle enters, at the same time.
0
1 7 :
Velocity
+ Position
xsL=xs (imposed)
Time
~
from the other side.
"
I A- - Vertical force ( FN ) = 3200N/m
- Footprint half length ( a ) = 0.08 m
a) Mk.d.d Time - Bristle stiffness ( K ) = 1.0 IO6 N/m3
t Velocity 0.4
0.3
0.2
0.1
I me
4 Position 0
-0.1
-0.2
-0.3 ..........
Force Time -0.4
-0.5
-10 -5 0 5 10
lis,,q.. ,q b) CK-Model Time Braking rest. Fig. 5.3-2 shows transient responses on
-
Velocity total longitudinal force Fx (per tyre width unit) to step
inputs in V,, with constant V ~ = 3 0m/s. The unusual
shapes of the transient responses are understandable if
~
0 I (Vs-Vsd
K(xs-xs~+C Time
on the intensity of the braking input: increases in Fv
values are due to the growth of bristle deformation in the
adhering part of the contact patch, while reductions in
F, values depend upon the increase in the extension of
the sliding zone. For small steps of sx ( bSr= 0.0303)
F, increases gradually to its steady state value because
the sliding zone is negligible. For greater steps, the Fv
time history shows an overshoot which becomes
Time
C) MCK-Model stronger as the steady state extension of the adhering part
of the contact patch decreases.
Fig. 5.1-3 Bristle responses From the steady state values of the time histories in
4-7
with 2Ka2
2,r=2a(l- vu); q = -
3p o FN
solving Eq.(I), the following function Fv= Fr(a)was
obtained:
-2( ~py]}.(3)
~=F,(3rov~(l-vu)~+p[3(11~)
In Fig. 5.3-4 the results of the analytical calculation for
V, = 5, IO. 1 5 , 2 0 , 2 5 , 3 0 m/s are shown. The 30 mls
nm.1.1
curves coincide with the ones obtained by numerical
Fig. 5.3-2 Transient response in longitudinal simulation. It is interesting that, when plotting fi as a
force to a braking step input function of slip velocity V,, at high slip values p
Fig. 5.3-2, p(u),~ ( 8 , )and p(V,) curves were built depends only upon V,.
for Vc=30 mls and are shown in Fig. 5.3-4. The shape
of these curves is fully understandable on the basis of . Transient response analysis. Fig. 5.3-6 shows the FI
responses for small steps of s, starting from different
steady state longitudinal force distributions along the
footprint shown in Fig. 5.3-5. The integral value of steady state braking conditions; tests have been carried
such distributions depends upon the extent of the sliding out at constant Vc=30 mls. These tests aim at
portion of the contact patch (where long. stress has a deepening system behaviour for small perturbations
parabolic shape) with respect to the non sliding one ahout steady state conditions, both at free rolling and in
(where long. stress has a linear shape). A maximum any braking stationary condition whatever.
Fig. 5.3-3Transient variation of longitudinal load distribution over the contact patch for a braking step input
4-8
Conclusions
The brush model presented, based on the bristle K
Model, is not able to adequately represent sliding
velocity VsL which, at the transition point, suddenly
0 1 2 3 4 5 6 7 8 changes from zero to the value of slip velocity V,y. This
unphysical behaviour does not fit in with experiments
Fig. 5.3-4 Friction coeff. versus slip c0eff.s that show (at least for small sx) a linear shape variation
and slip vel. in steady state cond. of VsL along the sliding zone [CLARK]. However all
the above does not seem to affect too much the system
The steady state values of F variations are of course response on the total longitudinal force which fits
different according to the local shape of the p(s,) curve experimental measurements quite well [ZANTI].
at the point corresponding to the initial steady state This may be due to the fact that for small s . ~values. the
condition. The shape of the transient responses on F, sliding zone is very small, while, for larger slip
are obviously linked to the time histories of load coefficients, the sliding zone is larger but the badly
modelled one is confined to a restricted region after the
transition point. In conclusion, the model seems to be
able to represent contact dynamics in a numerical
simulation of wheel behaviour.
Besides, the bristle K Model is simple enough to permit
the development of analytical models of contact
dynamics (see section 6).
Fig. 5.3-6 Transient variation of the longitudinal force for small braking steps
"
-= V , -RR, -RRe- at;
Jt ax
where c(x,t) is the longitudinal deformation of the
tread,(O; x,z)isaframeofreferencemoving with the
footprint (see Fig. 4-2) and the initial and boundary
conditions are:
c(.,o)=r,,(X,O); t;(O,t)'O. (6)
Linearization of Eq. (5) with conditions ( 6 ) (with the
assumption that R, is constant) leads to the equations:
de de
-
at
+
= v C - OJR,(I U,,) - R,, R,-
ax (7)
e(x,O) = 0 ; e(0,r) = 0
&. L a
7 1 0.k O.&
Drn.18,
o.& 0 . k o.&
where the symbols with subscript "0" indicate the
reference steady state condition quantities and the lower
case symbols indicate small variations around the
reference condition, i.e.:
e(x,t)iA{(x.t); Y C - A V C ; w * d R .
Close to the free rolling condition u,,=Oand the tread
slidingportion can be disregarded; so Eq.(7) is valid over
the whole footprint and can be solved, using the Laplace
transform technique, yielding the following result:
e(s.x) = Y S b){I-exp[-(xs)/(R, Q,,)]}/s (9)
where e(x,s) is the Laplace-transform of e(x.r) and
v~ (s) v,-(s) - R, w ( s ) the Laplace-transform of the
Tim,a,o.014 0.0?8 0.018 input. Now, with the same assumptions as in the K-
Model, by integrating the product K e ( s . r ) along the
Fig. 5.3-9 Transient variation of the long. force footprint, the following transfer function for the
for small steps in wheel speed longitudinal force (f,) can be obtained:
6 . AN ANALYTICAL MODEL OF
CONTACT DYNAMICS
with t9*2a/ReR,,. Fig. 6.2-1 shows the time
6 . 1 p(a) model response to a vs step obtained from Eq. (IO); it should
be noted that this response is identical to that shown in
If F, is expressed as a function of Vc and Vs by Fig. 5.3-2. evaluated numerically for an s,, variation
means of p curves, a linear model can be obtained as from zero to 0.0123.
4-1 1
1.2
I f ltlt I
1
IF
.,(.............
~
/ . I/
0.8 ................
.\ / \2 :
LANDING GEAR 2 = I
M~.S?+C~.~+K~;
WHEEL SPIN
DYNAMICS AQ=L(-AT-
SJ,, F,,,,,A~,.+AFJ~)
0 00 =o ROLLING
P= a,,Ks +(0"+l)K,
00 #O
FRICTION
0"
For the box of the contact dynamics, both the p(o)
4-12
model and the brush model (see sec. 6. I and 6.2) may be Dynamic Simulation and Load Prediction f o r
used. a Flexible Aircraft in Taxiing. Tuke-off and
In generic braking conditions, poles and zeroes of the Landing, AAAF - Proc. of the Int. Forum on
whole landing gear-wheel system, based on the brush Aeroelasticity and Struct. Dynamics,
model of the contact dynamics, have not been Strasbourg, v.2, pp. 839-856, 1993.
investigated in detail yet. When the p(o) model is [DOYLE] G . R . D o y l e , A Review <;f Computer
used, the presence of the wheel augments the period and Simulations f o r Aircraft-Surfcice Dynamics,
the dumping of the two complex poles of the landing AAIA Journal, v. 23, pp. 257-265, 1986
gear; a third real and stable pole, due to the wheel spin
dynamics, is also present and moves towards the origin [DUGOF] H . D u g o f f , P.S.Fancher, L.Sege1, A n
as (ap/ao)l0 increases, becoming unstable when Analysis of Tire Tractiori Properties arid Their
(a,u/ao)lo<O. When the brush model is used, at free Influence on Vehicle Dynamic Performance,
rolling one more pole, real and stable, due to the contact SAE Paper n o 700377, 1970, i n SAE
dynamics, is present. Transactions 1970.
It is interesting to note that, with the simplifying [ENRIG] J . J . E n r i g h t , Laboratory Siiuulation of
hypotheses assumed, it is possible to analytically Landing Gear Pitch-Plane Dynamics. SAE
express the poles of the wheel system as functions of Paper no 851937, 1985, in Selected Paper
system parameters. SAE PT-37, Ed. by J. A. Tanner, 1989.
[FANCH] P.S.Fancher, C.C.Macadam, Computer
8 . CONCLUSION Analysis of Antilock System Perfiormance in
the Braking of Commercial Vehicles, Braking
The review of the literature has shown that investigation
of Road Vehicles, I m e c h E C o n f .
into gear walk is difficult and few papers on specific
Publications, Published by Mechanical
problems of aircraft braking are available. Concerning
Engineering Pub. Ltd., 1976.
tyre modelling, a rigid ring model with residual
stiffness, for tyre structure, coupled with an unsteady [GERDE] J . C . G e r d e s , D.B.Maciuca, E.Devlin,
brush model, seems to be a good choice for numerical J.K.Hedrick, Brake Systeni Modelling ,for
investigation into gear walk. IVHS Longitudinal Control, Advances in
In the present paper a numerical unsteady brush model, Robust and Nonlinear Control Systems,
based on a Lagrangian approach, and an analytical model ASME, Dynamic Syst. & Control Division
of the tyre in contact dynamics have been developed. The DSC V. 53 pp. 119-126, 1993.
numerical model has been extensively tested providing a [GILES] N.A.Giles, Integrated Braking Systems, -
good qualitative understanding of contact phenomena. Aerospace Design Engineering, June 1993.
The analytical model is an extension of the ones in the
literature, as it is able to take velocity variations of the
[HIRZE] E.A.Hirze1, Antiskid and Modern Aircraft,
SAE Paper no 720868,1972.
wheel hub into account; this is essential when such a
hub is connected to an elastic landing gear leg. [HOHNO] H . 0 h n o , T . S u z u k i , K . A o k i . A.
Furthermore the poles of the brush subsystem in generic Takahasi, G.Sugimoto, Neural Netkrwrk
braking conditions are given. Cont ro 1 f o r A u torna tic B raking Con t ro I
System, Neural Networks, v. 7, n o 8, pp.
REFERENCES 1303-1312, 1994.
[HUANl] S . H u a n g , C.Su, In-Plane D~ncimics of
[BAKKl] E . B a k k e r , L.Nyborg, H.B.Pacejka,
Tyres o n the Road Based on a i l
Tyre Modelling for Use in Vehicle Dynamics
Experimentally Verified Rolling Ring Model,
Studies, SAE Paper no 870421, 1987, in SAE
Vehicle System Dynamics, v. 21, pp. 247-
Transactions 1987.
267, 1992.
[BAKK2] E . B a k k e r , H.B.Pacejka, L.Linder, A
New Tyre Model with an Application in [HUAN2] S . C . H u a n g , The Vibration of Rolling
Vehicle Dynamics Studies, SAE Paper no Tyres in Ground Contact, Int. J. of Vehicle
890087, 1989, in SAE Transactions 1989. Design, v. 13, no 1, pp. 78-95, 1992.
[BOWMA] J . E . B o w m a n , E.H.Law, A Feasibility [LINDO] W . C . L i n , D.J.Dobner, R.D.Fruechte,
Study of an Automotive Slip Control Braking Design and Analysis of an Antilock Brake
S y s t e m , SAE Paper no 930762, 1993, in Control System with Electric Brake Actuator, I
SAE Transactions 1993. Int. J. of Vehicle Design, v. 14, no I , pp. 13-
43, 1993.
[CLARK] S.K.Clark, Mechanics of Pneumatic Tires,
National Bureau of Standards, Depart. of [MAUER] G . F . M a u e r , G . F . G i s s i n g e r , Y.
"
Paper presented at the 81st Meeting of the AGARD SMP Panel on The Design, QualiJcation and Maintenance
of Vibration-Free Landing Gear, held in Ban& Canada from 4-5 October 199S, and published in R-800.
5-2
7 1 ProENGINEER
Level 1
A CAD I
I
NAsTRAN \ / MATLAB
ANDECS User
Code
Block
Level 2 SystemBuild
I I
Figure 1: SIMPACK and its Interfaces Figure 2: SIMPACK - MATRIXx Interface
One of the most important features of SIMPACK is the fact that Active control of vehicle suspensions is not only an issue in the
SIMPACK is an open system which possesses various inter- aircraft landing gear development. Cars, trucks and railway
faces to external standard software products (see figure 1) from vehicles are the major applications to profit from active suspen-
the domain of finite element (FEM), computer aided design sions, and great efforts are made among the respective manu-
(CAD) and control system analysis (CACE) programs. facturers to develop reliable and affordable systems.
The principles are largely independent of the specific vehicle . -
FEM: employment. Although passive suspensions, mostly consisting
To read in arbitrary flexible geometry a file interface to FEM- of a spring and a damping device, today have a high standard, -
programs was developed. This gives access to mass, stiffness they suffer from a disadvantage that lies in their principle.
and damping matrices and to the above mentioned terms of Optimized to isolate the vehicle body from the ground best at a
second order. Vice versa, loads computed by SIMPACK can be certain frequency range, their performance often diminishes at
transferred to the FEM-code. other frequencies. Additionally, the need to retain tire ground
contact at all times poses a design conflict with comfort
CAD: requirements.
For the purpose of incorporating physical and graphical CAD- For an aircraft the implications of the phenomena mentioned
data, SIMPACK can be linked by a function call interface to above is noteworthy: designed to absorb the energy of a hard
CAD-packages. This means an addition to the interactive SIM- landing impact, aircraft suspensions perform quite poorly in
PACK model setup tool. The link enables data consistency reduction of ground-induced loads during taxi and take-off.
between MBS- and CAD-data at each step of the model estab- Not surprisingly, supersonic aircraft and a new generation of
lishing process. In general this feature is not guaranteed by just stretched civil transport aircraft suffer the most under ground-
using file interfaces. induced structural vibration because of their increased struc-
tural flexibility inherent in their design with slender bodies and,
Control System Design and Analysis: at supersonic aircraft, their relatively thin wings.
Due to its ability to numerically linearize the system equations, Active and semi-active suspensions promise a solution to this
SIMPACK can be used as a simulation block within linear con- problem. An active suspension can be defined as a suspension
trol design tools. SIMPACK may also be linked as a fully non- layout which controls the forces acting in the shock absorber
linear block into nonlinear optimization and control tools like by control of energy dissipation, or, if required by the control
MATRIXx, MATLAB,ANDECS, etc. (figure 2 shows these law, by generating additional force. The usual means of control
links for the connection of SIMPACK to MATRIXx). is a closed-loop control with a control law acting with respect
This gives the user the possibility to use SIMPACK and its to measured states (often velocities or accelerations) at certain
parameter variation capability in a closed optimization loop to points of reference. A suspension termed semi-active has the
efficiently design dynamical systems. restriction that it cannot input energy into the system. Usually
realized as damping control, it is therefore only able to control
the amount of energy dissipation. Only forces of the same
direction as that of the instantaneous relative damper velocity
can be generated. Nevertheless, since fully active control sys-
5-3
tems usually require a heavy and costly force generation device A study of the properties of fully and semi-active oleos led to
(for aircraft landing gears mostly pressurized oil reservoirs the decision to take the semi-active damper as the system of
have been proposed) they are unlikely to be quickly introduced choice (see table 1) for the work presented here, see [IO].
in production aircraft. Semi-active suspension systems offer
considerable advantages of light weight and less complicated
mechanical requirements without suffering a great loss in per- passive semi-active active
formance. low weight low weight high weight
Summarizing, the main reasons for the introduction of semi- good performance good performance good performance
active landing gear control are: only for design over broad range over broad range
1, Minimize the load on the airframe structure, minimize case
force peak values and vibrations that can result in relative low com- medium com- high complexity
fatigue and reduce the life of the airframe. plexity plexity
2. Minimize accelerations acting on pilot and passengers
stable stable potentially un-
since the induced vertical and horizontal accelerations
stable
and vibrations can lead to passenger discomfort and
crew disorientation.
Table 1: System Features Comparison
1.3.1 Activities
There have been several preceding efforts to the improvement The semi-active damper seems to provide the best compromise
of aircraft of ground ride. A typical measure to adopt the land- between performance and additional weight handicap. The
ing gear force-deflection curve toward a lower slope is the use investigations of ground-load reduction are concentrated on the
of a two-chamber-oleo. In 1977, Somm/Straub/Kilner [4] pro- nose gear, since results show that this has the most potential for
posed an adaptive landing gear system for several military ride quality improvement.
transport aircraft to improve taxi performance on rough run-
ways. They worked with a secondary air chamber that could be 1.3.2 Control LAWS
pressurized shortly after touchdown in order to generate the Different control laws are investigated in the literature, mostly
desired softer pneumatic spring rate as a function of aircraft for cars and trucks. During aircraft operations, however, associ-
weight. Today, with reliable, inexpensive and powerful elec- ated parameters change during a wide range. Take-off weight
tronic signal processing available, research has moved toward will always differ drastically from landing weight, aircraft
. computerized closed-loop control. In 1984 a feasibility study ground velocities vary between very low speed at taxi and 60 to
for a series-hydraulic active control landing gear intended for 80 m/s at take-off and landing. Oleo and tire properties may
supersonic military aircraft was published by McGehee and vary depending on age and maintenance. Finally, runway qual-
- - Moms, [ 5 ] . Additionally to theoretical analysis a number of ity differs from airport to airport.
tests had been performed on a test-rig to permit experimental Suitable control laws have to be robust against these changes.
verification of the concept. The setup is typical for a fully Several solutions to the problem exist and have been discussed
active control system and has before and since been investi-
I gated several times. The gear force applied to the airframe is
regulated by the hydraulic pressure in the piston of the oleo
thoroughly, ranging from sky-hook damping to nonlinear and
adaptive control laws (for a state-of-the-art review see [ 1I]).
For many problems the algorithms have to go further than to
which is used as an actuator. solely observe center of gravity motion. Other points influenc-
The results obtained led to the modification of a F106 nose ing the results are maximum applicable actuator force and sus-
landing gear which was also tested in test-rig setups, [6]. pension travel limitations.
Another investigation that combined an analytical and an For the problem discussed in this paper, the sky-hook damping
experimental approach was conducted by Freymann, [7]. approach has been chosen.
Object of research was the simulation of a nose gear oleo.
Again the force generation was realized by pressurized oil, 1.3.3 Hardware
even though the valve layout was different. A laboratory setup
Since the idea of active and semi-active landing gear is rela-
was used to experimentally verify the basic control concept.
tively new, and no system for aircraft is as yet in the production
Even though the results obtained in both cases might be overly
stage, the hardware is still in the phase of concept evaluation.
optimistic when it comes to application to a realistic aircraft,
In the field of fully active oleo control so far heavy laboratory
the investigations showed nevertheless that a fully active land-
setups have been used to validate simulation results.
ing gear is feasible and can lead to significantly reduced air-
For semi-active control systems, some manufacturers, espe-
frame loads. Another somewhat more recent approach is the
cially in the domain of heavy-duty trucks, have proposed and
analysis of improvements gained with closed loop semi-active
built technically feasible layouts. Other suggestions include dry
oleo control. Studies by Kamopp [8] for automotive applica-
friction damping (which, for large aircraft, is likely to be unre-
tions show that the performance of a semi-active damper is
alistic due to the large energy dissipation). A relatively new
only marginally smaller than that of a fully active system, pro-
proposal is the use of electro-rheological fluids, i.e. fluids that
vided that an adequate control law is used. Catt/Cowling/Shep-
change their viscosity in response to an electrical field.
herd [9] come to a similar result in a simulation study of
Objectives for the final choice of a system will surely be led by
aircraft suspensions. The semi-active damping device is usu-
a consideration of performance against penalty. Surely low
ally proposed to be an oleo with a modulated damping orifice
weight is desirable, and high safety against failure is essential.
cross-section, completely omitting the metering pin.
A good performance also requires a high bandwidth, low reac-
5-4
tion time actuator. The possible achievable actuator peak force Hertz band above and below the critical range. A list of the
might also be a decisive criterion. rigid and flexible eigenmodes can be found in [lo].
An active controlled landing gear introduced to improve the
ride quality of a large aircraft would therefore have to concen-
2. Motivation trate on reduction of vertical accelerations due to elastic body
oscillations in the 4 to 8 Hertz frequency range. The excitation
2.1 The Problem of these oscillations is caused by runway roughness either due
to wear and tear or to other unavoidable reasons. Qpical exci-
2. I . I Ride Comfort Assessment tations used in the simulations are presented in chapter 3.1.4.
In chapter 1 it could be seen that the improvement of crew and
passenger ride comfort will be of growing concern for manu- 2.2 A Proposed Solution
facturers of large aircraft.
As base measurement for ride comfort assessment the vertical Recent research indicates (see chapter 1) that active gear tech-
accelerations are widely used. Since the individual perception nology may pose a possible solution to the problems of aircraft
is additionally influenced by other factors of the biomechanical ground loads and that the semi-active controlled front landing
human system, scales have been developed to weight the fact gear promises to be the most likely solution to be incorporated
that certain frequencies are perceived to be more uncomfort- in a production aircraft by the aircraft industry.
able than others for a given amplitude. In IS0 2631, [12] the The proposal of this paper to solve the ride comfort problem is
frequencies between 4 and 8 Hertz are denoted as the most cru- therefore: design the controller for a front landing gear, modify
cial for comfort (figure 3). Ride comfort improvement has a conventional landing gear to incorporate controlled damping
(thus creating a semi-active front landing gear), analyze it and
compare it with the conventional gear.
Fed with this data, the BEAM prepmcessor calculates the * load variations small enough to allow linearization of defleo
desired eigenmodes of the four airframe members. These tion curve.
eigenmodes are then used by SIMPACK to evaluate the aircraft This tire model was used for most of the taxi simulations.
eigenmodes implemented in the simulation. The selection of
the number of eigenmodes plays an important role. High struc- Oko-Pneumatic Shock Absorber Force Imv:
tural vibrations result in small integration stepsize and thus Conventional (Passive) System:
increase CPU work. On the other hand the SHANNON theo- As primary force generating element each aircraft landing gear
rem determines the maximum eigenfrequency implemented in is equipped with one central oleo, acting on the principle of a
a cenain model to be at least twice the maximum frequency one stage or single slope oleo like those used in the landing
under consideration ( 1151 and 1161). The eigenmodes for the gears of the A300 aircraft. The passive oleo displays two sepa-
taxi analysis were chosen in a way to guarantee a safe analysis rate force generating elements, the gas chamber that generates
up to 25Hz. the spring force and the set of orifices responsible for the
damping force. The total oleo force is composed of a gas spring
part (calculated using the general gas law). a damping part and
a friction part.
The damping force was calculated using the squared value of
the relative velocity of the sliding oleo member and the main
fitting attached to the fuselage. It is dependent on the direction
of the relative motion because of different dampingcoefficients
for compression and for extension. The values for theses coeffi-
cients, depending on the oleo orifice ROSS sections, were deter-
mined in tests by the landing gear manufacturer.
Friction forces are difficult to assess and depend on a multitude
of parameters. Especially for the landing impact case, longitu-
dinal forces increase. the friction and can in fact block the gear
and prevent spring deflection. leading to high structural loads
and even structural failure. The gas chamber pressure influ-
ences the fitting pressure on the outer oleo wall, resulting also
in higher friction forces. For the taxi condition the gas pressure
variation and applied x-forces remain small, so a constant fric-
tion term supplied by the landing gear manufacturer was incor-
poIated.
Semi-Active System:
The main difference between a serm-active and a conventional
Figure 6:First and Second Fuselage Bending Eigenmode oleo is found in the oil orifices responsible for the damping
force. To provide control over the damping force, the oil Bow
and thereby the damping coefficient is regulated. This can be
In figure 6, the first two fuselage bending eigenmodes are. achieved by variable diameter orifices or by implementation of
shown. a servo valve. Since the actual design may exhibit quite differ-
ent oil flow regulators, but all of them control the damping
3.1.3 Force Laws coefficient, this coefficient was chosen as control variable. The
To simulate the dynamicalbehavior of the A300 airplane, SIM- control range, i.e. the maximum and nunimum damping coeffi-
PACK had to be given a set of force laws. With the aid of the cient values achievable. was deliberately chosen to lie between
relative displacements and velocities at the specified node 20 and 500% of the original value used in the conventional
pints of the multi body system, the according forces were A300 oleo because this range was considered to be technically
derived. The force laws are rcalizcd as user routines within feasible.
SIMPACK and are explained in detail in [IO].
Aemdynantics Force LAW:
lire force law: Stationary motion of an aircraft through the air results in an
The tire force model has to yield the longitudinal and lateral aerodynamic force acting on the aircraft. Taken into account
forces and the torque moment acting behveen the road surface the small influence of the aetodynamicfaces and moments on
and the wheel and depending on the slip between tire and the overall dynamicsof the system, the d y n a m i c force user
ground Nevertbeless. even with a relatively simple model, routines were limited to the longitudinalmotion. To account for
computational effort has proven to be substantial for an aircraft non-stationary moments, the pitch damping coefficient was cal-
with ten or more wheels. Lkpending on the simulation case it culated.
proved to be helpful to linearize the tire behavior if possible.
For an aircraft taxiing on ground without corneringor braking, 3.1.4 Runway sufme
a linear spring model was incorporated. The following supposi- To analyze the dynamic behavior of a general nonlinear MBS
tions apply: model, time histories have to be evaluated for a broad range of
* negligible slip (no braking). frequencies. The excitations used on the model should include
* no side slip, a wide frequencies range. Additionally, single frequencies of
* no loss of ground contact for any tire, special concern might be agitated separately.
5-7
The following existing runway Nrbances have been taken A comparison between four measured runway ZDs (thin
into account in the present investigation: lioes) and the PSD of the implemented model runway (solid
General roughness: Especially on older ailports and on line) is shown in figure 9.
airports in the countries of the former Warsaw Pact, the
run- and taxiways display undulations of all frequencies
and with higher amplitudes than those measured on air-
ports in western countries, 1171. These conditions have Power spectru~Density [ft/radian/ftl
reportedly led to dangerous situations. Pilots of aircraft
equipped with Cathode Ray Displays (so called glass
cockpits) experienced difficulties to read the instru-
ments due to vertical accelerations during the take-off
Nn.
Concrete Plate Deforming: A widely used method to
construct fortified runways is the casting of large plates
using Liquid concrete. These plates are separated from
each other by gaps filled with rubber. Aging of concrete
runways causes the plates to settle unevenly, leading to
long wavelength bumps and steps at the gaps.
Center Line Lights: All run- and taxiways are equipped
with lights indicating the middle line during night time
opcsation. These lamps extend a few centimeters from
the ground (figure 7) and exert a shock when hit by a
1 - ModelRunway
-Measured Runway
tire. They are spaced regularly and can induce oscilla-
tions. 10-2 I U iff IO
Reduced Frequency [radianlftl
a: Housing Height (2...5 cm) Figure 9 Power Spectral Density of Measured and Model
Runway
b Length (I5...30cm)
c: Spacing 6.. 15 m)
4
Multiple Pulse Runway Model:
Step and pulse model inputs are especially suitable to induce a
broad band of frequencies into a model. In practical ground
. _ operation. pulse inputs can result from concrete plate gaps or
from center line light housings. In figure 7, a light housing built
Figure 7 Center Line Light Housings in a Runway into the ninway is sketched. In reality. the tire hitting the hous-
ing will deform, so that the impact of the heave motion of the
bottom of the tire spring will be muffled. The multiple pulse
To simulate these runway surface conditions, the following runway model uses this kind of pulse shape. Housing dimen-
three distinct models have been programmed. sions and spacing are programmed as variables and are defined
in the experiment model setup.
Qwi-stochastic Runwny Model
Thompson [17] gave a list of a number of measured Power Harmonic Wave Runway Model:
Spectral Densities (PSD) of runways in NATO countries. These The harmonic wave as excitation was included in the system
were used to build up a deflection curve using a set of harmonic analysis to perform a sort of worst case evaluation. With this
functions. The runway elevation curve resulting frnm this p m runway model, a distinct frequency acts on the landing gear to
cess is displayed in figure 8. induce potentially problematic oscillations like the main fuse-
lage bending eigenfrequencies. A simple sinusoidal wave was
used for the model. In figure IO, an example runway surface is
I
40. 80. im. 160. 200
Runway x-mrdinate [ml I
40. 80. im. 160. zw
Figure 8: Quasi-stochastic Model Runway Surface Runway x-coordinate [ml
given used for the excitation of the 3.6 Hertz fuselage bending using the Sensor Signal Processing and the control algorithm
eigenmode of the A3M) at a speed of 20 ds. implemented in the Contmller.
The velocity feedback used as controller input makes the cho-
sen control concept a variation of the sky-hook approach.
4. Control Concept Analysis
4.2 Optimization Process
4.1 Control System Lay-Out
4.2.I Criteria Selection
During the selection process of a control law to be imple- The goal of the implementation of a controlled front gear is to
mented, the rigorous safety requirements in aircraft design increase passenger and pilot comfort by reduction of fuselage
have to be taken into consideration. Thus, low system com- vibration. The criterion vector c must take into account the nec-
plexity is one of the design drivers for the control system. A essary states. The aircraft state vector Y
single input single output (SISO) control system requires the
minimum of one sensor signal and one actuator element and
was therefore chosen.
Aircraft currently in commercial operation are in part already
equipped with acceleration sensors. To keep additional system r= (4.1)
requirements low, a single acceleration sensor at the cockpit
position was assumed and used for the feedback signal. Figure
I 1 shows the resulting system layout.
contains the states describing the aircraft motion. Here,
It/ / I
stands for the velocity states of the joints.
The block diagram in figure 12 shows the general principle of expresses the deformation coordinates in the velocity domain.
the feedback system. The Oleo Closure Rate remains the main For calculation and analysis the states of several distinct loca-
tions on the elastic bodies are needed. Unlike on rigid bodies.
where states of any location can easily be derived if the states
of the reference system of the body are known, the states of
arbitrary locations on the elastic body depend on the deforma-
Olrn tion position and velocity cwrdinates.
Closun Aircraft To observe all interesting eigenmode shapes it is necessary to
introduce measurementsof suitable locations on all elastic bod-
Aircraft swcrurc
ies (fuselage and wings). For example, the first fuselage bend-
ing mode can be observed by measurements in the far front or
far rear end of the fuselage, whereas for the second eigenmode
additional information is needed from stations between the
Sensor Signal front and the center of gravity of the aircraft. As a compromise
Conmlkr Recessing between the resolution of high frequency eigenmodes and sim-
ulation time. 16 locations have been selected on the A300 air-
frame.
Figure 12: Control System Block Diagram The semi-active landing gear influences primarily the vertical
motion of the aircraft. Vertical acceleration was therefore
input for the Landing Gear Oleo. It determines the forces picked as criterion state to be measured. The Root Mean
applied on the Aircrafr Sfrucrure depending on the Runway Square (RMS) value is calculated using
Pmfile inputs and the Damping Orifice Area Actuator. The Air-
craft Response coming from the Aircraft Strucfure is fed back
5-9
with
CLRMS =
Jj(iz-Xm)
dt (4.6)
Yith
(4.13)
Besides the Root Mean Square value (RMS value), which rates
amplitudes and their duration time, peak values are includes as
well:
cl,pe*k= m a x ( a b s ( x , - x m ) ) (4.7)
where cLwk denotes the peak-eriterion value. SIM'ACK Simulation and Analysis I
With this double criterion strategy, a trade-in of low R M S val-
ues for very high but short peaks is avoided. RMS values and
pealr values of all sensors were weighted equally, with slight cdtcriaq
3 * 16 RMS values of vertical accelerations
emphasis on the fmnt fuselage part 3 * 16 maximum amplituaes of vertical accelerations
Cost function and constraints have to be continuous. hut no The optimization process described above was used to derive a
specific structure is requid. For the acceleration minimiza- set of controller parametersbest suited to minimize the design
tion, the cost function can now be stated as parameters 4. On the following pages the system performance
using those controller parameters is shown for the three exdta-
tion cases. The performance is subsequently compared with
that of the original conventionaltype landing gear. A complete
presentation and discussion of the results can be found in [lo].
5-10
5.2 Quasistochastic Excitation Performance gear can especially be seen for the first fuselage bending eigen-
mode at 3.6 Hertz and for the second symmetric and antimetric
In figure 14 the vertical cockpit acceleration response of both wing bending eigenmode at 5.6 and 5.7 Hertz, respectively.
Vertical Cockpit Acceleration [mls21 Figure 16 shows the system behavior for the repeated pulse
excitation as exerted by center line light cases or concrete plate
gaps.
lo 7
-1 . -Conventional Nose Gear
- Semi-active Nose Gear
0. 1. 2. 3. 4. 5.
lime [SI
-Conventional Nose Gear
Figure 1 4 lime Domain Response for Quasistochastic
Excitation
10-31
- Semi-active Nose Gear
the conventional gear (thin line) and the semi-active gear @old 0 IO 20 30
line) in the time domain are shown for a total simulation time Frequency [Hertz]
of five seconds. The system state at the beginning of the simu- Figure 1 6 Power Speceill Density of Center Line Light
lation was determined by a pre-analysis time integration to Excitation Response
guarantee a steady state response throughout the whole analy-
sis period. It can be seen that no fading oscillations are present.
The controlled gear performs well, the amplitudes of the accel- While the system response of the semi-active gear remains
eration are diminished by a factor of five. below the one of the passive gear over the whole frequency
The power spectral density (PSD) of the cockpit accelerations range, two distinct peaks can be seen below lOHz for the air-
of this time history. shown in figure15 for the frequency range craft structural eigenmodes. At those points the system perfor-
mance of the semi-active oleo retains only a small advantage
over the conventional system. This behavior differs from the
one seen for quasistochastic excitation and may be amplified
Power Specaal Density of Cockpit Acceleration [m2/s31 by the total absence of inputs between the bumps and on the
main gears.
Power Spectral Density of Cockpit Acceleration Im2/ s31 Power Spectral Density of Cockpit Acceleration [m2/ s3]
i
0 10 20 30
10-3
0 10
-Single Gear Excitation
-All
20
Gear Fkcitation
30
Frequency [Hertz] Frequency [Hertz]
Figure 1 7 Power Spectral Density of Harmonic Excitation Figure 1s:Comparison of Power Spectral Densities of Nose-
Response (3.6Hz) Gear-Only and All-Gear Excitation Responses
. -
[6] W.E. Howell. J.R. Mffiehee, R.H. Daugherty, W.A. puter codes in vehicle system ~ p a m i c sVehicle
.
Vogler: F-106 airplane active control landing gear System Dynamics, Supplementto Vol. 22. Swets &
drop test performance. Landing Gear Design Loads Zeitlinger, 1993
(CP W ) ,AGARD.1990
171R. Freymann: An active control landing gear for the
alleviation of aircrafttaxi ground loads. Zeitschrift fiir
Flugwissenschaft und Weltraumforschung, 1987
[8] D. Kamopp: Active Damping in Road Vehicle Suspen-
sion Systems. Vehicle System Dynamics, 12(6),Swets
& Zeithgex, Lisse, 1983
[9] m n e Catt, David Cowling and Alan Shepherd:
Active landing gear control for improved ride quality
during ground roll. Smart Stluctures for Aircraft and
Spacecraft (AGARD CP 531), Stirling Dynamics Ltd.
Bristol, 1993
[IO]H. Wentscher: Design and Analysis of Semi- Active
Landing Gears for Transport Aircraft. h.D. Thesis,
TU MUnchen, 1995
[ll]R.V. Dukkipati, S.S. Vallurupalli,M.O.M. Osman:
-
Adaptive Control of Active Suspension A State of
the Put Renew. The Archives of Transport. Vol. IV.
1992
[I21 International Organization for Standardization: Guide
for the Evaluation of Human Exposure to Whole Body
Vibration. IS0 2631-1974(E). July 1974
6- 1
Abstract
Dynamic behaviour of motorcycles front wheel are not different to similar wheel suspension
motion of airplanes and road vehicles. In every case self excited oscillations caused by
the tire occur. Only the different design of wheel suspension produces some different
coupling. The deviation of a moving wheel on a horizontal plane from straight rolling is
caused by unstable lateral contact forces. The kinematical generation of such forces can
be imaginated looking on the elastic tire. A first order time-depending behaviour together
with inertia mass and damping of the wheel using the HURWITZ stability criterion gives
important insight into the main parameters of the problem.
Nevertheless the nonlinearities of the system overshadow this clear insight. Cornering force
and cornering moment are nonlinear functions of cornering angle, so the use of the well
known Magic Formula together with the black box method of control theory is a logical
extension. The problem of instability arises again when thh wheel is excited from wavy
surface or similar from mass unbalances of the tire. A direct time depending of the lateral
tire forces can be used to show parametric instability of the system (MATHIEU equation).
The rolling of the wheel on a non-smooth surface needs then to introduce methods of
discrete system approach. The paper presents an overview of the problem of shimmy.
1 Introduction
The pneumatic tire as a rolling suspension is intended, on the one hand, t o ensure that track
unevennesses are traversed with as little vibration as possible, and, on the other, t o guarantee
the driving direction a n d stability of the vehicle. In sharp bends, the pneumatic tire should
also withstand the centrifugal forces, and it should generate the appropriate longitudinal
forces during braking and acceleration. These manifold requirements, however, can only be
met by means of a constructive compromise, wherein a n unavoidable coupling arises, by
virtue of the elasticity of the tire in all three dimensions, between the longitudinal, lateral
a n d vertical forces a n d between their associated dynamic deformations. These couplings lead
t o non-linearities, which complicate mathematical treatment.
A first attempt t o solve the problem starts out from the premise that the track is even a n d
the tire ideally round, a n d also assumes, with regard t o the contact time of a tire particle with
the ground, that the associated slippage between tire and ground occurs in linear fashion,
from the entry point of the contact surface t o the exit point thereof. If one additionally
takes into consideration an average deformation g of the contact surface, due t o the effective
lateral force S, then a linear elastic deformation equation can in principle be written down
for the lateral deformation as a function of lateral force:
s = cs g . (1)
Paper presented at the 81st Meeting of the AGARD SMP Panel on The Design, Qualification and Maintenance
of Vibration-Free Landing Gear, held in Ban8 Canada from 4-5 October 1995, and published in R-800.
6-2
The associated lateral spring coefficient cs can be determined experimentally. The linearly
increasing deformation from the entry point to the exit point, by contrast, is assigned to
the travelling surface (shear deformation) and is superposed onto the deformation g. This
enables one to write down a kinematic equation for the slip angle of the tire to which the
lateral force S is assigned, of the form
This effective slip angle 6 , j j , which is a function of the elastic lateral deformation of the tire
carcase, is equal to the stationary slip angle for the case S = constant. The function S is non-
linear, since it describes the partial sliding of the tire, and even progressively decreases overall
once it has passed beyond the sticking state. If one only takes into consideration small slip
angles in the region of zero (driving straight ahead, equation (2)), then the increase Po c' in
the cornering force function at the zero point yields the following linear differential equation
of first order, whose eigenvalue is negative and therefore exhibits damping behaviour. Any
disturbance of the stationary rolling state is therefore damped along the direction of rolling,
wherein the damping distance is approximately equal to the tire radius. Since, however, the
stationary slip angle characteristic curve is largely a function of the length of the contact
area, which in turn is a function of the possibly dynamic wheel load, there is automatically
a coupling between the horizontal and vertical dynamics of the tire.
axle is always in front of the contact patch, see figure 2. Let the wheel be elastic mounted to
the straight-ahead position by means of an elastic attachment (torsion spring coefficient C L ) .
There is a principal moment of inertia I about the vertical axis. For the pneumatic tracking
we designate n p = !$. In the case of damping proportional to velocity (steering damping),
a coefficient k~ is used. The characteristic curves for the tire are displayed in figure 3.
Analytically, the use of n p has the advantage of a simpler mathematical description. The
equations of motion are as follows:
(5)
6-3
Figure 1: Ezomples of wheel suspension: a) motorcycle, 4) motor car and c) aeroplane (after [Volk93])
S n
200
Figure 2: kinetic and dynamic quantities of the Figure 3 Lateral guiding force and pneumafic
steered wheel irackinglength
= 0. (8)
The expansion of this gives:
In the case of a wheel free to rotate around the steering axle CL = 0, k~ = 0, the ROUTH-
HURWITZstability conditions for the bow wheel of an aeroplane or for the steering of a
motorcycle give
a cs -Po c' > 0. (11)
holds:
pi ""'b
/.
d 5 6 i- K
a 2 1: c; = CL (1 - k(a))
The shape of the function k ( a ) is given in figure 5. From the polynomial of third degree for
the characteristic value X
6-5
from which it can be gathered that there is an additional term proportional to the factor c;.
From these formulae it can be seen that
1. without steering damping there can be no stabilisation and
2. if there is steering damping, sufficient stabilisation can also be provided with a dimi-
nished steering stiffness c:.
3. Since in general the geometric tracking length is significantly smaller than the relaxa-
tion distance, the additional terms in equation (9) must compensate for the negative
effects of the stability condition (equation (15)).
Thus steering tolerance is frequently not immediately recognised, since the graph of the
function k(cr) only slightly diminishes the effect of the steering damping. Thus a maintenance
problem arises from the fact that, even with a very small amount of steering tolerance,
vibrations are amplified.
Figure 6: Dependency of
steering restisiance on fric-
tion
4 MR
';= .IrA(-J
6-6
where MR represents the frictional moment in the steering. For the case where CL = 0, but
ki # 0, the characteristic polynomial is now:
x3 + (-Pov + f )x2 +
cs
c
[(U + no) -
U cs v cs k;
Ia
+ PO Ia
+ + no) I
(U ,- 0.
v cs (18)
The HURWITZdeterminant is given by:
Here one can also see that irrespective of whether A,,, i.e. the HURWITZdeterminant for
the case CL = k~ = 0, is positive or not, a suitably dimensioned frictional moment or ki
compels st ability.
This simple analytic investigation makes it clear that every rotating but elastically mounted
wheel is capable of exhibiting self-excited rotational vibrations about the axis of steering.
Even if wheel load variations are not taken into consideration, and thus an even ground is
assumed, there is a tendency towards so-called wobble vibrations in case of a motorcycle. The
effect of the non-linearity of the system on the amplitudes of vibration alone is also exhibited,
wherein more exhaustive efforts with the aid of harmonic balances also take into consideration
the non-linearity of the cornering curve. Unfortunately this effect is in principle non-linear, I
also with respect to the vertical dynamics of the wheel, leading to states of motion at several
frequencies. The method of harmonic balances is thereby rendered useless for dealing with
the problem. Notably the tire width and curvature of the contact area have a stabilising 4
effect [Boe77].
3 Steer-Shimmy of Motorcycles - i
Instabilities of the steering system are known for several types of vehicles, e. g. the tail- or I
mainwheel of airplanes and the (steer-)shimmy of automobiles. This instability also accurs
at single-track vehicles as bicycles, scooter, and motorcycle. It describes the resonance
movements of the steering-system around the steering axis.
Several authors examined this phenomenon of motorcycles (e. g. [Pac66], [Sha71], [Sha77],
and [Koch80]).
The modelling of the motorcycle by means of a linear differential equation system of 20.
order provides the main (eigen) figures of the driver - single-track vehicle system. [Sha71]
indicated the terms for three eigen-values:
0 capsize mode, as vibration-free movement of tilting (tipping) of the whole motorcycle
to one side, which is stable in the lower velocity range, and constantly instable in the
medium and higher velocity range.
0 wobble mode, which is called shimmy at other vehicles, and is steer motion which
occurs at motorcycles as a frequency between 7 and 10 Hz almost independent from
the vehicle velocity, and is good damped in the lower and medium velocity range.
0 weave mode, as vibration of the whole motorcycle around the yaw axis with an incre-
asing frequency at lower driving speed and in the medium and higher velocity range
movement tending towards a limit value between 2 and 3 Hz.
6-7
The linear theory, however, does in comparison with measurements on the motorcycle not
come up with real-world condition results for the wobble mode, since the whole calculation
model for the medium velocity range calculates a well damped vibration then, the mea-
surements show that especially here just very low, respectively no damping exists. Other
mechanical couplings are, therefore, effective as they are enclosed in the linear vehicle model.
From the observation of a real vehicle showing that shimmy movements are very much cou-
pled with the balancing condition of the wheel, [KochSO] developed a model, which contains
the steering system of the motorcycle alone, without the rear frame. As a fundamental addi-
tion compared to previous modelling, he encludes the periodic wheel load oscillation, which
comprise a force component around the steering axis, supplying the energy for exciting the
wobble mode. These wheel load oscillations occur through pitching and vertical motions
of the vehicle and front wheel which are connected with the radial deformation of the tire.
These oscillations are produced by road irregularities, as well as by the periodic forces of a
poor balanced wheel.
[Koch80], herewith, got a mechanical backup system, which encludes the radial deformation
of the tire in the wobble mode. Result of this modelling is a non-linear differential equation
system, which can be transformed to a MATHIEU differential equation, so that the kinetic
ranges of instability of the wobble mode can be demonstrated.
(mo:Mass of wheel, CR: Tire radial stiffness, y: vertical motion of wheel axis, z:
vertical motion of road surface). With y = yo cos(fl t) one gets:
6-8
The reaction force coming from the road excitation results as radial force to
FR= CR (Y - 2)
at the center of the tire contact.
2. Solution of the pitching motion
A road irregularity at the front wheel produces a pitching motion of the vehicle. Fig.
8 may be used as model. The following differential equation for the center of gravity
~~
4
(q: pitch angle, cRV/R: front/rear wheel stiffness, z: vertical motion of road surface,
u ~Eigenfrequency
: of pitch motion) so finally the solution is
q = qo cos(n t ) = Qv aG 20 ms(R t ) .
0, (U; - 521)
The following superposed movement results from the pitching and vertical motion of
the whole vehicle for the front wheel
(33)
(34)
~ ~~
.- ..
*..+ n Fsv + Mzv + FR n
C O S ~ sinX = 0, (35)
8 = Ivz + mv 1, + IGz + mo 2.
2
(36)
Linearization for small steering angles and with
results in
l+2w~D;\+w: (39)
If the radial force from the pitching and vertical motion are introduced into the quations of
motion around the steering axis the following nonlinear equation results out of it:
Although the solution of the MATHIEU differential equation demonstrates especially at dou-
ble excitation frequency its main range of instability, the chape of the non-linear function,
however, dictates what happens.
The shape of this function is shown in fig. 10 using the data of a motorcycle as function of
the driving speed U, respectively wheel frequency R. The shape of the front axle amplitude
Irl 4
as function of the driving speed, respectively wheel frequency is shown in fig. 11 showing two
resonances of the vertical and pitching eigenfrequency. Since the resonance of the pitching-
I
Figure 11: Amplitude rafe of fhe vertical fmnt-arle track (afler [Koch80])
eigenfrequency (wn = 7,78 Hz) approximately coincides with the medium range of instability,
also instable wobble motion is to be expected within this velocity range.
The measured frequency and velocity during the wobbling of this vehicle agrees with theo-
retically gained values of the nonlinear stability analysis. These results show beautifully
that the wobble mode can also be determined as a characteristic value with the help of the
linear theory, though it does not agree with measurements at the motorcycle, with respect
to either frequency or damping. Only by including in the modelling the vertical dynamics
6-1 I
of the wheel together with the lateral dynamic characteristics, which leads to a non-linear
differential equation, results in regions of instability that can be confirmed by measurements.
A first- or higher-order theory of tires for the lateral characteristics and a radial tire stiffness
dependent on the tire speed may further refine the results.
Even if [Koch801only used the balanced masses in its model as an excitation, the broad spec-
trum of road unevennesses also provides excitations. This is also confirmed by measurements
wherein greater instability arises with an unbalanced front wheel.
Also G. Nybakken is engaged in work in his dissertation [Nyb73]on investigating the stability
of rheolinear aeroplane landing gear.
6-12
brake force and self aligning torque with great accurac. That mathematical representation
is limited to steady-state conditions.
Figure 12 shows the basic form of each characteristic G :he tire. Bakker et al. took over the
I I I
Figure 12: Steady-state tire c h a r a c i e n x s (after [Bak87])
4
idea by Ruf [Ruf83] to approximate these tire behaviour by special angular functions.
The demands for that formula was expressed - -
0 to be able to describe all steady-state tire characzistics,
0 being easily obtainable from measured dates, - \
,
0 being physically meaningful; its parameters shoii characterize in some way the ty- ,I
pifying quantities of the tire,
0 being accurate.
These demands were found for all tire characteristics i: :he expression
with Y(X)standing for either side force, self aligning tcque or brake force and X denoting
slip angle ( C Y ) or longitudinal slip (6). I
D in equation (42) is the peak value and the produc- D C B equals the slip stiffness at
zero slip. The coefficient C governs the shape of the crve and the coefficient E makes it I
6-13
I
I possible to accomplish a local extra stretch or compression while stiffness and peak value
remain uneffected.
Due to ply steer, conicity, rolling resistance and camber, the characteristics will be shifted in
horizontal and/or vertical directions. These shifts are represented by Sh and S, respectively.
By these approximating functions Bakker et al. got in parts excellent fitting results of
measured and calculated tire behaviour as shown in Fig. 13 and 14.
B R A K E FORCE -Fa
SIDEFORCE - F t
(NI
8,000
6,000
c.000 UCASUREMEHTS:
= (7, = ODCO.
b 0, - 2 DEO.
2.000
0
0 10 20 30 10 60 80 10 80 90 100
-LONGITUDINAL
SLIP - K 111
SELF ALlONlNG
TORQUE -Yz
IWm)
2oofMEASUREMENTS
.---
-2001
For simulation respectively pre-calculation of the dynamic driving behaviour of single and
double track vehicles, the knowledge of the dynamic behaviour of the tire is necessary,
without time consuming sub-routines for a tire model. [Schu87] reports in detail on the
latter.
Easier and without time consuming appliance the transfer properties on a test stand with
time varying exitation signals can be evaluated. This leads to no explanation, and realization
on the physical processes inside the tire respectively in the contact area though, but enables
the vehicle engineer to easily integrate the tire dynamic into a model of vehicle dynamic or
driving dynamic. When measuring the tire with time varying exitation, the simultaneously
occurring inertia forces and gyroscopic torques, which do not occur at static measurements,
have to be corrected carefully. On this procedure, and the following description of the
transfer function by means of control equations reported [Wi1187], and [NastSl] applied to 4
different tires, two motorcycle, and two passenger car tires. The time varying side slip angle
a(t)was described with a frequency limited white noise up to approximately 17 Hz with a
variance of 0,21" (see fig. 15).
AMPLITUDE (degree'degree)
1.E41
1
1 10
I FREQUENCY (Hz)
~~~~
I
Figure 15: Power spectral density of side slip angle (after [NastSl])
If one applies the corrected dynamic self aligning torque, and the corrected dynamic side
force to the stationary values, that are these values, which occur at side slip angles (f = 0 Hz)
constant with time, so one will come to non-dimensional self aligning torques and side forces.
By computing the power spectral densities, one finally gets the transfer function Malaof the
self aligning torque, and Fa/aof the side force. The fig. 16, and 17 show this exemplarily
as Bodeplot.
6-15
LMPUTUDE(lldegree)
ia
U. I
1 10 20
FREQUENCY (Hz)
45 -
v3
-45- v*
I - -
-90- vi = 5 m/s
. v2 = 20 m/s
-135 -
- -
. v3 =50m/s
-1 80
1 10 20
FREQUENCY(Hz)
Figure 16: Measured transfer function of self aligning torque relative t o side slip angle
(after [NastSl])
6-16
LMPLITUDE(l/degree)
10
I v,J
1
v*
0.1 1 FREQUENCY(Hz) 10 20
45
0
"3
-45 v*
vl
-90
1 10 20
FREQUENCY (Hz)
Figure 17: Measured transfer function of side slip force relative t o side slip angle (after
[Nast 911)
A PD2-behaviour (proportional plus double derivative action controller) is assumed as ap-
proximation equation for the transfer function, which reproduces both amplitude channels
appropriately. Merely the phase has to be adapted by another approximation term - a
dead time term. That means that the dynamic behaviour of self aligning torque, and side
force has to be adapted by higher terms as 2nd order.
The approximation equation in the time domain for the self aligning torque are
QO
4#
Q4
QZ
U 2 2s 3
dQ0
Y-
I.J 0 S IO I5 I 30 U 40 4J W
Figure 18: Stiflness factors CR and Cs as function of wheel load; damping factors DR and Ds as func-
tion of wheel speed: 1. Bridgestone 120/80 V 16 (motorcycle); 2. Bridgestone 150/80 V 16 (motorcycle);
3. Conti Contact 155 R 13 (automobile); 4. Conti Contact 175/70 R 13 (automobile) (after [NastSl])
6-18
Massfactor MR [Nms2/degr.]
0,015
0,Ol
0,005
0
0 5 10 15 20 25 30 35 4 45 50
Velocity [m/s]
Velocity [mls]
0.01
0
0 5 10 1 5 2 0 2 5 3 4 3 5 4 0 4 5 5 0
verocrrvtw
0 the stiffness factor increases on increasing wheel load for the side force as well as for
the self aligning torque.
0 CR is constant for motorcycle tires (cross bias type) over the speed, for passenger car
tires (radial tires) the stiffness factor and CR increases degressively on increasing speed.
On further measurements at time varying wheel load F. Sun [Sun901 also reports in his
Ph.D. thesis, and reached similar results, in which he additionally determined the dynamic
transfer characteristics for time varying wheel load.
Since the whole transfer behaviour is still very much depending on the wheel velocity, except
the stiffness factor CR for the self aligning torque of the motorcycle tire (cross bias type),
one can assume that non linear effects, i.e. because of gyroscopic effects can not fully be
recorded by means of the linear approximation formulation.
Assuming black-box-formulas are available as break-, sideforces and torques are described in
their dynamic behaviour a straight forward integration procedure is possible which will be
shown here.
6-20
A stabilisation of the numerical or piecewise linear integration in the case > 0 can be
achieved through an extension of the equation by a damping term
and possibly a mass compensation term. These, however, does not render impossible the
correct assessment of short-wave surface unevennesses, with regard to their effect on P ( t )
and on S(t). As a result, the high-frequency rolling process cannot be treated further by
means of analytical mechanics.
Since further investigations of an analytic kind are not possible, one is compelled to resort to
numerical methods. In these, the contact surface is modelled as discrete, and the sequence
of contact points is investigated in discrete calculating steps. A unilateral contact condition
enables the numerical representation of partial lifting and traversing of track unevennesses.
In the process, the tire structure is resolved into many degrees of freedom, taking into
consideration the tire mass and damping. Such a structure is geometrically non-linear in
principle and has a relatively high frequency. It is guided by means of a rigid body system
of low frequency, consisting of a rim and shock absorber. If, as is customary for motorcycles,
one introduces a head lug angle A, then one obtains the following system of equations of
motion for the steering angle p and the spring deformation s of the wheel axle (see figure 1).
In addition, there are also the lateral motion and longitudinal motion of the motor vehicle
mass itself. The contact forces generated while traversing the ground unevennesses follow
from the equations of motion for the tire structure, which are not given here [Boe93]. The
sample calculation given in figure 21 shows the dynamic effect of wheel load variations on
a stable steering system that is nevertheless capable of damped wobble. The non-linearities
taken into account in the numerical integration are the deformations of finite magnitude, the
unilateral contact, i.e. local lifting and settling of the running surface on rough obstacles,
and the effects of traction and sliding friction. The model does not take into account the
effects of tire width or of cross-sectional bending. Wobble damping, however, is markedly
influenced by enlarging the tire width and increasing the length of the contact surfaces. The
resultant moment of pivoting has the effect of diminishing wobble. Thus, not surprisingly,
for small contact surfaces, which arise,,under low wheel loads, have an especially critical effect
on wobble. A surface-oriented theory of non-stationary rolling contact has been analytically
developed for perfect roadholding, taking into account the bending of the trajectories of the
6-21
II
1
i
f,<
R
. _ y
= long., v = vertical
I w-'
-_
- -
Breaker Tension Force
1
r
5
2x
Figure 21: Sample calculation for the steering system of the froni wheel of a moiorcycle
contact particles of a flattened tire. Wobble vibrations produced on the testing stand show
that the phase angle between the amplitude of wobble and the restoring torque of the tire can
certainly pass beyond go", which cannot be represented by tire behaviour of the first order.
The tendency of the system to auto-excitation is thereby markedly reduced when there is a
sufficient area of contact between the tire and the ground. If the ground is uneven, however,
I so that the contact surface is at times considerably reduced, one speaks of a transition to
I hard auto-excitation, since the stabilisation is now only produced by the kinematic terms
of f i s t order. At high speeds, therefore, the reduction of load on the front wheel due to
I the air resistance of the vehicle necessarily leads to critical states of wobble, while the axial
load is sufficiently large in the normal driving state. The non-linear coupling between the
horizontal and vertical dynamics of the wheel also calls for numerical treatment here, since
an analytic procedure only permits calculation in the case of perfect roadholding. Such a
6-22
numerical procedure, however, is very costly in term of computing time, and has only been
carried out until now with coarse substitute models with few degrees of freedom. If, however,
one starts out with an analytically based strategy, then the destabilisation effect can also
be demonstrated here by means of numerical examples for integrating the rolling action on
rough ground.
6 Conclusion
Starting from that given overview furture theoretical investigation should include. It is
possible to stabilise the system even in case of small contact area using the main influencing
parameters found by linear theory. If the parameters are near stability border the nonlinear
effects may suddenly disturb the rolling wheel and self-excitation occurs due to vanishing of
vertical spin effects in contact area. In any case also for nonlinear computation the linear
theory can be used as a guideline. All possible non-smooth surfaces could be taken into
account using discret system modeling for the tire.
Lit eratur
[Bak87] E. Bakker, L. Nyborg, and H. B. Pacejka:
Tyre Modelling for use in vehicle dynamics studies.
SAEP 87 04 21, Society of Automotive Studies, Inc., 1987.
[Bak89] E. Bakker, H. B. Pacejka: ..
A new tire model wih an application in vehicle dynamics studies.
SAEP 89 00 87, Society of Automotive Studies, Inc., 1989.
[Boe77] F. Biihm:
Computing and measurements of the handling qualities of the belted Tyre.
Proc. 5th VSD 2nd INTAM Symposium, Vienna, pp 85-103, Swets and Zeitlin-
gen, Amsterdam and Lisse, 1978.
[Boe93] F. Biihm:
Reifenmodell f i r hochfrequente Rollvorgfmge auf kurzwelligen Fahrbahnen.
(Tire model for high-frequency rolling processes on short-wave tracks.)
(Bericht Sonderforschungsbereich 181, Teilprojekt B1)
VDI-Berichte Nr. 1088, pp 65-81, VDI-Verlag, Diisseldorf/Germany 1991
[KochSO] J. Koch
Experimentelle und analytische Untersuchungen des Motorrad-Fahrer-Systems
(Experimental and analytical assessment of the Motorcycle-Driver-System;PhD-
Thesis 1978).
Fortschr.-Ber. VDI-Z., Reihe 12, Nr. 40, 1980.
6-23
[NYb731 G. Nybakken:
Investigation of Tire Parameter Variations in Wheel Shimmy.
PhD-Thesis, The University of Michigan, 1973.
[Pac66] H. B. Pacejka:
The wheel-Shimmy phenomenon.
PhD-Thesis, Technische Hogeschool, Delft, 1966
[Ruf831 E. Ruf
Theoretische Untersuchungen des Federungsverhaltens von Vierrad-StraBenfahr-
zeugen.
(Theoretical investigations into the suspension behaviour of four-wheel road ve-
hicles. )
Fortschr.-Ber. VDI-Z., Reihe 12, Nr. 44, S. 104-105, Dusseldorf, 1983.
[Schu87] D. H. Schulze:
Instationare Modelle des Luftreifens als Bindungselement in Mehrkorpersystemen
fur fahrdynamische Untersuchungen.
(Non-stationary models of the pneumatic tire as a linkage element in many-body
systems for investigations into the dynamics of driving.)
VDI-Fortschritt-Berichte, Reihe 12, Nr. 88, VDI-Verlag, Dusseldorf, 1987.
[Sha71 R. S. Sharp:
The Stability and Control of Motorcycles.
J. Mech. Engin. Sci., 13 No. 5, pp 316 - 329, 1971.
[Sha771 R. S. Sharp and C. J. Jones:
The straight running stability of single track vehicles.
Proc. 5th VSD 2nd INTAM Symposium, Vienna, pp 334-342, Swets and Zeitlin-
gen, Amsterdam and Lisse, 1978.
[Sun901 F. Sun:
Dynamic Properties of Motorcycle Tires.
Dissertation, Beijing Institute of Technology, 1990.
[Valk93] R. van der Valk, H. B. Pacejka:
An analysis of a civil aircraft main gear shimmy failure.
Vehicle System Dynamics, 22, Swets and Zeitlingen, BV Lisse, Netherlands, 1993.
6-24
1
7-1
Paper presented at the 81st Meeting of the AGARD SMP Panel on The Design, Qualification and Maintenance
of Vibration-Free Landing Gear, held in Ban& Canada from 4-5 October 1995, and published in R-800.
1-2
These examples occurred over I O years ago. While this may be actuator and the second, the more classical process addressed
considered ancient history by some standards - the computer during airplane development.
industry, for example - it is the authors opinion that there has
been little progress in the state of the art of shimmy design in the In the first case, shimmy occurred at relatively low speed and
intervening years. produced alarmingly large lateral oscillations of the airplane
nose. One report tells of tire marks on the runway indicating
The authors havereason to believe that the shimmy problem is angular deflections of up to 40 from centre in one direction and
perceived by many to be a black art with no basis in sound up to 35degrees in the other. This rather violent form of
engineering methods. The authors disagree and are concerned shimmy was caused by backlash at the tiller inpur potentiometer
that this black art perception has hindered progress in shimmy to the steering actuator and a modification was introduced to
design unnecessarily. The shimmy problem is no more myste- eliminated the backlash. This type of system induced shimmy
rious than, say, aircraft flutter with which the authors have expe- problem will not be discussed further.
rience. What is lacking is the accumulation, integration and
dissemination of currently available analytical and experimental The second type of shimmy, unlike the first, occurred at rela-
experience to support the development of practical design pro- tively high speed and with fairly modest amplitude oscillations
cesses. These, in turn, would yield structurally efficient, shimmy- of the gear assembly at a relatively high frequency (approxi-
free designs. The authors are optimistic that practical design pro- mately 20 Hz). It usually occurred with the steering system ON
cesses are feasible and believe that the aerospace industry would (ie., the mode with lower damping), but there were reports of
benefit if those processes were formulated and made generally shimmy with the steering system OFF as well.
available.
This latter shimmy problem stimulated a broad range of immedi-
Presented in this report are reviews of the DASH 7 and DASH 8 ate and longer term actions, all intended to achieve an expedient
shimmy problems experienced by de Havilland. The report out- and lasting solution. The actions taken were:
lines the analysis models and methods used and presents the sta-
bility results obtained from those methods. The report concludes (1) Introduce (large) unpaved runway tires across the
with a recommendation for the development of advisory material fleet because there were fewer reported incidents of
on shimmy design. shimmy involving these tires.
the purpose for presenting the above list is to indicate the effort As indicated in that figure, effective values for rotation arm
required to resolve a shimmy problem when it is experienced in length, inertia, damping and frequency are dependent on the
service. Following are some experimental and analytical results length of the shock strut. Coulomb damping or free play is not
which were important in helping us develop a physical under- included in the side displacement model.
standing of our nose gear shimmy process and in establishing our
analysis methods. The equations of motion for shimmy analysis, including the von
Schlippe and Moreland tire models, are presented in Appendix
4.2 Shimmy Test Results A. The tire coefficients used in the equations derive directly from
The test results presented in Figures 2, 3 and 4 are instructive. In Reference 2, except for Morelands tire time constant, 7, , which
all cases, the ground speed is about 55 knots. Figures 2 and 3 was selected to conform with test results. The structural damping
apply to a tight gear for which free play in the scissors - the ratio, 5 , . was assumed to lie in the range zero to 0.03.
only significant source of free play - was removed using shims.
In the case of Figure 2, the weight-on-wheels switch was deacti- 4.4 Results
vated to allow excitation by the spin-up transient and steering Stability analysis results for both the von Schlippe and Moreland
was ON at the outset. The result clearly shows the landing gear is tire models are presented in Figure 6 for the case of steering sys-
fundamentally unstable at 55 kts in this low-damped configura- tem ON and no free play. As noted, all curves except one apply
tion. for a single point model of the nose gear. Also plotted in Figure 6
is the damping result obtained from the test data given in
In Figure 3, the landing was made with the time delay circuitry Figure 2, a single point at 55 knots ground speed with a damping
operating in order to eliminate excitation by the spin-up tran- value of 0.007.
sient. With steering OFF at touchdown, the large transient oscil-
lations due to spin-up are observed to attenuate quickly, due to The predicted frequencies are not shown but lie in the 26 Hz
the relatively high damping in the castering mode. While small range. This is about 20% higher than the value measured and is
random vibrations persist, there is no evidence of shimmy after attributed to mass and stiffness modelling inaccuracies. The
the time delay period when steering is ON, showing that the time effect of this higher torsion frequency was determined by analy-
delay was effective in circumventing the excitation of shimmy sis to increase the predicted stability by 0.01; i.e., to lower the
by the wheel spin-up transient. curves in Figure 6 by that amount.
Figure 4a shows results for which the shims in the torque scissors With respect to damping, three curves are indicated to give rea-
were removed to produce approximately 1/8 inch of free play in sonably good agreement to the test result in the 55 knot speed
the scissor apex. Again, steering was OFF at touchdown and the range: the von Schlippe single point result for 6, = 0 , the von
transient oscillations are observed to decay to limit cycle oscilla- Schlippe distributed mass model result for 6, = 0.03 and the
tions within the free play range wherein, steering actuator damp- Moreland single point result for 5, = 0.03 and T~* - 0.002.
ing has negligible effect. With steering ON - activated after the Clearly, much better agreement with this single test point could
0.75 second time delay - and with torsion damping thus have been obtained both in frequency and damping by applying
decreased, the oscillations within the free-play band cause the arbitrary adjustments to the parameter values initially chosen.
gear to go unstable, that is, to shimmy (Figure 4b). However, this was not pursued because the problem remained
that no conclusion could be drawn as to which structural and tire
4.3 Analysis Model models would give the most reliable result due to the lack of test
The analysis model used to investigate shimmy follows the data.
approach used by Lockheed in Reference 1. It consists of a tor-
sion model defined relative to the inclined leg axis and a side dis- The analysis was extended to investigate the effect of free play
placement model (Figure 5). The torsion model is made up of 4 using the distributed mass model (Figure5a) and the von
concentrated inertias to give the freedoms required to represent Schlippe tire model. These results demonstrate why free play
free play (backlash), bearing friction (coulomb damping) and must be avoided in landing gear design. Its presence can result in
torque scissor stiffness and free play (both functions of stroke). limit cycle vibrations at speeds lower than the classical shimmy
Fuselage torsion stiffness and trunnion attachment free play are speed - the speed obtained from the solution for a tight gear.
also included in this model. Initially, the steering actuator stiff-
ness was defined solely by the value of bulk modulus for hydrau- The addition of mass balance to the shock strut, was investi-
lic fluid which was assumed to have a value of approximately gated as a means for improving nose gear stability. The results
100,000 psi (with typical air entrainment). This value was subse- presented in Figure 8 show that the forward placement of weight
quently reduced to about 34,000 psi to account for the flexibility has a stabilizing effect and that the most effective arm length, at
of the steering actuator structure determined from an isolated least for this landing gear, is equal to the tire radius. This
shake test of that unit. approach ultimately led to a modification - using a 25 pound bob
weight - which was made available to operators as an option.
The side displacement model used eliminates the y independent While the increase in stability is predicted to be only 0.01, it is
degree of freedom by assuming an infinitely stiff leg free to roll the authors opinion, based on reported performance, that the
( @) about an effective rotation point as illustrated in Figure 5b. mass balance modification is more effective than predicted.
5. DASH 8 MAIN LANDING GEAR SHIMMY EXPERI- Increasing the side bending stiffness was also determined to be I
ENCE important but not as important as increasing the torsion stiffness.
5.1 The Occurrence of Shimmy It was also determined that shimmy stability could be improved
by increasing the mechanical trail so the design was changed to
Early in the design phase, the vendor for the main landing gear,
increase the mechanical trail to the maximum extent possible.
using analysis methods which they considered reliable, reported
that shimmy would occur within the operating speed range of the The results from more elaborate analysis models indicated that a
airplane and that torsion and side bending stiffnesses would have
torsion stiffness of 2 . 2 ~ 1 Ib-idrad
0~ and a side bending stiffness
to be increased to achieve stability. Because of the weight
in the range of 6,000 to 8,000 Ibhn would be required to preclude
increase involved, many months of intense analysis effort and shimmy. By mid 1981, these stiffness requirements grew to a
debate were required before de Havilland and the landing gear
high of 2 . 4 ~ 1 0Ib-idrad
~ and 14,000 I b h . As noted above, it
vendor agreed on final design values. was only after a long, protracted period that the final design val-
ues of approximately 1 . 8 ~ 1 0Ib-idrad
~ and 7,600 I b h . with
5.2 Analysis Models
nacelle flexibility included, were agreed to by both the vendor
The main landing gear for the DASH 8 airplane is shown in Fig- and de Havilland.
ure 9. The DASH 8 is a high wing aircraft with the main landing
gear supported from within the nacelle. As a result, the landing Final analysis results are presented in Figures 1 1 and 12. Stabil-
gear is relatively long and flexible and therefore may be more ity in the side bending mode is of little interest because the
prone to shimmy. damping values in this mode are high and increasing at speeds
above the no-slip speed. Regarding torsion stability, Figure 11
During the course of the shimmy investigation, a variety of gear shows the rigid wing giving generally conservative results com-
structural models were developed ranging from simple stick pared to the full finite element analysis. Figure 12 shows stability
(Figure IO), to full finite element. In all cases, except for the in the torsion mode increases with stroke. The results indicate the
modal analysis discussed below, these structural models were landing gear to be stable for the complete operating speed range
used to define a 3x3 stiffness (or flexibility) matrix defined at the of interest with a damping margin of 0.03.
axle centreline in the single point analysis method.
In the modal analysis method, the finite element model for the
landing gear and nacelle structure was extended to include a stick
model for the wing with distributed mass and stiffness properties.
The mode shapes and frequencies forthis system were calculated
and the shimmy equations were then set up and solved in modal
form.
Summarized in Figure 13 are the results of a non-linear analysis
to establish the influence of bearing friction and free play on
shimmy. These results contributed to the specification of free
play limits in the main landing gear design.
6. DISCUSSION
I
I
In the case of the DASH 7, it appears that the nose landing gear - -
The single point analysis method was extended to include free problem has been contained provided the tires are matched and
play. Its effect was simulated (in a global sense) by replacing the in good condition, the wheels are kept in good dynamic balance,
linear stiffness matrix defined at the wheel axle centreline with a free play is limited and nose gear touchdown and taxi speeds are -
non-linear stiffness matrix having low angular stiffness, kyw , kept low. Nose gear mass balance is generally considered to be
within the free play range and full predicted stiffness outside that beneficial by those operators who elected to have that option.
range. These equations were solved in the time domain. Obviously, the containment of nose gear shimmy has increased
the maintenance cost, pilot workload and in some instances, the
Both Moreland and von Schlippe tire models were used to repre- airplane weight.
sent tire dynamic properties. However, the Moreland model was
considered less reliable for absolute stability results because In the case of the DASH 8, with over 400 airplanes in service, all
there was no empirical method for defining the tire time constant, shimmy free, it is safe to say that the main landing gear design
T,, (see comments in Appendix A). has now been proven. However, the main concerns in this case
are the time and effort required to establish the final design
For structural damping, a typical value of 6, = 0.03 was requirements and the degree of over-design (weight penalty) that
assumed in the computations. However, no allowance for this may have been included in the final product.
damping was given in defining stability; that is, stability curves
had to have a margin of 0.03 to be considered shimmy free.
7. RECOMMENDATIONS
5.3 Analysis Results At the time the DASH 8 was being designed, there was no well
established consensus on what constituted a reliable shimmy
In January 1981, using the simpler stick models shown in
design methodology. This still appears to be the case today. For
Figure IO, de havilland confirmed that the landing gear would be
this reason. it is recommended that some guidance information
unstable within the operating speed range and that the most
be provided - similar to an advisory circular issued by aircraft
effective way to increase the shimmy speed was to increase the
airworthiness authorities. This advisory material should:
torsion stiffness of the gear (at the time 1.13~10~ Ib-in/rad).
7-5
Establish standardized terminology lined herein. At the very least, their resolution would have been
more expeditious.
Give guidelines on structural modeling
REFERENCES
Define tire models, including the assumptions and
limitations (1)Active Shimmy Control System, AFFDL-TR-75-136.
Provide information on related vibration processes (10) von Schlippe, B., Shimmy of a Pneumatic Wheel, NACA
and problems - sometimes confused with shimmy - TM-1365,August 1954
such as wheel unbalance, tire out-of-round. judder.
anti-skid processes, steering feedback processes,
APPENDIX A. SHIMMY ANALYSIS MODELS
gear walk, etc. Some of these processes may warrant
separate advisory material of their own.
This Appendix describes the shimmy analysis methods devel-
oped at de Havilland.
Provide landing gear maintenance guidelines, based
on analytical and experimental data, for minimizing
A.l Equations of Motion
landing gear vibrations
A much simplified model for a single or dual wheel landing gear
Provide a bibliography is illustrated in Figure AI which shows an inclined shock strut
with mechanical trail, L, measured parallel to the torsion axis and
geometric trail, LG = R x sineg, where Os is the inclination
8. CONCLUDING REMARKS angle and R is the deflected tire radius.
As a consequence of the DASH 7 and DASH 8 shimmy prob-
lems experienced at de Havilland, the authors were required to Also shown in Figure AI are two equally valid coordinate analy-
develop analytical tools for problem solving and design. This sis systems, one located at the apex of the shock strut and the
work was accomplished successfully despite the lack of a stan- mechanical trail arm and aligned with the strut, and the other
dardized industry accepted methodology at the time. It is the located at the wheel axle with axes orthogonal to the ground
authors contention that reliable shimmy design procedures are plane. The former system simplifies the mechanical terms in the
required and that it is feasible to do so. Had advisory material equations of motion while the latter simplifies the tire force and
been available during the design of the DASH 7 and DASH 8 moment terms. In order to expose more clearly the less familiar
landing gears, de Havilland may not have had the problems out- tire terms, the latter coordinate system is used in this report to
develop the equations of motion. These equations are written in A.3 TireTerms
terms of the lateral ( y), roll ( 4 ) and yaw (y)degrees of free- During shimmy, the motions of the gear axle are coupled not
dom defined at the wheel axle; only through the gear structure as expressed by the terms on the
left hand side of equation (AI) but also through the interaction
between the tire and ground, as represented by the tire forces and
moments appearing on the right hand side of that equation.
where U = [i] the ground plane. This was done in order to develop the simplest
possible expressions for the tire forces and moments. Formula-
tions for single and dual wheel assemblies are included.
E;]
force, Fc , and pneumatic trail, E ~ as, follows (Figure A2):
K = stiffness matrix
A.2 Structural Terms
Formulations for the left hand side terms are not given in this = (A3)
report but they can be found in References 1 and 3. It is impera-
tive that utmost care be taken in defining and verifying the
assumptions made in establishing the analysis model and in where
developing the data used. For instance, if a landing gear is Fc = K T ~ A +C T ~ A (A41
assumed to be a uniform structure, its mass and inertia effects
can be lumped at the wheel axle and a single point model ana- and where K T =
~ tire lateral spring constant
lyzed. In this case the mass and inertias of the wheel rolling C T=~ tire lateral damping coefficient
assembly should be supplemented with appropriate effective = pneumatic trail parallel to ground
masses and inertias of the supporting gear structure. The 3x3 A = tire lateral deflection (elastic)
mass matrix, M, will contain (off-diagonal) terms that couple the R = tire (deflected) rolling radius
three degrees of freedom, y, 4 and y. If it is important to
include detailed representations of bearing friction or free play To enable a solution of the shimmy equations, Moreland postu- -
(backlash), then a distributed mass model will be required to pro- lates an alternate expression for the tire cornering force in terms
vide the additional degrees of freedom needed to represent these of a tire cornering power, C y , and time constant, T,, , as follows:
characteristics.
Fc = CwWt + CtyTu$t (45)
The damping matrix, D, should include both structural and bear-
ing friction damping terms, if appropriate. As indicated ahove, where y, = the tire elastic twist deflection
this may require the development of a more complex gear
dynamic model than can be expressed by a single point formula- Finally, Moreland assumes the absence of slipping or scrubbing I
tion with 3x3 matrices. between the tire and ground by imposing the kinematic relation-
exception - the tire time constantrv. Choosing a reliable value in terms of the yawed-rolling relaxation length, 2.0. This geo-
for rY is generally a problem because no consistency can be metric relationship can also be deduced from Figure A4.
found in the literature for this parameter as illustrated by the
compiled data of Figure A3. Since r, has a strong influence on
the results, it is recommended that the value used should be vali-
dated in some way if the results are to be accepted in absolute
terms.
Empirical values for all of the parameters required by the von
Also. while empirical values for the pneumatic trail, Ep. are Schlipp tire model can be obtained from Reference 2. It might
available. some analysts prefer to base the tire torsion moment, be noted that some analysts also extend the expression for M Iby
.
Mr directly on the tire torsional stiffness K T ~viz;
, including the tire torsion damping term, C T ~viz;
,
A 5 The von SchUppe Stretched String Model (SingleWheel) A.6 Dual Wheel G e m
The von Schlippe tire model (References IO and I) assumes that The Moreland and von Schlippe tire models given above apply
the tire equator ISrepresented by a "stretched string" that takes for a landing gcar with a single wheel. The following steps
on the deflected shape shown in Figure A4 due to out-of-plane extend those formulationsto uual wheel landing gears:
E] E]
mouons In this model the tire forces and moments are:
(1) Factor the force and moment expressions by 2
F,R = F.
2 +L ~ @ K ,
(A171
F
As before, h equals half the footprint length. F.L = ;-L~@K.
Z and 2 are defined in Figure A4 from which the following geo-
metric relationships can be deduced. and where LR = half axle length
Kz = the tire vertical stiffness
Z = Y-y-R@-hW
(3) Add the following yawing moment term to account
2 = ?-y-R$+hW for the ground friction force;
AM^ = L~~(F,L-F,R)
Finally, the kinematic relationship that assumes no slippage
between the tire and ground is given by the following expression where p = tire rolling friction coefficient
E9
7-9
7-10
.
7-1 I
.. .. ..
: : : :
:i: - I
; : t i ; -
: r
7-12
I 1%
3>
6 t
!Ii aa
5?
Y
8- 1
Unsteady Tire Dynamics and the Application Thereof to Shimmy and Landing
Load Computations
Klaus Koenig
Daimler-Benz Aerospace Airbus GmbH
HunefeldstraBe 1-5
28199 Bremen, Germany
Paper presented at the 81st Meeting of the AGARD SMP Panel on The Design, Qualification and Maintenance
of Vibration-Free Landing Gear, held in Banff, Canada from 4-5 October 1995, and published in R-800.
8-2
For larger deformations as they occur at landing impact theories use a relaxation length [m]. This leads to the
further investigations with larger amplitudes and with and fact that e.g. B2.. depend only on the rolling velo-
without rolling velocity are necessary. city v in the Mooreland theory but on v2 in the other
But this should be a second step to be defined in detail after theories.
successful execution of the first one with small perturba- - Moreland originally only studied I$. The extension
tions. by y was introduced by Rogers and Brewer, and the
Partially such studies have already been performed (s. ref. extension by $ by Edman (9).
5 ) but no complete test result is known to the author.
Nevertheless, since no principle obstacle to such an ap- 4. A COMPLETE SET OF TIRE EQUATIONS
proach to establish reliable equations of tire dynamic is ev- For the following, Smileys equations as given in chapter
ident, it can be hoped that perhaps in another 50 years such 3 will be used as a basis. They are modified by neglection
data will be available for the benefit of future engineers of minor unknown elements (t=O) inclusion of informa-
and of technical progress - if further progress is ever tion from ref. 2 (WO,,, WO,, WO,,), assumption of sym-
wanted at that time. metry (Wl,,, W2,,) and some other simplifications (WO?,
: analog to Wl,, and Wl,,, W2+, : reduction to the main
3. COMPARISON OF SOME PUBLISHED TIRE element). But this is not sufficient. Even for tire rolling in
EQUATIONS steady state one needs information for the x,z and 6 com-
The following publications will be compared: ponents. They are taken from ref. 2 and 6 and mainly based
- von Schlippe, and Dietrich, (s. ref. 1) on radial tire stiffness and respectively the assumption that
I - Smiley, (s. ref. 3) x loads depend on a friction coefficient and that this coeffi-
, - Pacejka, (s. ref. 4) cient itself depends on the slip ratio as shown in Fig. 4-1.
- Rogers, and Brewer, (s. ref. 5 ) with the Moreland
theory. For the not free rolling tire in the spin up phase of landing
If the equations of movement as given in these publica- the situation is more complicated. Nearly nothing is
tions are transformed to the coordinate system of Fig. 3-1 known. The only fact is, that for the pure sliding, all hori-
and linearized by series developments up to elements of zontal ground loads depend on the vertical load and a fric-
second order coefficients for matrices B and W as defined tion coefficient while their direction is opposite to the
in chapter 2 can be derived as given in table 3-1. movement.
I These equations are given in their original symbols with-
out further explanation of the symbols. To have at least something in hand a very crude assump-
Nevertheless some basic conclusions can be drawn: tion was made. It defines that during the spin up phase both
- There are really only diagonal elements in the B types of loads are acting but the one for pure rolling with
matrices. a factor between 0 and 1 increasing with the slip ratio and
- B.,, = B.,, = B.,, are equal for all theories . This the other one for pure sliding decreasing in the same way
may mean that the tire can vibrate in itself with a down from 1 to 0 until the slip ratio reaches its zero for the
single mode, showing at least 3 different load compo- first time.
nents More details of the derivation should not be given here to
I
- The different theories dont place coefficients at the avoid any boredom. But the final result is presented in
same places. table 4-1.
- A larger portion of the matrices is non-symmetric. These are rather complicated non linear equations but they
, - Von Schlippe and Dietrich do not study the tilt. should be valid for the whole landing process of horizontal
Pacejka does not give tilt moments . landings.
- Remarkable are the missing coefficients Wl,,
W2,, in the theories of Smiley and Pacejka, the mis- 5. SHIMMY ANALYSES
sing coefficient W2w, WO,,, WO,,, WO,, in More- 5.1. Tire Equations
lands theory and the missing coefficient W q , in First linearized tire equations should be derived from the
I
Pacejkas theory. nonlinear set given in table 4-1.
- The absence of the elements WO,,, WO,, WO,, is Details are given in table 5.1-1. This leads to extended
caused by the relaxation mechanism. Loadings equations of unsteady tire loads for all the 6 degrees of
caused by sideward deflections will disappear with wheel freedom.
time.
- The theories of von Schlippe and Dietrich, Smiley and 5.2. Influence of Tires on Shimmy
Pacejka look rather similar, they are based on equal To study this it is reasonable to look at a rigid undercar-
basic thougths. The Moreland theory is funda- riage to avoid coupling of structure and tire elasticity. For
mentally different the most simple rigid model of Fig. 5.2-1, after some rear-
- Moreland introduces a time constant C1 [sec] for de- rangements, one obtains the equations of motion
scription of the relaxation mechanism while the other
8-3
approach to Variant S.
Used are (Variant E) :
Type I11
o = z0 wi,, -
(2 d, q,, + rr e)-Ip e PR(O)= 1765.8 N
zp = 0.02 m
ro =0.13m
~0 = 0.12 m
+ is,(Wlqq-LR
* * (Wlqy + W1,) + L i * Wl,,) + dR = 0.085 m
pi = 225630 N/m2
+ qst (w2qq-LR (w2qy + w2,) + L i w2,,) Pr =Pi
with
=o
pk =0.8
* p = pq(R1) + pq(R2) - LR (py(R1) + py(R2)) sk Z0.2
m, = 2.4525 kg
8 = 8,(R2) - 8,(R1) IR = 0.05 1993 kg m2
Iw = I, +
m, (di Lk) + I, = 0.0981 kg m2
m,=O
I, = radial moment of inertia of wheel including tire
I, = polar moment of inertia of wheel including tire For the single wheel case the estimation did not lead to ex-
actly the same stability areas as Variant S.
d, = swivel damper constant For the double wheel first only the wheel distance dRwas
introduced to show the influence of the increased swivel
An eigenvalue analysis of these equations delivers fre-
inertia 2.rnw.d~z.Then the polar inertia I, of the wheel was
quencies and modal dampings. Different cases were ana-
introduced and the influence of the swivel damper shown.
lyzed and compared.
~
Finally the tire coefficients W(x,e) are added.
As one can see there are some influences. The shimmy ve-
Variant S :
I Ki = 44145 N/m
K, = 298.224 Nm/rad
LI = 0.142222 m
locity of the double wheel is smaller than that of the single
wheel. Both the swivel and polar inertia contribute to this
effect.
i -
I
!
-
m, = 2.4525 kg
K,= 0
WI,, = 0
and increases the shimmy velocity.
craft centre of gravity. pLis the total lift. The tire loads are 5. Rogers, L.C. and Brewer, H.K.
PX,R,p y l ~P, Z ~ Rand the undercarriage loads at the wheel Synthesis of Tire Equations for Use in Shimmy and
centre are pxuu,pyUu,pzuu.The P.~Rare measured parallel Other Dynamic Studies
to the ground, the peUu are measured in a coordinate system Journal of Aircraft VoI.8 No.9, Sept.1971, p.689-697
of the undercarriage. The undercarriage loads include the
tire loads and the wheel inertia loads.
In Fig. 6-1 one can see first of all how different the loads
6. Horne, W.B. and Leyland, T.J.W.
INFLUENCE OF TIRE TREAD PATTERN AND
RUNWAY SURFACE CONDITION ON BRAKING
c
of the two undercarriages are. FRICTION AND ROLLING RESISTANCE OF
Further, one can clearly discern that the loads oscillate and MODERN AIRCRAFT TIRE
that this oscillation was obviously induced by the wheel NASA-TN-D 1376
spin up / spring back effect.
The large sideloads are remarkable. They reach their max- 7. Koenig, K.
imum after spin up - on the rolling tire. They are nearly Die Lasten des Landesstoaes
doubled on the undercarriage if compared with the loads ZFW 3 (1979) H6, p.344-360
on the tire and they dont always decay with the vertical
loads. They are mainly caused by the tire dynamics and 8. Krabacher, W.E.
can not be predicted without such a complete set of equa- Aircraft Landing Gear Dynamics Present and
tions as presented above in table 4-1. Since the whole Future
movement of the aircraft in this unsymmetric landing case SAE-paper 93 1400 (1993, Dayton Conference)
is determined by the sideloads one can conclude that no
prediction of the unsymmetric impact loads which is in the 9. Edman,J.L.
least reliable can be computed without inclusion of the Experimental Study of Morelands Theory of
complete tire dynamics of the rolling tire under varying Shimmy
vertical deflections. WADC Techn. Rep. 56-197 (July 1958)
A measurement on the flying aircraft instead comes com-
paratively late and is tremendously more expensive.
cn
W
z
0
W
r
k
n
w
I
cn
i i
m
I- 3
e
U,
0
2
0
cn
z
8-6
wheel centre
A (axle, non rotating)
R (wheel, rotating)
definition of working point:
Pr = O (free rolling)
a =O (free rolling)
T- s(0) = 0 (free rolling)
y,(O) = Z,(O) = $JO) = $,,(O) = 0 (straight on)
zt 20 zrn
conditions of linearization:
< zpc < Zpk (no bottoming)
V
ro = zpx= zpy= zpv = 0 (simplification)
plane with ground plane
b, = 0 (linear vertical stiffness)
w.... = w....(0) (as would be measured)
B .... = B ....(0) (as would be measured)
l/X, E (1 - AX,/V)/V
abbreviations:
XO
(, ) = v
ztA(0) - 0,25 * pR(0)/Kx(O) = zo
bottom dew
rr = ro - zp(0)/3
additional tire coefficients for the extended
unsteady tire loads
Fig. 3-1 WHEEL COORDINATES
P
wlxx = - Pk pR(o)/(sk *
'
A = - pk PR(o)/(3 sk * rr)
wlxO = - pk ' pR(O) * rr/(sk
WO,, = - K,,(O)
WI,, = - 0,1 ,/E
wlOx = - V k pR(O) ' z(o)/(sk *
= - Pk PR(O) 2(0)/(3 sk * rr)
wlOO = - P k * pR(O) * rr zo/(sk *
slip ratio
1 Bl(x,z,8) = B2(x,z,8) = 0
Pi tire pressure
P' rated pressure
Pk
ClS
maximal friction coefficient
full gliding friction coefficient
Y \
0.20
I
-E
sE 30
3
-
g 20 -._
2 0.16
10
0.12
unstable . .'
. .
0.08
I 1, 51 shimmy mode (flutter critical yaw)
+ f2' 52 stable yaw (coupled tire mode)
: . -.____ _ _ _ - _-_
e l3 EJ __ decoupled tire mode
~.
0.04
0
1 2 3 4 5 6 7 8 9
velocity v [m/sec]
Fig. 5.2-4
INFLUENCE OF x AND 6 MOVEMENT
Fig. 5.2-2 variations :
0 single wheel W(x.8) = dR = ds = 0
EIGENFREQUENCIES AND 0 double wheel W(x.8) = Ip = d. = 0
I '8 W(x.8) = ds= 0
MODAL DAMPINGS +A 0)
"
'8
W(x,8) = Ip = 0; ds= 0.06
variant E ; ds=O.oS
single wheel ; rigid u/c
Smiley case I approx. B
d, = O.OGNm/(rad/sec) ; LR= 0.08 m
I I - .
unstable stable
- area 11 I I area I
= s s s s s s s
0;s a b b b 0;s 0;s
Fig. 5.2-5
INFLUENCE OF y AND q~COEFFICIENTS
Fig. 5.2-3 STABILITY DIAGRAM single wheel o test results
single wheel - d, = 0.6 Nm/(rad/sec)
I-
8-9
V
Fig. 6-la LANDING IMPACT, LEFT MAIN UNDERCARRIAGE
I.
P. WOERNER, 0. NOEL
Messier-Dowty S.A.
Zone ACronautique Louis Breguet
BP 10, 78142 VClizy Cedex
France
It is therefore necessary to evaluate and model landing as damped vibratory : because of damping influence, an
gear non-linearities in order to improve the pmiictability imaginary quantity appears, from which the denomination
of the shimmy phenomenon; an accurate shimmy analysis of complex modes. This damping term can be
must take into account a range of non-linearities values as suppressed in equation (2) to obtain the equation of the
wide as necessary, in order to reveal the most unstable conservating system associated :
configurations. Mx +m=o(3)
Paper presented at the 81st Meeting of the AGARD SMP Panel on The Design, Qualification and Maintenance
of Vibration-Free Landing Gear, held in Banfi Canada from 4-5 October 1995, and published in R-800.
9-2
11 -
-
The structure is loaded on one or many excitation nodes.
Tbe input signal is sinusoidal,transient, or random, and is
mass-spring system, a funuion of position but not of
velocity; on the other hand, with no contact between such
peaks,slipping occurs, and this force depends only on the
rnulhfxquency : all the modes can be simultaneously velocity. A good model must take both aspects into
- excited. account :cf reference [l].
I
I The modal analysis is performed by aaos&r function Elementary representation of friction as a function of
measurement and processing : the theoretical explanation relative velocity alone is shown below.
of this method follows from solving equation (4) and (5)
with whatever load form :
X = # q (4)
m q + c q + k q = f = # T F (5)
According to 2.1,this can be explained as a transfer matrix
FnFl-
between input F and output X,whose terms are of the form 4
N
# 4%
~ i (0)
j = C,:-
r-1
0 2 +j x, U 0,
Standard stick-slipfriction modcl
There is a strong nonlinearity around &=O, with a
multiple evaluation of friction between -Fh and Fh : this
The term Hij(U) links the dynamic -LUX in the is unacceptable for numeric computations,which impose
accelerometer localizsd in d.0.f. i , to the excitation the linearizationof such a zone.
applied in d0.f. j .
When the transfer maaix between input and output is
measured, we have to find an analytical expression such
as mentioned above, which best fits the experimental
result (identification process), then to deduce eigen
pulsations q and damping coefficients tr . The
determination of #, t e r n can be ma& easier by
measuring responses on tfK same points as excitation :we
have i=j, thus :
9-4
The friction nonlinear model must have a good physical Fex t(t )
interpretation of its parameters, as well as an easy
implementation in numerical software : a satisfactory e
representation for friction is the DUM model (1%8),
described in reference [2], in which friction force is the
solution of following first order differential equation : I f m is the inertia moment relative to torsion, described by
wheel r n m s with respect to vertical, we have a system
of 2 equations with 2 unknown values F and 8, with 2
valuesf, and Kj to identify :
I
OO 0.1 1 1.6 2 28
is:fi=kJG K +k
,therefore: fi=1.81Hz frequency is : fa = &&
K +k
,therefore : 4 = 1.81
HZ
- the second one, associated with higher excitation
levels (cf c w e with 100 Nm), when sliding can appear The system gets more stiff when excitation level grows;
; the stiffness Kj does not occur any more. This the frequency suddenly moves fromfi t o 6 when the
excitation level grows (4 > fi ) , which is often called the
frequencyis: fa = &&,therefore: A =0.87Hz "jmp"phenomenon
The system gets less stiff when excitation level grows: it This freeplay behaviour is not the same as friction
progressively slides between fr ami when excitation behaviour, whicb helps to characterize these two kinds of
level grows (12 fr ) . non-linearity during sweeping frequencies tests.
A second test campaign was performed to obtain a more - tbe second zone is governed by freeplay : dqlacement
complete dynamic characterization of the landing gear, in increases a lot, which shows a strong non-linearity ;
order to validate the numerical model on many modes and friction influence is less sensitive : frequency drop noted
to make it more accurate. Tests with transfer function in first zone fades.
measurements were well-suited, to measure many - the third zone approaches linear behaviour :
frequencies initially unknown. displacement increases lineary, and frequency reaches an
asymptotic value.
4.2 Identification of non hearities with respect to These tests with appropriation have shown that,when the
tests with appropriation studied landing gear has this linear behaviour, torsional
frequency is about 16 Hz and lateral flexion frequency
4.2.1 Experimental characterization of landing gear with about 14 Hz, which are close values. This result will be
appropriation recalled further in non-liaear simulation analyses (cf 5.3)
Furthermore, this test permits the deduction of damping C
Messier-Dowty SA. has performed a modal analysis on a
clamped main landing gear with no contaa between
in linear behaviour from force and displacement
measurements,according to relation shown in 23.1 when
wheels and ground, by means of appropriation of the
thewkmodeisexcited: = j q ~ k Ca
excitation on a wheel center, with 10 levels of excitation
force. Landing gear experimental curves show a 4.2.2 Friction parameter identification
non-linear behaviour in torsion, with dependance of
frequencies and displacement on excitation amplitude :in In addition to linear characterizationof landing gear, this
particular, for low values of excitation, frequency kind of test permits the identification of n o n - b a r
decreases, but reaches an asymptote for higher values of parameters to be introduced in a numeric model to fit the
excitation, for which the phenomenon becomes three zones previously identified.
quasi-linear. In panicular, we have to determine& and Kj values so that
the solution of above equations(6)and (7)(cf 33) fits tests
measurements with !iquency sollicitation : one
possibility is to solve this equation analytically and
Kf values , and looping on those 2 values to miaimize the
I gap between numerical solution and test result :this leads
I to relatively long computations. - -
i 1st ZOlK : A better way consists of finding a straightforward .
I approximated solution of these equations, insofar as the
e x t e d loading is sinusoidal : the assumption is that - -
response will be h o m c too, with the same pulsation :
this is the "harmonic linearization method", or "Rib
averaging criterion" : cf references [3] and [8].
If Fd (t) = FuI sin ut,we have F and 8 in the form :
F = Fccosot+Fssinot
e = u c c o s u t + ussincut
m i j + c ( e @ + + e p + ~ - ~ , s i n ~= t o (7)
These expressions of F and 8 are only approximate
9-7
(% 42
.2
.4 Limitations of tests with appropriztion
Despite the relative ensiness of these test proadure. many
limitations exist :
-in 0lde.r to reach asymptotic Values forfrequencyand
These 4 integdons can be computed numerically by damping, &on level most be sufficiently
s u m F and 0 by theii appmrimw cxpnssions : high, all the more as the landing gear has a more
they cancel out for @mi& v a l w of weighting factom n 0 n - h behavioolr : as long as frequencies depend
U,,U ~ , Fand~.,wttichdcfine~andesol~tio~.
, on excitationkvel, hear behsviour is not reached
- if excitation level is not sufficimt.identification of
4.2.3 Erample of experimentalfiningof nonlinear frictionparameters becomes less auauate.
landing gear model -duringthistestcampaigu,somemodeswerediffidt
to appropriate : their meammd fquencies atd
Theexperimeatallandinggeartorsionikquencyrespoase damping wuc not PCCOVsteenough.
m e is shown as a solid line: this frequency varies
between 16Hz (in the linear mne) and 24 Hz Besides, the identi6cation of fdction parameters shown
above may not be accmk enough when coupling rppears
For ea& &on kvel, friction parameters weLe between modes if coarespondhg friaion is not
adjusted to fit the experimental m e : stars indicate the neglipabk,iasofarasitshypomesisisthebehavimlike
frequencylevel readwd from W e values. a single 1 &of. model,llILcollpledwith other d0.f.
28
I
4.3.1 Experimental characterization of non-lineanties
with banfer%tion measmments
'b
I&m&n
.b M M . A &-a
P ~ ( W
3- In0
e'
9-8
Insteadofusingafrequencyrrsponaefrmction,liLemany
freqnency drmain algorithms, this method ia formulated
fromthe impulse functioo (due to Dirac impulse), which
may tc obtained fmm the invemeFwriertnwform of the
Sepoency nsponse fuuuiow the tbeonticalideotification
, expesedin the
process is the same 88 ~ w n i n 2 . 2 . 2but
time domain. One of the advaatagea of this methodis the
possibility to take many excitation mums mon easily
into account.
Transfer function meanurmmts are today in progress, so
we have not yet 0-d all thart ChSnaLnsa' 'cs. We
sbow below an example of mch identiEcarioIb in
agmmeat with the first 3 modes of an experimentalFRP
-thesecondfrequencypeakdependsontheexcitation
level, with a "jump" between two frequencies, whic4
also shows a noolioear behaviour, identical to the
freeplay 1 d0.f. model seen in 3.42.
-
Exprlnunt.IdidenU&dFRFVUL3&atndr.
As shown in 43.1. wch identificationmust bc made from
other tests on this system showednodhearbehaviour: with &&rmt &tation and diffaent
-wbenexcitingthesystemwithasinglefkxpency sweeping dinctions, asspr pry for &te&g
equal to this first peak, the response showed harmonic nO D - h W i t i G S .
cornponeMS of that frequency, which should mt
appiarina liaearsystan. 5. SAIMMY NUI)IERICALMODEL
- when exciting the system with a swept sine, some
peak frequencies depend on the direxion of sweeping 5.1 shtnmy lndd hypothesis
(haeasing or decreesing she),which shows -lay
,,* beha,+,u, with a ,*jump**pbenarmenOn 5.1.1 Degrees offncdom of shimmy mcdel
between the= frequencyvalues. A wmvlete shimmy numerical ' modelhasbeendeveloved
from +Vious &chnncteristics, including l&ar
43.2 Identification of non-lineanties fmm transfer
function measurements and non-linear behavioor of the main landing gear, in
a p m e n t with test measmments mentioned above.
The identitication of mn-lhdties consists of linding a This model takes into account 3 d0.f. of the landing gear,
model whose FRP (Fkqnency Response Fundom) agree located on the center of wheel axle :lateral displaQmeM
with all the experimental F W for each accelerometer : Uy,rotation Mx about ox longhdinal axis, whicb
this identitication is 0btaine.d with qe&c software, ~ple~en t laanl aexion, and rotation M. h u t m
both
integratedintothemeamwnentchaia vertical axis, which rrpresenrs tomion
Them are several ways to examine the measured This model must take into Dccount the following
frequency respome function to extract these data, which elemem:
are descrited in refezence 141. The numerical method used - s k l l d descriprion of landing gear
hen is called 'Tolyrefemce llme Domain Technique " -gymwpicefkcm
(cf nfezence [SI), and allows the extraction of modal - dylurmic tire
parameters, namely modal damping and Wuencies, -landing gear non-linearitiea
fromexperimentaldata. -appliedexciratioll
9-9
5.1.2 Structural description of landing gear This nowlinear system is solved iteratively. with a
dimetisation of the derived temporal terms (Newmark
Landing gear s t r u a m is taken into account with : schema) and 6iction f o m (Crank-Nicholsonschma).
- 2 manices 3x3 [M]and [K]obtainedfmm a dynamic The solution depends on airnaft velocity, which appears
condeosation (Guyan method) on the center of wheel
axle of the landing gear finite element model, whose in the gyroscopic metrix G and in the fn?quential tire
model P(s).
calculated modes agne with modal tests with
appropriation when linear behaviour is reached (cf
53 Readtn of model compntation
43.1)
- a damping matrix 3x3 [C] , directly built fmm From the previous temporal equation, we can obtain
m u w n d damping coefficients with the same modal torsion respooses of landing gear, for different non-linear
tests (cf 4.2.1) parameter values and airplane velocities :
5.1.3 Description of gyroscopic effects -torsional freeplay valuea :correspond to a madation
freeplay measured of a staadard point of measure,
We must take into Bccount gyroscopic effects due to whose valws are 0 , 0 . 5 , 1 , 1 5 mm.
mlling wheels of the gear,inthe form of a3x3 matrix [C]
which contains a term P V I R :Jis the moment of inenia - torsional fiction values : maximum forcefo values
for the system wheel + tyre, R is the tire radius, Vis the ive 0,25 ,50 Nm, with a constant friction slope Kj in
longibldioal airplane velocity. agreement with modal tests (cf4.2.2).
- airplane velocity values :45,55,65,75,85 mls.
5.1.4 Description of dynamic tire
The torsional damping can be deduced from logarithmic
A dynamic tire frequencial model was conveyed by its deament of temporal torsionresponse.
manufacturer, in the form :
The results can be summarized in diagrams dlowing
Fpnru(s)= P(s) U 0 = P(s) /UyCs)MAS) Mz(sJ1 torsional damping v e m airplane velocity :each diagram
P(s) is a 3x3 transfer matrix (s is Laplace variable), made corresponds to one freeplay value, with the 3 fiaion
up with 9 elementary transfer functions, whose terms A
values : x : ONm ; : 25 Nm ; :50Nm+
agree with manufactum dynamic tire tests, and depend 8,
explicitly on airplane velocity.
5.1.5 Landing gear non-linearities
Experimental results show more frequency variations in
torsion tests than in flexion : comequently landing gear
non-linearities, which are principdly fiiaicm and
&play, are only locaked in torsional d.0.f. MI, with the
expression described in 3.3 :cfequatiom (6)and (7) . -2O B 3 3 3 3
5.1.6 Applied excitation 40 50 60 70 80 90
Fn.lthyOnun
An excitation in moment along the oz axis is introduced
in the model,to load the landing gear in torsion :this load I I I I I
will be applied in fume shimmy tests.
8, I I I I I
6
4
2
0
-2-
40 50 60 70 80 90
F&Um 25
:. . .. , ... ,
.....
Tbis is cllncntly in pm- at Mem'er4bvly SA., with
amoreacauatenodinearchnnaarzstton * . of the llllldied
main Larding gear hnn teats with traosser fuaaion
meammeas, and witb the development of a complete
Y- lllmmld' model,insgreementwithmodaltan
d t a .
such unpredi- caused by mn-lioearities call also
be fouad in thco~tical~ramples,such as thc "D@g Fbhm ddmmy tern on a flywheel will enable the
spring", w h behaviour is highly d t i v e to initial validation of thb non-linear numeric model, cmd thc teats
conditions, so far as to become chaotic, in spite of a of Consprations able to improve landing gear stability.
comparstivelysimple behaviouraleapressim :
Tbisimpmvedmodclizatiwofthclnndinggcarwillallow
x*+o.asx"+=7.z5c~t to aaalyse ita dynamic behaviour in its devdopment. At
4
I this stage, it wiU be therefom convenient to delhe and
qurlifypoasible~~mtofthede~ifnecessrryto
a c b i ~the rrquestcd dynamic stability of me landing
gear in its range of operaiim colrdittons.
7. REFERENCES
[l] Computer Simulation of stick-slip friction in
Mechanical Dynamic Systems- D.Kamopp - Joumal of
Dynamic Systems, Mwsurements and Conml - Marcb
1985. VoL 107
[2] On the modelling and simulation of HcIim - D.A.
Hnessing, Jr. B.Riedlrod- Joumal of Dynamic Systems,
i -4
a S l O l s 2 0 a o S O I O * m
Measmementa and coonol- TMMctl'om of the ASME
-September 1991,W.113
- -
D u m # xplinr reax&iv@ to Idd*leodiUowx
x & h :X(0) -3 X'(0) - 4
l :x(o) 3.01 xyo) 4.01
h r ~ .linr
[3] Vi- -J ~ c o ~ S ~AI m~ ,
[4] Modal testiDg :theory and praCtire -DJ.E h -Jdm
Mo~eover,such chaotic behaviour has been expressed in w*y, 1984
pmrious landing gear nonlinear studies, in particular in
refe- [61 pamgraph 4. [SI The polydea'em b e domain teclmique - F.
Deblauwe, D.L.Bmwn, RJ. Allemang - MAC 1987 p.
However.it is imponam to notice rhat this qualitative 832,845
mente of mn-linearities CMOot apply to all landing
gem, insofar as toIsiw and latael fluioll iiquencies [6l Modeling and Analysis of a Dual-Whel Nostgear :
canbe more or less dose to each other; foriastance, if Shimmy Iostrbility and hnpaa Motions - G. Li
Sreqoenciesin linear behaviour are 10 Hz in t d o n and 15
Hz io lateral flexion, increasing fidm haeases tasim m Fl5 Nose Landing Gear Shimmy Taxi Test and
flqnency, which gas closer to lateral fquency: the Conelptive Adysea - D.GnwMolprm
system may bewme more uaatable, in opposition to this
studiedlanding gear. [SI The behDvioor Of N o I I - ~ ~Wbrptioa S m W. -
SzcmpIiasLesntpnca ' - K l m A ~ & d ~ h b b &
Jeff Baumann
Senior Engineer, Structural Dynamics and Loads
McDonnell Douglas Aerospace
P.O. Box 516
St. Louis, MO 63166-0516
USA
-
SUMMARY Fo = force the strut must transmit to the wheel at
This paper discusses a nonlinear shimmy model that point 0 for dynamic equilibrium, lb
includes components and features typically found in
the nose landing gear of fighter aircraft. The gear
FA = force of the strut on the effective rigid body at
structural model includes inertia, coulomb and viscous pointA, 16
damping, stiffness, and torsional freeplay. Single and FL = resultant force in the lower piston bearing,
dual non-corotating wheels are considered. The tires lb
are represented with the Moreland point-contact tire FMS = foreign military sales
model. FP = freeplay, rad (deg)
Qualitative analytical results from modeling the nose Fu = resultant force in the upper piston bearing,
, landing gear of the McDonnell Douglas FIA-18A are lb
- presented. This landing gear, in the normal US Navy
H = overall length of the strut, in
carrier-based configuration, is stable. When the
- launch bar and its power unit are removed for land- E = angular momentum, in lb sec
based versions of the aircraft, the nose gear may, [Z]bod,, = effective inertia tensor for the rotating lower
- - under certain adverse conditions, experience shimmy.
This condition can be aggravated by rotating forces strut, inlbsec 2
due to tire unbalance or other irregularities. A number = diametral moment of inertia of the wheel,
of corrective modification have been attempted with ZD
I varying degrees of success. These range from in lb sec2
relatively simple measures, such as increasing =polar moment of inertia of the wheel,
attention to tire maintenance, to more complicated
Zp
schemes, such as adding a friction collar, adding mass in lb sec2
dampers or even adding a dummy launch bar. For the JCOLuR = colladlaunch bar moment of inertia,
most severe cases, the dummy launch bar seems to be
required, however, the friction collar with mass in 16 sec2
dampers does show a significant improvement.
KO = tire lateral stiffness, lb / in
KSD = steering damper torsional stiffness, in lb
NOMENCLATURE KTL =torque link torsional stiffness, in lb
- Kw = tire radial stiffness, lb / in
A = acceleration, in / sec2
L = trailing arm length, in
C = tire coefficient of yaw, 1/ lb
Lp = length of the piston, in
Ci = tire yaw time constant, sec
ii?, = moment the strut must transmit to the wheel at
CD = tire lateral damping, lb sec/ in
point 0 for dynamic equilibrium, in 16
CPC = Cleveland Pneumatic Company
= exponential steering damper coefficient,
C ~ D @,, = moment of the strut on the effective rigid
body point A, in 16
in lb secMSD
MDA = McDonnell Douglas Aerospace
CTL = torque link structural damping, in lb sec
me = effective mass of the strut, lb sec2 / in
D = distance between piston bearings, in
DEAD = deadband function MR =rotating unbalance, lb sec2
Paper presented at the 81st Meeting of the AGARD SMP Panel on The Design, Qualification and Maintenance
of Vibration-Free Landing Gear, held in Ban& Canada from 4-5 October 1995, and published in R-800.
10-2
STEERING MOTOR /
Fig. 1a. Model Geometry for Cantilevered Landing Gear Fig. lb. Cantilevered Landing Gear
The second phase of Grossmans study involved information as the simulation ran. This made the
repeating the CPC analysis with more realistic values simulation modular in that new tire or structural
for the parameters. An attempt was also made to models could be added later.
linearize the nonlinear terms in the equations.
Since the details of this analytical model have been
The analysis correlated well with the shimmy speed - presented elsewhere (Baumann, Barker and Koval[7]
freeplay relations found in the taxi tests. Absolute and Baumann 181, a more general description will
stability was obtained for torsional freeplay less than follow. Table 1 lists the degrees of freedodstate
0.75 degrees. Also, the linearized models were shown variables used for this model.
to agree with the nonlinear model for the range of
parameters tested. 2.1. Kinematics
2.1.1. Coordinate systems
2. NONLINEAR LANDING GEAR MODEL The inertial coordinate system XoYoZo is attached to
the runway (see Fig. la). The Zo axis point upwards,
The objective of this work was to develop software
that would aid landing gear engineers in assessing the the Yo axis is aligned with the runway axis and points
shimmy stability of a more general class of landing aft, the X o axis completes the right-handed triad and
gear designs. With this analysis capability, potential points to the port side of the aircraft.
problems could be identified early in the design
process. It was also desired to be able to predict how A moving X I Y,Zl is fixed to the top of the aircraft
changes to existing design would affect the shimmy strut. It is rotated away from the XoYoZo about the
stability.
negative X o axis by the constant caster angle y .
Special care was taken to realistically model the
torsional components that are typically found in The X z Y 2 Z 2 system is fixed to the bottom of the strut
fighter aircraft nose landing gears. Three different and accounts for bending of the strut. It is rotated
landing gear geometries were studied: cantilevered, away from the XIY,Zl about the Y, axis by the
fully-levered and semilevered. Only the cantilevered
variable strut bending anglea .
type, as pictured in Fig. 1, will be discussed here.
The strategy employed was to analyze the landing gear The X3Y3Z3 system is fmed to the bottom of the strut
structure and the tire separately from each other. The and accounts for rotations of the strut. It is rotated
structural model and tire model were programmed as away from the X z Y 2 Z , about the Zz axis by the
separate computer subroutines that would share variable steering angle0 .
10-4
TABLE 1 - Degrees of Freedom The motion of this point will also be subject-to the
variable strut bending angle, a ,and steering angle,8 .
Strut Model
The center of mass of the effective rigid body, Point
x = lateral strut tip deflection, in
CG, of the effective rigid body which represents the
dx
x =- swiveling mass of the landing gear is located with
dt respect to point A by the arbitrary vector
a = strut tip bending angle, rad (deg)
da
a =-
dt Thus, motion of this point will also be subjected to the
8 = steering angle, rad (deg) three flexible degrees of freedom described above.
e =-
de
dt
8c = collar angle, rad (deg) 2.2. Tire Effects
Based on a complete kinematic description of the axle,
a separate tire model is used to determine the force
= damper angle, rad (deg) and moment that the strut must exert on the tire(s) to
OD
maintain dynamic equilibrium. This force may
include suspension, traction, rotating unbalance, and
Tire Model inertia effects. One such tire model is described in
p = position of rotating unbalance, rad (deg) detail in Section 3.
A = lateral tire deflection, in
For convenience, it is assumed that the force and
w = tire drift angle, rad (deg)
moment act at point 0, no matter where the tire is
mounted on the axle. Any transformations that are
away from the X 2 Y2Z, about the Z2 axis by the required are included in the tire model.
variable steering angle8 . The force and the moment will generally be of the
form
2.1.2. Rotational Kinematics
The angular velocity of the bottom of the strut is equal ig = FO(XO,E3)
to the angular velocity of the X3Y323 coordinate
system. This angular velocity accounts for the strut
Po = Po(Xo,,oL3)
bending rotation a plus changes in the steering angle which requires that previously defined expressions for
8 and is given by
To and E3 be substituted into the moment
z3= (asine >5+ (a case )j3+hi3. expressions.
TIRE MODEL
2.4.1. Translational Dynamics
From Newton's second law applied to the strut rigid I GROUND
body,
-
G = G + m , ACG Fig. 2. Model for Torques
where m, is the effective translational mass of the
strut. any torsional model. The modeling of the tire is
discussed in Section 3.
The lower shock strut can translate along the strut axis
2.4.2. Rotational Dynamics
with motion restricted by oleo-pneumatic forces
From Newton's second law for rigid body rotations within the strut. Coulomb friction will be present in
about the mass center of the strut rigid body, the upper and lower piston bearings. This friction will
oppose translation as well as rotation of the lower
shock strut. Gross rotations of the lower shock strut is
prevented with internal splines or with an external
torque linkage, as shown in Fig. lb.
The upper arm of the torque linkage is connected to
the steering collar. For carrier based aircraft, such as
Solving for ZA, the F/A-l8A described below, a launch bar would also
be attached to this steering collar and its effects would
have to be included in the torque model.
The steering collar is driven by a steering motor.
When not actively steering, this motor will serve as a
shimmy damper. This damper is modeled with a
spring, a Coulomb damper and an exponential
damper. The exponential damper will generate a
damping torque proportional to the relative angular
velocity across the damper raised to some exponent. If
2.4.3. Torque Models this exponent were 1.O, the damper would be a linear
Since shimmy involves rotations about the strut axis viscous damper. If this exponent were 2.0, the damper
and since most corrective measures involve torsional would be a "velocity-squared'' damper.
components (often torsional dampers), a good The connection between the steering collar and the
torsional model is very important. Figure 2 shows how damper will often have some freeplay. This freeplay,
the system of Fig. l b is modeled. caused by gear backlash and clearance between mating
parts, will allow some relative motion across the
Starting at the bottom of Fig. lb, the interaction
connection between the collar and the damper.
between the tire and the ground plays a large role in
10-6
A total of three degrees of freedom are required to set to zero, 8 is no longer considered to be a degree
describe the full torsional model described above. of freedom and one of the other torque models would
These degrees of freedom are the collar angle 8 , its be used.
derivative6 , and the damper angle 8 . The torque Applying Newtons second law of motion for rotations
model will contain the derivatives of these three to the collar yields
quantities as well as the torque TA which the above
components will apply to the lower strut rigid body.
Since some of these components may not be present or
may have only minimal effect, provisions are made to
remove one or more of the degrees of freedom from JCOLZAR
the model. Six separate cases involving various
Note that this equation requires that collar rotary
combinations of degrees of freedom were analyzed.
inertia be present ( j C O L L A#Ro ). If the rotary inertia
Each of these special cases are treated separately in
numerical computations. Only the first of these cases, is zero, 8 is not considered to be a degree of freedom
the full model as pictured in Figs. l b and 2, is and again one of the torque models would be used.
discussed here.
The final state variable derivative
It is first necessary to determine if the steering damper
is engaged to the collar. This depends on the absolute
difference between the collar angular position, 8 and
the damper angle, 8 D . is already defined since it is also a state variable.
Summing torques in the disengaged steering damper Because the steering damper is connected to the collar,
requires
and 0, = $ec)
T2 =ITrI-Tso.
If T2 is greater than zero, then the spring torque (q )
will overcome the Coulomb friction and need to be
reacted by the exponential damper. If T2 is less than
-cSDp C y sgn(6 C)] JCOLLAR
zero, the Coulomb friction alone will be sufficient to
react the spring torque. Summarizing:
Again notice that JCoLZAR must be nonzero for this
for T2 I O model to apply. This expression uses the deadband
function to account for freeplay. This function is given
by
Note that i. = li is the unit vector aligned with the The angular velocity of the axle is specified by the
strut model as
axle. Since liis specified in the strut model, it is
possible to solve for the angles g and h :
-S
O R =~ o x i o + o y j O +o,kO.
where
FmD = W - K w ( R - R s ) = tire radial force,
where h is the time rate of change of the wheel spin
rate. FN = K D A + CDA = Moreland tire side force
and
3.2.2. Vertical Constraint
FD = p (FmD cosh) = tire drag force.
This constraint sets the geometric radius, R ,of the
tire by requiring that the tire always maintains contact The final element of Moreland's model relates yaw to
with the ground at the contact point. The appropriate side force:
equation is
\ir = - -1( c F ~+w).
Cl
Differentiating,
3.3.2. Tire Unbalance
Tire unbalance forces are due to centrifugal effects of
These two equations can be solved for R and R a rotating unbalance. If the product of the unbalance
mass and radius is MR and the wheel spins at a
rate , then the unbalance force is
3.2.3. No Side Slip -
F~~~~~= M R R (-sin
~ pj2 + cos pi2)
A tire under the action of a side force will drift or yaw
(as it rolls) from the direction in which it is pointed where p is the rotation of the unbalance mass and is
(see Fig. 3). This angle is called the side-slip angle or
equal to ( a t + p ; p is the initial wheel position.
the drift angle and is denoted by y . A standard
condition is to require that as the wheel rolls and side
slips, the component of the velocity normal to the path
of travel becomes zero, i.e.,
10-9
UTAW
FIXED ADAPTER
Fig. 7.Steering Manifold
when shimmy is reprted. As this reduces the tire Life, While the wheels m.not balanced, the tires are
this makes the aircraft more costly to operate. marked by the manufactnm with a balance mark that
mbe aligned with the valve stem on the wheel to
Two specific fmrshave been observed to be minimize static unbalance. Additionally, the tires
important. These factors lead to cyclic forcing of the should be inflated with dry air or nitrogen (to
nose gear and. while this M not shimmy, it can drive minimizeinternal moisbm) and be allowed to grow
similar vibrations. Ibe first effect is rotating for atleast 12 homs before loading. This ensures an
unbalance. T%e wheels and tires are checked to insure even stretch of the rubber. An effort should also be
specification compliance but the combined assembly made to match tire diameters and pressures for the two -
is not. Apmduction tire-wheel assembly was tested by nose gear tires. The tire predsure should be regularly
MDA on a commercial automotive wheel balancing checked.
machine and the unbalance was measured at 12 inch-
ounces.While this is a severe unbalance, it is roughly To prevent formation of flat spots, the aiman should
the sum of the two unbalances allowed by the be moved evay 48 hours 01 jacked if unused for long
specifications.Since the tires wear qnickly with use, it periods of time.
would seem impractical to balance the assembly on a
routine hapis.
4.5.2. Addition of a Friction Collar to Replace the
The second factor is cold set of the rubbar redulting in
Launch Bar Power Unit Manfold
flat spots. This can reanlt from loading the tire before
the rubber is fully stretched by the inflation pres sur^ The seals of the launch bar power unit manifold
or by leaving the aircraft parked on the tires for represent about 20%of the breakout torque and about
s e v d days. This was observed during a winter flight 10% of the running torquerequired to swivel the
test in St. Louis.The tires were inadvertently allowed wheels. This damping (along with some stabilizing
to cold set and shimmy resulted. mass)can be regained by adding a collar where the
manifold is usually situated.
45. Corrections
4.5.3. Addition of Simple Masses in Conjmnion with
Several modifications to the FMS landing gear
the CollarAssnnbly
contlluration have been pmposed and/or attempted to
correct shimmy. Ibeee will be discussed in d e r of Additional weights can be added to the above collar,
increasing cost and complexity. as shown in Pig. 8, to replace some or all of the mass
removed with the launch bar. Since this mass has 110
4.5.1. Increased Attention to Tire Maintenance mechanical function it can be placed where it will not
As earliu discussed. the condition of the tires has a interfexe with gear stowage. Due to the malluradiw
big effect on shimmy and increased e m p h i s on of gyration, this mass adds up to about 30% of the
maintenance can do a lot to mitigate shimmy. original rotational inertia The weights can be
10-11
O.W)OE+OO
. *
- -o.5ooE+oo
-
0
REFERENCES
[l] vonschlippe, B., and Dietrich, R., Shimmying of
a Pneumatic Wheel, Papers on Shimmy and
Rolling Behavior of Landing Gears presented at
the Stuttgart Conference, October 16-17, 1941.
NACA TM 1365, August 1954.
[2] Moreland, W. J., The Story of Shimmy, Journal
of the Aeronautical Sciences, December 1954,
pp. 793-808.
[3] Collins, R. L., Theories on the Mechanics of Tires
and their Applications to Shimmy Analysis,
AIAA Journal of Aircraft, Vol. 8, No. 4,
April 1971, pp. 271-277.
[4] Black, R. J., Application of Tire Dynamics to
Aircraft Landing Gear Design Analysis, NASA
CP-2264, Workshop on Tire Modeling, at
Langley Research Center, September 7 - 9, 1982,
pp. 71-93.
REPORT DOCUMENTATION PAGE
1. Recipients Reference 2. Originators Reference 3. Further Reference 4. Security Classification
of Document
AGARD-R-800 ISBN 92-836-1032-6 UNCLASSIFIED/
UNLIMITED
I I ____ ____ I~_______
8. Author(s)/Editor(s) 9. Date
Multiple March 1996
12. Distribution Statement There are no restrictions on the distribution of this document.
Information about the availability of this and other AGARD
unclassified publications is given on the back cover.
13. KeywordsKlescriptors
Landing gear Stability
Aircraft landing Landing simulation
Aircraft Design
Vibration Maintenance
Damage Vibration analysis
14. Abstract
The Structures and Materials Panel of AGARD sponsored a Workshop focusing on the various
vibrational and stability problems (e.g. shimmy, antiskid induced vibrations) that must be
considered in the early design phase of landing gear systems, especially problems which are
related to vibrations of the combined structural system formed by the landing gear, its tyres and
the flexible aircraft structure. The intention was to indicate the impact of (combined) landing
gear/aircraft vibration problems on aircraft design and to discuss the state-of-the-art technology
in this area and to define possible future steps of development.
.
a@-
NATO @ OTAN
7 RUE ANCELLE 92200 NEUILLY-SUR-SEINE DIFFUSION DES PUBLICATIONS
FRANCE AGARD NON CLASSIFIEEX
TBlkopie (1)47.38.57.99 TBlex 810 176
Aucun staclr de publications n'a existh AGARD. A partir de 1993. AGARD dtiendraun stack limit6 des publications a w o c i h aux cycles
de conf&ences et corn sgciaux ainsi que les AGARDographies et les rapports drs groupes de travad. o r d s et publies & panir de 1993
inclus. Les demandes de mseignements doivent etre adresseeS i AGARD par lettre ou par fax i I'adresse indiqu& ci-dessus. VeuiUez ne
pes tdiphonor. La diffusion initiale de toutes les publications de I'AGARD est effectuk aup& des pays membrcs de I'OTAN par
I'intWnediaire des centres de distribution nationaux indiqnh cidessous. Des exemplaires suppl~menmirespeuvent pmfois @e o b h u s
aup& de ces centra (& I'exception des Etats-UNs). Si vous souhaitez recevoir toutes les publications de I'AGARD. ou simplement c e l b
qui Concement cenaiaS Panels, vous pouvez demander & &e inch sur Ln liste d'envoi de I'un de ces centres. Les publications de I'AGARD
son1 en vente au@s des agenw indiqUees ci-dessous. sow forme de pholocopie ou de microfiche.
CENTRES DE DIFFUSION NATIONAUX
ALLEMAGNE ISLANDE
Fachinformationszentrum Karlsruhe Director of Aviation
D-76344 Eeeenstein-Leowldshafen 2 cl0 Ruerad
BELoIQUE Reykjacik
Coordonnateur AGARD-VSL ITALIE
Etat-major de la Force atnenne Aeronautica Militare
@artier Reine Elisabeth Uffcio del Delegato Nazionale all'AGARD
CANADA
Rue d'Evere, 1140 Bmelles
LUXEMBOURG
Pomezia (Roma)
AGENCES DE VENTE
NASA Center for ESMnformation Retrieval Service The British Library
Aerospace Information (CASI) European Space Agency Document Supply Division
800 Elkridge Landing Road 10. rue Mario Nibs Boston Spa. Wetherhy
Lintbicum Heights. MD 21090-2934 75015 Paris West Yorkshin U 2 3 7BQ
Etats-Unis France Royaume-Uni
Les demandes de microfiches ou de photocopies de documents AGARD (y compris les demandes faites aup& du C A S ) doivent
comporter la denomination AGARD, ainsi que le numero de s&ie d'AGARD (par exemple AGARD-AG-315). Des informations
analogues, telles que le titre et la date de publication son1 souhaitables. Veuiller noter qu'il y a lieu de spkiiier AGARD-R-nnn et
AGARD-AR-nnn Ion de la commande des rapports AGARD et des rapports consultatifs AGARD respectivement. Des dftnnces
bihliographiques complhtes ainsi que des r&sumh des publications AGARD figurent dans les journaux suivants:
Scientific and Technical Aerospace Reports (STAR) Government Reports Announcements and Index (GRA&I)
publid par la NASA ScienWic and Technical publi&par le National Technical Information Service
Information Division Springfield
NASA Headquarters (J'IT) Virginia 22161
Washington D.C. 20546 Etats-Unis
Etats-Unis (accessible kgalement en mode interactif dam la base de
donnks bihliographiques en ligne du NTIS, et sur CD-ROM)
WCRWAY
NCWCplkPaWR
Apn: Biblioteket
P.O.Box 25
N-ZOOI Kjeller
PORTUGAL
EstndoMniordaFoq
-
SDFA Ccntro de Do
WAIN
INTA (AGARD Rtbli lions)
Pintor Roslles 34
28008 Madrid
TURKEY
MillS Ssvunma Ea@
ARGE Dairesi Bqkat
06650 Bakanhkh-An
UNlTED KINGDOM
Defence RLKaxh Inf<
Kcntigm House
..Bmvm
6s ~
Smct
. ~.~
Glasgow 0 2 8Ex
UNlTED STATES
NASA Hedquatbm
Coee JOB-1
Wadngton, D.C. 205
The UaHsd States Natlo0.l DLstribntiol Cmtre dom NOT hold rtodrs Ot ,GARI) pObtif.nolla
.
Applications for topics should kc made d i m t to the NASA Coara for AemSpnee lnfm tim (CAS0at the address below.
Change of address q u e s t s should also go to CASI.
Y'
SALES AGENCIES
NASA Canter for -
s
B . RetrievalseNice ThC British Libmy
Documcatsuppl Centre ,
Bmton Spa,We&y
West Yorlcrhin LS23 784
unimi states ~ FmEe United Kingdom
should indude the word 'ACARD'
rb a8 title and publication d.te in
ARD-R-nnn and AGARD-AR-mn,