Structural Mechanics 145

Post on 29-Nov-2015

251 Views

Category:

Documents

8 Downloads

Preview:

Click to see full reader

Transcript

© 2011 ANSYS, Inc. February 8, 2013 1

Advances in ANSYS R14.5 Structural Mechanics Solutions

Mai DOAN

Jiaping ZHANG

© 2011 ANSYS, Inc. February 8, 2013 2

• Large and Complex Models – Tree Filtering / Tagging Tree Objects / Connections Worksheet /

Object Generator / Sub modeling / Performance enhancement

• Customization – Acoustics ACT

• Linear Dynamics

• Material models

• Cyclic & Linear Periodic symmetry

• Fracture Mechanics

• Modeling Composite

• Additional Control for Contact Modeling

• Improved External Data Mapping

• Miscellaneous features in ANSYS Mechanical

Agenda

© 2011 ANSYS, Inc. February 8, 2013 3

ANSYS Release 14.5 Structural Mechanics Highlights: Part 1

Mai Doan ANSYS Senior Application Engineer, Austin

Mai.Doan@ansys.com

© 2011 ANSYS, Inc. February 8, 2013 4

Analysis of Large and Complex Models

© 2011 ANSYS, Inc. February 8, 2013 5

Motivation

Many users deal with always larger models. Increased performance and usability is required at each step of an analysis in order to reduce the time from geometry to results. Various tools and techniques for efficiently managing large models.

• Tree Filtering • Tagging Tree Objects • Connections Worksheet • Object Generator • Submodeling • Performance enhancement

© 2011 ANSYS, Inc. February 8, 2013 6

Tree Filtering

The Tree Filter:

• To streamline the management of tree objects, especially for large models, a tree filter can be used to find objects or to reduce the length of the tree that is displayed.

• Filter by: Name, Tag (discussed later), Type or State.

• Examples:

Filter by state = underdefined to find only those objects in the

tree.

Filter by type = results to only display result items.

© 2011 ANSYS, Inc. February 8, 2013 7

Tagging

Tags are designations that can be assigned to any branch in the tree and used with the previously discussed filtering capability.

• Activate the tags icon from the “View” menu

With the tags view active, a new “tags” window is

displayed below the details window

© 2011 ANSYS, Inc. February 8, 2013 8

Connections Worksheet The “Worksheet” for the connections branches (e.g. “Connections”, “Contacts”, “Joints”, etc.), now adds a preferences section to allow users to configure the worksheet display.

© 2011 ANSYS, Inc. February 8, 2013 9

Connections Worksheet

New at v 14.5 is the Connection

Matrix

© 2011 ANSYS, Inc. February 8, 2013 10

Object Generator The Object Generator uses an existing object in the tree as a template for replication.

Example: mesh size control on 1 of the holes.

• Highlight the tree object to be replicated (the Face Sizing)

and activate the Object Generator icon

• The object generator window opens and indicates objects will be

generated from the “Current Selection” (graphical selection)

© 2011 ANSYS, Inc. February 8, 2013 11

Object Generator

After generating, the new objects are displayed graphically and they appear in the tree. The object names reflect that of the parent object along with the prefix chosen (if any).

In this example the object to be generated (the mesh size control) required only a single geometry selection, a face. Object generator can be used to apply BC where multiple scoping selections are required (e.g. contacts, bolt pretensions, beam connections, etc.).

© 2011 ANSYS, Inc. February 8, 2013 12

Global Model Sub Model

Submodeling is a technique where a coarsely meshed model can be solved followed by a subsequent solution using only a portion of the coarse model with a more refined mesh. Submodeling is available for structural and thermal analysis types with solid geometry.

• Data can be transferred from a 2D coarse model [Full Model] to a 3D submodel.

Submodeling for local refinement

© 2011 ANSYS, Inc. February 8, 2013 13

Performance enhancements. HPC and GPU

•Linux workstation : Two Intel Xeon X5677 processors (3.47 GHz, 8 cores total), 48 GB RAM, 2 NVIDIA Tesla C2075, Red Hat EL 6.1

2.1M DOF, Nonlinear Static Analysis

The sparse solver can take advantage of multiple GPUs to reduce solution time

© 2011 ANSYS, Inc. February 8, 2013 14

Result file size reduction

Sample model: 778,400 10-node tet, 16,800 contact elements, 1,429,000 nodes (4,3 MDOF)

To reduce storage needs,

R14.5 stores single

precision element results.

This leads to result files up to 50% smaller

© 2011 ANSYS, Inc. February 8, 2013 15

Faster post-processing

To save time when creating animations, the post-processing of models containing large number of bodies and elements has been improved to be up to 40% faster

© 2011 ANSYS, Inc. February 8, 2013 16

Customization for Structural Mechanics ACT

© 2011 ANSYS, Inc. February 8, 2013 17

Motivation

Companies want to include best practices in their simulation environment so as to provide the ability to use and re-use expert knowledge, beyond standard capabilities

© 2011 ANSYS, Inc. February 8, 2013 18

© 2011 ANSYS, Inc. February 8, 2013 19

Adding loads and boundary conditions

New loads and boundary conditions can be added to ANSYS Mechanical – for example acoustics boundary conditions

© 2011 ANSYS, Inc. February 8, 2013 20

Add custom results

Create your own custom results such as a maximum admissible criteria based on a ratio of stresses vs a given material property (Rm…)

© 2011 ANSYS, Inc. February 8, 2013 21

Encapsulate APDL scripts

APDL

! APDL_script_for_convection.inp ! Input parameters: esel,s,type,,10 cm,component,ELEM thickness = 0.005 film_coefficient = 200. temperature = 226.85 ! Treatment: /prep7 et,100,152 keyop,100,8,2. et,1001,131 keyo,1001,3,2 sectype,1001,shell secdata,thickness,10 secoff,mid cmsel,s,component emodif,all,type,1001 emodif,all,secnum,1001 type,100 esurf fini alls /solu esel,s,type,,100 nsle sf,all,conv,film_coefficient,temperature alls

© 2011 ANSYS, Inc. February 8, 2013 22

Embed 3rd party solutions

External solvers or tools can be plugged into the Mechanical Interface

© 2011 ANSYS, Inc. February 8, 2013 23

XML/Python based interfaces

XML File to create buttons and dialogs

Python scripts takes data from the tree and sends appropriate information to the solver

© 2011 ANSYS, Inc. February 8, 2013 24

Now available!

ACT-based Extensions library : Customization forum enabling customers to upload/share their own Extensions ACT-based extensions posted on ANSYS customer portal:

→Acoustics →FE Info →Morphing →Advanced post-processing →Post-processing for beams…

© 2011 ANSYS, Inc. February 8, 2013 25

ACT Acoustics extension can be used to create acoustic boundary conditions and defining fluid bodies (elements & material properties):

• Pressure

• Mass Flow

• Impedance

• FSI Interface

• …many more

• Fluid body (FLUID30, 220 or 221

elements according to the mesh)

ACT Acoustics Extension

© 2011 ANSYS, Inc. February 8, 2013 26

Acoustic Body:

• Transform SOLID185 in FLUID30, SOLID186 in FLUID220 and SOLID187 in FLUID221 (EMODIF)

• Define Sound Velocity, Density and Dynamic Viscosity… (MP)

– Material properties can be temperature dependent

• Define Perforated material

• Define Reference Pressure (R)

• Define acoustic-structural coupling

• Define PML and PML options

ACT Acoustics Extension

© 2011 ANSYS, Inc. February 8, 2013 27

ANSYS Linear Dynamics

© 2011 ANSYS, Inc. February 8, 2013 28

Structural and Material-Dependent Damping

• Damping, both structural and material-dependent, includes: Alpha Damping, Beta Damping, and Constant Damping Ratio

Enhancements include: – Damping parameterizations

– Better analysis environment state checks

– Better warning/error messages alerting users whenever appropriate

– Better visibility handling of applicable damping with respect to analysis type and solution type

© 2011 ANSYS, Inc. February 8, 2013 29

Harmonic solution time can be improved by using frequency sweep capability of VT

VT is only supported for FULL harmonic analysis, only available when FULL is selected in solution method

Harmonic solutions at different frequencies of VT is calculated based on one FULL harmonic solution

VT Harmonic Analysis

© 2011 ANSYS, Inc. February 8, 2013 30 © 2012 ANSYS, Inc. February 8, 2013

Material Models added in Workbench R14.5

© 2011 ANSYS, Inc. February 8, 2013 31 © 2012 ANSYS, Inc.

• Enthalpy input in Transient Thermal analysis (Heat Transfer analysis)

February 8, 2013

• Shape Memory Alloy (SMA) Material Model a) This model is always available for Static Structural and Transient Structural analyses. b) Available for Modal, Linear Buckling, Random Vibration and Response Spectrum analyses when performed using linear perturbation. c) Two properties are added here i) Superelasticity ii) Shape Memory Effect

© 2011 ANSYS, Inc. February 8, 2013 32 © 2012 ANSYS, Inc.

February 8, 2013

Hyperelastic Material Models These new models are available now for Static Structural and Transient Structural analyses.

• Ogden Compressible Foam Hyperelastic Material model i) Ogden Foam 1st Order

ii) Ogden Foam 2nd Order and iii) Ogden Foam 3rd Order are available.

• Extended Tube Hyperelastic Material model • Special Hyperelastic Material Model : Mullins Effect Mullins Effect is an extension to existing nearly and fully incompressible

hyperelastic models to model softening or damage.

© 2011 ANSYS, Inc. February 8, 2013 33 © 2012 ANSYS, Inc.

February 8, 2013

Damage Model Two properties are supported for Static Structural and Transient Structural analyses. i) Damage Initiation Criteria ii) Damage Evolution Law

© 2011 ANSYS, Inc. February 8, 2013 34 © 2012 ANSYS, Inc. February 8, 2013

Cyclic and Linear Periodic Symmetry

© 2011 ANSYS, Inc. February 8, 2013 35

• Through R14, large cyclic symmetry models would either post very slowly due to excessive use of memory and computation cycles during result expansion

• R14.5 enables the user to expand a fraction of the model

Cyclic symmetry

n = Controls the extent of expansion n Є [1, N] n = 0 for Program Controlled.

i = Controls display range (i, i + n – 1) i Є [1, N] i = 0 for Program Controlled

© 2011 ANSYS, Inc. February 8, 2013 36

Example Expansions for N=13

4 sectors 7 sectors 13 sectors

You can now display and animate results over a fraction of the total symmetry expansions and incur a corresponding fraction of the memory and computational expense.

© 2011 ANSYS, Inc. February 8, 2013 37

New Symmetry type created to simulate linear periodic models with structural translational symmetry in Mechanical using Symmetry Region object. The Low and High Boundaries of that object assumed to be identical, and they will be automatically meshed using Mesh Matched controls. During solution those matched nodes are coupled together using MAPDL CP command.

Linear Periodic Symmetry

The existed Symmetry Region object with new Symmetry type Linear Periodic shown together with old Symmetric and Anti- Symmetric types used before.

© 2011 ANSYS, Inc. February 8, 2013 38

Linear Periodic Symmetry

Simple model with translational Symmetry which is described by Symmetry Region object with new Linear Periodic type is defined below.

© 2011 ANSYS, Inc. February 8, 2013 39

Total Deformation Solution Results

The Total Deformation results are presented for two cases: first solution was done with new Symmetry object and the second one – without it. Symmetry exists No Symmetry

© 2011 ANSYS, Inc. February 8, 2013 40

Post-processing with Symmetry Expansion

© 2011 ANSYS, Inc. February 8, 2013 41

ANSYS Release 14.5 Structural Mechanics Highlights: Part 2

Jiaping Zhang ANSYS Structural Technical Support, Houston

Jiaping.Zhang@ansys.com

© 2011 ANSYS, Inc. February 8, 2013 42

Outline

o Fracture Mechanics

o Modeling Composite

o Additional Control in Contact Modeling

o Improved External Data Mapping

o Miscellaneous features in ANSYS Mechanical

© 2011 ANSYS, Inc. February 8, 2013 43

Outline

o Fracture Mechanics

© 2011 ANSYS, Inc. February 8, 2013 44

Motivation Cracks and flaws can occur for several reasons:

• The material may be inherently defective

• Long term operation

• Environmental conditions

Static fracture mechanics calculation could be used to

• Determine maximum allowable flaw size

• Basic for fatigue crack growth

Crack in gas pipeline due to corrosion

Crack in welded joint

© 2011 ANSYS, Inc. February 8, 2013 45

A crack inside ANSYS Mechanical is defined using a

• Pre-Meshed Crack object

• Crack object (Does not require geometry editing)

These objects can be inserted under the Fracture folder.

• Fracture Tool

o Stress Intensity Factors (K1, K2 and K3)

o J-Integral (JINT)

o Energy Release Rate (G1, G2 and G3) (VCCT)

Fracture Mechanics in ANSYS mechanical

© 2011 ANSYS, Inc. February 8, 2013 46

Pre-Meshed Crack object

A Pre-Meshed Crack is based on a previously-generated mesh (Mesh object in the tree outline). • 2D analyses(could still used for 3D) • Meshes imported with FE Modeler Crack front is defined using a nodal named selection. Coordinate system is required : • Y axis of the coordinate system defines the crack

surface normal • origin of the coordinate system represents the open

side of the crack

Solution Contours: Number of contours for which you want to compute the fracture result parameters.

© 2011 ANSYS, Inc. February 8, 2013 47

Crack Object Defines the semi-elliptical shape of the crack in three dimensional static analyses.

Base Mesh (Tetrahedrons quadratic elements) will then be modified to create a crack mesh.

The crack mesh generation is performed after the creation of the base mesh.

© 2011 ANSYS, Inc. February 8, 2013 48

…Crack Object

© 2011 ANSYS, Inc. February 8, 2013 49

•Coordinate system •Crack definition

•Major and Minor radius •Largest contour radius •Circumferential divisions •Mesh Contours •Crack front divisions •Solution contours

•Buffer Zone

…Crack object

© 2011 ANSYS, Inc. February 8, 2013 50

…Crack object

Buffer Zone: It represents the zone in which the elements will be removed to insert the crack template. The transition between the crack and the base mesh will be filled with elements.

Small “buffer zone”

Rapid mesh transition

Larger “buffer zone”

Smooth mesh transition

© 2011 ANSYS, Inc. February 8, 2013 51

Post Processing for crack

Path Plot and Graphic visualization of:

o Stress Intensity Factors (K1, K2 and K3) o J-Integral (JINT) o Energy Release Rate (G1, G2 and G3)

© 2011 ANSYS, Inc. February 8, 2013 52

Outline

o Modeling Composite

© 2011 ANSYS, Inc. February 8, 2013 53

Motivation

• Light weight • Corrosion resist • Fire protection • Impact resistant

Carbon fiber umbilical (Courtesy from Aker Solutions)

Phenolic gratings used on a tension leg platform (Courtesy from Strongwell)

Multi-Layer composite pipe (Courtesy from MLC Pipes)

Excellent attributes for composites:

High pressure gas accumulator bottles used in Riser tension assembly (Courtesy of Lincoln composites)

© 2011 ANSYS, Inc. February 8, 2013 54

Composite Integration in Workbench

Preprocessing Simulation Postprocessing

Evaluate the Composite Design • Evaluate Stresses in the Layer • Evaluate Failure Criteria

Build up the Composite Design • Define the Composite Layup • Define Fiber Directions

© 2011 ANSYS, Inc. February 8, 2013 55

Extrude Shell to Solid inside ACP

Launch Mechanical

•Material, Mesh from ACP is consumed in system B •Geometry is synthesized in B •Solid composite section from ACP is imported into B

In the case of thick composites, the layered shell theory can cause significant errors in the obtained results.

Extruded Solid in ACP

© 2011 ANSYS, Inc. February 8, 2013 56

Maintain topology between shell & solid model

Extruded Solid in Mechanical now has same face topology as shell

© 2011 ANSYS, Inc. February 8, 2013 57

Post Processing in ACP-Post

• ACP Post will share the material, geometry and Model with ACP Pre (System A) • Solution from Static Structural (System B) will be used by ACP Post for post processing

Post processing in ACP-Post

© 2011 ANSYS, Inc. February 8, 2013 58

Composite

Composite

Metal

Composite Non-composite inside Mechanical Assembling composite and other materials together in Mechanical

© 2011 ANSYS, Inc. February 8, 2013 59

Progressive Damage Analysis in Workbench

Result Support in Mechanical Material Support in Engineering Data

© 2011 ANSYS, Inc. February 8, 2013 60

Parameterization in ACP Pre-Post

Define parameters in ACP Pre and Post

© 2011 ANSYS, Inc. February 8, 2013 61

Outline

o Additional Control in Contact Modeling

© 2011 ANSYS, Inc. February 8, 2013 62

Motivation

Structural engineers need to go deeper and further in the detailed analysis of a given assembly, especially to better investigate connections

© 2011 ANSYS, Inc. February 8, 2013 63

“Trim Contact” automatically reduces the number of contact elements generated within each pair, thereby speeding up processer time.

“Program Controlled” will typically turn Trim Contact ON. However, no trimming is done for manually created pairs and/or when Large Deflection is turned ON.

WB Mechanical – Trim Contact

© 2011 ANSYS, Inc. February 8, 2013 64

…WB Mechanical – Penetration Tolerance

Penetration Tolerance: • Contact compatibility is satisfied in normal

direction if normal penetration is within allowable tolerance (TOLN )

• Can be defined as a Factor (of underlying element depth )or as a Value.

• Default =0.1*element depth (surf-surf) • Only exposed for Aug Lagrange

Formulation.

Fn

xp

npenetrationormalnormal xkFAugmented Lagrange:

© 2011 ANSYS, Inc. February 8, 2013 65

…WB Mechanical – Elastic Slip Tolerance

Ftangential

ELSI

ELSIkF tangentialtangential

Elastic Slip Tolerance: • Contact compatibility is satisfied in

tangential direction if Elastic Slip (ELSI) is within allowable tolerance (SLTO) • Defined as a Factor of average

underlying element length or as a value

• Applies to bonded, rough and frictional contact behaviors to enforce compatibility in tangential direction.

• Default SLTO=1% of average element length

© 2011 ANSYS, Inc. February 8, 2013 66

…WB Mechanical – Shell Thickness Effect

Shell Thickness Effect

• Applies only to surface bodies

Shell Thickness Effect = NO Shell Thickness Effect = Yes

© 2011 ANSYS, Inc. February 8, 2013 67

…WB Mechanical - Connection Matrix

New Connection Matrix in the Connections Worksheet

© 2011 ANSYS, Inc. February 8, 2013 68

Fluid Pressure Penetration Visualization

New contact result for fluid pressure penetration in the contact tool

Command snippets are still required to apply the loading to create this result

No Fluid zone

Fluid zone

© 2011 ANSYS, Inc. February 8, 2013 69

WB RBD – Joint Technology Enhancements Nonlinear Springs:

(RBD stands for Rigid Body Dynamics)

© 2011 ANSYS, Inc. February 8, 2013 70

Efficient modeling of small clearance between shaft and bearings are available for spherical, general and bushing joint and do not require full contact modeling

WB RBD – Joint Technology Enhancements

© 2011 ANSYS, Inc. February 8, 2013 71

Contact properties can be defined using User Programmable Functions or by tabular input for real constants (FKN, CNOF, TCC etc.) as function of time, temperature, pressure, gap/penetration etc.

MAPDL- Contact Tabular Real Constants

© 2011 ANSYS, Inc. February 8, 2013 72

• Example: User defined FKN vs Contact Pressure

… MAPDL- Contact Tabular Real Constants

*DIM,CNFKN_top,TABLE,6,1,1,PRESSURE CNFKN_top(1)=0.1,0.08,0.07,0.07,0.05,0.04 !FKN CNFKN_top(1,0)=1,5,10,50,100,150 !Contact Pressure RMODIF,cid,3,%CNFKN_top%

© 2011 ANSYS, Inc. February 8, 2013 73

MAPDL – User Routine Enhancements Enhancements to User-Defined Friction (USERFRIC):

• Supports structural and non-structural degrees of freedom

• The frictional stresses can be defined as a function of the slip increments and the contact/target temperatures.

New USERINTER Routine:

• Define complex interactions between contacting surfaces when the existing interaction models are not adequate. Supports Normal/Tangential Directions and Coupled multi-physics applications.

New USERCNPROP Routine :

• Program your own real constants for all current-technology contact elements (CONTA17x).

- Perform any nonlinear contact pressure/penetration operation instead of a linear operation (that is, use a varied contact stiffness instead of a constant stiffness). The contact stiffness can vary with pressure, penetration, temperature, and your own defined state variables.

© 2011 ANSYS, Inc. February 8, 2013 74

Outline

o Improved External Data Mapping

© 2011 ANSYS, Inc. February 8, 2013 75

Motivation

Exchange files are frequently used to transfer quantities from one simulation to another.

Efficient mapping of point cloud data is required to account for misalignment, non matching units or scaling issues.

© 2011 ANSYS, Inc. February 8, 2013 76

Users can import displacement and force data from external files into Mechanical using the External Data system

Support for Displacement and Force Transfer

Worksheet view to specify the components and scale results

Unavailable data can be ignored

© 2011 ANSYS, Inc. February 8, 2013 77

Users can import pressure data and apply in a downstream harmonic analysis

Ability to prescribe real and imaginary components (normal to and defined by components)

Import Pressure in a Harmonic Analysis

© 2011 ANSYS, Inc. February 8, 2013 78

Choice of coordinate system for vector data alignment

Default is “Source CS”

UY in Source CS

Any Local CS including Cylindrical

UY in Local Cylindrical CS

© 2011 ANSYS, Inc. February 8, 2013 79

2D-3D Mapping • Imported vector data are transformed for 2D-3D mapping.

• Unavailable data can be ignored in the definition.

Z Component ignored because imported data only available for X and Y components

© 2011 ANSYS, Inc. February 8, 2013 80

Further options to validate the mapping

Additional visual tools have been implemented to control how well the data has been mapped onto the target structure.

Isolines comparizon

Components display

Contours display

© 2011 ANSYS, Inc. February 8, 2013 81

Improved mapping performance

Kriging weighting calculations now use multiple cores to faster results. For larger meshes utilizing 8 cores, a 6 to 7 time’s speedup can be obtained.

# Source nodes # Target Nodes R14 (1 Core) R14.5 (8 cores) Speedup

207602 79152 41.46 8.27 5.01

387419 437816 366.94 48.18 7.62

2570789 526082 3510.19 462.53 7.59

240178 700911 336.25 50.11 6.71

378779 4457975 3289.4 424.96 7.74

© 2011 ANSYS, Inc. February 8, 2013 82

• The only file with the coordinate information is regarded as the master file. • The rest of them will have data to be mapped against the Node IDs. • This removes the requirement of storing coordinates in each file.

Mapping multiple data sets: Concept of a master file

© 2011 ANSYS, Inc. February 8, 2013 83

Outline

o Miscellaneous features in ANSYS Mechanical

© 2011 ANSYS, Inc. February 8, 2013 84

The user can now use the Named Selection worksheet to find nodes and geometry items using a “Distance” criteria. The distance is calculated from the origin of the selected coordinate system.

Note: The distance when searching for bodies, faces, edges, and vertices is calculated from the centroid of that item.

Distance Criteria For Named Selections

© 2011 ANSYS, Inc. February 8, 2013 85

Hide/Show options on Named Selections

© 2011 ANSYS, Inc. February 8, 2013 86

2D Pretension Bolt Loads Pretension Bolt loads can now be applied in 2D Analysis

© 2011 ANSYS, Inc. February 8, 2013 87

Surface Force(Moment) Reaction Probe

© 2011 ANSYS, Inc. February 8, 2013 88

Display annotation improvement Maximum annotation=50, with different colors

© 2011 ANSYS, Inc. February 8, 2013 89

From Left to right: 1. Create a View 2. Delete View(s) 3. Apply a View 4. Rename 5. Replace saved view based on current

graphics 6. Import Views from a formatted xml in disk 7. Export Views as a formatted xml to disk

Added Manage View

© 2011 ANSYS, Inc. February 8, 2013 90

Show Errors Button New Button in the toolbar to populate message window with explanation for under-defined objects. In prior release, a solution had to be invoked to populate the window.

© 2011 ANSYS, Inc. February 8, 2013 91

Check ANSYS 14.5 release notes

© 2011 ANSYS, Inc. February 8, 2013 92

Thank you!

top related