29/06/04, Page 1
SAMTECH Group SAMTECH Group Machine Tools Applications Machine Tools Applications
August 10, 2004 August 10, 2004
SL/04/SAM/MKG_ppt/28an_a
SAMTECH Group SAMTECH Group Machine Tools Applications Machine - - PowerPoint PPT Presentation
SAMTECH Group SAMTECH Group Machine Tools Applications Machine Tools Applications August 10, 2004 August 10, 2004 29/06/04, Page 1 SL/04/SAM/MKG_ppt/28an_a SAMTECH Mechatronic Mechatronic Solution Solution SAMTECH using Finite Element
29/06/04, Page 1
SL/04/SAM/MKG_ppt/28an_a
29/06/04, Page 2
29/06/04, Page 3
What is the problem of High Speed Machine Tools? Why is it difficult to analyze numerically Machine Tools? Methodology suited to compute dynamic behavior of vibrating mechanical systems SAMCEF for Machine Tools Presentation of examples Conclusions and future investigations&developments Questions/Answers
29/06/04, Page 4
In order to increase the productivity, it is more and more necessary to increase the velocity of Machine Tools moving parts If maximum velocities are higher and higher, accelerations are increasing significantly Accelerations prescribed by controllers are generating huge inertia forces due to the mass of moving parts If the mass of moving parts is reduced due to the removal of material, there can be a loss of stiffness and thus of precision Inertia effects are thus major excitations for vibrations of High Speed Machine Tools in working conditions Vibrations cannot be cut in pieces
29/06/04, Page 5
Need of a coherent global model describing the whole dynamic flexible behavior of the Machine Tool (mechanisms + structures + controllers) Sensitive components are the ones transmitting the vibrations (guideways, linear motors, ball screws…) These components have to be modeled using advanced numerical computation techniques able to manage the interactions between moving flexible bodies NL Finite Element coupled with kinematical constraints and large relative motions contact techniques are appropriate to achieve this goal SAMTECH has been working on this innovative approach since 15 years and brought recently this new powerful and integrated technology to the market, starting with the sector of Machine-Tools
29/06/04, Page 6
Examples of possible sources of “Non-Linearities” involved in FEM models of Machine Tools:
Need of an implicit unconditionally stable time integration scheme:
sense
29/06/04, Page 7
Augmented Lagrangian form of the equations of motion
⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = − − = = + + ) , ( ) , , ( ) (
inertia internal externa
t q k g g g t q q g p k q
l T
φ φ λ r r r r & r r r & & r B M
[ ]
3 / 1 , ) ( ) ( ) 1 ( ) , , ( ) , , ( ) 1 ( ) 1 (
1 1 1 1 1 1 1
∈ ⎭ ⎬ ⎫ ⎩ ⎨ ⎧ = + − + − = + − +
+ + + + + + +
α αφ φ α α α α α
n n n n n n n n n T n n T n n
q k q k t q q t q q q & & & & g g λ B λ B M HHT-Form of discretised equilibrium equation
Equations are solved iteratively each time integration step: “HHT-Form” Mixed system of equations in FEM-DOF and Lagrange Multipliers:
M: Mass matrix B: Gradient of the constraint matrix φ: Kinematical constraints
Classical non-linear Finite Element equations ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ + = − − = = ) K( ) C( with ) , , (
internal inertia internal externa
q q g g g g t q q g q
l
& r r r r & r & & r M
λ: Lagrange multiplier p: Penalty factor k: Scaling factor
29/06/04, Page 8
General integrated environment giving access to different levels of model fidelity and different analysis types from the same CAD based model
29/06/04, Page 9
Linear Motor Gantry Rotational axis Parallel Translation axis C-Frame Box-In-Box High speed dynamics Flexible Screw
29/06/04, Page 10
29/06/04, Page 11
Finite Element Analysis using SAMCEF
Structure/Contact FE Analysis
creation/restitution
and controller design:
MATLAB Simulink
Processing Environment
Optimization of parameterized mechatronic systems using BOSS quattro
29/06/04, Page 12
NL Finite Elements (linear motor…) and S.E. Flexible kinematical joints (Guideways, flexible ball screw…) Digital Control (Sensors and Actuators)
SIMULINK
…
Rigid kinematical joints
29/06/04, Page 13
Description of all the parts and the connection devices of the mechanical system Different Model Fidelity levels (rigid body, detailed meshed structure or Super-Element) Different types of connection devices (rigid or flexible kinematical joints, contact/friction conditions) “Sensors” and “actuators” to be connected to the digital control boxes Output of linearized mechanical models for control design
29/06/04, Page 14
Linear Static Analysis (SAMCEF Asef) Modal Analysis (SAMCEF Dynam) Super-Element creation and restitution (SAMCEF Dynam) Integrated environment with MBS and NL FEA
29/06/04, Page 15
SAMCEF Mecano Control Box
Next Sample Time Command Vector Updated State Variables Current Time Sensor Vector State Variable
Sensor Variables:
Positions, Displacements, Velocities, Accelerations and/or Reaction Forces
Command Vector:
Forces, Displacements, Velocities, Accelerations and/or functions
29/06/04, Page 16
Controllers can be imported from a Controller Design Software (for example MATLAB Simulink using Real Time Workshop (RTW) to generate C subroutine) They can also be fully written by the user (open to in-house controllers) Linked with SAMCEF Mecano
29/06/04, Page 17
Positions, displacements, velocities, accelerations as well as reaction forces are provided independently by SAMCEF Mecano to the controller Activation and deactivation times Several controllers can be used at the same time The same controller can be used several times Different schematic boxes (sensors, controllers and actuators) can be connected between them Prescribed values and gains can be defined through the SAMCEF data set
29/06/04, Page 18
The computation time step is bound to the time discretization of the Controller tComp <= t Controller Automatic time step choice between instants of “Rendez-vous”
29/06/04, Page 19
29/06/04, Page 20
Modeling defined directly on the CAD geometry Easy modifications of the modeling fidelity for bodies
FE package
Easy modifications of the modeling fidelity for connections
Easy switch from an analysis type to another one (SAMCEF Mecano, SAMCEF Dynam, SAMCEF Asef…) Import of controllers from MATLAB Simulink
29/06/04, Page 21
29/06/04, Page 22
29/06/04, Page 23
29/06/04, Page 24
29/06/04, Page 25
29/06/04, Page 26
29/06/04, Page 27
Linear motors Ball screws Guideways Rack and pinion Ball screw Guideways & rack and pinion Linear motor Definition on geometry Multiple models for each physical component Only additional parameters between different models
29/06/04, Page 28
29/06/04, Page 29
Driven hinge joint SCRF stiffening elements MEAN stiffening MEAN nut node
COUPLING JOINT
ROTOR STATOR BEAM STICK BEARING elements stiffening elements MEAN STICK STICK
29/06/04, Page 30
29/06/04, Page 31
29/06/04, Page 32
Efficient choice of the right controller Fine tuning of the initial values for the gains
29/06/04, Page 33
29/06/04, Page 34
Integration of various Digital Controllers into SAMCEF Non-Linear Simulations in order to get a mechatronic FE model combining mechanisms, structures and controllers in one single model
29/06/04, Page 35
Static stresses evaluation Vibration Modes calculation Super-Element reduction Time response (static, kinematical or dynamic)
29/06/04, Page 36
29/06/04, Page 37
0.1 0.2 0.3 0.4 0.5 0.01 0.02 0.03 0.04 0.05 0.06 time(s) position(m) step response: siemens controller step input position horizontal slider
Siemens Controller Position Controller Speed Controller Current Controller Controllers can be Imported from a Controller Design Software (for example MATLAB Simulink) Fully written by the user
29/06/04, Page 38
Bode Diagram Frequency (Hz) Phase (deg) Magnitude (dB)
10 10
110
2180 L = 34 cm L = 37.75 cm L = 41.5 cm L = 45.25 cm L = 49 cm
Bode Diagram Frequency (Hz) Phase (deg) Magnitude (dB)
20 40 60 80 10
110
2L = 34 cm L = 37.75 cm L = 41.5 cm L = 45.25 cm L = 49 cm
Identification of the system: FRFs from the motor force input to
Clear dependence of the eigenfrequencies on the arm length
29/06/04, Page 39
Movement of the linear motor and the end-tip of the beam
29/06/04, Page 40
0.5 1 1.5 2 2.5 3 3.5 4Pole placement controller
0.5 1 1.5 2 2.5 3 3.5 4LQR controller
0.5 1 1.5 2 2.5 3 3.5 4MIMO controller
0.5 1 1.5 2 2.5 3 3.5 4MISO controller
0.5 1 1.5 2 2.5 3 3.5 4HAC/LAC controller
0.5 1 1.5 2 2.5 3 3.5 4met regelaar zonder regelaar
PD controller
29/06/04, Page 41
0.2 0.4 0.6 0.8
0.05 0.1 0.15 0.2 time(s) position(m) sinus response: synchrone horizontal position input horizontal position vertical position input vertical position 0.2 0.4 0.6 0.8
0.05 0.1 0.15 0.2 time(s) position(m) sinus response: non synchrone horizontal position input horizontal position vertical position input vertical position
0.05 0.06 0.08 0.1 0.12 0.14 0.16 X(m ) Y(m) Circle test : syn ch ron e traject in pu t traject ou tpu t
0.05 0.06 0.08 0.1 0.12 0.14 0.16 X(m) Y(m) Circle test : non synchrone traject input traject output
Tuning individual controllers Concurrent Tuning of all controllers
29/06/04, Page 42
29/06/04, Page 43
29/06/04, Page 44
RamZ Fixed Part SliderX
9517
SliderY
9443 8343 9484
29/06/04, Page 45
29/06/04, Page 46
29/06/04, Page 47
29/06/04, Page 48
29/06/04, Page 49
29/06/04, Page 50
CAD Geometry defined in SAMCEF Field or recovered from CAD through STEP, IGES, CATIA exchange
X slide Y slide Portal Ram Electro spindle
Guides ways
SAMCEF Mecano assembly conditions (rigid bodies, rigid kinematical joints and flexible sliders, gluing) to connect moving flexible bodies
Super-Elements to manage large meshes
29/06/04, Page 51
0.5 1 1.5 2 2.5 3 3.5 4
2000 4000 6000 Tijd (s) Totale fout ( μ m)
met regelaar zonder regelaar
Performance Hac-Lac controller for ramp position trajectory
0.5 1 1.5 2 2.5 3 3.5 4
1000 2000 3000 4000 5000 Tijd (seconden) Volgfout ( ) 탆
Vibration controller (LAC controller) Linear Motor controller (HAC controller)
Time Tracking error Tracking error
Without and with controller
Time
29/06/04, Page 52
Same mechanical modeling principles
as previous examples
Linear motors implemented in the
generic non linear force element of SAMCEF Mecano
Digital controllers acting on axes
modeled in MATLAB Simulink and integrated in SAMCEF Mecano
Comparison
between numerical results and experiments: less than 8%
precision of execution
29/06/04, Page 53
Parallel machine using 3 linear motors and 1 rotational axis implemented in the non-linear pre- post processor SAMCEF Field Rigid and flexible behavior
29/06/04, Page 54
Rigid and Flexible behavior on parallel machine
29/06/04, Page 55
control structure mechanism
Three axes machine Based on a “delta” parallel structure Strokes : X 500 Y 500 Z 300 Linear speed: 100 m/min on all axes Acceleration: 35 m/s² on all axes Jerk: 5000 m/s3 on all axes Repeatability: 5 µm on all axes Spindle: 15 kW 40 000 rpm
29/06/04, Page 56
Mechanism Mechanism + Structure Mechanism + Control Mechanism + Structure + Control
Mechanism + Structure
29/06/04, Page 57
29/06/04, Page 58
29/06/04, Page 59
29/06/04, Page 60
Introduction of harmonic response (SAMCEF Repdyn) in SAMCEF Field Introduction of thermal transient analysis (SAMCEF Thermal) in SAMCEF Field Extension of the Mechatronic approach by Finite Element Analysis to other industrial sectors (robotics, textile machines, automotive, aeronautics, space, defense…)
29/06/04, Page 61
29/06/04, Page 62
Machine Tools manufacturers currently develop high speed machines were inertial effects are becoming critical A Rigid Multi-Body model is clearly not sufficient Super-Elements are not always sufficient (inefficient management of contact conditions between moving flexible bodies) SAMCEF for Machine Tools proposes an implicit NL Dynamic solution using NL Finite Element/Contact techniques SAMCEF for Machine Tools is also a global and integrated CAD based solution covering all the Machine Tools components (linear motor, guideways, ball screw, gears, rack-and-pinions…), the model fidelity levels and analysis types (linear or not) requested by this industrial sector