Title: Two Special Finite Elements for Modelling Modelling of Rolling Contact in a Multibody Environment
1Delft University of Technology Design Engineering
and Production Mechanical Engineering
Modelling of Rolling Contact in a Multibody
Environment
Arend L. Schwab Laboratory for Engineering
Mechanics Delft University of Technology The
Netherlands
Workshop on Multibody System Dynamics, University
of Illinois at Chicago , May 12, 2003
2Contents
- -FEM modelling
- Wheel Element
- Wheel-Rail Contact Element
- Example Single Wheelset
- Example Bicycle Dynamics
- Conclusions
3FEM modelling
2D Truss Element
4 Nodal Coordinates
3 Degrees of Freedom as a Rigid Body leaves 1
Generalized Strain
Rigid Body Motion
Constraint Equation
4Wheel Element
Nodes
Generalized Nodes
Position Wheel Centre
Euler parameters
Rotation Matrix R(q)
Contact Point
In total 10 generalized coordinates
Rigid body pure rolling 3 degrees of freedom
Impose 7 Constraints
5Wheel Element
Strains
Holonomic Constraints as zero generalized strains
Elongation
Lateral Bending
Contact point on the surface
Wheel perpendicular to the surface
Radius vector
Rotated wheel axle
Normalization condition on Euler par
Surface
Normal on surface
6Wheel Element
Slips
Non-Holonomic Constraints as zero generalized
slips
Velocity of material point of wheel at contact in
c
Generalized Slips
Longitudinal slip
Radius vector
Lateral slip
Two tangent vectors in c
Angular velocity wheel
7Wheel-Rail Contact Element
Nodes
Generalized Nodes
Position Wheel Centre
Euler parameters
Rotation Matrix R(q)
Contact Point
In total 10 generalized coordinates
Rigid body pure rolling 2 degrees of freedom
Impose 8 Constraints
8Wheel-Rail Contact Element
Strains
Holonomic Constraints as zero generalized strains
Distance from c to Wheel surface
Distance from c to Rail surface
Wheel and Rail in Point Contact
Wheel Rail surface
Normalization condition on Euler par
Local radius vector
Normal on Wheel surface
Two Tangents in c
9Wheel-Rail Contact Element
Slips
Non-Holonomic Constraints as zero generalized
slips
Velocity of material point of Wheel in contact
point c
Generalized Slips
Longitudinal slip
Lateral slip
Wheel Rail surface
Spin
Two Tangents in c
Normal on Rail Surface
Angular velocity wheel
10Single Wheelset
Example
Klingel Motion of a Wheelset
Wheel bands S1002 Rails UIC60 Gauge 1.435
m Rail Slant 1/40
FEM-model 2 Wheel-Rail, 2 Beams, 3 Hinges Pure
Rolling, Released Spin 1 DOF
11Single Wheelset
Profiles
Wheel band S1002
Rail profile UIC60
12Single Wheelset
Motion
Klingel Motion of a Wheelset
Wheel bands S1002 Rails UIC60 Gauge 1.435
m Rail Slant 1/40
Theoretical Wave Length
13Single Wheelset
Example
Critical Speed of a Single Wheelset
Wheel bands S1002, Rails UIC60 Gauge 1.435 m,
Rail Slant 1/20 m1887 kg, I1000,100,1000
kgm2 Vertical Load 173 226 N Yaw Spring Stiffness
816 kNm/rad
FEM-model 2 Wheel-Rail, 2 Beams, 3
Hinges Linear Creep Saturation 4 DOF
14Single Wheelset
Constitutive
Critical Speed of a Single Wheelset
Linear Creep Saturation according to Vermeulen
Johnson (1964)
Tangential Force Maximal Friction Force
Total Creep
15Single Wheelset
Limit Cycle
Limit Cycle Motion at v131 m/s
Critical Speed of a Single Wheelset
Vcr130 m/s
16Bicycle Dynamics
Example
Bicycle with Rigid Rider and No-Hands
Standard Dutch Bike
FEM-model 2 Wheels, 2 Beams, 6 Hinges Pure
Rolling 3 DOF
17Bicycle Dynamics
Root Loci
Stability of the Forward Upright Steady Motion
Root Loci from the Linearized Equations of
Motion. Parameter forward speed v
18Bicycle Dynamics
Motion
Full Non-Linear Forward Dynamic Analysis at
different speeds
Forward Speed v m/s
18
14
11
10
5
0
19Conclusions
- Proposed Contact Elements are Suitable for
Modelling Dynamic Behaviour of Road and Track
Guided Vehicles.
Further Investigation
- Curvature Jumps in Unworn Profiles, they Cause
Jumps in the Speed of and Forces in the Contact
Point. - Difficulty to take into account Closely Spaced
Double Point Contact.