abaqus new function in v6.11

121
New Functions and Enhancements in V6.11E March 2011

Upload: ds-alfonso

Post on 11-Mar-2015

3.060 views

Category:

Documents


5 download

TRANSCRIPT

Page 1: Abaqus New Function in V6.11

New Functions and Enhancements in V6.11E

March 2011

Page 2: Abaqus New Function in V6.11

Overview

Abaqus/CAE

Abaqus/Standard

Abaqus/Explicit

Abaqus/CFD

2

New functions and enhancements

Page 3: Abaqus New Function in V6.11

Abaqus/CAE

Page 4: Abaqus New Function in V6.11

6.11 Enhancements

Modeling & CAD Interfaces

Meshing

Attributes & Analysis Support

Predefined Field Support

Topology and Shape Optimization

Visualization

Page 5: Abaqus New Function in V6.11

Modeling & CAD Interfaces

Page 6: Abaqus New Function in V6.11

• CATIA V5 Bidirectional Associative Interface

• CATIA parameters can be modified from Abaqus/CAE

• Model updated automatically

• Support for CATIA V5 R20

6

CAD Interfaces

CAD geometry and parameters

export to Abaqus/CAE

Updated parameters

export to CATIA V5

Page 7: Abaqus New Function in V6.11

Substructures

7

• Continuation of 6.10-EF project

• Support for:

• Substructure load cases

• Substructure load

• Improved display of retained nodal dofs

• Translucency control in part/assembly display options

• Substructure statistics query

Page 8: Abaqus New Function in V6.11

• Enhanced spline feature

• Create spline wires through points

• Define using table or import points from file

• Option to create sets

8

Modeling

Page 9: Abaqus New Function in V6.11

Meshing

Page 10: Abaqus New Function in V6.11

• New tool for partitioning faces by edge

projection

• Option to extend edges at free ends

• Elements don‟t cross boundaries between

regions with different thickness

10

Partitioning

Page 11: Abaqus New Function in V6.11

• Reduce picking needed to create mid-surface

• Improved robustness

• Offset operation performance

• Feature regeneration

• Enhanced heuristics for Extend and Blend geometry tools

• Thickness data propagates correctly with virtual topology

11

Mid-surfacing enhancements

Page 12: Abaqus New Function in V6.11

Tet meshing

• Minimum element size specification

• Tetrahedral element size growth control for interior volume

• Improved quality and robustness

• Control deviation between boundary mesh and surface geometry

• Reduced likelihood of creating short element edges

• Better gradation on surface meshes

12

Page 13: Abaqus New Function in V6.11

• New mesh edit functions

• Merge/subdivide elements

• Grow/collapse short element edges

• Bottom-up meshing

• Now available for orphan meshes

• Generate elements by offsetting

• Additional options for extrude method

13

Mesh Editing

Page 14: Abaqus New Function in V6.11

• XFEM support 2nd order tetrahedral elements

• Visualization support

• Performance improvements

• Support global reservoir modeling workflows

• Support for new coupled displacement-pore pressure

elements (C3D4P & C3D6P)

14

Miscellaneous Meshing

Page 15: Abaqus New Function in V6.11

Attributes & Analysis Support

Page 16: Abaqus New Function in V6.11

Mapping Capability

• Interface for:

• Importing spatially varying point cloud field data

• Applying data sets as loads, predefined fields and interactions

• Examples:

• Pressure, temperature and

film coefficients

• Shell thickness, density

• Permits mapping for scalar values

• Mapping options & controls

• Default value, algorithm, search tolerance

• Visualization tools planned

Page 17: Abaqus New Function in V6.11

Mapping Capability

• Import data using

• Text files & spreadsheets

• Existing Abaqus output database

• Set field output, frame, results options

• Viewport snapshot

• Use with clients

• Apply scale

• Data formats

• X, Y, Z, value

• Grid

Page 18: Abaqus New Function in V6.11

Mapped Field Clients

Shell SectionsElement and Nodal Thickness

Homogeneous, Composite Shell sections

Conventional Shell Composite Layups

Loads -- Pressure

Predefined FieldsTemperature, Pore Pressure, Void Ratio, Saturation

InteractionsSurface Film Condition, Concentrated Film Condition

Film and Sink Temperature values

Surface Radiation

Materials -- Density

Page 19: Abaqus New Function in V6.11

• Capabilities for realistic modeling of fasteners

• Create Template model

• Separate from actual analysis model.

• Contains surfaces, constraints and connectors

• Assign to a region

• Attachment points, orientations, and surfaces

specified to create an “assembled fastener”.

• Allows specification of a calibration script.

19

Assembled Fasteners

3 plate template model

Page 20: Abaqus New Function in V6.11

Predefined Field Support

Page 21: Abaqus New Function in V6.11

Predefined Field Support

Added support for existing Initial Condition

keywords as predefined fields:

*INITIAL CONDITIONS, TYPE=STRESS

*INITIAL CONDITIONS, TYPE=STRESS, GEOSTATIC

*INITIAL CONDITIONS, TYPE=PORE PRESSURE

*INITIAL CONDITIONS, TYPE=RATIO

*INITIAL CONDITIONS, TYPE=SATURATION

Page 22: Abaqus New Function in V6.11

Stress Predefined Field

Supports Direct Specification and From File (ODB)

Depending on the type of region selected the data table for stress component will change:

Page 23: Abaqus New Function in V6.11

Geostatic Stress Predefined Field

Page 24: Abaqus New Function in V6.11

Pore Pressure Predefined Field

• Can be defined using a Uniform magnitude, From File, User defined and Expression, Mapped and Discrete fields

• Supports constant or linear pressure distributions

Page 25: Abaqus New Function in V6.11

Void Ratio Predefined Field

• Can be defined using a Uniform magnitude, From File, User defined and Expression, Mapped and Discrete fields

• Supports constant or linear void ratio distributions

• Supports different distributions for each supplied ratio.

Page 26: Abaqus New Function in V6.11

Saturation Predefined Field

• Can be defined using a Uniform magnitude, or Expression, Mapped and Discrete fields

Page 27: Abaqus New Function in V6.11

Miscellaneous Enhancements (Abaqus/CFD)

Distributions for Velocity on

Inlet/Outlet and Wall Condition

BCs

Support for analytical fields

– Values are calculated using a

simplified integration scheme

at element nodes to determine

final value for each element

Enabled Keyword editor for CFD

Models

Page 28: Abaqus New Function in V6.11

Miscellaneous Enhancements

Added support for

Expression, Mapped

and Discrete Fields

for Material Density

Write amplitudes

with 16 digits of

precision to input

files

Page 29: Abaqus New Function in V6.11

Miscellaneous Enhancements

Added Expression,

Mapped and Discrete

field support for Film

condition sink

temperatures

Surface and Concentrated

Film Conditions

Can be used with

Embedded Coefficients,

Property, Analytical or

Discrete field definitions of

the film condition

Page 30: Abaqus New Function in V6.11

Topology and Shape Optimization

Page 31: Abaqus New Function in V6.11

Topology and Shape Optimization (ATOM)

• Topology optimization

• Modify stiffness

• Good for evolving optimum shape

• Shape optimization

• Moves nodes

• Good for fine tweaking of shape

Both support:

Contact

Geometric non-linearity

Nonlinear materials

• Export smoothed shape to STL or INP

31

Page 32: Abaqus New Function in V6.11

Optimization Workflow

32

ATOM is a new module in Abaqus/CAE

Specify problem

Modify .inp file

Standard

Postprocess

Final Solution ?

No

Visualize Smooth output

Export to CAD

Write .inp file

Shape or Topology

Optimization

components

Page 33: Abaqus New Function in V6.11

Visualization

Page 34: Abaqus New Function in V6.11

• Contour plots on beam sections

• Available for Box, Rectangle, Circle, Pipe, I and L sections

• New „BEAM_STRESS‟ field output variable

• SF and SM required

• View cuts enabled with beam profile rendering

34

Visualization

Page 35: Abaqus New Function in V6.11

• FBD enhancements

• Section force/moment history output

35

Visualization

Page 36: Abaqus New Function in V6.11

• FBD enhancements

• Section force/moment display on multiple view cuts

• Multiple free bodies on a single view cut

36

Visualization

Page 37: Abaqus New Function in V6.11

• Multi-point constraints visualization

• Display probed node/element labels and values

• Particle (PC3D) elements display

37

Visualization

Page 38: Abaqus New Function in V6.11

Abaqus/Standard

Page 39: Abaqus New Function in V6.11

6.11 Enhancements

Parallel Cavity Radiation

XFEM V – Fracture and Failure Enhancements

Coupled Electrical-Thermal-Structural Analysis (ETS)

Electromagnetics I - Low frequency (Eddy current)

AMS Solver

GPGPU support in the direct solver

Contact

Page 40: Abaqus New Function in V6.11

Parallel Cavity Radiation

Page 41: Abaqus New Function in V6.11

Parallel Cavity Radiation

Goals:

Enable cavities larger than v6.10-EF limit 16,000 facets

Provide a parallel and scalable cavity radiation feature

Approach:

Solve same equations as old serial code (gray, diffuse, surface

radiation)

Avoid inversion of large matrix by solving problem with iterative

solver

(New algorithm - only available in 64-bit platforms)

Page 42: Abaqus New Function in V6.11

Parallel Cavity Radiation

Serial cavity radiation: 4 minutes(v6.11 , 8 cpu)

Parallel cavity radiation: 33 seconds (v6.11, 8 cpu)

Example problem: exhaust manifold (4,500 facets)

8x speed-up

Page 43: Abaqus New Function in V6.11

Parallel Cavity Radiation

Serial cavity radiation limit: 16,000 facets

Limitation due to internal 2GB limit for element storage.

Largest cavity radiation run to date: 128,000 facets

128,000 facets ran in 128 cpus.

Job completed in 63 iterations (38 minutes).

In serial mode (if possible) would require 131GB on a single

node to run.

Very large cavities

Page 44: Abaqus New Function in V6.11

XFEM V - Fracture & Failure Enhancements

Page 45: Abaqus New Function in V6.11

XFEM

Continued enhancements to XFEM:Functionality

Support 2nd order tets (C3D10 and C3D10H) with XFEM

Output strain energy release rate (ENRRTXFEM) for XFEM based LEFM approach

Use cases/drivers

1st order tets and 1st order bricks with XFEM were supported since 6.9

Most real engineering structures in Auto, Aerospace and medical applications are made of 2nd order tetrahedron elements.

Direct requests from Boeing, Daimler, Dana etc.

Go mainstream with XFEM

Crack path/surface is more stable with 2nd order tet elements

Usage

No user interface changed

Can be performed in a static procedure, implicit dynamic procedure and the low cycle fatigue analysis.

Page 46: Abaqus New Function in V6.11

XFEM – Quadratic Tets

Page 47: Abaqus New Function in V6.11

Coupled Electrical-Thermal-Structural

Procedure (ETS)

Page 48: Abaqus New Function in V6.11

ETS Overview

Fully coupled three-field analysis with the following fields

Electrical Potential (Steady-State)

Temperature (Transient or Steady-State)

Displacement (Steady-State)

Three-dimensional continuum elements

Q3D4 4-node tetrahedron, linear displacement, electric potential, and temperature.

Q3D6 6-node triangular prism, linear displacement, electric potential, and

temperature.

Q3D8(RH) 8-node hexahedral, tri-linear displacement, electric potential, and

temperature.

Q3D10M(H) 10-node tetrahedron, modified displacement, electric potential, and

temperature.

Q3D20(RH) 20-node hexahedral, tri-quadratic displacement, tri-linear electric

potential, and temperature.

New keyword interface

*COUPLED TEMPERATURE-DISPLACEMENT, ELECTRICAL

Page 49: Abaqus New Function in V6.11

Spot Weld Example

Q3D8 ElementsQuarter-symmetry model

Page 50: Abaqus New Function in V6.11

Spot Weld Example, cont

Page 51: Abaqus New Function in V6.11

Electromagnetics I

Low frequency (Eddy current)

Page 52: Abaqus New Function in V6.11

Low Frequency Time-Harmonic Procedure

Time-Harmonic response for a given cyclic current

excitation

Same as direct steady state dynamic procedure in structures

Compute EM wave response in both air and conducting media

Use cases/drivers

EM Time-harmonic followed by thermal / mechanical / thermo-

mechanical analyses by transferring results (applications: EM

induction heating and/or forming)

Theory

Neglect high frequency term

Page 53: Abaqus New Function in V6.11

Sample Results

Long annular cylinder in an oscillating uniform

magnetic fieldMagnetic induction in vertical direction compared against that of

benchmark results

Page 54: Abaqus New Function in V6.11

Sample Results

Multiple conductors in uniform oscillating magnetic

fieldTime-harmonic electric field in air and conductors

Time-harmonic magnetic induction (flux density) in conductors

Electric field Magnetic induction

Page 55: Abaqus New Function in V6.11

Sample Results

Eddy fields in a spherical conductor sitting inside a

magnetic field

In phase and out of phase (curling) electric fields

In phase and out of phase electric field magnitudes

Page 56: Abaqus New Function in V6.11

Magnetic Field in a Straight Conductor

Current flow to z direction in the conductor

Page 57: Abaqus New Function in V6.11

Magnetic Field of Two Straight Conductors

Current direction is the same in the two conductors.Zero D EM Potential is applied on the outer surfaces.

Magnetic FieldElectric Field

Page 58: Abaqus New Function in V6.11

Conducting infinite cylinder in a uniform time harmonic

magnetic field

Conductor

Surface current is applied on the outer surface

Magnetic Field

Electric FieldElectric Field

Page 59: Abaqus New Function in V6.11

AMS Solver

Page 60: Abaqus New Function in V6.11

AMS Performance Improvements

60

4.3M DOF Powertrain Model:

4500Hz cutoff frequency, 1709 modes, selective recovery (167,618 dofs)

on Intel Nehalem-EX with128GB memory

0

500

1,000

1,500

2,000

2,500

3,000

3,500

1 4 8 16

Elap

sed

Tim

e (

sec.

)

Number of Cores

FREQ (6.10-EF)

AMS (6.10-EF)

FREQ (6.11)

AMS (6.11)

Page 61: Abaqus New Function in V6.11

AMS Performance Improvements

61

4.3M DOF Powertrain Model:

4500Hz cutoff frequency, 1709 modes, selective recovery (167,618 dofs)

on Intel Nehalem-EX with128GB memory

0.00

1.00

2.00

3.00

4.00

5.00

6.00

7.00

8.00

9.00

10.00

1 4 8 16

Spe

ed

up

Number of Cores

FREQ (6.10-EF)

AMS (6.10-EF)

FREQ (6.11)

AMS (6.11)

Page 62: Abaqus New Function in V6.11

AMS Performance Improvements

62

Neon 2M DOF Vehicle Body Model:

600Hz cutoff frequency, 3064 eigenmodes, selective recovery, 2000

residual vectors on Intel Nehalem-EP with 32GB memory

74

5861

70

55 58

47

2014

44

1711

0

10

20

30

40

50

60

70

80

1 4 8

Elap

sed

Tim

e (

min

.)

Number of Cores

FREQ (6.10-EF)

AMS (6.10-EF)

FREQ (6.11)

AMS (6.11)

Page 63: Abaqus New Function in V6.11

GPGPU Accelerated Direct Solver

Page 64: Abaqus New Function in V6.11

Performance targets and required hardware

The basic target was to achieve an overall speedup

for Abaqus/Standard (standard.exe wall time) of 2x

versus a 4 core cpu only time for our benchmark

model s4b.

Requires compute specific GPGPU

NVIDIA Tesla C2050, C2070 GPU Computing Processor

Support for compute cards from AMD is expected by release.

The solver can be run on lesser cards, but performance

expectations must be reduced

Optimal performance requires the factorization to

remain in core

Page 65: Abaqus New Function in V6.11

Performance data

4.34E+11 4.45E+11 6.59E+11 9.90E+11 1.91E+12 2.19E+12 4.37E+12 5.76E+12 1.03E+13 1.68E+13 1.70E+13 2.63E+13 1.08E+14

Solver 1.83 1.37 1.69 1.82 2.32 2.15 2.20 2.75 2.90 3.38 3.69 3.36 1.96

Standard 1.14 1.02 1.16 0.83 1.52 1.52 1.49 2.19 2.01 2.00 2.29 2.66 1.79

0.00

0.50

1.00

1.50

2.00

2.50

3.00

3.50

4.00

4 c

ore

/ (

4 c

ore

+ g

pu

)

GPGPU speedup (4 core / (4 core + gpu))

Page 66: Abaqus New Function in V6.11

Contact

Page 67: Abaqus New Function in V6.11

Contact stress error indicators

Provide some perspective on accuracy of contact stresses

Hertz contact example

Analytical solution

Abaqussolutions

Error indicators

Maximum contact pressure

Position

Pre

ssu

re Maximum error

indicator

• Tend to be large where local variation of base variable is more complex than what can be captured by the mesh

• Not normalized; same units as base variable• Not conservative or precise estimates of error

• Points to remember for error indicators:

Page 68: Abaqus New Function in V6.11

Contact stress error indicators

Recall improvements for 2nd-order elements in 6.10EF

Less noise

Less noise

Prior versions

6.10EF & 6.11

Prior versions

6.10EF & 6.10

Error indicator

6.11

Error indicator

6.11

• Accurate prediction of maximum CPRESS• Some uncertainty where gradient is large but

pressure is low

• Need finer mesh to predict maximum contact pressure

Page 69: Abaqus New Function in V6.11

Contact stress error

indicatorsTwo deformable blocks

Error indicator does not show evidence of inaccuracy if mesh is too coarse!

p=1

p=1

Edge effects evident after mesh refinement

One element per block

Page 70: Abaqus New Function in V6.11

Edge-to-surface contact

Supplementary edge-to-surface formulation for general contact

Targets situations in which the active contact zone in a numerical

model corresponds to a line associated with a feature edge

Whereas the surface-to-surface formulation best treats contact

over a finite area

General contact with S-to-S formulationGeneral contact with S-to-S and

E-to-S formulations

• Future: add edge-to-edge formulation

• Diverges 25% into simulation • Penetration near feature edge• 36 increments; 317 iterations

• Runs to completion• Good resolution of contact• 28 increments; 130 iterations

Page 71: Abaqus New Function in V6.11

Edge-to-surface contact

Tests featuring edge contact

Two views of same analysis

Page 72: Abaqus New Function in V6.11

Abaqus/Explicit

Page 73: Abaqus New Function in V6.11

6.11 Enhancements

Eulerian Heat Transfer Element

Mass Adjust

Subcycling improvements

Additional AQUA Wave types

SPH

Page 74: Abaqus New Function in V6.11

Eulerian Heat Transfer Element

Page 75: Abaqus New Function in V6.11

Eulerian Heat Transfer Element

New Element Type: EC3D8RT

8-node thermally coupled linear multi-material Eulerian brick

Active degrees of freedom 1,2,3,11

Temperature calculated as part of the fully-coupled problem

Procedure:

*Dynamic Temperature-Displacement, Explicit

Mechanical loads valid for EC3D8R also supported for

EC3D8RT

Thermal loads that are supported:

*Dflux, *Film, *Radiate

*Dsflux, *Sfilm, *Sradiate

Page 76: Abaqus New Function in V6.11

Eulerian Heat Transfer Element

Example 1:Deep indentation problem (white area is initially void )

temperature is fixed at the bottom; the heat source is the plastic dissipation

ran with 2 cpus and 4 domains

Page 77: Abaqus New Function in V6.11

Eulerian Heat Transfer Element

Example 2:Eulerian heat transfer in a progressively filled block

Material flows in at 100, Film condition on the lefts side (1st step);

additional film condition at the top in the 2nd step

Page 78: Abaqus New Function in V6.11

Eulerian Heat Transfer Element

Example 3:Rivet Forming –changed existing example problem to use EC3D8RT instead of EC3D8R in a *Dynamic

Temperature-Displacement analysis

Initial temperature at 20, will increase due to plastic work effects

Run with 8 cpus and 16 domains

Page 79: Abaqus New Function in V6.11

Mass Adjust

Specify a target mass with optional target

time increment

Page 80: Abaqus New Function in V6.11

*MASS ADJUST

You can specify a target mass or a “trim level” for an

ELSET.

Abaqus/Explicit will adjust the mass up or down to meet

the target.

You can further redistribute the mass within the ELSET to

raise the stable time increment to a target value.

You can even redistribute just the “current” mass to

achieve the target time increment without adding extra

pounds!!!

Page 81: Abaqus New Function in V6.11

Verification

81

1. Verify specified mass: test and reference elements have different densities

but the element set masses are adjusted to be the same. The dynamic

responses are therefore similar.

Test Reference

Page 82: Abaqus New Function in V6.11

Subcycling improvements

Page 83: Abaqus New Function in V6.11

Subcycling robustness improvements

Allows different time increments to be used for different groups of elements

Reduces run time for an analysis when a small region of elements (subcycling zone) in the model controls the stable time increment

Keyword Interface to define a subcycling zone

*SUBCYCLING, ELSET=element_set_name

Subcycling feature available in Abaqus/Explicit version v6.8-EFHowever functionality not robust enough.

Analysis becomes unstable/shows unphysical behavior in the subcyclingzone

Energy balance not achieved when general contact is defined

Page 84: Abaqus New Function in V6.11

Subcycling robustness improvements

Page 85: Abaqus New Function in V6.11

Subcycling robustness improvements

No Subcycling 20 h 42 m (8 cpus)

Subcycling 9h 23 m (8 cpus)

Number of elements in subcycling zone = 29623in non-subcycling zone = 508977

Subcycling ratio = 6

Page 86: Abaqus New Function in V6.11

Additional AQUA Wave types

Page 87: Abaqus New Function in V6.11

AQUA Waves that generate loads on structures

Waves acting on structures generate loads such as

buoyancy, drag, and inertial loads.

For A/Explicit rel6-10ef, the wave definition was

limited to only the 5th Order Stokes wave

formulation.

For A/Explicit 6.11, You could also use the Airy wave

formulation or provide a more general definition

through VWAVE user-subroutine.

87

Additional wave formulations in A/Explicit

Page 88: Abaqus New Function in V6.11

SPH (Smoothed Particle Hydrodynamics)

in Abaqus/Explicit

Page 89: Abaqus New Function in V6.11

Water Splash In a Square Pan

89

100 K particles, 53 particles/element

86 mins on a PC

5770 incs, MSPEI 8.6

Support is more local in this case as the number of

particles per element is almost the same

Page 90: Abaqus New Function in V6.11

Bird Fan Blade Slashing

A cylindrical bird strikes an initially straight edge of a rotating turbofan blade

The blade deforms and the bird disintegrates

Contour plots of pressure shown

90

4.2 K particles

47 to10 particles/element

0:47 mins on a PC

2200 incs, MSPEI 5.6

EOS material with tensile failure

Elasto-plastic blade

Page 91: Abaqus New Function in V6.11

Wave Impact

A block of water falls under gravity (dam rupture)

Velocity vector plots on the left

91

220 K particles, 32 particles/element

140 hours (not sure which machine)

117K incs, MSPEI 20.2

Tabular EOS with tensile failure

Page 92: Abaqus New Function in V6.11

Wave Impact

A block of water falls under gravity (dam rupture) Simplified boulders are being

Velocity vector plots on the left

92

220 K particles, 32 particles/element

140 hours (on storm)

117K incs, MSPEI 14 (before the latest improvement); expect 7

Tabular EOS with tensile failure

Page 93: Abaqus New Function in V6.11

Water Splashing of a Figurehead

A block of water hits an object

USUP EOS

53K particles

93

*eos, type=usup1500e+3,0,0*tensile failure,element deletion=no,pressure=ductile,shear=ductile2.0*viscosity1.0e-8*density1.e-9

Page 94: Abaqus New Function in V6.11

Priming a Pump

94

182 K particles, ? particles/element

89 hours (not sure which machine)

150K incs, MSPEI 11

Tabular EOS with tensile failure

A block of water is pushed by a piston while the pump

is rotating

Page 95: Abaqus New Function in V6.11

Bottle Drop

A filled water bottle gets dropped on the floor

Comparison with CEL

Differences in sloshing

Similar deformed shapes for the plastic bottle

95

12 K particles

CEL: 2h52m; SPH: 1h48m

Page 96: Abaqus New Function in V6.11

Smashing of a Figurehead

Figurehead with initial velocities is smashed into a

wall

Tooth paste like material (from our example manual)

96

8.2 K particles, 40 particles/element

90 mins (not sure which machine)

120K incs, MSPEI 6

Tabular EOS with tensile failure

Page 97: Abaqus New Function in V6.11

Smashing of a Figurehead

With cohesive contact

97

53 K particles, 47 particles/element

90 mins (not sure which machine)

360K incs, MSPEI 6

Tabular EOS with tensile failure?

Page 98: Abaqus New Function in V6.11

Ball Drop in Water Tank

A relatively light ball falls into a tank of water

SPH vs CEL

Initial velocity

Hard to tell what‟s going on in SPH

98

Page 99: Abaqus New Function in V6.11

Taylor Test

A perfectly plastic cylindrical copper bar is impacting a rigid wall84K particles

Finger pattern develops and some particles fly offLikely due to the mesh being non-uniform to start with

Tensile instability could also be the issue

99

Page 100: Abaqus New Function in V6.11

Taylor Test – comparison with CEL and C3D8R

100

Stress contour plots match OK at various stages

during the analysis

Left: CEL

Center: C3D8R

Right: SPH

Page 101: Abaqus New Function in V6.11

Taylor Test – reaction forces

101

C3D8R and CEL match well

All curves are unfiltered

Default options in all three cases

Page 102: Abaqus New Function in V6.11

Garden hose: pressurization + spraying

Very high number of increments

Is the EBE DT excessively conservative?

102

27 K particles, 54 particles/element

2375K increments, 104 hours CPU time

MSPEI 5.3

Page 103: Abaqus New Function in V6.11

Projectile Impact on Plate – slow bullet

A cylindrical rigid projectile impacts a steel plate

Properties:rate dependent hardening + damage initiation (ductile and shear) with energy based evolution

Friction 0.3 between the bullet and the plate

The circular particle patch in the center is TIE-ed to the FE plate

103

V = 500m/sec, T=0.5 msec

Bullet gets stuck in the hole

103 K particles, 44 particles/element

9K increments, 1h48mins, MSPEI 6.3

Page 104: Abaqus New Function in V6.11

Projectile Impact on Plate – fast bullet

Whole analysis shown on the left

Slower motion of the perforation shown on the right

104

V = 1000m/sec, T=0.2 msec

Bullet perforates

103 K particles, 44 particles/element

4K increments, 47mins, MSPEI 6.1

Page 105: Abaqus New Function in V6.11

Performance – 3rd model

Taylor test: comparison with C3D8R and CELMaterial: Perfectly plastic copper, no damage

All analyses ran using 1 CPU on a lnx86_64 v6Intel machine (gladius)

Old results: for this size model the SPH analyses are probably 30% faster.

105

Nr ofElements

Nr of nodes per element

DT stable MSPEI Nr of increments

Total CPU time

CEL 585K 8 2e-08 constant

3 to 4 3864 146:25 mins

C3D8R 78K 8 1.7e-8to1.2e-9

1.15 33188 49:28 mins

SPH_A 84K 39 8.2e-9constant

6.7 9783 93:35 mins

SPH_B 84K 17 6.8e-9constant

4.2 61:15 mins

SPH_B

Page 106: Abaqus New Function in V6.11

Abaqus/CFD

Page 107: Abaqus New Function in V6.11

Overview of Enhancements in 6-11

Keyword support and documentation

Surface output variables

RNG k- model improvements

Improved robustness

“Resolution-insensitive” wall functions

Temperature-dependent properties

Improved co-simulation job submission

Page 108: Abaqus New Function in V6.11

Keyword support and documentation

Page 109: Abaqus New Function in V6.11

Input file usage with documented keywords are provided in 6-11. For example,

*CFD

*Momentum Equation Solver

*Transport Equation Solver

*Pressure Equation Solver

*Turbulence Model

*Fluid Boundary

*Surface Output

The required and optional parameters for all keywords are documented in 6-11 Abaqus Keywords Reference Manual.

Documented Keywords

109

Abaqus/CFD - 6.11 Enhancements

Page 110: Abaqus New Function in V6.11

Surface Output Variables

Page 111: Abaqus New Function in V6.11

Surface Output Quantities

111

Abaqus/CFD - 6.11 Enhancements

Scalar Quantities Vector Quantities

Field History Field History

y+

(YPLUS)Also defined for laminar flows

Mass Flow Rate

(MASSFLOW)

Surface Traction vector

(STRACTION)

Int. Traction (Forces)

(FORCE)

Surface Output Variables

Wall Shear Stress

(WALLSHEAR)

Normal Heat Flux

(HFLN)

Volume Flow Rate

(VOLFLOW)

Int. Heat Flux

(HEATFLOW) Optional Output:

Pressure, viscous forces

( PRESSFORCE,

VISCFORCE)

Heat Flux vector

(HFL)

Total Traction vector

(TRACTION)

Normal Traction vector

(NTRACTION)

y*

(YSTAR)Defined only for “k”-family models

Area

(SURFAREA)

Average Temperature

(AVGTEMP)

Area Average Velocity

(AVGVEL)

Average Pressure

(AVGPRESS)

Page 112: Abaqus New Function in V6.11

112

Abaqus/CFD - 6.11 Enhancements

Surface Output Examples

Wall shear stress contours Surface traction vectors superimposed on pressure contours

Aortic Aneurysm

Drag Force Lift Force

Velocity contours

Vortex Shedding behind a cylinder

Flow around a cylinder y+ plot

Page 113: Abaqus New Function in V6.11

RNG k- model improvements

Page 114: Abaqus New Function in V6.11

Abaqus/CFD - 6.11 Enhancements

k- model is subject to spurious overproduction of k

in highly strained flows („stagnation point anomaly‟)

Based on experimental evidence in shear layers and

on mathematical grounds, the turbulent eddy

viscosity is limited using an upper bound

This method has shown to significantly improve the

stability of the k- model for highly strained flows

114

Improved robustness – Time scale limiters

Page 115: Abaqus New Function in V6.11

Temperature-dependent properties

Page 116: Abaqus New Function in V6.11

Abaqus/CFD - 6.11 Enhancements

116

Temperature Dependent Viscosity

622.1780,,)( 2

9896.12

11

2

CeCeCT T

C

Channel with isothermal walls at 800 C and inlet fluid enters at a constant velocity and temperature of 200 C.

Page 117: Abaqus New Function in V6.11

Improved co-simulation job submission

Page 118: Abaqus New Function in V6.11

Abaqus/CFD - 6.11 Enhancements

Single command line job submissionNo port numbers necessary

abaqus -cosimulation cosim_job -job cfd_job,std_job -cpus 8,2

Queue submission of co-simulation jobs

Without restart

abaqus -cosimulation cosim_job -job cfd_job,std_job -cpus 8,2 -queue general_mem5_wall60

With restart

abaqus -cosimulation cosim_job -job cfd_job,std_job -oldjob cfd_job_old,std_job_old –cpus 8,2 -queue general_mem5_wall60

Higher cpu assignment flexibility : specification of cpu ratio

Specification of cpu ratio

abaqus -cosimulation cosim_job -job cfd_job,std_job –cpus 10 –cpuratios 0.8,0.2 -queue general_mem5_wall60

118

Improved co-simulation job submission

Page 119: Abaqus New Function in V6.11

Example – Electronic Cooling

119

Application Examples

Application

• Thermal performance of electronic

components and systems

Motivation

• Miniaturization of devices,

superior performance,

higher reliability and lower

cost

Approach

• Full system structural and

thermal analysis in conjunction

with natural/forced convection

cooling

Chip

Chips

Heat Sink

Capacitors

PCBPower Source

• Non-linear system-level simulation

• Linear simulation – Shock & vibration, Linear dynamics

• Fracture & failure (Cohesive, XFEM)

• Advanced materials and elements

• Complete set up in Abaqus/CAE

• CFD volume mesh from structural model

Extract skin using Shell From Solid feature

Remove faces and cover open faces to create a closed enclosure

Create volume using Solid From Shell feature and mesh the volume

• Conjugate heat transfer analysis

• Abaqus/Standard & Abaqus /CFD

Temperature contours

Temperature isosurfaces

Velocity vectors on intermediate plane

• Sequential thermal-stress analysis

• Abaqus/Standard

Mises Stress contours

Page 120: Abaqus New Function in V6.11

Top 3

Abaqus/CAE

ATOM

Abaqus/Standard

Electromagnetics

Abaqus/Explicit

SPH

Abaqus/CFD

120

Page 121: Abaqus New Function in V6.11

New Functions and Enhancements in V6.11E

March 2011