RD.14/185701.1
Ricardo Software Engine Dynamics Simulation
ENGDYN DOCUMENTATION/USER MANUAL VERSION 2014.1 May 2014
Contents
Contents WHAT’S NEW IN ENGDYN 2014.1? .............................................................................................. V A. 1 2
B. 1
2
3
4
5
6 7
KNOWLEDGE CENTRE ........................................................................................................... 1 W HAT IS ENGDYN? ................................................................................................................... 1 TUTORIALS .................................................................................................................................. 3 2.1 Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis .......................... 3 2.2 Tutorial 2: Crankshaft Dynamic Analysis ........................................................................43 2.3 Tutorial 3: NVH Analysis .................................................................................................79 2.4 Tutorial 4: Block Stress Analysis ...................................................................................125 2.5 Tutorial 5: Crankshaft Stress Analysis ..........................................................................149 2.6 Tutorial 6: EHL Big End Bearing Analysis .....................................................................183 HELP .....................................................................................................................................213 USING ENGDYN ....................................................................................................................213 1.1 Overview .......................................................................................................................213 1.2 Getting Started ..............................................................................................................213 1.3 Description of the Main Panel .......................................................................................214 ENGDYN MODELS .................................................................................................................223 2.1 Overview .......................................................................................................................223 2.2 ENGDYN Co-ordinate System ......................................................................................223 2.3 Finite Element Model Data ............................................................................................224 2.4 Crank Train Model .........................................................................................................225 2.5 Cylinder Block Model ....................................................................................................231 2.6 Journal Bearing Oil Film Model .....................................................................................232 2.7 In-Cylinder Model ..........................................................................................................235 2.8 Gas Cylinder Pressure ..................................................................................................236 MODEL GENERATION ...............................................................................................................237 3.1 General ..........................................................................................................................237 3.2 Configure Engine ..........................................................................................................237 3.3 Define Models ...............................................................................................................246 3.4 Editing the Crank Train Model .......................................................................................253 3.5 Editing the Cylinder Block Model ..................................................................................324 SOLUTION ...............................................................................................................................353 4.1 General ..........................................................................................................................353 4.2 Lubrication .....................................................................................................................354 4.3 Loading..........................................................................................................................359 4.4 Evaluate Solution ..........................................................................................................374 POST-PROCESSING .................................................................................................................391 5.1 General ..........................................................................................................................391 5.2 Selecting Loadcases .....................................................................................................391 5.3 Selecting Crank Train Modes ........................................................................................392 5.4 Plotting EHL Results .....................................................................................................393 5.5 Plotting Results .............................................................................................................402 5.6 Exporting Results ..........................................................................................................407 5.7 Backsubstitution ............................................................................................................413 5.8 Animate Results ............................................................................................................413 CRANK SHAFT STRESS ANALYSIS .............................................................................................419 6.1 Overview .......................................................................................................................419 6.2 Using the Graphical User Interface ...............................................................................434 CYLINDER BLOCK ANALYSIS .....................................................................................................463 7.1 Overview .......................................................................................................................463
i
KNOWLEDGE CENTRE What is ENGDYN?
7.2 Frequency Response and Acoustic Analysis ............................................................... 463 7.3 Quasi-Static Analysis .................................................................................................... 468 7.4 Using the Graphical User Interface .............................................................................. 476 8 ADDITIONAL POST-PROCESSING .............................................................................................. 499 8.1 Overview ....................................................................................................................... 499 8.2 Quasi-Static Analysis of a Piston .................................................................................. 499 8.3 Applying Loads to a Finite Element Model of a Connecting Rod ................................. 500 C.
COMMAND FILES ............................................................................................................... 505
1
PRE-PROCESSOR COMMANDS ................................................................................................. 506 1.1 AXES ............................................................................................................................ 507 1.2 BEARING ...................................................................................................................... 508 1.3 BLOCK .......................................................................................................................... 515 1.4 CHILD ........................................................................................................................... 519 1.5 CONROD ...................................................................................................................... 521 1.6 CRANK ......................................................................................................................... 525 1.7 CYLINDER .................................................................................................................... 530 1.8 DAMPER ...................................................................................................................... 534 1.9 DIRECTION .................................................................................................................. 539 1.10 DRIVE ........................................................................................................................... 540 1.11 ELEMENT ..................................................................................................................... 541 1.12 ENGINE ........................................................................................................................ 548 1.13 IN_CYLINDER .............................................................................................................. 549 1.14 LINK .............................................................................................................................. 550 1.15 LOADING ...................................................................................................................... 552 1.16 LUBRICANT ................................................................................................................. 557 1.17 MASS ............................................................................................................................ 559 1.18 MATERIAL .................................................................................................................... 564 1.19 MOUNT ......................................................................................................................... 565 1.20 NODE............................................................................................................................ 567 1.21 OPEN ............................................................................................................................ 568 1.22 PROFILE ...................................................................................................................... 569 1.23 TITLE ............................................................................................................................ 572 2 SOLVER COMMANDS ............................................................................................................... 573 2.1 BEARING ...................................................................................................................... 574 2.2 BLOCK .......................................................................................................................... 575 2.3 COUPLING ................................................................................................................... 576 2.4 DAMPING ..................................................................................................................... 578 2.5 LUBRICANT ................................................................................................................. 580 2.6 OPEN ............................................................................................................................ 581 2.7 SOLUTION ................................................................................................................... 582 3 POST-PROCESSOR COMMANDS ............................................................................................... 592 3.1 CRANK_ANALYSIS...................................................................................................... 592 3.2 BLOCK_ANALYSIS ...................................................................................................... 606 D.
THEORY ............................................................................................................................... 616
1 2
INTRODUCTION........................................................................................................................ 616 STATIC SOLUTIONS ................................................................................................................. 617 2.1 Determinate Solution .................................................................................................... 617 2.2 Indeterminate Solution .................................................................................................. 617 3 DYNAMIC SOLUTIONS .............................................................................................................. 618 4 BEARING MODEL ..................................................................................................................... 619 4.1 Mobility Method ............................................................................................................. 621 4.2 Solution of Reynolds Equation on a Computational Mesh ........................................... 629
ii
Contents
4.3 4.4 4.5 4.6 4.7 4.8 4.9 E.
Oil Viscosity ...................................................................................................................641 Boundary Lubrication Model .........................................................................................643 Integration of Bearing Forces and Moments .................................................................645 Bearing Thermal Balance Solution ...............................................................................645 Compliant Model for EHD Bearing ................................................................................650 Nomenclature ................................................................................................................653 References ....................................................................................................................655
APPENDICES .......................................................................................................................657
1
CALCULATION OF CRANKSHAFT STIFFNESS’ USING THE FINITE ELEMENT METHOD .......................657 1.1 Overview .......................................................................................................................657 1.2 Calculation of Crankshaft Web Stiffnesses ...................................................................657 1.3 Calculation of Crankshaft Element Stiffnesses .............................................................660 2 A METHOD OF EVALUATING THE MASS AND STIFFNESS PROPERTIES OF A FLYWHEEL .................663 2.1 Overview .......................................................................................................................663 2.2 Modelling Approach ......................................................................................................663 2.3 Calculation of Mass and Stiffness Properties ...............................................................663 3 CALCULATION OF MATERIAL STRENGTH PROPERTIES FOR CRANKSHAFT STRESS ANALYSIS .......669 3.1 Overview .......................................................................................................................669 3.2 Base Strengths ..............................................................................................................669 3.3 Elevated Strengths ........................................................................................................669 4 CALCULATION OF STRESS CONCENTRATION FACTORS FOR CRANKSHAFT STRESS ANALYSIS ......671 4.1 Overview .......................................................................................................................671 4.2 Nomenclature ................................................................................................................671 4.3 Journal Fillets ................................................................................................................671 4.4 Oil Hole Breakouts ........................................................................................................673 5 TREATMENT OF NOTCH SENSITIVITY IN CRANKSHAFT STRESS ANALYSIS ....................................677 5.1 Overview .......................................................................................................................677 5.2 Treatment in Finite Element Analysis ...........................................................................678 6 USING OUTPUT FROM THE RICARDO PROGRAM VALDYN AND TVFORCED IN CRANKSHAFT STRESS CALCULATIONS ...................................................................................................................681 6.1 Using VALDYN for Crankshaft Stress Calculations ......................................................681 6.2 Using TVFORCED for Crankshaft Stress Calculations ................................................684 7 SAFETY FACTOR CALCULATIONS ..............................................................................................687 7.1 The Goodman Diagram .................................................................................................687 7.2 Goodman Diagram Construction...................................................................................688 7.3 Equivalent Stress Options .............................................................................................689 7.4 Multi-axial Fatigue Safety Factor...................................................................................691 7.5 Dang Van Fatigue Safety Factor ...................................................................................692 7.6 LinearSWT Safety Factor ..............................................................................................695 7.7 Cycles to Failure Calculation.........................................................................................696 8 RADIATED NOISE CALCULATIONS .............................................................................................699 8.1 Overview .......................................................................................................................699 8.2 Acoustic Equations and the Rayleigh Integral ..............................................................699 8.3 Radiation Efficiency .......................................................................................................703 8.4 Multiple Face Sets .........................................................................................................703 9 MODAL FREQUENCY RESPONSE CALCULATIONS .......................................................................705 9.1 Overview .......................................................................................................................705 9.2 Calculation of Modal Contributions and Vibration Response ........................................705 10 APPLYING LOADS PREDICTED BY VALDYN TO THE CYLINDER BLOCK AND CRANKSHAFT MODELS 707 10.1 Overview .......................................................................................................................707 10.2 Writing Force Profile Data From VALDYN ....................................................................707 10.3 Applying Force Profile Data to the ENGDYN Model .....................................................710
iii
KNOWLEDGE CENTRE What is ENGDYN?
10.4 Applying Force Profile Data to the Cylinder Block Model ............................................. 710 11 ASCII DATA FILE FORMAT READ BY SDF_READ_ASCII() ........................................................ 715 11.1 Introduction ................................................................................................................... 715 11.2 Format........................................................................................................................... 715 12 MATUTIL .......................................................................................................................... 720 12.1 Overview ....................................................................................................................... 720 12.2 Using matutil ................................................................................................................. 720 12.3 Theory ........................................................................................................................... 727
iv
What’s New in ENGDYN 2014.1?
What’s New in ENGDYN 2014.1? Enhancements Feature
Documentation
Updated NVH Analysis Tutorial
Knowledge Tutorial 3
Centre
Updated Block Stress Analysis Tutorial
Knowledge Tutorial 4
Centre
Various bug fixes
Release Notes
v
KNOWLEDGE CENTRE What is ENGDYN?
A. KNOWLEDGE CENTRE 1 What is ENGDYN? ENGDYN is a computer program for analysing the dynamics of the engine and in particular the dynamics of the crank train and its interaction with the cylinder block. ENGDYN provides a number of different solution techniques for predicting engine dynamics using models of varying degrees of sophistication. The crank train and cylinder block models can either be defined as rigid, compliant and dynamic. In its simplest form ENGDYN can be used to perform a statically-determinate solution, whilst in its most sophisticated form it can be used to predict the time-domain response of the 3-dimensional vibration of the coupled crank train and cylinder block system with nonlinear oil films at each of the main journal bearings. This flexibility enables the user to generate an engine model and to perform a solution to meet his particular needs. ENGDYN provides a Graphical User Interface (GUI) to enable the user to perform the model generation, solution and results presentation phases within an easy to use graphical environment. The GUI contains a built-in unit’s converter and automatically converts parameters from their defined units to SI units. Alternatively, for the more experienced user and to provide compatibility with models that were generated before the development of the GUI, ENGDYN also provides a non-graphical environment that uses command data files for model generation, solution and results presentation. ENGDYN uses a Ricardo binary standard data file to store both the model and the results of the solution. When saving, ENGDYN stores the model at the current position and as the model is built the file is appended to. Similarly, once a simulation has been executed, the results are appended to the file. Crank train simulation using ENGDYN consists of three stages. Firstly, the engine model must be generated and consists of the crank train and cylinder block. The GUI is designed such that the user builds the model in a sequential order by using a series of forms. Once the key engine parameters have been defined ENGDYN draws the reduced model of the crank train that the user can pick to edit particular items. Finite element models of the crank train and cylinder block can be viewed with the reduced models. ENGDYN performs its own finite element solution, to evaluate for example crankshaft web stiffness’. Once the model has been generated, the GUI can be used to generate the input command file for performing a solution. This file contains the engine conditions to be simulated, the type of solution to be performed, the models to be included in the solution and the solution parameters for controlling the solution. The solver performs the simulation and appends the results to the binary standard data file. An ASCII file is also written summarising the results of the solution for each load case. Finally, once the simulation has finished running, the user is able within the GUI to select load cases for post-processing. The interface lists all the load cases within the file and for each load case gives the solution parameters defined for that case. The interface currently allows simulation results to be plotted, animated and exported to an ASCII file. 1
KNOWLEDGE CENTRE What is ENGDYN?
The interface has a built in library of pre-defined plots that enables the user to rapidly view the results. The plots can either be printed or saved to file. The user is also able to export results to an ASCII file to allow further data manipulation by the user not currently supported by ENGDYN. In addition, for solutions in which the displacement or vibratory motion of the engine is predicted, the results can be animated interactively in both the time and frequency domains. In addition, the loads calculated by ENGDYN can be used to evaluate the stresses and fatigue safety factors at the critical locations of the crankshaft and can also be applied to the cylinder block either in the frequency domain to perform modal frequency response and acoustic analyses or in the time domain to perform quasi-static analyses. The user can then perform further analyses if required. A new solution can be executed and on completion the results are appended to the binary standard data file.
2
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
2 Tutorials 2.1 Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft
Analysis 2.1.1
Introduction
Objective: Gain a general understanding of how to use ENGDYN GUI and carry out the simplest analysis using user defined (non-FE) data. The user will perform statically determinate analysis to calculate loads for both a Mobility and Hydrodynamic (HD) bearing analysis User will also gain experience in using plotting and animation features. Items covered:
Building the engine model Bearing analysis using the Mobility Method Bearing analysis using the Hydrodynamic Method (Finite Volume Solver) Plotting and Animation
Estimated duration: <0.5 days (overall including performing solutions) Engine:
Audi V6 TDI V6 High speed direct injection diesel 90 degree vee-angle 30 degree crank-pin offset 4 main bearings 6 big end bearings
3
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
The data for this tutorial was obtained from a benchmark study of this engine by Ricardo. Required Files: ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\v6\cp700.PRES ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\v6\cp1000.PRES ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\v6\cp2000.PRES ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\v6\cp4000.PRES
2.1.2
Getting Started
Copy the necessary files from the example directory to a working directory and ensure that you have write permissions for all the files. Start engdyn On Unix or Linux platforms simply type engdyn On windows click on the shortcut, otherwise go to Start>Programs>Ricardo Software>2014.1>Mechanical Suite>ENGDYN>ENGDYN
2.1.3
Building engine model
2.1.3.1 Configure Engine Select ‘Configure Engine’ from the buttons on the left side of the Main Panel
Complete the Engine Configuration Panel as shown
4
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Note that warning messages will pop up as long as all the needed data is not properly entered. If you close the Engine Configuration Window while warning messages still appear, all the data will be lost. Therefore, please make sure all the data is properly entered before closing this window.
5
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
6
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
To enter the main bearings non-uniform information, the user must use the mouse to click on the red dimensions in the diagram in the panel. On completion select ‘Apply’ which will display
Press OK button. The crankshaft model will appear in the Main Panel as shown Use CTRL + middle mouse in order to zoom in and zoom out Use SHIFT + middle mouse in order to rotate
2.1.3.2 Define Models Select ‘Define Models’ from the buttons on the left side of the GUI.
7
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
Complete the model definitions for the crankshaft, cylinder block, in-cylinder and connecting rod as shown
2.1.3.3 Edit Cranktrain Select ‘Edit Cranktrain’ from the buttons on the left side of Main Panel.
Items are edited by selecting from the list on the left of the panel Many items shown in the panel below need not be defined due to the type of model being edited. Highlight ‘Crank Web’ as shown This will then display each of the web elements in green Click on ‘Select All’
8
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Each of the web elements will then turn red
Click on ‘Edit Selected’ to display the Web Panel Complete the ‘Web’ panel as shown below
The web thickness and journal lengths are taken from the data entered in the ‘configure engine’ section The Web panel allows counterweight data to be entered for the crank webs The counterweight geometry is defined only for balancing the crankshaft at a later stage, NOT for defining the mass properties. As we shall be defining the mass properties explicitly in this example (rather than using an FE model), we do not require any counterweight data to be added Select ‘OK’ Highlight ‘Big End Bearing’ This will then display each of the pin journal nodes in green Click on ‘Select All’ Each of the main journal nodes will then turn red
9
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Click on ‘Edit Selected’ to display the Pin Bearing Panel Complete the ‘Bearing’ panel as shown below
The oil hole angular position can be defined either using the Height column (as in this case) or using the Position column. The value not supplied is calculated from the journal diameter. The remaining tabs (Mesh, Material, Profile) need not and cannot be edited because the model type is ‘Mobility’. These tabs are only required with Hydrodynamic or Elastohydrodynamic models Select ‘OK’ Highlight ‘Main Bearing’ Click on ‘Select All’ to highlight all main bearing nodes Click on ‘Edit Selected’ to display the Main Bearing Panel Complete the Main Bearing Panel as shown
10
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
The ‘Oilholes’ tab cannot be edited because the bearing is fed by a groove. The remaining tabs (Mesh, Material, Profile, Stiffness) need not and cannot be edited because the model type is ‘Mobility’. These tabs are only required with Hydrodynamic or Elastohydrodynamic models Select ‘OK’ Highlight ‘Connecting Rod’ Click on ‘Select All’ to select all pin journal nodes Click on ‘Edit Selected’ Complete the ‘Conrod’ panel as shown.
Select ‘OK’ Highlight ‘Piston’
11
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Click on ‘Select All’ to select all pin journal nodes Click on ‘Edit Selected’ Complete the ‘Piston’ panel as shown.
Highlight ‘Cranknose Assembly’ In this instance since there is only one cranknose there is no need to use Select All or to pick using the mouse. The program automatically selects the cranknose and displays it in red Select ‘Edit Selected’ to display the Cranknose Assembly Panel Complete the Cranknose Assembly Panel as shown
The ‘Element Length’ field does not have to be entered accurately because the crankshaft is not modelled using an FE model. This data will only be used for showing a representative model on the screen. Select ‘OK’ Highlight ‘Flywheel Assembly’
12
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Again there is no need to use Select All or to pick using the mouse. The program automatically selects the cranknose and displays it in red Select ‘Edit Selected’ to select the Flywheel Assembly Panel Complete the ‘Flywheel Assembly’ panel as shown
Select on ‘OK’ Highlight ‘Lumped Masses’ 13
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis The nodes are numbered from the front of the crankshaft, number one representing the damper hub For a rigid or compliant crankshaft model it is not necessary to define masses at all the nodes Select node 3 using the left mouse button. This is the web node of the first web on the crankshaft axis. Alternatively nodes can be picked by dragging the mouse with the left button held down To make multiple selections, hold the SHIFT key down whilst selecting nodes Select Edit Selected to display the Lumped Mass Panel as shown Complete the panel as shown
Data can either be defined with known properties or using a geometric shape and can be with respective to a Cartesian or Polar coordinate system. Select the Add or Update button to add the data Multiple masses can be added to each node Select OK The mass will be shown in white
14
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Repeat this process defining the masses as listed in the table below. Node
3 5 8 10 12 14 17 19 21 23 26 28
Description
Web 1 counterweight Pin 1 mass Pin 2 mass Web 2 counterweight Web 3 counterweight Pin 3 mass Pin 4 mass Web 4 counterweight Web 5 counterweight Pin 5 mass Pin 6 mass Web 6 counterweight
Mass [kg] 2.610 0.408 0.408 1.500 1.760 0.408 0.408 1.760 1.500 0.408 0.408 2.630
x [mm] 0.470 -0.875 0.875 -0.075 0.300 -0.875 0.875 -0.050 0.297 -0.875 0.875 -0.051
Offset y [mm] -23.220 0.000 0.000 -20.460 13.090 0.000 0.000 12.280 7.215 0.000 0.000 10.140
z [mm] 10.270 0.000 0.000 7.247 12.280 0.000 0.000 13.100 -20.510 0.000 0.000 -22.720
When all lumped masses have been entered the model should appear as shown
Click on ‘Define Material’ to display the Crankshaft Material Properties Panel
15
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Complete the Crankshaft Material Properties Panel as shown
Select on ‘OK’ Click on ‘Calculate Masses’ The message ‘Calculation of mass properties completed successfully, Balance Not Set’ should appear in the message box at the bottom of the Main Panel. Click on ‘Set Balance’ to display the Primary Balance Panel as shown
No additional balancing is required because the Lumped Masses include the effect of balancing the crankshaft. Click on ‘Assemble Model’ The message ‘Model assembly completed successfully, Block not assembled’ should appear in the message box at the bottom of the Main Panel Click on ‘OK’ Click on the ‘Dismiss’ button of the Edit Cranktrain panel
16
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis The message ‘Ready to Perform Solution’ should appear in the message box at the bottom of the Engdyn screen Save the model
2.1.4
Bearing analysis using the Mobility Solution
2.1.4.1 Define Lubricant Properties Click on ‘Lubrication’ button from the buttons on the left hand side of the Main Panel
Use the Browse button to select the lubricant SAE5W30 from the database By default the program should initially select the database directory at ..\Ricardo\2014.1\Common\Materials\Lubricants This database contains the most common lubricants Use either Add or Update to add the lubricant
Click on OK 2.1.4.2 Define Loading Conditions The cylinder pressure diagram and any additional loadings (e.g., loads from Valdyn and gravity forces) are entered in this step.
17
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Click on the ‘Loading’ button on the Main Panel
This will display the Loading Definition Panel as shown
A number of different loading maps can be defined, Full Load, Part Load and No Load. The solver will interpolate at speeds between those defined using this panel Type in a speed of 750 rev/min Position the mouse over the File Name column and use the right button to display the pop-up menu. Use ‘Select Pressure file’ to select the file cp700.PRES from your working directory. The panel will appear as shown.
You may wish to remove the pathname in front of the file.
18
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis The pressure file can contain just pressures (at equal intervals) or as in this case pressures and angles. You may wish to inspect and understand the file using an appropriate editor Define the Ambient and Crankcase Pressure as 1.0 bar The crankcase pressure is only currently used for Hydrodynamic and Elastohydrodynamic bearing models to define the boundary condition at the edge of the bearing Add an extra line to the table by positioning the mouse over the left column (line number) and using the right button to display a pop-up. Select ‘Insert Row After’ Complete the panel as shown
The remaining tabs, Force Profile, Force Equation and Distortion need not be completed for this tutorial. When a model has a piston pin offset or crank offset, the cylinder pressure diagram is often redefined by ENGDYN internally to have the interval angle reduced to 0.25 degrees (sometimes lower) by a process of interpolation. This is because the offsets cause the TDC angle to change slightly and it is important to have an accurate cylinder pressure definition at TDC to achieve good results. Use the Plot button to display the applied loads and to plot for example the Indicated Torque as shown
19
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
This diagram is not particular smooth. If the indicated torque curve is known then you can factor each diagram using the Factor column. Enter 400[N.m]/T at 3000 rev/min and see the graph change. Set back to 1.0 before proceeding Click on ‘OK’ Save the model using the File menu from the top of the Main Panel 2.1.4.3 Evaluate Solution We now have sufficient data to proceed with an analysis. Click on the ‘Evaluate Solution’ button on the Main Panel
20
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
The panel will appear as shown. No changes are required We can only perform a Determinate solution because we have only built a rigid model Select the Cases tab Click on ‘Select’ button to add an arithmetic speed sweep series of engine speeds to the loadcases
21
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Select the Bearing Model tab Complete the panel as shown
The remaining tabs do not need to be completed for a solution with a rigid block and crankshaft. Select Solve Directly.
The analysis should take a few seconds to run. On completion of the analysis the summary file
.EDSUM will be written. This file contains summary data for the solution. Open this file with an appropriate editor to view the results. Now that the analysis had been completed, the results can be plotted. The next two steps show how to plot the results of the completed analysis. 2.1.4.4 Select Loadcases
22
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Click on the ‘Select Loadcases’ button on the Main Panel
Select the loadcases as shown
This table displays the solution parameters defined for each loadcase Loadcases are grouped by Solution Type and by Loading 2.1.4.5 Plot Results Click on the ‘Plot results’ button on the Main Panel
Select ‘Journal Bearing’ from the Model list and highlight the Subset, Plot and Results.
23
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
Move the panel so that the crankshaft can be seen on the screen as shown As with editing the cranktrain each of the main journal nodes are displayed in green Select main bearings 1 and 2 by dragging the mouse whilst keeping the left button held down Each of the main journal nodes will then turn red Alternatively select the nodes by using the Select All button or by selecting each node individually and using the SHIFT key Click on the Apply button to display the Graph Panel. Use the Page Up and Page Down buttons to view the following graphs
24
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
25
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
26
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
27
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
28
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
29
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis 2.1.4.6 Exporting the data
In addition to the pre defined plots ENGDYN also has the option of exporting formatted data into ASCII files Click on the Export Results button on the left hand side of the software and the following panel appears
As an example, we are going to export the journal bearing load at the main bearing 1. Select JOURNAL_BEARING_LOADS in the dataset and Main bearings on the subset Select the bearing 1 in the 3D viewport Name the ASCII file with something suitable, such as ‘bearing_load’ Change format to block and select only Y & Z directions Click Apply
30
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Two new .EDRES files, for loadcase 1 & 4, were written to your directory. Their content can be reviewed through an appropriate text editor.
2.1.5
Bearing analysis using the Hydrodynamic Solution
Main bearing number 1 will be modified to solve using the Hydrodynamic model 2.1.5.1 Copy the ENGDYN Model to a new filename Use the File menu on the Main Panel to ‘Copy Design’ This will copy all model and loading data 2.1.5.2 Set up Main Bearing 1 as a Hydrodynamic Model Click on the ‘Edit Cranktrain’ button on the Main Panel to display the Edit Cranktrain Panel
Highlight ‘Main Bearing’ as shown
31
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
This will then display each of the main journal nodes in green Select main bearings 1 by dragging the mouse whilst keeping the left button held down or by picking the node. The main journal nodes will turn red Select ‘Edit Selected’ to display the Main Bearing Panel Change the Model Type to Hydrodynamic as shown
32
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis For a Hydrodynamic Model Type it is necessary to use the Mesh and Material tabs to define additional data for this model. In this tutorial we will assume the bearing and journal are circular, and therefore it is not necessary to define a profile using the Profile tab. Select the Mesh tab, and define a mesh 11 x 73 Select the Material tab
This material is used for the boundary lubrication model It is necessary to define the journal material and the material of the bearing lining Click on Define adjacent to Bearing Material, to display the Material Properties Panel. Complete the panel as shown
These data values are typical for a bearing surface. These data are only used by the boundary lubrication model. 33
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Asperity RMS Height h, Density and Asperity radius can be calculated from measured data using the Ricardo MATUTIL program supplied with the Ricardo Software installation which is described in Appendix 12. The Select button is used to select materials that have previously been defined or that are in the SFE file of finite element model defining the bearing. Select OK Click on Define adjacent to Journal Material, to display the Material Properties Panel as shown.
The journal material defaults to the crankshaft material STEEL which was defined in Section 2.1.3.3 Complete the panel as shown The asperity data can be calculated from surface profile measurements using the MATUTIL program supplied with the Ricardo Software installation which is described in Appendix 12.
Select OK and enter the wear and friction coefficients for the bearing surface as shown
34
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
Select OK Click on ‘Define Material’ to display the Crankshaft Material Properties Panel Select OK No edits are required since we have defined the material previously. Click on ‘Calculate Masses’ The message ‘Calculation of mass properties completed successfully, Balance Not Set’ should appear in the message box at the bottom of the Main Panel. Click on ‘Set Balance’ No additional balancing is required because the Lumped Masses include the effect of balancing the crankshaft. Click on ‘Assemble Model’ The message ‘Model assembly completed successfully, Block not assembled’ should appear in the message box at the bottom of the Main Panel Click on ‘OK’ Click on the ‘Dismiss’ button of the Edit Cranktrain panel The message ‘Ready to Perform Solution’ should appear in the message box at the bottom of the Engdyn screen The model should appear as shown.
35
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
Save the model using the File menu from the top of the Main Panel 2.1.5.3 Evaluate Solution We now have sufficient data to proceed with an analysis. Click on the ‘Evaluate Solution’ button on the Main Panel
The Evaluate Solution Panel will appear.
36
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
Set the End Angle to 1440 deg This is equivalent to two and a half cycles and should be sufficient to obtain a converged solution of the main bearing HD solution Select the Cases tab and define a single speed of 3000 rev/min as shown
Select the Bearing Model tab and complete the panel as shown
37
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
The remaining tabs do not need to be completed for this solution since the other parameters are related to indeterminate or dynamic analyses. Select ‘Define Oil Temps’ to display the Bearing Oil Temperatures Panel
We will assume that the temperature of the oil in the bearing is at the inlet temperature. Select OK Select Solve Directly on the Evaluate Solution Panel.
38
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis On completion of the analysis the summary file .EDSUM will be written. This file contains summary data for the solution. Open this file with an appropriate editor to view the results. The solution time for this solution on a Pentium M 1.4 GHz Laptop with 512 Mbytes RAM is approximately 9 minutes. Solution time will be dependent on how heavily loaded the bearing is and whether there is any boundary lubrication occurring. This model at this engine condition has some asperity contact which can be viewed by selecting Contact Pressure when animating (See 2.1.5.5) The memory usage is 67 Mbytes. 2.1.5.4 Select Loadcases Click on the ‘Select Loadcases’ button on the Main Panel
Highlight the only speed selectable (3000 revs/min) and then select OK 2.1.5.5 Animate Results Click on the ‘Animate Results’ button on the Main Panel
Complete the ‘Animation Results’ panel as shown
39
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis
Zoom onto main bearing 1 Hold down both the ‘Ctrl’ button on the keyboard and the middle mouse button whilst dragging the mouse. Change the angle of the view Hold down both the ‘Shift’ button on the keyboard and the middle mouse button whilst dragging the mouse. Click on the ‘Play’ button of the animation control toolbar to animate the results as shown.
40
KNOWLEDGE CENTRE Tutorials Tutorial 1: An Introduction to ENGDYN – Concept Crankshaft Analysis Experiment with looking at different results by selecting various parameters from the Animation Results Panel Close the animation panels and exit the program.
41
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
2.2 Tutorial 2: Crankshaft Dynamic Analysis 2.2.1
Introduction
Objective: To build an ENGDYN model suitable for performing a dynamic analysis of a crankshaft. Introducing FE models of the crankshaft The results of this tutorial are used as input to the NVH and Crankshaft Stress tutorials Items Covered: Building the engine model Introducing an FE model of the crankshaft Matrix Reduction Dynamic solution Estimated duration: 0.5 day (model preparation) 1 day (overall including performing solutions) Engine: Inline 4 gasoline engine 75.0 x 44.75 mm 1.6 litre Required Files: ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\il4_crank.inp ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\il4_2000.PRES ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\il4_3000.PRES ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\il4_4500.PRES ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\il4_5000.PRES ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\il4_5500.PRES ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\il4_6000.PRES Finite element model requirements:
43
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis The crankshaft FE model should be modelled as described in the “Geometry of the Finite Element Model” Section (Section B.2.4.5.2 in the “Help” part of the manual”)
2.2.2
Getting Started
Copy the necessary files from the example directory to a working directory and ensure that you have write permissions for all the files. Start engdyn On Unix or Linux platforms simply type engdyn On windows click on the shortcut, otherwise go to Start>Programs>Ricardo Software>2014.1>Mechanical Suite>ENGDYN>ENGDYN
2.2.3
Building the Model
2.2.3.1 Configure Engine In order to build an ENGDYN model it is first necessary to define the major dimensions and features of the engine. This is done using the Engine Configuration Panel. Select ‘Configure Engine’ from the buttons on the left side of the Main Panel.
Complete the Engine Configuration Panel as shown.
44
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
The Firing Order can be selected using the Options button or supplied as a hyphen separated list as shown Journal lengths are the lengths of the journals including the journal fillets NOT the bearing lengths which are defined later All data in this example are uniform. In cases where data are non-uniform is selected. This will display in red on the 2D graphic those items that can be edited by selecting with the mouse. Items that are not shown are edited using the Edit button. On completion select ‘Apply’ which will display
Select OK. The crankshaft model will appear in the Main Panel as shown
45
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
Note that the dimensions of the cranknose and flywheel assemblies on screen are only representations as no data has yet been entered for these. These are defined in 2.2.3.3. 2.2.3.2 Define Models An ENGDYN model consists of a number of sub-models which are defined using the Model Definitions Panel. Select ‘Define Models’ from the buttons on the left side of the Main Panel.
Complete the model definition for the Crankshaft
46
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
Use default values for Model, Data and Matrix Formulation The crankshaft in this example is already in the same co-ordinate system and orientation as defined in the 2.2.3.1. No transformation is therefore necessary Click on Browse button to display the ‘Model Translation’ Panel as shown
This panel is used to translate FE models to Ricardo SFE format Models can alternatively be translated outside the graphical interface using the appropriate FEARCE translator. If the user does this then Origin is left as Ricardo-SFE. The crankshaft in this example is already in the same co-ordinate system and orientation as defined in the 2.2.3.1. No transformation is therefore necessary
>
47
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis Set origin to ABAQUS and Units to mm Units is the units of the ABAQUS model During the translation the model will be translated so that the SFE file is in SI units. Use the Browse button to select the file il4_crank.inp from the working directory
The program will automatically set the output name to il4_crank.SFE although this can be changed Click Translate to translate the model from the ABAQUS format into SFE format Select OK to close the panel. On completion the message Translation successful will appear at the bottom of the Main Panel. Then there will be a SFE file named il4_crank.SFE in the current directory. Go to the current directory, find the file il4_crank.SFE Double-click on this file and the RICARDO FE viewing and interface tool (RDESK) will launch showing our crankshaft FE model
Note that on UNIX/LINUX systems, type rdesk into the shell at the working directory.
48
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis The R-Desk viewport can be controlled in a similar way to ENGDYN’s. To move the model in the panel, Middle mouse and move = translation + middle mouse and move = rotation middle mouse scroll wheel = zoom in/out Select the Cylinder Block tab and define the Model as Rigid with User Defined data as shown.
Click on OK. This will display the model as shown.
By default only the edges of the FE model of the crankshaft will be displayed as shown. To change the appearance of the model select the Model Appearance Panel from the View Menu. Note that the webs and journals align with the reduced model (green), but that the cranknose and flywheel assemblies are not yet correctly defined. This will be addressed in the next step
49
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis 2.2.3.3 Editing the Cranktrain The next step is to define the data related to the cranktrain. Select ‘Edit Cranktrain’ from the buttons on the left side of the Main Panel to display the Cranktrain Tool Panel
Given we defined the crankshaft as Dynamic (See 2.2.3.2) it will be necessary to define everything listed on the left hand side of the panel except Lumped Mass and Mechanical Links Lumped Mass is used to define any additional masses that are not included in the FE model or defined using Flywheel Assembly or Cranknose Assembly. Mechanical Links is used to define links for co-simulation. Highlight ‘Crank Web’ as shown This will then display each of the crank web numbers and elements in green. Select webs 1, 2, 7 and 8 by dragging the mouse with the left mouse button held down. Each of the selected webs will then turn red as shown
50
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
Select ‘Edit Selected’ to display the Web Panel Set Counterweight to Present and complete the counterweight data as shown
51
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
The counterweight geometry is defined only for the purposes of balancing the crankshaft in 2.2.3.4 NOT for defining its mass properties. The angles are with respect to the Engdyn coordinate system not with respect to the web. The journal and thickness data has been calculated from the data supplied in 2.2.3.1. The stiffness data is calculated in 2.2.3.5 Select OK Select webs 3, 4, 5 and 6 by dragging the mouse with the left mouse button held down. Select ‘Edit Selected’ to display the Web Panel Set Counterweight to Present and complete the counterweight data as shown
52
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
Select OK Highlight Pin Journal Element This will then display each of the pin journal elements in green Select ‘Select All’ Each of the pin journal elements will then turn red Select ‘Edit Selected’ to display the Element Panel as shown
53
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis This panel allows the internal diameter of a hollow journal to be defined. The pin journals of this crankshaft are solid. The Outside Diameter was defined during 2.2.3.1. Select OK Highlight Main Journal Element This will then display each of the main journal elements in green Select ‘Select All’ Each of the main journal elements will then turn red Select ‘Edit Selected’ to display the Element Panel as shown
As with the pin journals this panel allows the internal diameter of a hollow journal to be defined. The main journals of this crankshaft are solid. The Outside Diameter was defined during 2.2.3.1. Select OK Highlight Big End Bearing This will then display each of the cylinder numbers and pin journals in green. Select ‘Select All’ Each of the pin journals will then turn red Select ‘Edit Selected’ to display the Bearing Panel Complete the Bearing Panel as shown
54
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
Given the Model Type is Mobility and the bearings are plain it is not necessary to enter any other data than shown. These bearing are feed from the journal via a feed from the adjacent main bearings. The oil hole angular position can be defined either using the Position column (as in this case) or using the Height column. The height is calculated from the angle and vice-versa. The remaining tabs (Mesh, Material, Profile) need not and cannot be edited because the model type is ‘Mobility’. These tabs are only required with Hydrodynamic or Elastohydrodynamic models Select OK Highlight Main Bearing This will then display each of the bearing numbers and main journals in green. Select ‘Select All’ Each of the main journals will then turn red Select ‘Edit Selected’ to display the Main Bearing Panel Complete the Main Bearing Panel as shown
55
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
Given the Model Type is Mobility and the bearings are partial grooved it is not necessary to enter any other data than shown. Any journal oil holes of a grooved bearing are not currently considered. The remaining tabs (Mesh, Material, Profile) need not and cannot be edited because the model type is ‘Mobility’. These tabs are only required with Hydrodynamic or Elastohydrodynamic models Select OK Highlight Thrust Bearing This will then display the single main journal node in red corresponding to the thrust bearing defined in 2.2.3.1. There is therefore no need to ‘Select All’ Select Edit Selected to display the Thrust Bearing Panel Complete the panel as shown
56
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
Supply geometric extent of the axial thrust face and the bearing clearance ENGDYN will calculate appropriate stiffness/damping based upon the system material properties Select OK Highlight ‘Cylinder’ Click on ‘Select All’ to select all pin journal nodes Click on ‘Edit Selected’ Complete the ‘Cylinder’ panel as shown.
The height Head Height is used to derive a node for the cylinder block model in 2.2.3.6. Select ‘OK’ Highlight ‘Connecting Rod’ Click on ‘Select All’ to select all pin journal nodes Click on ‘Edit Selected’ Complete the ‘Conrod’ panel as shown.
57
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis Note the units for the connecting rod inertia! Select ‘OK’ Highlight ‘Piston’ Click on ‘Select All’ to select all pin journal nodes Click on ‘Edit Selected’ Complete the ‘Piston’ panel as shown.
The skirt and top ring data has been set to 0 since this data are only required for block stress analysis and is therefore not required for this tutorial. Select OK Highlight Cranknose Assembly This will then display the cranknose assembly in red. Therefore there is no need to ‘Select All’ Select Edit Selected to display the Cranknose Assembly Panel Complete the panel as shown
58
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
The Element Length effectively positions the first node of the reduced model of the crankshaft. The position of this node must correspond to a plane of nodes in the FE model. The natural frequency of the vibration damper is 123.3 Hz. The seismic mass of the damper defined here must not be included in the FE model of the crankshaft The hub of the vibration damper may be included in the FE model or as in this case defined using this panel. Any other additional masses, for example gears, sprockets and pulleys, not included in the FE model can be defined by selecting Lumped Mass on the Edit Cranktrain Panel. For the purposes of this exercise we will assume there is no additional masses. Select ‘Edit’ button adjacent to Hub Data to display the Lumped Mass Panel Complete the Panel as shown.
59
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
This defines the mass and inertia of the hub. Given X offset is 0 this lumped mass is assumed to be at the node (defined by the element length on the Cranknose Assembly Panel) Click on OK Click on OK on the Cranknose Assembly Panel Highlight Flywheel Assembly This will then display the flywheel assembly in red. Therefore there is no need to ‘Select All’ Select Edit Selected to display the Flywheel Assembly Panel Complete the panel as shown
60
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
The distance X1 effectively positions the last node of the reduced model of the crankshaft. The position of this node must correspond to a plane of nodes in the FE model. In this example, by setting Type to Conventional, we are assuming the flywheel to be rigid. The Offset is with respect to the node defined by X1. The inertia data is with respect to its centroid. Click on the Clutch tab and enter the mass and inertia data for the clutch.
The other tabs to not apply to a conventional flywheel assembly Select OK There are no additional lumped masses to be defined or mechanical links to be defined. Select Define Material to display the Crankshaft Material Properties Panel as shown
61
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
The material data will be read from il4_crank.SFE. Use select to list the materials in the SFE file and select the material STEEL from the list as shown.
This table lists all the material in the SFE file. The selected material is the material of the web overlap section. Materials used by ENGDYN and those stored in the .SFE are identified by a single unique name. Select OK and the material is added to the panel as shown
These data may be edited if required Select OK These data are now stored by ENGDYN and also written to the .SFE file. Previous data will be overwritten. The message “Mass Not Calculated’ will be written to the Main Panel.
62
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis 2.2.3.4 Calculating the Mass Properties and Balancing the Crankshaft Click on ‘Calculate Masses’ to calculate lumped mass properties at each reduced node The progress bar as shown will be displayed
The program derives an element set (as described 2.4.5.3 of the manual) for each reduced node. These sets can be viewed using FEVIEWER. The mass and inertia of each set is then assigned to the appropriate node. On completion the message ‘Balance Not Set’ will be written to the Main Panel. Click on ‘Set Balance’ to display the Primary Balance Panel as shown
63
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis This panel shows the primary imbalance of the model as given by ‘As modelled’. This is the imbalance of the FE model and any additional lumped masses. The imbalance will invariably be dependent on the accuracy of the finite element model. For an in-line 4 engine we would expect the crankshaft to have primary balance. This panel allows the user to simulate balancing the crankshaft by drilling the counterweights. Set Calculation to ‘Using Web data’ Use select adjacent to Web Numbers to display the Select Web Panel as shown
This displays all webs that have counterweights as defined in 2.2.3.3 Select webs 3, 4, 5, and 6 as shown
Select OK Use Define adjacent to Material Name to display the Counterweight Material Properties as shown
64
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
This by default will display the crankshaft material properties as defined in 2.2.3.3. Select can be used to select a different material for the counterweight material. Click on OK. The Primary Balance Panel will appear as shown.
For an in-line 4 engine the ‘Required’ balance is 0 Click on Apply to balance the crankshaft The message ‘Balance calculation completed successfully’ will be written to the Main Panel. Click on OK The message ‘Balance Set, Stiffness Not Calculated’ will be written to the Main Panel.
65
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis Scroll up on the messages on the Main Panel to see the output related to the balance drillings.
2.2.3.5 Matrix Reduction Click on ‘Calculate Stiffnesses’ to calculate the crank web and element stiffness using the finite element method. 2 solvers are available. The default Vectorized Sparse Solver (VSS) or the Symmetric Conjugent Gradient Click on ‘No’ to select the VSS solver The program splits the FE model of the crankshaft into a number of submodels. FE calculations are performed on each model. These calculations take circa 5 minutes on am Intel I7-2620M 2.70 GHz Laptop with 4 GB RAM. On completion of the finite element solutions the reduced mass and stiffness matrices of the crankshaft are derived and an eigenvalue solution performed. The message ‘Stiffness Calculated, Model Not Assembled’ will be written to the Main Panel on completion. Inspect the contents of your working directory to see the files generated. Click on ‘Calculate Stiffnesses’ again This time a query panel will be displayed as shown
The program checks whether finite element solutions for each web and element already exist. Click on No No finite element solutions will be performed The reduced mass and stiffness matrices of the crankshaft are derived and an eigenvalue solution performed. Click on Assemble Model This verifies the validity of the data The message ‘Model Assembly completed successfully. Block not assembled’ will be written to the Main Panel if assembly was successful. Click on Dismiss on the Cranktrain Tool Panel. The main panel will now appear as shown
66
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
2.2.3.6 Assembling the Cylinder Block Model The next step is to assembly the cylinder block model. Select ‘Edit Block’ from the buttons on the left side of the Main Panel to display the Cylinder Block Tool Panel The Main Panel will appear as shown
67
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
The program creates nodes at each main bearing and at the reversal positions of the small end of the connecting rod (Thrust and Anti-Thrust) and at the centre of the cylinder head for each cylinder. These nodal positions are derived from data defined in 2.2.3.1 and 2.2.3.3. Simply click on Assemble Model. Click on OK The model is now completed and should appear as shown
68
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis Save the model using the File menu from the top of Main Panel
2.2.4
Solution
2.2.4.1 Define Lubricant Properties Click on ‘Lubrication’ button from the buttons on the left hand side of the Main Panel
Use the Browse button to select the lubricant SAE5W30 from the database By default the program will initially select the database directory at ..\Ricardo\2014.1\Common\Materials\LubricantsThis database contains the most common lubricants Use either Add or Update to add the lubricant
Click on OK
2.2.4.2 Define Loading Conditions The cylinder pressure diagram and any additional loadings (for example loads from Valdyn and gravity forces) are entered in this step. Click on the ‘Loading’ button on the Main Panel
69
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis This will display the Loading Definition Panel as shown
A number of different loading maps can be defined, Full Load, Part Load and No Load. The solver will interpolate at speeds between those defined using this panel. Type in a speed of 2000 rev/min Position the mouse over the File Name column and use the right button to display the pop-up menu. Use ‘Select Pressure file’ to select the file il4_2000.PRES from your working directory. The panel will appear as shown.
You may wish to remove the pathname in front of the file. You may wish to inspect and understand the selected file using an appropriate editor Define the Ambient and Crankcase Pressure as 1.0 bar
70
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis The crankcase pressure is only currently used for Hydrodynamic and Elastohydrodynamic bearing models to define the boundary condition at the edge of the bearing Position the mouse over the column number and with the right button depressed display the Row Operations pop-up. Select Insert After. The Panel will appear as shown
Define the pressure diagram at 3000 rev/min as shown
Complete the panel as shown.
71
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
The remaining tabs, Force Profile, Force Equation and Distortion need not be completed for this tutorial. Use the Plot button to display the applied pressure loading as shown
Save the model using the File menu from the top of the Main Panel We now have sufficient data to proceed with an analysis 2.2.4.3 Evaluate Dynamic Solution Click on the ‘Evaluate Solution’ button on the Main Panel
72
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis The Evaluate Solution Panel will appear. Set the solution Type to Dynamic and edit the solution tolerances as shown
Select the Cases tab and click on Select to define a speed sweep at full load between 2000 and 6000 rev/min in steps of 250 rev/min as shown
73
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
Click on OK to populate the table on the Evaluate Solution Panel as shown
Select the Model Options tab and change the Cylinder Damping to 500 N.s/m and Frequency to 50 Hz as shown.
74
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
The cylinder damping parameter is used to tune the torsional response of the crankshaft The coupling is a torsional spring and damper to prevent rigid motion of the crankshaft about its axis. This is analogous to the stiffness and damping of a dynamometer connected to the crankshaft. The damping value of the coupling can be used to tune the torsional response of the crankshaft. The natural frequency of the coupling must be significantly lower than the first flexible torsion mode of the crankshaft. Select the Bearing Model tab and complete the panel as shown
75
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
Select ‘Define’ to display the Bearing Oil Temperatures Panel
We will assume that the temperature of the oil in the bearing is at the inlet temperature Select OK Select Solve Directly on the Evaluate Solution Panel.
76
KNOWLEDGE CENTRE Tutorials Tutorial 2: Crankshaft Dynamic Analysis
On completion of the analysis the summary file .EDSUM will be written. This file contains summary data for the solution. Open this file with an appropriate editor to view the results
77
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
2.3 Tutorial 3: NVH Analysis 2.3.1
Overview
Objective: To use ENGDYN to perform vibration and acoustic analysis on a powertrain The analysis uses finite element (FE) models of both the powertrain assembly and the crankshaft The tutorial requires the ENGDYN model from the Crankshaft Dynamic Crankshaft Tutorial Items Covered: Dynamic cylinder block model Engine mount modelling Component Mode Synthesis (CMS) matrix reduction Vibration Analysis Acoustic Analysis Rayleigh Integral Method Indirect Boundary Element Method (BEM) Estimated duration: 0.5 day (model preparation) 1 day (overall including performing solutions) Engine: 1.6 litre Inline 4 gasoline engine 75.0 x 44.75 mm Required Files: ..\Ricardo\2014.1\Products\ENGDYN/Tutorials\il4\IL4_BLOCK.SFE ..\Ricardo\2014.1\Products\ENGDYN/Tutorials\il4\il4_block_noise.VPT ..\Ricardo\2014.1\Products\ENGDYN/Tutorials\il4\il4_block_vibration.VPT ..\Ricardo\2014.1\Products\ENGDYN/Tutorials\il4\IL4_BLOCK_BoundaryElementM odel.SFE The .EDSF file from Crankshaft Dynamic Analysis Tutorial, so all the files used in tutorial 2 are needed here. Finite element model requirements: For dynamic and acoustic analysis of the powertrain it is necessary to have a fully dressed finite element model with engine mount brackets to support the engine
79
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis mounts. The model used in this tutorial is not representative but is sufficient for demonstrating the analysis steps.
2.3.2
Getting Started
Go to the working directory in which you completed the Crankshaft Dynamic Analysis Tutorial. Copy the file IL4_BLOCK.SFE from the Tutorials directory in the installation to this working directory and ensure that you have write permissions for this file. Start engdyn On Unix or Linux platforms simply type engdyn On windows click on the shortcut, otherwise go to Start>Programs>Ricardo Software>2014.1>Mechanical Suite>ENGDYN>ENGDYN In the Crankshaft Dynamic Analysis Tutorial the cylinder block (or powertrain) was modelled as rigid. Obviously for analysing the vibration and acoustic behaviour of the powertrain it is necessary to model the powertrain as dynamic. This section describes the steps in updating the model derived in the previous tutorial. Open the .EDSF file used for the Crankshaft Dynamic Analysis Tutorial using the File menu at the top of the Main Panel From the same menu ‘Copy Design’ to a new file
2.3.3
Updating the Model
2.3.3.1 Define Models In this step we need to redefine the cylinder block model as dynamic rather than rigid. Select ‘Define Models’ from the buttons on the left side of the Main Panel to display the Model Definitions Panel.
Select the ‘Cylinder Block’ tab and define the model as Dynamic as shown.
80
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Click on Browse button to display the ‘Model Translation’ Panel as shown
Use the Browse button to select the file IL4_BLOCK.SFE from the working directory.
The program will automatically set the output name to IL4_BLOCK.SFE although this can be changed Click on the Transformation tab to define the translation vector as shown to transform the finite element model of the cylinder block.
81
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Click on OK to dismiss the panel In this case no translation is performed since the model has already been translated into SFE format using the appropriate translator. Click on OK Note how the previously defined reduced model of the cylinder block has disappeared. By default the finite element model of the cylinder block is not displayed. The model may be displayed using the Model Appearance Panel from the View Menu at the top of the Main Panel. 2.3.3.2 Reassembling the Crankshaft Model The program forces you to reassemble the crankshaft model because the cylinder block model has changed. Select ‘Edit Cranktrain’ from the buttons on the left side of the Main Panel to display the Cranktrain Tool Panel
82
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis We do not need to redefine any data. However the program will force us to re-edit the main and thrust bearing data. Click on Assemble Model An error message will be displayed in the pop-up as shown.
Dismiss the error message Highlight ‘Main Bearing’, click on ‘Select All’ followed by ‘Edit Selected’ to display the Main Bearing Panel. Simply select OK since we do not want to re-edit the data. Click on Assemble Model again An error message will be displayed in the pop-up as shown.
Dismiss the error message Highlight ‘Thrust Bearing’ and click on ‘Edit Selected’ to display the Thrust Bearing Panel Simply select OK since we have no data to edit. This is an error in the current release of the program since it should not force you to edit the thrust bearing since the stiffness is later calculated from the cylinder block model. Click on Define Material and OK the panel. Follow Step 4 and Step 5 of the Crankshaft Dynamic Analysis Tutorial On completion the message ‘Stiffness Calculated, Block Not Assembled’ should be displayed in the message area of the Main Panel.
2.3.3.3 Edit Cylinder Block - Model Definition
83
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis It is necessary to define a reduced model of the cylinder block. The reduced model is a number of nodes and degrees of freedom that are a subset of the complete model. We need to define a number of sets. A constrained node set for each main bearing. For this set the bearing is defined by a single node with 6 degrees of freedom whose movement is the average of the nodes on the bearing surface. Four node sets for each cylinder to define the upper and lower reversal points of the small end of the connecting rod on the thrust and anti-thrust sides of the cylinder bore. Each node set will contain a single structural node. A constrained node set for the head gas face of each cylinder. For this set the gas face is defined by a single node with 3 translational degrees of freedom whose movement is the average of the nodes on the gas face. A constrained node set that defines the attachment point of each engine mount. For this set the mount is defined by a single node with 3 translational degrees of freedom whose movement is the average of the nodes on the bracket. Select ‘Edit Block’ from the buttons on the left side of the Main Panel to display the Cylinder Block Tool Panel The following Query Panels will be displayed in succession.
These are displayed because we have already defined a cylinder block model and in doing so have defined a number of sets Click on Yes in each case In most cases you want to answer No to the question. You only want to answer Yes if either the block model type has changed (as in this case) or the main bearing or cylinder geometry has changed. Messages will appear in the bottom of the Main Panel and the Cylinder Block Tool Panel will appear as shown
84
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
In this action the program has transformed the finite element model of the cylinder block, calculated its mass properties and attempted to automatically define the sets of the reduced model using the data supplied previously The reduced model is shown in orange. The sets defining each main bearing are incomplete and nothing has been defined for the cylinders. This is because the nodes are outside the geometric tolerance. It is necessary to define these sets using the Define Model Panel. Click on ‘Define Model’ to display the Define Model Panel as shown
Consider firstly the constrained node sets defining each main bearing Change Type to Constrained Node and Name to Main Bearing
85
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis The Cyl and Val fields will now be ghosted and the ID field will be editable. Click on ‘Select’ adjacent to the ID field and select 1 from the list so that the panel appears as shown
ID 1 denotes main bearing 1 Alternatively you can simply type 1 in the field Select the Definition tab
Select each tab and understand the default values Each set is ‘clipped’ based on geometric shape Each shape is defined by a Centre, Axis, Extent and Diameter Select the Tolerance tab
86
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
The default linear tolerance is set to half the minimum distance between any two adjacent finite element nodes. Each set has its own tolerance values Change the Linear Tolerance to 0.5 mm and press ‘Clip Set’.
The ‘clipped’ mesh is shown in pink The program only clips the external surfaces of the model. This set is OK so click on ‘Add Set’ to add the set to the reduced model as shown
87
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Repeat this procedure for each of the remaining main bearing by changing the ID on the Name tab and the linear tolerance on the Tolerance tab. On completion the panel should appear as shown.
Secondly consider the node sets defining the thrust and anti-thrust nodes for each cylinder.
88
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis Select the Name tab Change Type to Node and Name to cylinder as shown
Click on ‘Select’ adjacent to the Cyl field and select all the cylinders Click on ‘Select’ adjacent to the ID field and select lowerAntiThrust so that the panel appears as shown
Select the Definition tab and again understand the default data The data for cylinder, valve seat and valve spring are with respect to a cylinder local coordinate system with the origin at the top of the cylinder. This enables sets to be defined for multiple cylinders in a single action. The shape is set to Sphere in this case to find the nearest relative to the specified centre. Select the Tolerance tab and again change the linear tolerance to 0.5 mm and click on ‘Clip Set’
89
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
The node numbers are displayed The coordinates of each node are displayed in the message area with respect to the cylinder and the global system. Select ‘Add Set’ to add the nodes to the reduced model. Repeat this procedure for each of the remaining cylinder nodes lowerThrust, upperAntiThrust and upperThrust by changing the ID on the Name tab and the linear tolerance on the Tolerance tab. On completion the panel should appear as shown.
90
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Finally add sets defining the attachment points of the engine mounts. Again these sets will be from the SFE file. In this example there are no engine mount brackets. Arbitrary sets have been defined to attach the engine mounts Select the engine mount sets from the list and click on Read Set.
91
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Select ‘Add Set to obtain the following
92
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Use the mouse and cursor keys to rotate and zoom the model if you haven’t already done so. OK the Define Model Panel 2.3.3.4 Edit Cylinder Block - Engine Mount Definition Highlight Elastomeric Mounts in the list on the Cylinder Block Tool Panel This will then display each of the engine mount nodes in green Select the front engine mount using the mouse The node will turn red Select ‘Edit Selected’ to display the Elastomeric Engine Mount Panel Complete the properties for the three degrees of freedom as shown
The stiffness’ are with respect to the local coordinate system The default orientation is with X vertical. These stiffness values are typical values for an automotive application. However the rigid body modes using these values will be higher than you would normally expect given the mass of the block assembly in this case is lower given it is not a fully dressed model. Click on OK Repeat for the right engine mount
93
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Repeat for the left engine mount
Repeat for the rear engine mount
94
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
2.3.3.5 Edit Cylinder Block – Model Reduction The model is reduced using Component Mode Synthesis. This technique ensures that the modal behaviour is the same as the complete model for a specified number of modes. The model is run free-free to give 6 rigid body modes which are latter restrained during model assembly by the engine mounts. Select ‘Matrix Reduction’ to display the Matrix Reduction Panel
For this tutorial we will use the FEARCE Vectorized Sparse Solver (VSS) Alternatively MSC Nastran can be used If you have limited memory then this option can be used to limit the memory during forward elimination and backsubstitution. This will result in a longer solution time. Set the Number of Modes to 50 For this example this will obtain modes up to 3.5 kHz. The number of modes will determine the size of reduced matrix and therefore the speed of the subsequent dynamic solution. In this case the model has 102 reduced freedoms, so the size of the reduced matrices will be 152 square.
Select the Default output name The output name specifies the name of the FEARCE .FRC command file that will be written. Selecting the Default button will give the output name the same name as the .SFE file
Switch on the Limited Memory toggle and specify a memory requirement of 1000 Mb or less if your machine has less available memory
95
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Select Solve. The following Query Panel will be displayed.
Select OK. Another Query Panel will be displayed.
Select ‘Direct’. This will execute the command file and start the solution. A Progress Panel will be displayed. This model has approximately 190000 degrees of freedom which are reduced to 102 plus 50 modes for the reduced model to be used by ENGDYN. On completion the message ‘Run completed successfully’ is written to the Main Panel and the following Information Panel is displayed.
96
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis 2.3.3.6 Edit Cylinder Block – Model Assembly Select ‘Assemble Model’
This will assemble the cylinder block model An eigenvalue solution of the reduced model is performed and the axial stiffness of the thrust bearing housing is calculated Click on OK to dismiss the Cylinder Block Tool Panel Save the model using the File menu from the top of the Main Panel
2.3.4
Solution
Repeat the dynamic analysis as described in the Crankshaft Dynamic Analysis Tutorial (Section 2.2.4.3) but for a single speed 2000 rev/min and with the block set to Dynamic.
97
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis 2.3.5 Vibration Analysis The dynamic analysis will predict the forces and displacements (if the block is defined as dynamic) at each of the nodes of the reduced model. This section describes how we can backsubstitute using a forced response analysis using the predicted forces at the reduced nodes to predict the displacements at any point on the cylinder block model. We will predict the response on one side of the block and at each engine mount at a single speed up to 1 kHz. 2.3.5.1 Selecting the Required Loadcase Click on ‘Select Loadcases’ from the buttons on the left side of the Main Panel. Select the loadcase at 2000 rev/min full load as shown
Click on OK to dismiss the panel
2.3.5.2 Performing the Forced Response Analysis Click on ‘Block Analysis’ from the buttons on the left side of the Main Panel to display the Cylinder Block FEA Panel
Click on ‘Select Output’ to display the Cylinder Block FEA Output Panel Highlight ‘Nodal Vibrations’ from the list as shown
98
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Frequency Response Loading will automatically be highlighted. Click on OK to dismiss the panel. Click on ‘Select Model’ to display the Cylinder Block Model Panel
The default will be the FE model used to derive the reduced model. This defines the model we wish to excite. The selected file must have modal data and must have reduced nodes at the same location as the reduced model. Click on OK to dismiss the panel. Click on ‘Define Output’ to display the Frequency Response and Acoustic Analysis Panel. Set the Maximum Frequency to 1 kHz as shown and click on Select to list the node sets in the model.
99
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis This will calculate the response at all the harmonic frequencies up to 1kHz for each selected engine speed. In this case, given we have selected a case at 2000 rev/min and this is a 4 stroke engine we will have results at 16.67, 33.33, 50.0 …Hz Click on ‘Select’ and select the sets LEFTSIDE, engineMount:ID_left, engineMount:ID_right, engineMount:ID_front and engineMount:ID_rear from the node set list as shown
The set LEFTSIDE was defined using FEVIEWER Select the Output tab Change Export to “Sets Only” and change the default output name to that shown
The output name is the base name of all the output files.
100
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis Setting Export to Sets will create an SFE file that contains only the geometry of LEFTSIDE and each of the engine mounts and the reduced nodes. The results will also be written to this file.. Click on OK to dismiss the panel. Click on ‘Generate Output’ to perform the frequency response analysis A message report will be displayed giving the output from the FEARCE program.
On completion the message ‘Run completed successfully’ is written to the Message Panel and the following Information Panel is displayed
OK this panel and the Message Panel On completion the following files will be written IL4_BLOCK_VIBRATION.FRC command file IL4_BLOCK_VIBRATION.FRL containing the Frequency Response Loading LC001.ALPHA containing the Modal Contribution factors for load case 1 LC001.TOTALS containing the total loads applied to the model for load case 1 IL4_BLOCK_VIBRATION_SUBSET.SFE containing the results for LEFTSIDE and each of the engine mounts and all the reduced nodes. 2.3.5.3 Plotting Results Go in the working directory and double-click on IL4_BLOCK_VIBRATION_SUBSET.SFE in order to open the .SFE file in R-Desk. The model should appear as shown in the screen shot below.
101
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Go to the FE Graph tab as shown on the screen shot below:
102
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
103
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis In the FE Graph section, click on New -> NVH -> Vibration graph, and the following panel will apppear on the right.
Select Displacement in the left panel, as shown below:
104
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Then select the mount sets in the Edit panel on the right (use CTRL + left mouse), and then click “Finish” at the bottom at the edit panel.
105
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Finally, click “Plot” in the left hand panel, and you will obtain the R-Plot Graph shown below.
106
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
You can explore the process with other quantities, for example velocity. You can plot several plot at the same time, by highlighting the wanted plots simultaneously (CTRL + left mouse), and the click on “Plot”.
107
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
2.3.5.4 Viewing Results Go to the “Data” tab at the bottom of the left panel as shown in the screen shot below.
108
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
109
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
In the “Deformation” tab, select “Displacement” for a frequency of 133.3 Hz, as shown below. This corresponds to the 4th order.
Select “Displacement” in the “Contour” tab. The screen will then look as shown below.
110
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Go to the Animation tab that can be selected at the bottom of the middle panel, and animate the results by clicking on the green run bottom.
Experiment with other quantities, like for example velocity and acceleration. To rotate: SHIFT + move mouse with middle bottom pressed To zoom in/out: CTRL + move mouse with middle bottom pressed
111
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
2.3.6
Noise Prediction
This section describes how we perform an Acoustic Analysis using the Rayleigh Integral Method. The Rayleigh calculation uses a simple method treating each node as a piston in a baffle assuming that the vibration is the same over the surface of the piston. This is normally acceptable for models where the mesh density is sufficient to give a good representation of the mode shapes, and the wavelength is large compared to the element size. We will predict the radiated sound power and sound intensity from each side of the block at a single speed up to 3 kHz. 2.3.6.1 Performing the Acoustic Analysis Return to the Cylinder Block FEA Panel Click on ‘Select Output’ to display the Cylinder Block FEA Output Panel Highlight ‘Radiated Noise’ from the list as shown
Frequency Response Loading and Nodal Vibrations will be automatically highlighted. Click on OK to dismiss the panel. Click on ‘Select Model’ to display the Cylinder Block Model Panel
Click on OK to dismiss the panel. Click on ‘Define Output’ to display the Frequency Response and Acoustic Analysis Panel.
112
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis Set the Maximum Frequency to 3 kHz as shown
Select the Acoustic tab and Options tab and switch off Delete Vibrations as shown.
If we are performing this calculation for a large number of speeds and frequency it is usually prudent to leave Delete Vibrations on in order to minimise the size of the resultant files. Click on the Boundary tab and ‘Radiating Sets’ and select the sets LEFTSIDE and RIGHTSIDE from the face set list as shown
113
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
These sets were generated using R-DESK VIEWER. Select the Output and change the output name to that shown
Click on OK to dismiss the panel. Click on ‘Generate Output’ to perform the frequency response and acoustic analysis A message panel will be displayed giving the output from the FEARCE program.
114
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
On completion the message ‘Run completed successfully’ is written to the Message Panel and the following Information Panel is displayed
OK this panel and the Message Panel On completion the following files will be written IL4_BLOCK_NOISE.FRC command file IL4_BLOCK_NOISE.FRL containing the Frequency Response Loading LC001.ALPHA containing the Modal Contribution factors for load case 1 LC001.TOTALS containing the total loads applied to the model for load case 1 IL4_BLOCK_NOISE.L001.RES is an ASCII file containing the sound power for load case 1 IL4_BLOCK_NOISE_SUBSET.SFE containing the results for LEFTSIDE and RIGHTSIDE. The size of this file is 63.5 Mbytes.
2.3.6.2 Viewing and plotting Results Open in R-Desk Viewer IL_BLOCK_NOISE.SFE. The process for viewing the results is the same as that described in the Vibration Analysis section (2.3.5). 1. In the Deformation Tab, select Displacement 2. Select “All” 3. In the Frequency tab, select Case 1, and then 1016.67 Hz 4. In the Contour Tab (Data panel on the left), select Sound Intensity.
115
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis 5. Select dBA below Your screen should look similar to the screen shot below.
Note that you can modify the background colour from the tools many in the top left corner. Click
Then Select Options -> Preferences -> Background Colours. You can animate the results using the Animation Tab (at the bottom of the middle panel) In order to plot the results, use the FE Graph tool similarly to what was done in the Vibration analysis section. 1. In the FE Graph panel, when you click “New” and the select “NVH”, chose the “Acoustic Graph” option. 2. In the “Edit” panel on the right, chose “Sound Power – Spectra – XY” 3. Select the desired sets (for example LEFTSIDE and RIGHTSIDE) 4. Click finish 5. Finally click “Plot” in the FE Graph panel. You will obtain a plot similar to the one below.
116
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
117
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
2.3.7
Noise Prediction Using Boundary Element Method
This section describes how we perform an Acoustic Analysis using the indirect Boundary Element Method. This method is a more rigorous method and calculates the noise generated by a closed volume. The need for the solution of a set of simultaneous equations, of order N (where N is the number of boundary elements), at each frequency, means that the BEM is substantially slower than the Rayleigh method. 2.3.7.1 Performing the Acoustic Analysis Using BEM Return to the Cylinder Block FEA Panel Click on ‘Define Output’, select Acoustic tab and choose indirect BEM in the Method field
Choose the Boundary tab, Click on Browse button to load SFE file named IL4_BLOCK_BoundaryElementModel.SFE
Choose Acoustic Mesh tab and select Defined in the Mesh field
118
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Click on the Define button that has become active, with the Acoustic Mesh panel
Enter a suitable filename for the new mesh, such as Acoustic_Mesh, and select Sphere in the Shape field
119
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis We shall define the centre of the sphere at the centre of the gravity of the engine model. We can calculate the centre of gravity of the model using the viewer tool Double click on the IL4_BLOCK.SFE model. Right click the Plot item under the model name in the Plots panel Select Calculate from options list, then select Centre of Gravity
The result will be displayed in the ‘Information’ panel at the bottom of the canvas. The values are returned in SI units and are (176.19, 109.85, 2.81) mm. We have to account for the ENGDYN system. We have translated the model in X direction by 172.0 mm when importing. This gives a Centre_of_Gravity in the ENGDYN system as approximately (4.0, 110.0, 3.0) mm. Apply the Centre_of_Gravity to the panel as shown
120
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Click on the Radius tab, define the sphere to be 1 meter in radius. Set the Number of Elements to be 40, which gives a suitable refined mesh
Click Generate when the panel is complete A new SFE file named Acoustic_Mesh.SFE is created in the current working directory. Double click on the SFE file to look at the model in R-Desk Viewer.
121
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
The mesh density can be easily changed by changing the number of elements in the Acoustic Mesh tab. Click OK to confirm the use of acoustic model that has been defined, then the File Name field has been updated as shown below.
Click on the Fluid tab, define the Fluid Density and Speed of Sound as below
122
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Click on the Options tab, set the Integration Increment to be zero. This field allows the user to increase the order of integration to improve the calculation accuracy, but at the expense of run time. We leave it as zero, meaning the solver will automatically use the optimum integration order.
Return to the Boundary tab, and the acoustic mesh is fully defined. Select the Frequency Response tab and set the maximum frequency to 2000 Hz Click on the Output tab and set the panel as below
123
KNOWLEDGE CENTRE Tutorials Tutorial 3: NVH Analysis
Click OK Click Generate Output to run the noise calculation in the Cylinder Block FEA panel Note that the solution will take approximately 8 hours on a laptop – a 2.5 GHz processor and 4.0 GB of RAM. Rather than running through the FEARCE NVH solver automatically, we can run the calculation through a batch file. On completion the following files will be written IL4_BLOCK_NOISE_BEM.FRC – FEARCE command file IL4_BLOCK_NOISE_BEM.FRL – containing the frequency response loading IL4_BLOCK_NOISE_BEM.L001.RES – containing the sound power for loadcase 1 IL4_BLOCK_NOISE_BEM.OUT – an output file from FEARCE IL4_BLOCK_NOISE_BEM_EXTERNAL.SFE – containing the vibration results which are interpolated onto the BE model IL4_BLOCK_BoundaryElementModel.SFE – our BE model containing sound power results Acoustic_Mesh.SFE – containing sound pressure results LC001.ALPHA containing the modal contribution factors for loadcase 1 LC001.TOTALS containing the total loads applied to the model for loadcase 1 2.3.7.2 Plotting and viewing results Open the IL4_BLOCK_BoundaryElementModel.SFE by double-clicking this file. Also open Acoustic_Mesh.SFE. In order to open it is the same session, right click in the top left Plots panel and then select “Open Model”. The results can be viewed and analyzed following the same procedure as that described earlier for the Rayleigh method. The sound intensity results are available in IL4_BLOCK_BoundaryElementModel.SFE and the sound pressure results are available in Acoustic_Mesh.SFE.
124
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
2.4 Tutorial 4: Block Stress Analysis
2.4.1
Overview
Objective: Build an ENGDYN model suitable for performing a stress analysis of the cylinder block The analysis uses finite element (FE) models of both the powertrain assembly and the crankshaft This tutorial requires the ENGDYN model from Tutorial 2 (Crankshaft Dynamic Crankshaft) Items Covered: Compliant cylinder block model Restraining the system Static condensation reduction Quasi static stress analysis Estimated duration: 0.5 day Required Files: ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\COMPLIANT_BLOCK.SFE Finite element model requirements: For a stress analysis of the cylinder block it is recommended to have a powertrain assembly that includes all of the major components contributing to the stiffness and hence load paths of the system Method: In many cases a quasi-static stress analysis of the powertrain is a suitable method for analysing the in service behaviour of the cylinder block The quasi-static method has a number of advantages over a fully dynamic analysis The method is computationally less intense allowing iterations to be carried out swiftly The loads are fully resolved statically and can be applied to both linear and non-linear FE solutions
125
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis Thermal loads can be combined with the mechanical loads generated by ENGDYN providing a full thermal mechanical assessment For cases where a known dynamic stress problem exists, ENGDYN can also be used to provide a solution and this is covered in a separate tutorial
2.4.2
Getting Started
The model generated in tutorial 2 (dynamic crankshaft analysis) is used in this tutorial, so we will make a copy of the model created in tutorial 2. Start engdyn On Unix or Linux platforms simply type engdyn On windows click on the shortcut, otherwise go to Start>Programs>Ricardo Software>2014.1>Mechanical Suite>ENGDYN>ENGDYN Open the model from tutorial 2 Use the File > Open command from the top menu bar Select Copy Design from the File menu on the top bar Select a new working directory in the browser and save the model under a new name
Note that saving the file in this way keeps the references to the component files. In this case, it refers to the crank FE model (and reduced components) and the lubrication and pressure data.
126
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis 2.4.3
Updating the Model
2.4.3.1 Define Models In this step, we need to redefine the level of our component cylinder block model Select ‘Define Models’ from the buttons on the left side of the Main Panel to display the Model Definitions Panel
Select the ‘Cranktrain’ tab and define the parameters as follows
Select the ‘Cylinder Block’ tab and set out model level to compliant as shown
Click on the Select button in the File name field and refer to the cylinder block FE model COMPLANT_BLOCK.SFE
127
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis Click on Transformation tab, and translate the model by 172 mm in the X direction
Click OK 2.4.3.2 Reassembling the Crankshaft Model Click on the Edit Cranktrain button Note that we do not need to redefine any data as we haven’t changed the crankshaft. However the program will force us to re-edit the main and thrust bearing data. This is because the thrust bearing stiffness will now be calculated from the FE model. Select the Trust Bearing item from the Cranktrain tool panel Click Edit Selected Click on Stiffness tab, and set 1 in the Housing field Click OK
Select the Main Bearing item from the menu Click on Select All, and then Edit Selected Click Ok in the bearing panel as we do not wish to alter the values
128
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis Run through the crank assembly steps, including ‘Define Material’, ‘Calculate Masses’, and ‘Set Balance’ as we did in tutorial 2 (dynamic crankshaft analysis)
Click on the Calculate Stiffness button
Select No as we haven’t changed the crankshaft model in any way 2.4.3.3 Defining the Reduced Model We now need to define a ‘reduced model’ of the FE cylinder block that we have imported. The reduced model is a number of nodes and degrees of freedom that are a subset of the complete FE model. In order to still be able to interact with our FE model once we have performed the reduction, we need to define the degrees of freedom (nodes) whose details we wish to retain. At a minimum this needs the definition at each main bearing, the definition for each cylinder bore and the definition for each cylinder head gas face. We will need to restrain the FE model in space before performing the model reduction. Suitable restraint points need to be selected - using engine mount locations would provide sensible load paths. For consistency of solution these restraints should be the same as used in the final FE stress analysis. Click on the Edit Block button, and the following query appears
129
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
Click Yes as we have changed the block model type A secondary query panel will appear as below relating to the cylinder definitions
Click YES as we have changed our block model type, and then the Cylinder Block Tool panel appears
Click on the Define Model button in the Cylinder Block Tool panel, and the Define Model panel appears as below
The three main tools are: 1) Clip: This allows us to adjust the geometric definitions of the required regions and tune the search of the FE model for the required nodes 2) Read: This allows us to read pre-defined ‘sets’ of data directly from the FE model to use in our reduced node definition
130
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis 3) Delete: This allows us to delete any definitions which are not wanted Click on the Delete tab and select the Constrained Sets tab
Highlight all five items in the window and click Delete
Click on the Clip tab (we are going to tune the definition of the main bearings to try and find the complete faces.) Under the name tab, define the parameters as shown below
131
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
Click on Select button next to the ID field and choose main bearing 1 from the list of options.
Click OK Click on Definition tab This section reveals the geometric definition that ENGDYN considers to be the inner surface of main bearing 1
132
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis Click on the Tolerance tab ENGDYN has calculated this tolerance by reading the FE model and determining the smallest length of an element side. The tolerance is half of that smallest length. But if some elements are very small, then this tolerance will be too small to correctly locate the surface that we require. Change the tolerance to 0.5 mm and click on Clip Set button. ENGDYN finds the complete bearing surface shown in the viewport.
Click Add Set and the viewport will update to show the complete bearing defined in orange
Repeat this procedure for each of the five main bearings In each case, change the bearing ID in the Name tab, then change the geometric tolerance to 0.5mm, Click Clip set and Click Add Set On completion the viewport should like as shown below. The next job is to correctly define the nodes required to represent the four cylinders.
133
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
Return to the Name tab in the Clip section of our Define Model panel Select Node in the Type field and select Cylinder in the Name field, as shown
Click on the Select button in the ID field and select upperThrust as shown
134
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis Click on the Select button in the Cyl field and click Select All button in the Cylinders field
Click OK and the Define Model panel is shown as below
Click on the Tolerance Tab and set 0.5 in the Linear field
Click on Clip Set and click Add Set
135
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
Repeat this step for each of the remaining cylinder nodes. You should be able to clip for all four cylinders at once each time upperAntiThrust lowerThrust lowerAntiThrust Set the tolerance to 0.5 mm Remember to clip the set to locate the nodes in the FE model before adding them
We will define the next sets using the Read option. This is used to select regions already included as named sets on the FE model. Because of the complexity of the shape of the flame face it is often easier to define them in this way than to try and create a shape to automatically clip from. So we can prepare our models using an FE package (either a third party package or Ricardo’s R-Desk VIEWER tool) by ‘painting’ required regions and naming these as face sets.
136
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis Select the Read tab on the Define Model panel and select type as Constrained Sets You will see a list of sets that exist on the model which we can choose from.
To add these all we need to do is highlight the ones we want, click Read Set, then Click Add Set. Highlight all of the sets and add them to the model Click on OK in the Define Model panel. The viewport should update to look similar to that shown below
Next, we can perform the matrix reduction. Go back to the Cylinder Block Tool panel
137
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
Select the Matrix Reduction option Select FEARCE option from the Solver list, as shown below
The Solver field allows us to choose which solver we would like to use for the reduction. Currently ENGDYN can support MSC/NASTRAN, ABAQUS or its own internal solver FEARCE. A translator can be supplied to export an appropriate analysis deck to NASTRAN or ABAQUS. If FEARCE is selected, no external licenses are required. The Output Name field determines what these files will be called and where they will be located on your system. Click the Default button next to Output Name field, which means the files will be placed in the working directory and named after the powertrain FE model The Element Check allows a check of the FE model to be performed. Click Element Check and limited memory. Set a memory limit approximate to your machine as shown below
138
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
Press Solve and the following dialogue panel appears
This dialogue is telling us that the command file to set up the reduction has been written. Choosing to continue will execute this, causing the following actions. 1) The reduced nodes will be added to the FE model; 2) The FEARCE solver will then perform the reduction. If NASTRAN was chosen for the solver, the NASTRAN input deck will be written out at this stage. Choose Ok You will then be asked whether you want to perform the reduction immediately or later (in batch) Select the Direct option The command file will execute and the reduction will begin This powertrain model has about 190000, and this is to be reduced to about 100. Click Assemble Model when the reduction is complete Close the Cylinder Block Tool panel 2.4.3.4 Performing ENGDYN solution We are now ready to perform the ENGDYN solution. The lubrication and cylinder pressure was defined when the original model was set up for the dynamic crank analysis. Hence we can go straight to the Evaluate Solution option. Select the Indeterminate option from the solution type list
139
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis As we are not considering the dynamics of the system, the solution convergence is usually straight forward. Make sure that the cylinder block model is set to Compliant as shown Set up other parameters as follows
Move to the Cases tab Add a single loadcase of 5000 rev/min in the speed table as shown
Click on the Bearing Model tab
140
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis Select the thermal balance option This will calculate the rise in temperature across the bearings due to work This will not slow an indeterminate solution too much so is worthwhile Finally select the required oil properties file The solution should only take a few minutes to complete with the message shown below
Close the message window 2.4.3.5 Block Stress Analysis First step is to select what loading condition we require. Click on the Select Loadcase button Make sure that the Type of results is Statically Indeterminate Dynamic results will not be fully resolved statically and are inappropriate We only have one set of results, so select this in the list
141
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
Click Ok when done Select Block Analysis form the post processing menu The cylinder block FEA panel will appear. This provides us with the workflow required to set up the FE analysis.
Click on Select Output Choose the Quasi-Static loading option from the list, shown as below
Click OK
142
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis Click on the Select Model option and the Cylinder Block Model appears
Select the block required for the stress analysis shown as follows
The model transformation used in the ENGDYN solution will be loaded by default. We can change this as before in using the panel under the Transformation tab.
When done, click OK Click on OK on the Cylinder Block Model panel Click on the define the output
143
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
We shall select which FE solver we want to use. Choose Ricardo-SFE. This will automatically apply the loads to the FE model. Click ‘Select’ next to the Crank Angle field
Select the ‘Use Interval’ option and set the interval to 10 degrees This will give us 72 load cases from ENGDYN
144
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
Click OK when complete The final main field that we can edit is the ‘Start Loadcase’ section. This defines what the first loadcase number that the ENGDYN loads are written to is. For a complete analysis we would usually add assembly and thermal loads to the system in addition to the ENGDYN loads. These would often be added as loadcases 1 and 2 and be added to the model as part of the preparation before ENGDYN. Hence we could start the ENGDYN loads from loadcase 3. In this example we will use the default as we do not have other loads The Output Name field determines where the files are written that control the FE solution.
Keep this as default and click OK click on Generate Output
145
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis ENGDYN will create a FEARCE run script that will apply the loads onto the FE model and run the solution or set up an appropriate third party solver deck. With the SFE option chosen for loading, this will be executed automatically. This will take several minutes to complete. The STRESS_BLOCK.SFE model is now loaded and ready for solution.
These loads can be viewed in R-Desk viewer Use the Data > Contours > Force option
To implement the FEARCE solver a small script file needs to be created. The following file has been included named Run_stress.FRC
146
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
This can be executed from a command line in the appropriate directory. Simply type FEARCE and the file name shown as follows.
Finally, open the STRESS_BLOCK.SFE file in order to view the results. The process is similar to the one described in the NVH tutorial for viewing the vibration results (section 2.3.5.4). Please note that this tutorial model in only provided for the purpose of illustrating the process and the used elements are not appropriate for stress analysis. A plot similar to the one below can be obtained by adjusting the scales.
147
KNOWLEDGE CENTRE Tutorials Tutorial 4: Block Stress Analysis
148
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
2.5 Tutorial 5: Crankshaft Stress Analysis 2.5.1
Overview
Objective: To use ENGDYN to perform a stress and durability analysis of a crankshaft The crankshaft utilises a finite element (FE) model. The cylinder block model is assumed rigid in this example. Items Covered:
Preparing the Engdyn model Geometry Definition Assigning Material Properties Stress Concentration and Fatigue Notch Factors Classical Stress Method Finite Element Method Plotting and Viewing Results
Estimated duration: 0.5 day (model preparation) 1 day (overall including performing solutions) Engine: Inline 4 gasoline engine Required Files: ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\il4_crank.SFE ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\IL4_STRESS_CRANK.SFE ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\il4*.PRES ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\il4\dynamic_crank.EDSF Finite element model requirements: Crankshaft FE model should be split into its individual web sections as described in Finite Element Model section of the manual
149
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis 2.5.2 Preparing the Engdyn Model 2.5.2.1 Performing the Solution A pre prepared .EDSF file is provided for this tutorial. This is a copy of the Engdyn model as prepared in the ‘Crankshaft Dynamic Analysis’ tutorial. The reader is invited to review this tutorial to familiarise themselves with the methods of preparation. We shall begin this tutorial by performing a solution suitable for the crankshaft stress analysis Open Engdyn and click on the File button on the main toolbar Click on Open from the drop down menu Select the ‘dynamic_crank’ .EDSF from the tutorial folder The display should look like the figure below
The lubrication should have been defined, but check by clicking on the Lubrication button and ensuring that the file SAE5W30.MAT is being referenced
150
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis If not click on the Browse button and reference the file so that the lubrication panel looks like the figure below
Click on OK in the Lubrication Definitions panel Click on the Loading button Ensure that the In-Cylinder panel is set as shown
Click on OK Click on the Evaluate Solution button Configure the Solution panel as shown
151
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
We shall carry out a fully dynamic solution, but will also save the statically indeterminate results – ensure that the ‘Save Static Results’ box is checked Configure the Cases panel as shown – setting a speed range of 2500-5500 rev/min in 500 rev/min intervals
We will use the default settings for the cylinder block, crankshaft, thrust bearing and coupling panels Click on the Bearing Model panel and configure it as shown below
152
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
To save time in this example, we shall explicitly define the oil temperatures at each bearing – check that these are all 110 OC using the Define Oil Temps button Click on Solve Directly.
When the solution has completed, check that it has converged for every loadcase We are now ready to carry out the crank stress analysis 2.5.2.2 Geometry Definition To obtain stress results at specific points on the crankshaft, we must define the geometry at these locations in more detail. This is done through the Crankshaft Stress Analysis panel, accessed by clicking on the Crank Analysis button. Click on the Crank Analysis button to display the following panel
153
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
We shall calculate stresses at the main and pin journal fillets at the rear of the crank To do this we must also define the geometry of the surrounding webs We shall regard all the webs as being uniform and so define all components at once Highlight ‘Crankshaft webs’ in the left of the panel Click on Select All The webs will show in red on the display as shown
Click on Edit Selected to display the Crankshaft Webs panel
154
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Fill in the ‘white’ boxes as shown Click on Apply The equivalent width box will update Click on OK Highlight ‘Pin journal fillets’ in the left hand of the panel Click on Select All Click on Edit Selected to display the Pin journal Fillets panel Set the Radius to 1.5 Click on Apply The panel should update as shown
155
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Click on OK Highlight ‘Main journal fillets’ in the left hand of the panel Click on Select All Click on Edit Selected to display the Main journal Fillets panel Set the Radius to 1.5 Click on Apply The panel should update as shown
156
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Click on OK The geometric data required for this tutorial has now been applied. The next step is to define the material properties of the crankshaft. 2.5.2.3 Defining the Material Properties At this stage we are not referencing an FE model, so we will define the material properties from the Engdyn database. The procedure that will be shown can also be used to enter the value explicitly. Click on the Define Material button as shown to display the Crankshaft Material Properties panel
157
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
In the Base Properties tab, ensure that the panel is configured as shown below
158
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis Only be concerned with the upper half of the panel (shown above) for the moment. The data in the lower half will be updated when all of the properties have been added Click on the Base Strengths tab
We will use properties for Steel 38MnS6 from the Engdyn database Click on Select (circled) to access the Engdyn materials database Highlight the STEEL38MnS6 properties as shown then click OK The base properties will be updated in the upper part of the Base Strengths panel
159
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Now click on the Elevated Strengths tab
160
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Select Calculated in the Strengths field Click on Define next to Size Factor to display the following panel
Ensure Default is selected as the Definition of the size factor
161
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis Click Apply, then click OK Back in the Crankshaft Material Properties panel, select Cold Rolled in the pin and main journal fillet fields The panel should appear as shown below
Click on OK 2.5.2.4 Defining Stress Concentration and Fatigue Notch Factors Engdyn will take into account any stress increase or relief due to geometric effects by applying stress concentration and fatigue notch factors to the journal fillets. These factors are applied by altering the material strengths at the locations. To apply these factors we must access the fatigue notch factors panel. Click on the Notch Factors button as shown
162
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
We will apply the Modified Lowell equation to calculate the required factors Ensure that Modified Lowell is selected in the Characteristic field Click on Apply, so that the panel appears as shown
Click on OK 2.5.2.5 Classical Evaluation of Stress and Durability Now that we have defined the geometry and material properties of the regions of interest on the crankshaft, we are in a position to carry out a classical evaluation of the stresses and durability. First we must select which loadcases we are interested in for the solution.
163
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis Click on the Select Loadcases button on the main panel to bring up the following sub panel
We shall perform a speed sweep across all defined engine speeds which have dynamic results Display the dynamic results by setting the Solution Type field to dynamic Click on the Select All button to highlight all speeds Click on OK Click on the Define Output button on the Crankshaft Stress Analysis panel Set the Classical Stress Analysis sub panel as shown Note: the output name should reference the current Engdyn model
164
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Click OK Click on Solution on the Crankshaft Stress Analysis panel Engdyn will now perform the stress solutions at each part of the crankshaft for which geometry has been defined and for each loadcase When this is complete we can plot results 2.5.2.6 Plotting Classical Stress Results We can now view results at any part of the crank for which the geometry has been defined. As an example we shall look at results for the rear pin journal fillet on the fourth pin. Highlight the pin journal fillets on the left hand of the panel
165
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
The pin journal fillets that have been defined will turn green on the reduced model Highlight the rear fillet on pin journal four by clicking on it with the left mouse button
The fillet will turn red Now click on the Plot Results button
166
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
The following panel will appear
We will look first at the stress history at the fillet Highlight History, Stress and Goodman-Standard-General options as shown
Click on Apply The graph shown overleaf will appear
167
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis The plot can be printed or closed through the File tab on the top menu of the plot
Now we will look at the predict durability safety factors at the fillet Highlight Summary, Safety Factor and Goodman-Standard-General as shown
Note: options that are not applicable disappear from the list of choices Click on Apply to display the plot
168
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
The top graph shows the maximum and minimum (alternating) stresses acting on the fillet The lower graph shows the lowest predicted durability safety factors on the filler across the speed range The red line is the design factor which can be changed in the Plot Results panel
169
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis 2.5.2.7 Performing a Finite Element Solution Now we have carried out a classical analysis, we are ready to move to the next level of complexity and perform an FE analysis on the crankshaft. The first step is to create the boundary conditions so that we can perform the solution. These are applied to specific sets that should already have been added to the crankshaft FE model. Details of these sets can be found in section 7.1.2.2 of the manual. A brief description of the sets required for loading follows. A layer of low stiffness (approx. 1/10th crankshaft stiffness) and massless elements need to be added to each end of the crankshaft to ameliorate the stresses applied to each end (shown in red below). The length of these should be about ¼ of a typical bearing for the crankshaft
A layer of massless shells (all other properties the same as the crankshaft) should be overlaid onto the low stiffness elements on the cranknose end Face sets need to be added to all of the Pin and Main journals. These represent the oil contact areas of the bearings The sets need to be named with the convention shown in the figures below
170
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Once the crankshaft has been properly prepared we are ready to generate the boundary conditions Click on the Select Loadcases from the main panel For an FE solution we must select a static loadcase – hence the statically indeterminate solutions were saved Click on the Select Loadcases from the main panel
171
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis For an FE solution we must select a static loadcase – hence the statically indeterminate solutions were saved.
Select 5500 rev/min from the statically indeterminate results Click on OK Now click on the Crank Analysis button in the main panel Set the Stress Analysis sub panel to FE Quasi Static method and Unit Loads output as shown
172
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Reference the stress quality model by clicking on the Select Model button on the right of the panel
Ensure that the model is in the same orientation as the Engdyn model by supplying the correct transformation from the appropriate tab Warnings may appear in the message box of the Engdyn panel if the defined geometry does not match up to the FE model Care must be taken to ensure that the FE model aligns correctly with the defined geometry Click on Apply, then OK Now click on the Define Output button to bring up the panel shown below Set the desired Solver and ensure that the Centrifugal Load and Unit Moments boxes are checked The output name will define what the .FRC file is called
Click on OK
173
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis Click on Solution Engdyn will create the .FRC file to load the crankshaft. When complete the following box will appear prompting as to whether the .FRC file should be run
Click Yes The .FRC file will run, and where appropriate the input deck for the solver will be written out If FEARCE(VSS) solver was chosen, then the solution will begin automatically The FE solution should now be performed. If a solver other than the FEARCE(VSS) solver was used, then the results need to be appended onto the stress FE model so that they can be post-processed and viewed. 2.5.2.8 Post Processing the FE Results Once the FE stress analysis has run, ensure that the results have been appended to the stress model .SFE where appropriate. We can now post-process and view the results. Click on Select Loadcases Ensure that the statically indeterminate loadcase that has been run is selected >
Click OK
174
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis Click on the Crank Analysis button and select FE Quasi Static method and Stress and Safety Factors as output (as shown below)
From the Select Model button, ensure that the correct model is referenced (stress model with results) and that it is in the correct orientation We shall again just consider the rear pin fillet of pin journal four Highlight Pin Journal Fillets in the left hand side of the panel (the pin fillets on the reduced model turn green) Select the rear journal fillet on pin four by clicking on it with the left mouse button (it will turn red)
Click on Edit Selected
175
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis As we have previously defined the geometry we do not need to change the entry
Click OK on the Pin Journal Fillets panel Click on the Define Material button if the material defined in the FE model is not identical to that in the Engdyn model (for example hydrostatic fatigue factors may not be defined in the FE model), the following panel will appear
Click on Select to display the panel below
Highlight the ENGDYN properties (in this case) and click on OK Click on the Notch Factors button and check that the definitions are correct
176
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Click OK in the Notch Factors panel We now need to define the type of algorithm for the durability safety factor calculation Click on Define Output
Click on the Select button next to the Node Sets field Highlight REAR_PIN_FILLET_04 from the sub panel and click OK
177
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis The Node Set sub panel will display all of the relevant sets defined on the FE model – we can only select sets that are also have their geometry defined in the Crankshaft Stress Analysis panel Click on the Safety tab For this example we shall only look at the Alternative Goodman criteria Highlight as shown and click OK Check that the rest of the FE Stress Analysis panel is configured as above and click OK The Export field will determine how much of the FE model is exported for post-processing. The alternatives are ALL, EXTERNAL or SETS. The smallest model will be just the sets The post-processing is now ready to run. Engdyn will first create a new .SFE of the chosen subset of the model (in this case just the selected sets). It will then combine and factor the unit loadcase stress results before finally performing the durability calculation. The results can then be plotted in Engdyn, or colour contours viewed with R-Desk. Click on Solution Engdyn will write out an .FPO file to perform the post-processing. When this is done the following panel will appear summarising the forces in each separate loadcase >
Click on OK Click on OK in the Query panel (below) to perform the post-processing
178
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
When the post-processing is complete the following panel will appear
Click on OK
2.5.2.9 Viewing and Plotting the FE Results If the post-processing completed successfully a new .SFE file would have been written to the directory containing the original stress quality FE model. This will have the word ‘STRESSED’ appended to the original model name, i.e. _STRESSED.SFE If we click on the Plot Results button in Engdyn, we can generate graphs from this new .SFE (referenced automatically) in exactly the same way as detailed in the classic analysis section. In addition, we can also view 3 dimensional contour plots of the areas that were selected using R-Desk. Click Crank Analysis on the left side of the software
179
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Select the pin journal fillets item from the menu list and select the rear pin fillet on the crank model in the 3D viewport We have calculated FE results for this item so can now plot the results Click Plot results We get the same panel as we used in the classic analysis
Have a look at a selection of the plots
180
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis
Next let us view one of the FE models
181
KNOWLEDGE CENTRE Tutorials Tutorial 5: Crankshaft Stress Analysis Double click on the file crank_stress_L0009.SFE The sub model is just the section that we exported There is now an extensive amount of data stored on the model A lot of this is resulting from sub-steps in the post processing and can be discarded Select the Fatigue tab on the data panel You will see three data sets, all labelled 09 for the load case number 9 QS09, Durability for quasi-static loads only QST09, Durability for quasi-static loads plus torsion vibration QSV09, Durability for quasi-static loads plus torsion and bending vibration
Explore the different results and views available
182
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
2.6 Tutorial 6: EHL Big End Bearing Analysis 2.6.1
Introduction
Objective: To build an ENGDYN model suitable for performing a static EHL analysis of a big end bearing. The crankshaft and cylinder block models are both assumed to be rigid. This tutorial can be applied to a single cylinder of a multi-cylinder engine Items Covered:
Building the engine model Connecting rod model definition Static Matrix Reduction Static analysis of a big end bearing Plotting and animating results
Estimated duration: 0.5 day (model preparation) 1 day (overall including performing solutions) Engine: Single Cylinder HYDRA Research Engine 80.25 x 88.9 mm Required Files: ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\hydra\ROD.SFE ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\hydra\SMALLEND.SFE ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\hydra\BIGEND.SFE ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\hydra\pcyl2000.PRES ..\Ricardo\2014.1\Products\ENGDYN\Tutorials\hydra\RotellaX15W40.MAT Finite element model requirements: Model should contain bearing shells and big end cap bolts A material should be defined for each bearing (if different) to define the lining of the bearing shell
183
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
2.6.2
Getting Started
Copy the necessary files from the example directory to a working directory and ensure that you have write permissions for all the files. Start engdyn On Unix or Linux platforms simply type engdyn On windows click on the shortcut, otherwise go to Start>Programs>Ricardo Software>2014.1>Mechanical Suite>ENGDYN>ENGDYN
2.6.3
Building the Model
2.6.3.1 Configure Engine In order to build an ENGDYN model it is first necessary to define the major dimensions and features of the engine. This is done using the ‘Configure Engine’ Panel. Select ‘Configure Engine’ from the buttons on the left side of the Main Panel.
Complete the Engine Configuration Panel as shown.
On completion select ‘Apply’ which will display
184
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
Select OK. The crankshaft model will appear in the Main Panel as shown.
2.6.3.2 Define Models An ENGDYN model consists of a number of sub-models which are defined using the Model Definitions Panel. Select ‘Define Models’ from the buttons on the left side of the Main Panel.
Complete the model definitions for the Crankshaft, Cylinder Block and In-Cylinder as shown.
185
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
The crankshaft can be defined as Rigid (Massless) since for the purposes of this static calculation the loads at the big end bearing will not be influenced by the crankshaft mass. This will mean the minimum amount of data will be required subsequently. Similarly the Cylinder Block can be defined as rigid. Define the Connecting Rod model type by positioning the mouse over the centre column and selecting the right mouse button. This will display the pop-up as shown.
Different models can be defined for each cylinder (when there are multiple cylinders) Compliant Big End, Compliant Small End and Compliant are applicable to static calculations only Rigid and Dynamic are applicable to dynamic calculations only. If nothing is defined or None is selected then the solution for that cylinder is equivalent to v3.0 and earlier Complaint Big End and Compliant Small End are used where a complete model of the connecting rod is not available. Examples of these models are SMALLEND.SFE and BIGEND.SFE. (You may wish to look at these models using FEVIEWER) The Complaint Big End and Compliant Small End require a cut-plane somewhere along the shank of the rod.
186
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
Select Compliant
187
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis Select the finite model by positioning the mouse over the centre column and selecting the right mouse button. This will display the pop-up as shown. Releasing the button will pop-up the model translation panel as shown.
Use the Browse button to select the file ROD.SFE from the working directory Define the transformation vectors to transform the model from the finite element model coordinate system to the ENGDYN coordinate system
These vectors are with respect to the ENGDYN coordinate system The connecting rod of each cylinder has its own local coordinate system such that the origin is at the centre of the big end bearing, Y is up the connecting rod, and X runs along the bearing axis Select OK on both the Model Translation and the Define Models Panels. This will display the connecting rod model as shown. The program determines the number of unique connecting rod models and writes a message at the bottom of the Main Panel as shown
188
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
2.6.3.3 Edit Cranktrain Select ‘Edit Cranktrain’ from the buttons on the left side of the Main Panel to display the Cranktrain Tool Panel
Given we defined the crankshaft as Rigid (Massless) (See 2.6.3.2) it will only be necessary to define Main Bearings, Big End Bearings, Small End Bearings, Cylinder, Connecting Rod and Piston.
189
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
Highlight ‘Main Bearing’ as shown This will then display each of the main journal nodes in green
Select ‘Select All’ Each of the main journal nodes will then turn red Select ‘Edit Selected’ to display the Bearing Panel Complete the Bearing Panel as shown
Given the Model Type is Mobility and the bearings are fully grooved it is not necessary to enter any other data than shown Select Apply to display the panel showing the bearing grooves
190
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
Select OK Highlight ‘Big End Bearing’ This will then display the big end journal node in red since there is only one cylinder. There is therefore no need to select ‘Select All’ Select ‘Edit Selected’ to display the Bearing Panel Complete the Bearing Panel as shown
This bearing is feed from the journal via a feed from the adjacent second main bearing The oil hole angular position can be defined either using the Position column (as in this case) or using the Height column For an EHD Model Type it is necessary to use the Mesh and Material tabs to define additional data for this model. In this first exercise we will assume the bearing and journal are circular, and therefore it is not necessary to define a profile using the Profile tab. Select the Mesh tab, and define a mesh 9 x 55
191
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
This mesh is the computational mesh for the bearing model Note how the mesh is refined to resolve the oil hole Select the Material tab
This material data is used for the boundary lubrication model. It is necessary to define the journal material and the material of the bearing lining. Click on Define adjacent to Bearing Material, to display the Material Properties Panel
192
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
In this case the material property data will be read from ROD.SFE It is necessary to have defined a material in the finite element model that corresponds to the bearing lining, although this material will not be assigned to any finite elements Use Select to list the materials and select the material BEARING from the list as shown
This table shows all the materials in the SFE file, relevant to the material being edited, together with any materials stored by ENGDYN. (In this case there will be no ENGDYN materials in since we have not yet added any materials) Materials used by ENGDYN and those stored in the .SFE are identified by a single unique name. These data have been defined previously outside ENGDYN, although they may be and can be zero. Select OK and the material is added to the panel as shown.
193
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
These data values are typical for a bearing surface and can be edited using these panel These data are only used by the boundary lubrication model Asperity RMS Height h, Density and Asperity radius can be calculated from measured data using the Ricardo MATUTIL program which is described in Appendix 12. Select OK These data are now stored by ENGDYN and also written to the .SFE file. Previous data will be overwritten. Click on Define adjacent to Journal Material, to display the Material Properties Panel There is no FE model of the crankshaft so this time don’t use the Select button (since it will only list the previously defined material called BEARING.) Complete the Material Properties Panel as shown
Select OK and enter the wear and friction coefficients for the bearing interface
194
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
These data are only used by the boundary lubrication model The wear coefficient is not used for the solution, but only as a scalar for postprocessing. Select OK to complete entering data for the big end bearing. Highlight ‘Small End Bearing’ This will then display the big end journal node in red since there is only one cylinder. There is therefore no need to select ‘Select All’ It is not strictly necessary to define data for the small end bearing but by doing so helps defining the connecting rod reduced model in 2.6.3.4. A bearing solution of this bearing will be performed. Select ‘Edit Selected’ to display the Bearing Panel Complete the Bearing Panel as shown
195
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
Given the Model Type is Mobility and the bearing has a single oil hole feed it is not necessary to enter any other data than shown Select Apply to display the panel showing the bearing oil hole Select OK Highlight ‘Connecting Rod’ as shown
This will then display the big end journal node in red since there is only one cylinder. There is therefore no need to select ‘Select All’ Note that the Define Model and Matrix Reduction buttons can now be selected. These will be used in 2.6.3.4 and 2.6.3.5 respectively.
196
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis >
Select ‘Edit Selected’ to display the Conrod Panel Complete the Conrod Panel as shown
Note that we have changed the units of Inertia to kg.mm^2 In this case because we have a complete model of the connecting rod this data will be overwritten with that of the FE when we perform 2.6.3.6. Select OK Highlight ‘Piston’ This will then display the big end journal node in red since there is only one cylinder. There is therefore no need to select ‘Select All’ Select ‘Edit Selected’ to display the Piston Panel Complete the Piston Panel as shown
Select OK.
197
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
This completes the process of defining the cranktrain data other than defining the connecting rod reduced model and performing the matrix reduction for that model. This is dealt with in the next two steps. 2.6.3.4 Edit Cranktrain – Connecting Rod Model Definition It is necessary to define a reduced model of the connecting rod. The reduced model is a number of nodes and degrees of freedom that are a subset of the complete model. We need to define two sets. A face set that defines the big end bearing, since this is bearing is an EHD model. For this set all the degrees of freedom of the bearing surface are included in the reduced model. In addition we require a constrained node set that defines the small end bearing. For this set the bearing is defined by a single node with 6 degrees of freedom whose movement is the average of the nodes on the bearing surface. Highlight ‘Connecting Rod’ as before
Select ‘Define Model’. Messages will appear in the bottom of the Main Panel and the Define Model Panel will appear as shown
198
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
In this action the program has attempted to automatically define the two sets using the data supplied previously – bearing position, diameter and length. The reduced model is shown in orange The set defining the big end bearing is incomplete and nothing has been defined for the small end bearing. This is because the nodes are outside tolerance. It is necessary to define these sets using the Define Model Panel shown on the right. Consider firstly the face set defining the big end bearing Select the Definition tab
Select each tab and understand the default values Each set is ‘clipped’ based on a geometric shape
199
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis Each shape is defined by a Centre, Axis, Extent and Diameter Select the Tolerance tab
The default linear tolerance is set to half the minimum distance between any two adjacent finite element nodes. Each set has its own tolerance values Change the Linear Tolerance to 0.5 mm and press ‘Clip Set’.
Select ‘Add Set’ to add the set to the reduced model
200
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
Secondly consider the constrained node set defining the small end bearing Select the Name tab Change Type to Constrained Node and Name to Small End Bearing
Select the Definition tab and again understand the default data Select the Tolerance tab and again change the linear tolerance to 0.5 mm Select ‘Clip Set’ followed by ‘Add Set’ to obtain the following
201
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
Use the mouse and cursor keys to rotate and zoom the model if you haven’t already done so OK the Define Model Panel This completes the process of defining the reduced model of the connecting rod
202
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
2.6.3.5 Edit Cranktrain – Connecting Rod Model Reduction Select ‘Matrix Reduction’ to display the Matrix Reduction Panel
For this tutorial we will use the FEARCE Vectorized Sparse Solver (VSS) If you have limited memory then this option can be used to limit the memory during forward elimination and backsubstitution. This will result in a longer solution time.
Select the Default output name The output name specifies the name of the FEARCE .FRC command file that will be written. Selecting the Default button will give the output name the same name as the .SFE file Note: the include field is for adding extra commands directly to the run script – we will not need to use this feature in this example
Switch on the Limited Memory toggle and specify a memory requirement of 250 Mb or less if your machine has less available memory
Select Solve. The following Query Panel will be displayed.
203
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
Loadcases 1 to 6 are reserved for 6 body loads in each of the six global directions. Any additional loads the user may which to apply, such as bolt loads, must be applied using loadcase 7 and higher. Any number of loadcases may be used. Select OK. Another Query Panel will be displayed.
Select ‘Direct’. This will execute the command file and start the solution. A Progress Panel will be displayed.
This model has approximately 100000 degrees of freedom which are reduced to 372 for the reduced model to be used by ENGDYN. On completion the message ‘Run completed successfully’ is written to the Main Panel and the following Information Panel is displayed.
204
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis OK this panel and the Matrix Reduction Panel 2.6.3.6 Edit Cranktrain – Model Assembly Select ‘Assemble Model’ The following panel will appear as we are using an FE model for the connecting rod
Click no as the FE model does not contain bolts etc. and so the values entered in the setting up of the model are more accurate
This will assemble the connecting rod model. Messages will be written to the Main Panel as shown Save the model using the File menu from the top of the Main Panel 2.6.4
Solution
2.6.4.1 Define Lubricant Properties
205
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis Click on ‘Lubrication’ button from the buttons on the left hand side of the Main Panel
Use the Browse button to select the lubricant RotellaX15W40.MAT from the working directory By default on Unix the program will initially select the database directory at ../Ricardo/engdyn/3.0/database/Fluid This database contains the most common lubricants Use either Add or Update to add the lubricant
Click on OK
2.6.4.2 Define Loading Conditions The cylinder pressure diagram and any additional loadings (for example loads from Valdyn and gravity forces) are entered in this step. Click on the ‘Loading’ button on the Main Panel
This will display the Loading Definition Panel as shown
206
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
A number of different loading maps can be defined, Full Load, Part Load and No Load. The solver will interpolate at speeds between those defined using this panel. For this tutorial we will define a cylinder pressure at a single speed. Any solutions performed at other speeds will use this pressure diagram. Type in a speed of 2000 rev/min Set the Interval to 10 deg The pressure file can contain just pressures at equal intervals as in this case at an interval of 10 deg or as an array of pressures and angles (which may be at unequal intervals) Position the mouse over the File Name column and use the right button to display the pop-up menu. Use ‘Select Pressure file’ to select the file pcyl2000.PRES from your working directory. The panel will appear as shown.
You may wish to remove the pathname in front of the file. You may wish to inspect and understand the selected file using an appropriate editor Define the Ambient and Crankcase Pressure as 1.0 bar The crankcase pressure is only currently used for Hydrodynamic and Elastohydrodynamic bearing models to define the boundary condition at the edge of the bearing
207
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
The completed panel will appear as shown
The remaining tabs, Force Profile, Force Equation and Distortion need not be completed for this tutorial. Use the Plot button to display the applied pressure loading as shown
Save the model using the File menu from the top of the Main Panel We now have sufficient data to proceed with an analysis
208
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
2.6.4.3 Evaluate Solution Click on the ‘Evaluate Solution’ button on the Main Panel
The Evaluate Solution Panel will appear.
Set the End Angle to 1080 deg This is equivalent to three cycles and should be sufficient to obtain a converged solution of the big end bearing EHL solution Define the Connecting Rod Model as Compliant This ensures we run an EHL solution rather than a rigid bearing solution Select the Cases tab and define a single speed of 2000 rev/min as shown
209
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
Select the Bearing Model tab and complete the panel as shown
The Half-Sommerfeld, Mass Conserving (Recommended) model is recommended for EHL calculations This solver is mass conserving but without oil film history via the Reynolds cavitation condition. The remaining tabs do not need to be completed for this solution since the other parameters are related to indeterminate or dynamic analyses. Select ‘Define Oil Temps’ to display the Bearing Oil Temperatures Panel
210
KNOWLEDGE CENTRE Tutorials Tutorial 6: EHL Big End Bearing Analysis
We will assume that the temperature of the oil in the bearing is at the inlet temperature Select OK Select Solve Directly on the Evaluate Solution Panel.
On completion of the analysis the summary file .EDSUM will be written. This file contains summary data for the solution. Open this file with an appropriate editor to view the results.
211
HELP Using ENGDYN
B. HELP 1 Using ENGDYN 1.1 Overview The ENGDYN Graphical User Interface provides the user with complete control of the model generation, solution and results presentation phases from within an easy to use graphical environment. Complex models can be rapidly generated, solved and the results viewed, without the need to generate input files by hand. The ENGDYN model can be saved at any point during construction. ENGDYN uses a Ricardo binary standard data file to store both the model and the results of the solution and has an .EDSF suffix. When saving, ENGDYN stores the model at the current position and as the model is built up, so the filename.EDSF file is appended to. Similarly, once the simulation has been run, so the results are appended to this file.
1.2 Getting Started To start ENGDYN enter the following command: engdyn [-V ] [-debug