DIANA 10 Finite Element Analysis. Civil Engineering Geotechnical Engineering Petroleum Engineering

Size: px

Start display at page:

Download "DIANA 10 Finite Element Analysis. Civil Engineering Geotechnical Engineering Petroleum Engineering"

Kenneth Wright
5 years ago
Views:

1 DIANA 10 Finite Element Analysis Civil Engineering Geotechnical Engineering Petroleum Engineering

History of TNO DIANA 1972 TNO Department of Computational Mechanics 1975 Start of DIANA (Displacement ANAlyzer) as a service 1980 DIANA commercially released in the Netherlands 1984 Establishment of

2 History of TNO DIANA 1972 TNO Department of Computational Mechanics 1975 Start of DIANA (Displacement ANAlyzer) as a service 1980 DIANA commercially released in the Netherlands 1984 Establishment of DIANA User Association 1990 Partnership set up with JIP Techno Science Corporation in Japan 1994 First International DIANA Conference 1998 TNO acquires FemSys Ltd. Its core product FEMGV becomes the new pre/post processor for DIANA 1999 Geomec development project with Shell International Exploration & Production 2000 Well-Life development project with Halliburton 2002 Fempal development project with EU store-rack manufacturing companies 2003 Establishment of TNO DIANA BV, subsidiary company of TNO 2004 Opening of US office TNO DIANA North America inc Release of DIANA Cooperation started with Midas IT 2006 Opening of UK office TNO DIANA UK Ltd 2006 Release of DIANA 9.2 with new pre-post-processor FX+ for DIANA 2007 Development of Hybrid Frequency-Time Domain Analysis module 2008 Release of DIANA Joint development with Midas IT of midasgts kernel with DIANA solver inside 2009 Release of DIANA Start of development of new state-of-the-art graphical environment for DIANA 2012 Recognition of DIANA as the optimized solution for hydropower applications worldwide 2014 Release of DIANA 9.5 including first release of new MeshEdit graphical user interface 2014 Establishment of DIANA Engineering office Arnhem 2015 Release of DIANA Partnership set up with NSG (Brazil) and DUNPU (China) 2016 Landmark Release of DIANA 10 with integrated and state-of-the-art pre-/post-processor Now Advancing in new numerical analysis techniques and providing cutting-edge technologies for engineering applications

3 idiana Integrated pre/post processor Working Window Works-Tree with full access to model s components Model Navigator Command line Interactive dialogue boxes

4 FX4D Task oriented pre/post processor Main Menu Tabbed Toolbar Works Tree Property Window Table Window Output Window

Performing Analysis Output to FX4D DIANA FILOS File

5 Work Flow Geometry Modeling Mesh Generation Load & B.C. s Post-processing FX4D Post-Neutral File Performing Analysis Output to FX4D DIANA FILOS File Pre-Neutral File MeshEditor Check & Final Editing Analysis Options

FX4D Advantages Graphical tools Intuitive

Advanced selection methods CAD-like modeling

& repairing tools Meshing tools Various easy &

elements Embedded reinforcements generation

boundary conditions Full support of load and

geometry entities Function based definition

6 FX4D Advantages Graphical tools Intuitive graphical user interface Various display options Advanced selection methods CAD-like modeling functions Efficient and robust geometry checking & repairing tools Meshing tools Various easy & strong meshing algorithms for 1st- and 2nd-order elements Embedded reinforcements generation Automatic DIANA element type selection Loads and boundary conditions Full support of load and boundary conditions Both applicable on mesh and geometry entities Function based definition MS-Excel compatible tables Post-processing Complete solution to result interpretation

7 DIANA Interactive Environment Toolbar Model window with full access to model s components Working Window Analysis Window Output Window Properties window Python Window Message Window

8 DIANA Interactive Environment Workflow(s): Geometry definition Assign properties Generate mesh Setup analysis Change in physical properties Run analysis Check results

DIANA Interactive Environment Graphical tools Intuitive graphical user

modeling functions Automated and robust geometry checking & repairing tools

Hybrid mesher) Embedded reinforcements generation Automatic DIANA element type

conditions Applicable on geometry entities with some flexibilities at mesh

9 DIANA Interactive Environment Graphical tools Intuitive graphical user interface Various display options Advanced selection methods Parasolid modeling functions Automated and robust geometry checking & repairing tools Meshing tools Powerful and automated mesh engines 2D / 3D (incl. Hybrid mesher) Embedded reinforcements generation Automatic DIANA element type selection Loads and boundary conditions Full support of load and boundary conditions Applicable on geometry entities with some flexibilities at mesh level Function based definition MS-Excel compatible tables Post-processing Complete solution to result interpretation

DIANA Interactive Environment Material models classifications Separate dialog boxes with various

functions and variable dependencies Load combination Generates design load combinations Generates

space) functions Properties Window Full control on parameterizing the selected operations

10 DIANA Interactive Environment Material models classifications Separate dialog boxes with various input fields Unique aspects to be activated in the material models Interactive tables for functions and variable dependencies Load combination Generates design load combinations Generates distinctive load combinations for analyses Functions Generates various (Spatial/non spatial and space) functions Properties Window Full control on parameterizing the selected operations Analysis Window Intuitive environment to setup any types of analysis including interactive phased analysis Many more

11 Background and development direction idiana Scripting Integrated environment FXD Geometry tools Mesh engines Application oriented International design codes Dedicated wizards User-friendliness 11

12 After DIANA 10.0 Continuous extension of pre-post functions Import/export of CAD/BIM formats Graphs for results Report generation Application functions Solver Performance, faster solving larger models State of art non-linear models Design checks

simulations for specific engineering applications Can be used by

13 After DIANA 10.0 Fit for purpose applications: State of the art nonlinear numerical simulations for specific engineering applications Can be used by application engineers, who are not analysis specialists Do analysis in short time (< 1 hour) Engineering input and engineering output

14 Modeling concepts in DIANA Interactive Environment

15 Geometry modelling Geometrical modelling functions: Import STEP/IGES file formats Primitives, basic shapes Move, Rotate, Scale Boolean operations Selection of points/lines/faces/bodies and operations Automatic clash detection Convert bodies and sheets Imprint points/lines on faces Rotate and move faces of a shape Extrude shapes Extract a sub-shape Align faces Many more 15

16 Modeling primitives 3D Import STEP/IGES file formats Import of CAD models Geometrical modelling based on Parasolid 16

17 Modeling primitives 3D solids Block Cylinder Cone Prism Torus Sphere 17

18 Modeling primitives 2D primitives 1D Polygonal sheet Circle sheet Line Circle Polyline Bezier curve Point 18

19 Modeling operations Boolean operations Unite Subtract Intersect Transformations Move Scale Rotate 19

20 Selection - Geometry parts Graphical selection of points, lines, surfaces and bodies 20

21 Modeling operations Automatic clash detection Edges of adjacent bodies and lines are automatically imprinted and considered in meshing 21

22 Modeling operations Sew two or more sheets Convert bodies and sheets 22

23 Modeling operations Imprint points and lines on surface 23

24 Modeling operations Rotate and move faces of a shape 24

25 Modeling operations Extract a sub-shape Extrude shapes Lines Faces 25

26 Modeling operations Align faces Coincident Co-planar Parallel Co-centric 26

Modeling operations - Spatial functions Spatial functions to geometry parameters Shell model with spatial functions to thickness and eccentricity of elements. 'FUNCTI' NAME "Girders" X -9.00000E+00-7.

00000E+00 4.00000E+00 3.50000E+00 3.25000E+00 3.10000E+00 3.00000E+00 3.00000E+00 3.10000E+00 3.25000E+00 3.50000E+00 4.00000E+00 5.00000E+00 5.00000E+00 5.00000E+00 / NAME "Flange 1" Y 2.25000E+00 3.75000E+00 / FACTOR 2.

00000E+00 0.00000E+00 2 NAME "Girders" THICK 2.00000E-01 FUNCTI "Girders" THICK ECCENZ -1.00000E-01 FUNCTI "Girders" ECCENZ XAXIS 1.00000E+00 0.00000E+00 0.00000E+00 3 NAME "Flange 1" THICK 2.

27 Modeling operations - Spatial functions Spatial functions to geometry parameters Shell model with spatial functions to thickness and eccentricity of elements. 'FUNCTI' NAME "Girders" X E E E E E E E E E E E E E E E+00 / FACTOR E E E E E E E E E E E E E E E+00 / NAME "Flange 1" Y E E+00 / FACTOR E E+00 / NAME "Flange 2" Y E E+00 / FACTOR E E+00 / 'GEOMET' 1 NAME "Central slab" THICK E-01 ECCENZ E-01 XAXIS E E E+00 2 NAME "Girders" THICK E-01 FUNCTI "Girders" THICK ECCENZ E-01 FUNCTI "Girders" ECCENZ XAXIS E E E+00 3 NAME "Flange 1" THICK E-01 FUNCTI "Flange 1" THICK ECCENZ E-01 FUNCTI "Flange 1" ECCENZ XAXIS E E E+00 4 NAME "Flange 2" THICK E-01 FUNCTI "Flange 2" THICK ECCENZ E-01 FUNCTI "Flange 2" ECCENZ XAXIS E E E+00

28 Model setup Property assignment Property assignment: Element-classes (Material, Geometry, Data) Reinforcements (Material, Geometry, Data) Interfaces (Material, Geometry, Data) Boundaries for heat-flow, groundwater analysis (Material) Supports (structural analysis) Loads acting on parts Loads acting on model Initial fields acting on parts 28

29 Model setup Property assignment Element classes 29

30 Model setup Property assignment Material classes 30

31 Model setup Property assignment Geometry classes 31

32 Model setup Property assignment Reinforcement classes Bodies Faces Faces Lines Lines Points 32

Model setup Property assignment Interface classes Defined on sub-shapes Faces (2D) Edges (1D) Connected (Interface element between continuum elements)

33 Model setup Property assignment Interface classes Defined on sub-shapes Faces (2D) Edges (1D) Connected (Interface element between continuum elements) Boundary (Interface element between continuum element and supported edge/line) Not connected (Deactivate clash detection continuum elements not connected) 33

External temperature Radiative temperature Thermal flux

34 Model setup Boundary assignment Boundaries for thermal and groundwater flow analysis Heat Flow Prescribed temperature External temperature Radiative temperature Thermal flux Groundwater Flow Prescribed head External head Groundwater flux 34

35 Model setup Boundary assignment Supports for structural analysis On points, lines, faces and bodies In global and local coordinate-system Defined for all nodes related to selected part 35

36 Model setup Load assignment Loads acting on parts Type of loads Points Lines Faces Bodies Force X X X Moment X X X Prescribed Deformation Prescribed Acceleration Distributed Force X Distributed Moment X X Hydrostatic Pressure X X Pore Pressure X X Reinforcement Pre-stress X Reinforcement Post-tensioning X X X 36

37 Model setup Load assignment Loads acting on parts 37

38 Model setup Load assignment Loads acting on parts 38

39 Model setup Load assignment Loads acting on model Dead weight Equivalent acceleration Centrifugal load Base acceleration

40 Model setup Load assignment Initial fields Displacements Rotations Velocities Temperatures Heads 40

41 Meshing options Generate the mesh Seedings per body/line (divisions/element-size) In clashes the smallest element-size is applied Preview of seeding Linear or quadratic elements Automatic quad/hexa mesher or triangle/tetrahedron mesher Surface meshes are defined as structured meshes when possible Body meshes are not structured Reinforcements are embedded in mesh Interfaces are embedded in mesh Boundary elements are embedded in mesh 41

42 Meshing options Generate the mesh 42

43 Analysis setup Setup analysis Define analysis sequence Define analysis details via dialogues or property box Options for loading or saving.dcf files Run analysis Progress logging in GUI Option for running analysis in batch mode Error message refer graphically to model Logging selected analysis results in GUI Automatic import of results at the end of analysis Option for loading results manually 43

Post-processing features Checking results: Deformed models Contour-plots Iso-surfaces Clipping planes Diagrams Vectors Stress and strain tensor rosettes Cracks

44 Post-processing features Checking results: Deformed models Contour-plots Iso-surfaces Clipping planes Diagrams Vectors Stress and strain tensor rosettes Cracks Tables, export to Excel Result view setting management (inc. saving the post-processing settings) Optional output in FX4D, idiana, original tabular format (.txt) 44

45 Post-processing features Deformed shape 45

46 Post-processing features Contour plot settings (iso-surfaces)

47 Post-processing features Contour plot settings (Clipping planes)

48 Post-processing features Contour plot settings (Clipping planes)

49 Post-processing features Contour plot settings (Clipping planes)

50 Post-processing features Contour plot settings (Clipping planes)

51 Post-processing features Contour plot settings

52 Post-processing features Reinforcements contours in grids

53 Post-processing features Vector plot settings

54 Post-processing features Tensor plot settings

55 Post-processing features Tensor plot settings (clipping plane)

56 Post-processing features Stress/strain normal to cutting plane

57 Post-processing features Crack plot settings

58 Post-processing features Crack plot settings (clipping plane)

59 Post-processing features Crack plot settings (clipping plane)

60 Post-processing features Tables Copy and paste to MS-Excel

61 Post-processing features Diagram plot settings

62 Post-processing features Layers in shell and beam elements

$"HEX_QUAD" ) n = 4 m = 6 icnt = 0 for i in range( 0, n ): for j in range( 0, m ): icnt = icnt+1 name = "Pile " numb = '{0:6d}'.$

63 Post-processing features Python scripting newproject( "Untitled", 100 ) setmodelanalysisaspects( [ "STRUCT" ] ) setmodeldimension( "3D" ) setdefaultmeshorder( "LINEAR" ) setdefaultmeshertype( "HEX_QUAD" ) n = 4 m = 6 icnt = 0 for i in range( 0, n ): for j in range( 0, m ): icnt = icnt+1 name = "Pile " numb = '{0:6d}'.format(icnt) pilename = name + numb createblock( pilename, [ 2*i, 2*j, 1 ], [ 1, 1, 5 ] ) createblock( "Plate", [ 0, 0, 0 ], [ 2*(n-1)+1, 2*(m-1)+1, 1 ] ) > newproject( "Untitled", 100 ) > setmodelanalysisaspects( [ "STRUCT" ] ) > setmodeldimension( "3D" ) > setdefaultmeshorder( "LINEAR" ) > setdefaultmeshertype( "HEX_QUAD" ) > createblock( "Pile 1", [ 0, 0, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 2", [ 0, 2, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 3", [ 0, 4, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 4", [ 0, 6, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 5", [ 0, 8, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 6", [ 0, 10, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 7", [ 2, 0, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 8", [ 2, 2, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 9", [ 2, 4, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 10", [ 2, 6, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 11", [ 2, 8, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 12", [ 2, 10, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 13", [ 4, 0, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 14", [ 4, 2, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 15", [ 4, 4, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 16", [ 4, 6, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 17", [ 4, 8, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 18", [ 4, 10, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 19", [ 6, 0, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 20", [ 6, 2, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 21", [ 6, 4, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 22", [ 6, 6, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 23", [ 6, 8, 1 ], [ 1, 1, 5 ] ) > createblock( "Pile 24", [ 6, 10, 1 ], [ 1, 1, 5 ] ) > createblock( "Plate", [ 0, 0, 0 ], [ 7, 11, 1 ] )

64 More new features in DIANA

65 New in DIANA 10 Element Library Element-library: 2-node class III (Timoshenko) beams Layered material input for flat shell elements Averaged normals in plane interface elements Matrix reinforcements in 3-dimensional solids and axi-symmetric solids

66 New in DIANA 10 Element Library Pyramid elements Dominant hexahedron mesher (Five-node isoparametric solid pyramid element and a Thirteen-node isoparametric solid pyramid element)

67 New in DIANA 10 Element Library Matrix reinforcements in 3-dimensional and axi-symmetric solids

New in DIANA 10 Material Library Material-library: Modified Maekawa model replaced by Maekawa-Fukuura model (non-orthogonal cracks) More options in Total Strain

68 New in DIANA 10 Material Library Material-library: Modified Maekawa model replaced by Maekawa-Fukuura model (non-orthogonal cracks) More options in Total Strain Crack model Shear-stiffness curves for simple soil models JCSS Probabilistic Model Code for concrete Predefined materials with reference to international design codes.

$cracking, dominant crack Elasto-plastic fracture model$

69 New in DIANA 10 Material Library Modified Maekawa model replaced by Maekawa-Fukuura model Non-orthogonal cracking, dominant crack Elasto-plastic fracture model (compressive zone) Non-orthogonal total-strain based crack model

70 New in DIANA 10 Material Library Shear-stiffness curves for simple soil models Hardin-Drnevich Ramberg-Osgood Diagram ( ) Diagram G/G0( )

71 New in DIANA 10 Material Library JCSS Probabilistic model code: Concrete Properties Reference property is basic compressive strength f c0. Other parameters are related to f c0 : In-situ compressive strength f c = α f 0.96 c0 [ MPa ] Tensile strength f t = 0.3 f 2/3 c [ MPa ] Modulus of Elasticity E c = 10.5 f 1/3 c ( 1/(1+ẞϕ)) [ GPa ] Ultimate compression strain ε u = f 1/6 c (1+ẞϕ) [ - ] Average and standard deviation values provided for concrete grades. Random field generated for f c0 Log-normal distributions Covariance length = 5 m, threshold-value = 0.5.

72 New in DIANA 10 Material Library JCSS Probabilistic model code: Concrete Properties Example of crack pattern in concrete floor on elastic bedding with probabilistic concrete properties defined according JCSS model code.

distributed properties (random fields) - Functions are defined in global coordinate system, independent from FE

73 New in DIANA 10 Spatial functions Spatial functions to material parameters - For geotechnical applications with parameters being dependent on depth - For plate structures with varying thickness and eccentricities - Random distributed properties (random fields) - Functions are defined in global coordinate system, independent from FE mesh - Properties at internally interpolated at integration point level - Functions can be defined in 1, 2 or 3 global directions

New in DIANA 10 Analysis Procedures Analysis procedures:

(spreading the results) Output of model parameters Input of

Eigen-modes Output reinforcement moment/forces for composed

74 New in DIANA 10 Analysis Procedures Analysis procedures: Parallel processing of element loops Spatial functions (spreading the results) Output of model parameters Input of frequency loads in terms of periods Extension of logging Eigen-modes Output reinforcement moment/forces for composed surface Eurocode 8 EN Dutch NEN NPR 1998 (Dynamic Response Spectrum Input)

DIANA 10 allows to make use of Parallel processing of Eigen Solver DIANA 10 allows to make use of parallel processing for the following element loops: Setting up linear elastic stiffness matrix

75 New in DIANA 10 Analysis Procedures Parallel processing of element loops In DIANA 9.6 only the linear solvers (PARDIS, ITERAT) divides the analysis job in parts that can be run at different processors at the some time. DIANA 10 allows to make use of Parallel processing of Eigen Solver DIANA 10 allows to make use of parallel processing for the following element loops: Setting up linear elastic stiffness matrix Calculating nonlinear stresses and internal forces Setting up tangent stiffness matrix Parallel processing will be available in standard version Number of threads can/should be managed by user For single and multiple processing a new procedure was introduced to calculate the integration point contribution which gives significant speed up at non-virtual machines.

76 New in DIANA 10 Analysis Procedures Parallel processing of element loops Examples of performance gains in setting up linear stiffness matrices

77 New in DIANA 10 Analysis Procedures Hidden feature. Random fields Following field generation procedures are considered: Covariance method Fast Fourier transformation method Local average subdivision *FILOS INITIA *INPUT *RANFLD BEGIN FIELD FLDTYP=COVARI BEGIN COVARI NX=40 NY=20 END COVARI COVTYP=EXPONE CORLEN=10. STDDEV= 20. MATIDX 1 MATPAR YOUNG END FIELD *LINSTA BEGIN OUTPUT NDIANA DISPLA STRESS INTPNT PARAME YOUNG INTPNT END OUTPUT *END

78 New in DIANA 10 New Modules Worthwhile Special projects Affordable Reasonable More detailed analysis Daily Engineering Projects

79 TNO DIANA BV Delftechpark 19a 2628 XJ Delft The Netherlands T +31 (0) F +31 (0) TNO DIANA BV (Engineering Division) Engineering (TD/E) Vlamoven TN Arnhem The Netherlands T +31 (0) TNO DIANA UK Ltd Ground Floor Building 1000 Lakeside North Harbour Western Road Portsmouth PO6 3EZ United Kingdom T +44 (0) F +44 (0) TNO DIANA North America Inc N.Laurel Park Drive, Suite 205 MI Livonia T +1 (0) F +1 (0)

DIANA. Finite Element Analysis. Civil Engineering Geotechnical Engineering Petroleum Engineering

DIANA. Finite Element Analysis. Civil Engineering Geotechnical Engineering Petroleum Engineering DIANA Finite Element Analysis Civil Engineering Geotechnical Engineering Petroleum Engineering NOW Advancing in new numerical analysis techniques Developing state-of-the-art solution for engineering applications