Skip to content

FEM Simulation of Deep Hole Drilling Processes

Deep hole drilling is the most optically blind operation in machining. The cutting edge is buried at the bottom of a deep, narrow bore, submerged in opaque coolant, and loaded with pressures that destroy any sensor placed near the cutting zone. No camera sees the chip form. No thermometer reads the 800°C at the tool tip. No probe measures the guide pad pressure against the bore wall. Finite element simulation is not an academic exercise in deep hole drilling — it is the primary tool for understanding what happens inside the bore. A validated FE model of a gun drill or BTA tool predicts cutting forces within 10–15% of experimental measurement, temperatures within ±50°C, and residual stress profiles that correlate with X-ray diffraction data. The model reveals whether the chip will break cleanly, whether the guide pads will overheat, and whether the tool will deflect off-centre — all before a single hole is drilled.

Why FEM for Deep Hole Drilling

Deep hole drilling presents unique challenges that make simulation particularly valuable:

ChallengeWhy It MattersHow FEM Helps
Inaccessibility of cutting zoneNo direct observation possibleProvides full-field temperature, stress, strain data
Long, slender toolDeflection and vibration affect bore qualityPredicts straightness and dynamic behaviour
Guide pad contactBurnishing dominates bore surface integrityModels contact pressure and friction at pads
High pressure coolantCools and evacuates chips — hard to modelCoupled CFD-FEM captures thermal effects
Residual stressDetermines fatigue life of componentPredicts subsurface stress distribution
Tool geometry complexitySingle-lip design with multiple clearance anglesOptimises geometry before manufacturing

What FEM Can Predict

OutputTypical AccuracyValidation Method
Feed force±10–15% of measuredDynamometer
Torque±10–15% of measuredRotary torque sensor
Cutting temperature±50°CThermocouple, thermography
Chip morphologyQualitative agreementMicroscope imaging
Residual stressMatches XRD profileX-ray diffraction
Guide pad pressureIndirect validationWear pattern analysis
Tool deflection±20%Bore straightness measurement

Material Constitutive Models

Johnson-Cook Model

The Johnson-Cook (J-C) plasticity and damage model is the most widely used constitutive model for machining simulation:

[ \sigma = (A + B\epsilon^n)(1 + C\ln\dot{\epsilon}^*)(1 - T^{*m}) ]

Where:

  • (A) = yield strength at reference strain rate
  • (B) = strain hardening coefficient
  • (n) = strain hardening exponent
  • (C) = strain rate sensitivity coefficient
  • (m) = thermal softening exponent
  • (\dot{\epsilon}^*) = dimensionless plastic strain rate
  • (T^*) = homologous temperature

J-C Parameters for Common Deep Hole Drilling Materials

MaterialA (MPa)B (MPa)nCmMelting Temp (°C)
AISI 4140 (low alloy steel)6124360.150.0081.461,520
AISI 4150 (gun barrel steel)6504800.120.0101.301,500
SA-5083 (nuclear grade)5804200.180.0121.351,500
Ti-6Al-4V (titanium)8623310.340.0120.801,660
Inconel 7184501,7000.650.0171.301,300
Al 6061-T63241140.420.0021.34582

Johnson-Cook Damage Parameters

The J-C damage model defines the equivalent strain at failure:

[ \epsilon_f = (d_1 + d_2\exp(d_3\sigma^))(1 + d_4\ln\dot{\epsilon}^)(1 + d_5T^*) ]

Materiald1d2d3d4d5
AISI 41400.151.25−0.750.0120.55
Ti-6Al-4V−0.090.25−0.500.0143.87
Inconel 7180.040.65−1.500.0100.00

Alternative Material Models

ModelWhen to UseAdvantage
Johnson-CookMost metals, high strain rateStandard, well-characterised
Zerilli-ArmstrongBCC/FCC metals separatelyPhysically based on dislocation mechanics
Bammann-Chiesa-JohnsonLarge deformation, high temperatureIncludes damage evolution
Oxley's modelSteels in machiningAnalytical, fast computation

Solver Selection

SoftwareFormulationBest ForMesh Type
Abaqus/ExplicitCEL or LagrangianBTA, guide pad contact, large deformationEC3D8RT (hex), C3D8RT
DEFORM 3DLagrangian with continuous remeshingGun drilling, chip formation, tool wearTetrahedral, automatic remeshing
Thirdwave AdvantEdgeLagrangian with adaptive meshingQuick parameter studies, cutting forcesAutomatic tetrahedral
LS-DYNAExplicit, various formulationsDynamic effects, impact, high speedHex, tetrahedral

Coupled Eulerian-Lagrangian (CEL) Method

The CEL method, implemented in Abaqus 2024, is the most advanced approach for BTA deep hole drilling simulation:

FeatureCEL Advantage
Handles extreme deformationEulerian mesh allows material to flow through fixed grid
No remeshing requiredAvoids distortion-induced termination
Combined cutting + guide padSingle model simulates both actions simultaneously
Mesh size~375,000 linear hexahedral elements (EC3D8RT)
Computation timeDays (older hardware) — reduced with mass scaling (factor 20)

Lagrangian with Continuous Remeshing

For gun drilling simulation, DEFORM 3D with continuous remeshing is the established approach:

FeatureLagrangian Advantage
Sharp material interfaceLagrangian mesh tracks free surfaces
Multiple operationsCan simulate sequential drilling steps
Tool wear modelsBuilt-in wear models (Usui, Takeyama-Murata)
Computational efficiencyLower element count than CEL

Chip Formation Simulation

Element Elimination Technique (EET)

Most FE models of deep hole drilling use the Element Elimination Technique to simulate chip formation:

StepDescription
1Elements reach J-C damage criterion
2Stress-carrying capacity reduced to near zero
3Element removed from mesh
4Chip separates from workpiece
5New elements exposed on cutting path

Chip Formation Parameters

ParameterTypical SettingEffect
Damage initiation criterionJ-C or constant strainDetermines when chip separates
Element size at cutting zone1–5 µm (graded mesh)Smaller = sharper chip geometry
Maximum element distortion0.2–0.5 (element aspect ratio)Prevents premature deletion
Friction coefficient (tool-chip)0.3–0.6 (Coulomb)Affects chip curl and temperature
Heat fraction (inelastic)0.9 (90% of plastic work → heat)Standard machining assumption

TIP

The mesh at the cutting zone must be at least 3–5× finer than the expected chip thickness. For gun drilling of steel with 0.05 mm/rev feed and a 2:1 chip thickness ratio, the minimum element size at the primary shear zone should be 3–8 µm. An overly coarse mesh will not capture the shear localisation that drives chip formation, and the simulation will produce a blocky, unrealistic chip shape.

Guide Pad Contact Modeling

The guide pads are unique to deep hole drilling and critical to bore quality. Their contact conditions must be modelled correctly:

ParameterTypical ValueWhy It Matters
Number of guide pads2 (BTA) or 1 pad + 1 support (gun drill)Determines bore support
Guide pad materialTungsten carbide or PCDHigh modulus, low friction
Contact pressure50–500 MPaDrives burnishing and bore diameter
Friction coefficient (pad-bore)0.05–0.15Burnishing + lubrication effect
Extrusion deformation10–50 µmDetermines bore size and surface finish
Contact angle90–135° (BTA)Affects self-guiding stability

Contact Modeling Approaches

MethodDescriptionApplication
Hertzian contact theoryAnalytical elastic-plastic contactQuick sizing of pad pressures
FE contact (Lagrangian)Surface-to-surface contact in Abaqus/DEFORMDetailed pressure distribution
CEL contactEulerian material interacts with pad surfaceFull process simulation

Thermal Modeling and Coolant Effects

Heat Generation and Distribution

Heat SourceFraction of Total HeatLocation
Primary shear zone (chip formation)60–70%Shear plane
Secondary shear zone (tool-chip friction)20–30%Rake face
Tertiary zone (tool-workpiece friction)5–15%Flank face, guide pads

Heat Partition

In deep hole drilling, the coolant carries away a significant fraction of the heat:

MediumHeat Carried AwayTypical Temperature
Chip60–75%500–900°C (depends on material)
Coolant15–30%30–60°C (exit temperature rise)
Workpiece5–10%50–200°C (near bore surface)
Tool2–5%200–600°C (cutting edge)

Coupled CFD-FEM Thermal Modeling

For accurate temperature prediction, CFD coolant simulation must be coupled with the FE model:

ApproachCoolant ModelingAccuracyComplexity
Convection coefficient (constant)Simple — single h valueLow (±100°C)Minimal
Convection coefficient (correlation)Nusselt-based (Re, Pr)Moderate (±50°C)Moderate
Full CFD-FEM couplingANSYS CFX / Fluent + AbaqusHigh (±20°C)High (expertise needed)

Validation Methods

Experimental Validation Techniques

Measured QuantitySensor/MethodCalibration Required
Cutting forces (Fx, Fy, Fz)Kistler dynamometer (9257B or similar)Static calibration before each test
TorqueRotary dynamometer or strain gaugeDynamic calibration
TemperatureEmbedded thermocouple (K-type) at cutting edgeIn-situ melting point verification
Temperature (field)Infrared thermography (open access zones)Emissivity calibration
Residual stressX-ray diffraction (sin²ψ method)Stress-free reference sample
Chip geometryOptical microscopyPixel-to-length calibration
Surface roughnessContact profilometerReference standard

Typical Validation Results

PaperMaterialModel TypeForce ErrorTemp Error
MDPI Simulation Study (2021)C45 steelAbaqus/Explicit EET±12%±45°C
TU Dortmund BTA CEL (2024)AISI 4140Abaqus CEL±15%±50°C
DEFORM Gun Drilling (2012)AISI 4150DEFORM Lagrangian±10%±40°C
BTA Analytical + FEM (2022)SA-5083MPL + Hertzian contact±18%N/A

Practical Simulation Workflow

StepActivityTime Required
1Define geometry (tool, workpiece) in CAD2–4 hours
2Select material model with verified J-C parameters1–2 hours
3Create mesh with graded refinement at cutting zone2–6 hours
4Define contact conditions (tool-chip, pad-bore)1–2 hours
5Set boundary conditions and loads1 hour
6Run simulation6–72 hours (depends on model size)
7Validate against experimental data4–8 hours
8Iterate parametersVariable

Mesh Configuration Guidelines

FeatureGun Drilling (DEFORM)BTA Drilling (Abaqus CEL)
Element typeTetrahedral (Lagrangian)EC3D8RT hex (Eulerian)
Element count50,000–150,000300,000–500,000
Minimum element size3–8 µm5–15 µm
Grading ratio5:1 (fine to coarse)10:1
Mass scalingNot typically usedFactor 10–20
Typical run time (24-core)6–12 hours24–72 hours

FAQ

Q: What software is best for simulating deep hole drilling? Abaqus/Explicit with CEL formulation is best for BTA deep hole drilling because it handles the large deformation of the chip and combined cutting/burnishing action. DEFORM 3D is better for gun drilling where tool wear prediction and chip morphology are the primary interests.

Q: What material model is used for deep hole drilling simulation? The Johnson-Cook constitutive model is the standard. It captures strain hardening, strain rate sensitivity, and thermal softening — the three essential effects in machining. Johnson-Cook damage parameters define when the chip separates from the workpiece.

Q: How accurate are FEM predictions for deep hole drilling? Validated FE models predict cutting forces within 10–15% of measured values, temperatures within ±50°C, and residual stress profiles that correlate well with XRD measurements. Chip morphology is predicted qualitatively — the general shape and segmentation pattern are correct, but exact dimensions may vary.

Q: Can FEM predict bore straightness in deep hole drilling? Indirectly. FEM predicts tool deflection and guide pad contact pressures, which are the primary drivers of bore straightness. The deflection data can be integrated to estimate bore deviation, but the accuracy is limited (±20%) because the drilled pilot hole and machine alignment are not typically modelled.

Q: How long does an FEM simulation of deep hole drilling take? A gun drilling simulation in DEFORM typically runs 6–12 hours on a 24-core workstation. A BTA drilling simulation with CEL in Abaqus runs 24–72 hours. Mass scaling can reduce run times by 5–10× with acceptable accuracy loss.

Q: How are guide pads modelled in FEM for deep hole drilling? Guide pads are modelled as rigid or elastic bodies in contact with the bore wall. Contact pressure, friction coefficient (0.05–0.15), and extrusion deformation (10–50 µm) are the key parameters. Modern CEL models in Abaqus can simulate the combined cutting edge and guide pad action simultaneously.

Q: Can coolant effects be included in the FEM simulation? Yes, but typically through coupled CFD-FEM analysis rather than a single model. The FE simulation predicts heat generation, and a separate CFD model predicts coolant heat transfer. The results are iterated until convergence. Simple convection coefficient approximations are faster but less accurate.

Q: What mesh size is needed at the cutting zone? Minimum element size of 3–8 µm at the primary shear zone for gun drilling, and 5–15 µm for BTA. The mesh must be 3–5× finer than the expected chip thickness to capture shear localisation. Graded meshes transition from fine at the cutting zone to coarse far from it.

Q: Is FEM simulation worth the effort for a job shop? For most job shops: no. Setting up and validating an FE model requires specialist expertise and is justified for high-value components (nuclear, aerospace), process development for difficult materials, or tool geometry optimisation. For routine production, empirical parameter tables are more practical.

Q: What is the Coupled Eulerian-Lagrangian (CEL) method? CEL uses a fixed Eulerian mesh through which material flows, combined with a Lagrangian mesh for the tool. This avoids mesh distortion problems that terminate conventional Lagrangian simulations when the chip deforms excessively. CEL is the most advanced method for BTA deep hole drilling simulation.

Deep Hole Drilling Hub — Your Trusted Third-Party Industry Resource