72_olm.pdf

IN FIELD MEASUREMENT OF FORCES AND DEFORMATIONS AT
THE REAR PART OF A MOTORCYCLE AND STRUCTURAL
OPTIMIZATION
G. Olmi, A. Freddi and D. Croccolo
DIEM Department, Engineering Faculty
University of Bologna
Viale del Risorgimento, 2, 40136 Bologna, Italy
ABSTRACT
The present work is about the analysis of the asymmetric displacement behaviour of the rear part of a motorcycle. Stylistic
reasons led to the design of a vehicle with only one suspension located at the left side. Experimental tests performed on a
circuit with seven obstacles along a straight line confirmed that bending displacement is higher on the right side than on the left
one. This work aims to perform a structural optimization of the components involved in the rear side of the motorcycle, in order
to find a solution to the problem. A hybrid approach is applied: the force acting on the suspension and bending displacements
at the rear side were simultaneously measured in working conditions, then a FEM model was set up, validated and applied for
design optimization purposes. Both methodological aspects and applicative results are presented and discussed. Finally a
solution in accordance with design specifics is proposed.
Introduction
Nowadays in the sector of sport motorcycles, designers aim to the development of products with better characteristics than
competitors’ not only from the mechanical point of view, but which can be attractive also for other features. Stylistic design is
becoming more and more important. Aiming to more original designs, motorbike models with only one suspension were
developed. In this case mechanical connection between the frame and the fork is performed on just one side (usually the left
one), which implies an asymmetric design for the fork itself and for the whole rear part of the motorbike. Experimental tests
performed on a prototype led to the result that, due to the asymmetric layout at the rear side, also the structural response was
different on two sides. Displacement sensors (resistance potentiometers) measured the relative movement between the frame
and the fork on both sides. In particular work conditions it resulted some millimetres higher on the right side than on the left
one (where the suspension is applied). This unacceptable behaviour was observed during an in field experimental campaign:
the motorbike was driven at constant speed (50 km/h) on a straight line with seven obstacles: asymmetric behaviour was quite
evident when overcoming each obstacle. This behaviour is graphically explained in Fig. 1: the difference between relative
approaching on left and right sides, sum of the two contributions ∆t (related to the frame) and ∆f (related to the fork) is later
indicated by ∆.
The aim of this work is to perform a structural optimization of the motorbike rear side, in order to reduce asymmetry effects. A
hybrid approach was applied: firstly the force acting on the suspension was measured during in field tests. Such load is
transmitted between the frame and the fork by the suspension. The force R (Fig. 1), is the result of the superposition of several
effects due to suspension mechanical response. Under a compressive load the reaction is to the sum of elastic, oleo dynamic
and pneumatic contributions and is not easy to estimate by a theoretical or numerical approach, especially in real conditions.
On the basis of references [1-4] it can be argued that performing in field tests, mainly with the application of strain gage
sensors on the most loaded locations is often the only way to obtain a rigorous measurement of forces. For example Petrone
et al. [5-9] used strain gage devices to measure load spectra on the wheels of a scooter or on main components of a
motorbike. Strain gages were directly applied on components or on particularly shaped transducers. In all the papers the
numerical results always need an experimental validation, because it is very difficult to simulate real working conditions with
numerical analysis [10]. Experimental results were the first step for a further analysis conducted by finite element modelling of
the fork and the frame. Such models, validated by comparing simulated deformations at the rear part to experimental results
concerning displacements, were then used for structural optimization purposes.
Left side
Right side
Left side
Right side
FRAME
∆t
R
R
∆f
FORK
Figure 1. Scheme of the motorbike rear side
Figure 2. Suspension at the left side of the motorbike
Experimental tests and results
For the force measurement strain gages rosettes were applied on a location where strain intensity was sufficiently high and
acceptable with reference to the gage sensitivity.
In order to measure the force acting on the suspension, two-grid (0°/90°) strain gage rosettes were used. They were applied
closely to the hinge hole. A little extension of the threaded bar near the hole was machined so that a cylindrical shape was
made available for rosettes application and strain were adjusted to increase measurement sensitivity. In Fig. 3 a detail of
rosettes collocation and of Wheatstone bridge connection are shown.
Figure 3. Application of (0°/90°) strain gage rosettes on the suspension
The first step regarded the calibration of the sensing system. It was required to determine an equation for the conversion of
Wheatstone bridge output into the measurement of the force acting on the suspension. A testing machine was used for such
operation. The machine (Fig. 3) is a two-column type, with the loading cell at its higher side. On the lower side the lower grip is
connected to the cam device transmitting up & down motion. The suspension is subjected to a compression cyclic load, which
is applied in displacement control conditions (closed loop control on displacements). In order to determine the force-bridge
output relationship on a sufficiently wide load range, several tests were performed for maximum strokes varying between 25
and 80 mm. A null value of the stroke is associated to a null value of force acting on the suspension. During the tests the force
acting on the suspension was measured by a strain gage loading cell, meanwhile the strain signal was acquired by a Leane
MCDR-M-128 device. The load cycle was repeated for ten times. Finally peak values of force distribution, corresponding to
maximum stroke condition, were related to the highest values of strain output (an example is shown in Fig. 5, above, right
2
hand). The result is shown in Fig. 5 (below, left hand), emphasizing the high linearity (R = 0.99) between force and strain.
Both force and strain are considered with negative values due to compression.
Strain output
[microstrain]
Lettura [microepsilo
Adjustable transverse element
Loading cell
Higher grip
0
-200
-400
-600
-800
-1000
-1200
0
10
20
30
40
50
-200
0
Time [s]
Force on the
suspension
[N]
Forza
[N
Suspension
Lower grip
Cam device
0
-1200
-1000
-1000
-800
-600
-400
-2000
-3000
-4000
-5000
Strain output [microstrain]
Figure 5. Calibration tests outputs
Figure 4. Testing machine used during calibration tests
1000
N.L. (spring pre-load)
0
Force
Forza [N]
[N
0
20
40
60
80
-1000
2
R = 0.951
-2000
-3000
-4000
N.L. rubber disk
-5000
Stroke [mm]
Figure 6. Force – Stroke diagram
In the calibration phase the diagram relating the suspension mechanical response vs. its stroke (Fig. 6) was also determined.
At low value non linearity is due to spring pre-load, while at high values it is due to the compression of a rubber disk when the
2
theoretical maximum value of the stroke (80 mm) is overcome. In the central zone the linearity is very high (R = 0.95).
Figure 7. The application of potentiometers on both motorcycle sides
The calibration was essential for further experimental phase aimed to the measurement of the force on the suspension in
working conditions. The suspension, equipped by strain gages was mounted on a motorbike, then potentiometers were
applied both on the left and on the right sides. They measured relative movement between the frame and the fork, to which
they were connected by hinges, as shown in Fig. 7. Also the output signals from potentiometers were acquired by the same
equipment for simultaneous measurement of force and displacements at left and right rear sides.
1
2
3
4
5
6
7
4000
6
3000
4
2000
2
1000
0
-1000362
363
364
365
366
0
368
-2
367
-2000
-4
-3000
-4000
-6
Force [N]
Diff.
∆ [mm]
Time [s]
Force on the Forza
susupension
[N]
[N
Force on the
susupension
[N]
Forza
[N
Figure 8. Circuit for in field tests with seven obstacles
Figure 9 a. Results of in field tests (40 km/h)
1000
500
-1
R2 = 0.9365
0
-500 0
1
2
3
4
5
6
-1000
-1500
-2000
-2500
-3000
-3500
[mm]
∆D[mm]
Figure 9 b. Stiffness diagram for in field tests (40 km/h)
4000
6
3000
4
2000
2
1000
0
-1000445
446
447
448
449
0
451
-2
450
-2000
-4
-3000
-4000
-6
Force [N]
Diff.
∆ [mm]
Time [s]
[N
Force on theForza
susupension
[N]
Force on the
susupension
[N]
Forza
[N
The in field tests were performed on a circuit with seven obstacles on a straight line, as shown in Fig. 8. Four obstacles out of
seven (referenced with numbers 4, 5, 6, 7) were larger sized than the others. The obstacles were overcome at three different
speeds (40, 50, 60 km/h), with three replications: passages through each obstacle were totally nine. The acquiring system is
particularly designed for in field tests. The forces and the differential displacement ∆ were related and compared by plotting
diagrams as shown in Figg. 9 a ÷ 11 a. See also [11].
6
4000
3000
4
2000
2
1000
0
-1000631
631.5
632
632.5
633
633.5
634
634.5
0
635
-2
-2000
-4
-3000
-4000
-6
Time [s]
Figure 11 a. Results of in field tests (60 km/h)
500
-1
R2 = 0.9346
0
-500 0
1
2
3
4
5
6
-1000
-1500
-2000
-2500
-3000
-3500
∆D[mm]
[mm]
Figure 10 b. Stiffness diagram for in field tests (50 km/h)
Force [N]
Diff.
∆ [mm]
Force on theForza
susupension
[N]
[N
Force on the
susupension
[N]
Forza
[N
Figure 10 a. Results of in field tests (50 km/h)
1000
500
-2
R2 = 0.9202
-1
0
-500 0
1
2
3
4
5
-1000
-1500
-2000
-2500
-3000
-3500
-4000
[mm]
∆D[mm]
Figure 11 b. Stiffness diagram for in field tests (60 km/h)
6
Discussion
The force in Figg. 9 a ÷ 11 a is indicated with negative values, due to suspension compression: at the first impact with the
obstacle, force module rapidly increases up to a maximum. At the same time the differential displacement ∆ increases, and the
asymmetric behaviour of the rear side is particularly evident: the frame and the fork tend to approach more at the right side
than at the left one (Fig. 12 a). Then, after a drop-kick, the load on the suspension suddenly decreases. Positive values in the
force diagram indicate that the suspension becomes slightly in tension. When the force sign changes, ∆ becomes negative,
meaning that the distance between the frame and the fork is higher at the right side than at the left one (Fig. 12 b). Then the
force trend is inverted and the suspension is again subjected to compression, with maximum positive value for ∆. The
described trend is repeated for seven times, as seven are the obstacles. Plotting results together shows that ∆ distribution
appears to be qualitatively mirrored with respect to force’s. When load module has a maximum value, also ∆ becomes
maximum, with the same trend also for minimum values. The sign changes of ∆ and of the force take place at the same times.
Such observations suggest that the relationship between the force acting on the suspension and ∆ is linear: to verify this
hypothesis, stiffness diagrams were plotted in Figg. 9 b ÷ 11 b with reference to results shown on the left side. For sake of
clearness only data points were plotted: in all the cases response can be well interpolated by linear regression lines with high
2
values of R (greater than 0.9). Some non linearity effects can be observed when force and ∆ assume high values with
impulsive increase and following decrease. Such effects appear to be more evident at higher speeds, at 50 and 60 km/h: when
the suspension stroke is very high, the rubber disk is compressed and the structural response becomes non linear as shown in
Fig. 6.
On circuit straight line with no obstacles, the measured force on the suspension is about 1,700 N with no strong variation. This
value is due to the pilot and motorbike weights transmitted to the rear part of the vehicle. Experimental results emphasize that
such value is highly increased by dynamic effects due to the impacts against obstacles. The highest value for the compressive
force is 3,600 N, i.e. more than 100% incremented with respect to load in static condition. As previously observed, as the
motorcycle overcomes obstacles, the force has two peak values, with in-between decrease. Such values are summarized in
histograms below (Figg. 13 ÷ 15) for tests at different speeds. Values of force peaks [N] were processed in an ANOVA test [12]
(Table 1): two factors were considered: the mean speed during the test and the type of the obstacle indicated by progressive
number (see Fig. 8). Two close peaks were interpreted as two replications for each treatment combination. The result of the
statistical analysis was that no significant differences exist among force peaks at the 5% significance level. It implies that force
peaks are not influenced by the speed, ranging from 40 and 60 km/h, and by the obstacle size (obstacles 4, 5, 6, 7 were larger
sized). Consequently the grand mean value for the increased force due to dynamic effects could be calculated: such value is
2,850 N, corresponding to a 60% increase.
Left side
Right side
Left side
Right side
FRAME
3500
FRAME
∆>0
Force peak [N]
3000
∆<0
a)
b)
FORK
2500
2000
1500
1000
FORK
500
0
1
3500
Force peak [N]
Force peak [N]
4000
3000
1500
1000
500
0
1
2
3
4
5
6
7
Obstacle
Figure 14. Force peaks for in field test at 50 km/h
4
5
6
7
Figure 13. Force peaks for in field test at 40 km/h
3500
2000
3
Obstacle
Obstacle
Figure 12. Deformation at the rear side (cases for ∆ > 0 and ∆ < 0)
2500
2
3000
2500
2000
1500
1000
500
0
1
2
3
4
5
6
Obstacle
Figure 15. Force peaks for in field test at 60 km/h
7
SPEED
40
km/h
50
km/h
60
km/h
1
2831
2895
2501
2914
2377
2792
Table 1. Factorial plan and two-factor ANOVA statistical analysis (Force peaks in N)
OBSTACLE NUMBER
2
3
4
5
6
7
2600 2903 3042 3118 3058 2851
SSQ
DOF
MSQ
2819 3011 2792 2823 2764 2768 Speed
64,910
2
32,455
2485 2585 2660 2788 2792 3106 Obstacle
865,371
6 144,229
2823 2995 3082 2887 2935 2684 Interaction
671,883
12
55,990
2214 2389 3365 3194 3019 3636 Error
1,296,757
21
61,750
2716 2847 3238 3027 2923 2764 Total
2,898,920
41
Fcalc.
0.526
2.336
0.907
p-v.
59.9
6.94
55.6
Frame and fork FEM modelling with the use of shell elements
During experimental tests the force on the suspension and relative displacements between the fork and the frame were
simultaneously measured. FEM structural analysis both of the frame and of the fork was then performed considering a mean
value of the measured force peaks. Frame and fork were modelled by the software I-DEAS (release 11) in the elastic field, due
to the linear relationship between force and displacement (Figg. 9 b ÷ 11 b). 3D CAD models with shell elements [13] were
adopted. The algorithm assumes a linear variation for stress along the thickness.
FEM models were validated by experimental results, concerning the elements choice, models meshing, constraining and
loading systems in order to use them for structural optimization.
In the frame model the engine, usually supposed to have infinite stiffness, was modelled by the application of beam elements,
having square sections with finite stiffness value. The cavities of holes for connections to the engine were modelled by radial
elements with infinite stiffness: central points were connected to the engine by “couple degree of freedom” constraints [13] .
Only the rotation along the horizontal axis is left free. Hinges for connections between the fork and the anterior wheel were
meshed by beam elements with the real circular section. The ∆t (Fig. 1) was evaluated for the nodes located where
potentiometers had been mounted. Suspension reactive force was applied to the node at the centre of radial elements with
infinite stiffness located at the hinge hole.
In Fig. 16 a sketch of the structure is shown: “beam elements” are indicated by lines with related transverse sections. Arrows
represent constraints preventing a translational motion, and double arrows are located on the hinge linking the frame to the
fork. In this case just one rotation, around the hinge axis, is permitted, while constraints prevent rotating motions around other
perpendicular axes.
A negative displacement for the right side and a simultaneous positive one for the left one can be observed. Such a behaviour
is quite reasonable because the reactive force, transmitted by the suspension, is applied at the left side. By supposing the
mean value of 2,850 N for the force, the difference between two sides displacements ∆t assumed the value of 0.93 mm.
Sspension
connection
Bearings for
connection to
the frame
Backward
wheel hinge
Suspension
connection
Hinge for
connection
to the fork
"beam"
elements
Forward
wheel hinge
Figure 16. Constraining and loading system of the frame.
"Couple degree of
freedom" element
Figure 17. Constraining and loading system of the fork
Triangular parabolic shell elements were used for fork modelling, while elements named “constraints” [13] were applied for the
welding connections between rods. The hinge for connection to the wheel was meshed by “beam” elements. An annular
section reproducing the real one was considered. Radial elements with infinite stiffness were created at the slotted hole
internal cavity. The central node was then connected to the hinge terminal nodes by “couple degree of freedom” constraints.
The only rotation around backward wheel axis was left. Suspension reaction was applied at the centre of radial elements,
considering the mean peak value of the force.
Fig. 17 shows a loading and constraints systems . The “couple degree of freedom” constraint is indicated by a blue circle at the
centre of radial elements. Only vertical displacements are considered for data processing: for force value of 2,850 N the
difference in displacements at the two fork sides (term ∆f, Fig. 1) was estimated as 3.05 mm.
On the basis of results in Figg. 9 b ÷ 11 b, a linear model was adopted in structural analysis. Estimated displacements under
the action of 2850 N force were 3.05 mm for the fork (∆f) and 0.93 for the frame (∆t), which implies a total value ∆ of 3.98 mm.
By assuming a linear relationship between the force and differential displacement ∆, each of the measured force values was
processed for the calculation of the corresponding value of ∆. Values of ∆, simultaneously measured with forces and
numerically simulated, were compared in diagrams for all the performed tests. Some of the results, related to tests at 40, 50
and 60 km/h, are shown on the left in Figg. 18 ÷ 20. “Truth Diagrams” (Figg. 18 ÷ 20 on the right) show a direct comparison
between experimental and numerical results. Mean absolute error is about 0.3 mm, while mean relative errors range between
12 and 16%.
6
6
5
5
∆ [mm]
4
4
∆ [mm] (simulated)
3
Experimental
3
Simulated
2
2
1
1
0
362
-1
363
364
365
366
367
368
0
-1
Time [s]
0
1
-1
2
3
4
5
6
(experimental)
D [mm]
(experimental)
∆ [mm]
Figure 18. Validation of the FEM model for the test at 40 km/h
6
6
5
5
∆ [mm]
∆ [mm] (simulated)
4
4
3
Experimental
3
Simulated
2
2
1
1
0
445
-1
446
447
448
449
450
451
0
-1
Time [s]
0
1
2
3
4
5
6
5
6
-1
D [mm] (experimental)
∆ [mm]
(experimental)
Figure 19. Validation of the FEM model for the test at 50 km/h
6
6
5
5
∆ [mm]
∆ [mm] (simulated)
4
4
3
3
Experimental
2
Simulated
2
1
0
631
-1
631.5
632
632.5
633
633.5
634
634.5
1
635
0
-2
-1
Time [s]
0
-1
Figure 20. Validation of the FEM model for the test at 60 km/h
1
2
3
4
D [mm](experimental)
(experimental)
∆ [mm]
Structural optimization
After validation, the FEM model was applied even to modified versions of the fork and the frame for several simulation tasks.
The aim was to determine the best structural solution with low values for the differential displacement ∆. A key issue in design
improvement was not to modify the aesthetics and the stylistic line of the motorbike. Consequently strong modifications
involving the structural layouts of the frame and the fork could not be acceptable, and the main feature of the motorbike, i.e.
the adoption of just one suspension at the left side had to be maintained.
By comparing the values of ∆t and ∆f, it can be guessed that the greatest contribution is due to fork deformations: term ∆
depends on fork displacements for the 80% of its total value. Such a result, with the term ∆f strongly higher than ∆t, suggested
that any modification of the frame would imply just a little improvement in displacement reduction with the risk of altering the
stylistic line. Consequently further simulating tasks involved the fork, with the aim of optimizing its design and to reduce ∆f.
Beams of the fork structure are tubes with annular section (external diameter: 60 mm, internal diameter: 56mm). According to
specifics, proposed modifications did not involve the external shape or dimensions, but only the thickness. Modified versions of
the fork, with tubes having different internal diameters were simulated with two goals: firstly to investigate the sensitivity of ∆f to
thickness variations, then to optimize the value of thickness to reduce ∆, without strongly modifying the original project.
Simulations were conducted, by supposing to modify the thickness of just one tube: the tube at the left side, or the tube at the
right side or the central one (Fig. 21). The diagram in Fig. 22 shows the percentage variation of ∆: sensitivity is doubled when
the variation involves the central beam. This result suggests the best design strategy to reduce displacements ∆f and ∆.
Values of ∆ are plotted vs. thickness: it can be observed that the relationship is non linear: the improvement on displacement
reduction becomes smaller as the thickness is increased. Such a trend is well interpolated by a quadratic curve which has a
plateau for a thickness value of 4.5 mm: for this value the reduction on ∆ is a bit greater than 20%. On the basis of this result, it
seems to be not convenient to increase the thickness any more, because no further appreciable reductions could be achieved.
The only way to obtain a stronger reduction would be to change the whole project, by increasing the thickness not only of the
central tube, but also of the lateral ones. However this solution would not be acceptable because it would imply a strong
increase in the fork weight, which would be doubled. The performed numerical analysis shows that acceptable improvements
can be obtained by adopting a quite easy modification of the original design. Reduction of ∆ to very low values would require
more strong modifications, which would alter the stylistic solution and imply a too high increase in weight.
Central tube
30
Right side tube
25
20
2
∆ reduction [%]
R = 0.9979
Central tube
Left or right tube
15
10
Interpolation
5
Left side tube
0
2
Figure 21. FEM simulations on modified versions of the fork
2.5
3
3.5
4
4.5
Thickness [mm]
Figure 22. Sensitivity analysis and design optimization
Conclusions
Some in field tests performed on a motorbike prototype with only one suspension on just left side had shown an asymmetric
mechanical response of the vehicle rear part. Relative bending displacements between the frame and the fork were quite
different at the left and the right side. The object of this research was to develop a modified version of original design, with a
reduction of the effect and by maintaining the original stylistic line.
A hybrid approach was applied aiming to structural optimization of the rear part, by integration of experimental testing and
numerical FEM analysis. Main steps of the work and related results can be so summarized:
•
•
•
•
Two (0°/90°) strain gage rosettes were applied on the suspension for measuring force in working conditions.
The whole device for force measurement was calibrated in lab. The linear relationship between the force on the
2
2
suspension and Wheatstone bridge output reading was determined with a high value of R (R = 0.99).
In field tests were performed on a circuit with seven obstacles at different speeds ranging from 40 to 60 km/h. Both
force on the suspension (by strain gage rosettes) and differential displacement ∆ (by potentiometers) were
simultaneously measured.
Experimental results show a good repeatability: force is about 2 kN in normal driving conditions and is strongly
increased when overcoming obstacles. An ANOVA test confirmed that differences among force peaks are non
•
•
•
•
•
significant for speed variation and obstacle size: mean value is about 2850 N (60% more with respect to static
conditions).
A further analysis was conducted to investigate the linearity of the relationship between measured force and ∆.
2
Linearity hypothesis was acceptable, as confirmed by high values of R of interpolating lines (greater than 0.9%)
FEM models with shell meshing elements were developed to investigate displacements ∆t and ∆f due to the frame
and the fork deformations.
Numerical models were validated by comparing simulated displacements to experimental measured ones. A good
agreement was observed, with mean relative error of about 13 %.
One of main results was that ∆ (=∆t+∆f) strongly (80%) depends on fork deformations, which led to the determination
of the component which mostly needed modifications. Consequently further analysis was focused to optimize its
design.
Validated fork model was finally simulated to investigate the sensitivity of the differential displacement to tubes
thickness variations. Different design versions were considered and finally a solution was proposed with a reduction
of ∆ up to 20%. An optimization was performed on the thickness variation to be applied to the original design. This
solution does not require a radical modification of the original project, does not imply strong weight variation and
maintains the same stylistic line.
References
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
Ladda V., Zaschel J.M., “Langzeituntersuchungen mit Dehnungsmeßstreifen zur Ermittlung der Fahrwerkskräfte am
Airbus A 300”, MTB, 15, 3 (1979).
Tu J.F., “Strain field analysis and sensor design for monitoring machine tool spindle bearing force”, Int. J. Math. Tools
Manufact., 36, 2, 203-216 (1996).
Freddi A., Olmi G., Croccolo D., “Ein neuer Aufnehmer zur direkten Messung der von einem Skifahrer auf die Skier
übertragenen Kräfte und Momente”, MTB, 1, 3-10 (2002).
Olmi G., Freddi A., “Experimental stress analysis of an asphalt milling machine in the assembly and working phases”, In
Proceedings of XXIII Danubia-Adria Symposium, Podbanské, Slovakia (2006).
Forni A., Meneghetti G., Petrone N., “Structural behaviour of a sport motorcycle in typical road events”, In Proceedings of
“Power two wheels”, ATA International Conference, Pisa, Italy (1998).
Petrone N., Biliato N., Toselli G., Fazi F., “Acquisizione ed analisi dei carichi su scooter 125 in prove su pista e su banco a
rulli”, In Proceedings of the XXIX AIAS National Conference, Lucca, Italy (2000).
Petrone N., Saraceni M., Cecchetto M., Pegoraro D., “Sviluppo di un banco prova per l’analisi strutturale e a fatica di telai
motociclistici”, In Proceedings of the XXX AIAS National Conference, Alghero, Italy (2001).
Petrone N., Saraceni M., Prugnola F., “Sviluppo di un simulatore di strada per prove complesse su forcelle
motociclistiche”, In Proceedings of the XXXIV AIAS National Conference, Milan, Italy (2005).
Saraceni M., Petrone N., “Sviluppo di un sistema di celle dinamometriche per la misura completa dei carichi agli assi ruota
di uno scooter 150 cc.”, In Proceedings of the XXXII AIAS National Conference, Salerno, Italy (2003).
Cossalter V., Cinematica e dinamica della motocicletta, Ed. Progetto, Padova, Italy (1999).
Olmi G., Freddi A., Croccolo D., Ferrazzo G., “Un approccio ibrido per la valutazione delle inflessioni in esercizio del
retrotreno di un motociclo”, In Proceedings of the XXXV AIAS National Conference, Ancona, Italy (2006).
Berger P.D., Maurer R.E., Experimental Design, Duxbury-Thomson Learning, Boston, USA (2002).
I-DEAS V11 On Line Tutorial