PAM-CRASH/PAM-SAFE V2006 Examples

Size: px

Start display at page:

Download "PAM-CRASH/PAM-SAFE V2006 Examples"

Loren Griffith
5 years ago
Views:

1 PAM-CRASH/PAM-SAFE V2006 Examples Author(s): Support Team 1

2 Multi Model Coupling 2

3 Multi-Scale features Multi-Scale Modelling Main Assembly File Module Definitions (limited to 2 Modules) Inter-Model Interface through CNTAC Inter-Model Connection through MGRID Running MMC job with pamworld 3

4 Commented Example Truck Frontal Impact vs. ODB barrier MODULE 100 MODULE 200 4

5 Commented Example Module Input Cards In-line, INCLU or both 5

6 Commented Example Inter-Model Interface through CNTAC MOD, General Entity Selection, END_MOD SLAVE SELECTION MASTER SELECTION 6

7 Commented Example Run Coupling Job Command line: -pamarg N=n1+n2 7

8 Commented Example Run Coupling Job Cluster file supported: -cf cfile - Virtual Machine of N processors - Same rule : N=n1+n2 Resulting Files : *.DSY, *.THP, *.out, *.msg - FILE Names Ignored - replaced by MODULE Names 8

9 Material type 304 Tied interface 9

10 Material type 304 for TIED interface Main Features Material type dedicated to model Glue Orthotropic properties for tension and compression defined by curves (Young modulus) Orthotropic properties for shearing direction (directions are uncoupled) Strain limit for all directions for element elimination Damage function for all directions 10

11 Material 304 definition hcont Master segment Slave node Element normal Ideal position Actual position dε n dx dx n t = ; dγ t = ; dγ u hcont hcont = dx h cont u Displacement decomposition hcont - is the user imposed tied thickness, which permits to correct modelisation error between the master and slave part (used by solver for glue response calculation). 11

12 Material 304 definition Yield curve for tension σ Yield curve for compression σ ε Positive strain +ε ε Negative strain Positive strain +ε Compression unloading Tensile loading E tension E compression Compression loading Tensile unloading 12

13 Max. Strain limit in tension MATER / Definition of strain limit for element elimination (Material 304 definition) Imposed displacement Example-1 Imposed Displacement NAME TIED mater-304 $ $ SDMP1 hcont D $ EC_NN ICUC_NN 0.6CURVE 30 $ ET_NN ICUT_NN 0.6CURVE 30 $ G_NT ICUT_NT 5 $ G_NU ICUT_NU 5 $ LKT_N INT_N INT_NT INT_NU $ LKT_N INT_N INT_NT INT_NU $ LKT_N INT_N INT_NT INT_NU $ LKT_N INT_N INT_NT INT_NU $ CELC_NN ELIC_NN CELT_NN ELIT_NN 4.5 CELT_NT ELIT_NT CELT_NU ELIT_NU Strain limit for element elimination Result Strain limit Yielding curve in tension 13

14 Damage function for tension MATER / NAME TIED mater-304 $ $ SDMP1 hcont D $ EC_NN ICUC_NN 0.6CURVE 30 $ ET_NN ICUT_NN 0.6CURVE 30 $ G_NT ICUT_NT 5 $ G_NU ICUT_NU 5 $ LKT_N INT_N INT_NT INT_NU Imposed displacement LOOKU / 1 NAME Damage table $ ARGUMENT END FUNCTION END Definition of damage function (Material 304 definition) Example-2 $ LKT_N INT_N INT_NT INT_NU 1 1 $ LKT_N INT_N INT_NT INT_NU $ LKT_N INT_N INT_NT INT_NU $ Imposed Displacement Complete Damage and element elimination follow 14

15 Material Type 42 Eroding Contact 15

16 Impact of a Crash box with Barrier ODB Material type 42 and eroding contact Material type 41 replaced by material type 42 with equivalent properties Application of eroding contact 16

17 Material type 42 and eroding contact Contact treatment improvement in case of solid element elimination Useful in case of crash barriers Takes into account internal faces of solid element in the initial list of contact segment CPU increase of 15% 17

18 MPC Plink Material 224 6DOF Penalty Spring-beam element 18

19 Main Features: Correspond to Penalty spring-beam element which connects two nodes at distance 0 Penalty Stiffness Calculation and Damping Material 224 displacement/rotation and forces/moments are internally computed Element can be used to automatic definition of the mechanical properties in MPC PLINK interfaces (in user friendly way) equivalent mass is automatically computed at the initialization step XMASS and INERTIA are ignored 2 different ways of constrained DOFs are available: only translational degrees of freedom are constrained both rotational and translational degrees of freedom are constrained 2 types of rupture model are implemented User defined Type definition 5 T N 1 S N 2 M R 19

20 Material 224 definition Penalty scale factor for translation Penalty scale factor for rotation Stiffness proportional damping ratio (<1) MATER / e NAME mat-spring224 $ $ SLFACMT SLFACMR SDMP1 XMASS INERTIA I3DOF IDRUPT $ $ PART / 5 PLINK 5 NAME MPC-Plink 224 Flag for releasing rotational degrees of freedom Rupture model identifier END_PART $ RUPMO / 10 5 NAME Rupture for mater 224 $ $ Torsion Force Bending Force Normal Force Shear Force 20

21 Rupture model (type 5): ( ) a ( ) a ( ) a ( ) a NORMALFORC E + SHEARFORCE + TORSIONFOR CE + BENDINGFOR CE 1. 0 AFAILN 1 AFAILS 2 AFAILT 3 Material 224 AFAILB 4 Shear force rupture limit Normal force rupture limit Torsion force rupture limit 21

22 Monitoring for RUPMO 22

23 Monitoring for RUPMO Activation of monitoring option in case of a «T» section welded with PLinks with rupture model and loaded in the normal direction No monitoring (default) Monitoring activated No split visualized after rupture/contour display 23

24 Monitoring for RUPMO Activation of monitoring option in case of a welded boxbeam crushed by a rigid body 24

25 Monitoring for RUPMO Activation of monitoring option in case of a welded boxbeam crushed by a rigid body No monitoring (default) Monitoring activated 25

26 Monitoring for RUPMO Activation of monitoring option in case of a welded boxbeam crushed by a rigid body No monitoring Monitoring activated: No split visualized Contour displayed 26

27 Slipring friction - locking 27

28 Slip-Ring Friction Reminder Friction force β Pulley Rim Formulation: F1 = F2 exp ( µβ ) F1 F2 µ = FRICPR, pulley rim friction β = contact angle of the belt on the slip-ring F1>F2 28

29 Slip-Ring Locking Example α Bar β = PI-α Bar

30 Slip-Ring Friction Example (2/2)

31 Slip-Ring Locking Example (1/2) 0: Locking Enabled (Dflt.) 1: Locking Disabled (old behavior) Reminder : L < L 0 x FRACT 31

32 Slip-Ring Locking Example (2/2) 32

33 Nodal Face Damping 33

34 Nodal Face Damping Example (1/2) Perturbation Drag force coeff. For various bodies Face Damping Air Density (MM,KG,MS) 34

35 Nodal Face Damping Example (2/2) Perturbation Without face damping With face damping 35

36 3D Boundary Condition Reaction Forces 36

37 Reaction Forces in the Nodes with 3DBC 3D Boundary Condition Prescribed Displacement Rigid body COG Reaction forces can by easily evaluated for all Nodes with applied 3D boundary condition (DIS3D, VEL3D, ACC3D, RAN3D, RAV3D and RAC3D) as /SECFO -> SUPPORT Comparison between Reaction Force in loaded end and Contact Force (Rigid-plate & Beam) Group of Nodes with 3D displacement 37

38 Stiffness scale test rigid ball impacts a thick plate Stiffness scale factor contour plot OCTRL / GLBTHP DFLT DMSC SHLPLOT STSC SOLPLOT STSC END_OCTRL TCTRL / INITIAL 0. NODAL YES DYNA_MASS_SCALE STIFFNESS_SCALE END_TCTRL Stiffness scale factor contour Time step affected via Stiffness and Dynamic Mass Scaling Pamview shows stifness scale factor for each element (not affected elements has value 1). In time 0 s File: ballimpact2_05.pc In time 1 s Applicable for Shell and Solid materials: 102, 103, 105, 106, 107, 108, 109 nad 1, 2, 7, 16 38

39 3-nodes TRSFM Transformation is defined by 3 sets of 3 nodes 39

40 TRSFM test simple impact test 3 nodes TRSFM TRSFM / NAME Impactor_Transformation PART 4000 END NPOS MOVE MOVE ROTA END Task: It is necessary to transform impactor to the right position and correct its size (via rotation, translation and scaling). Initial position source nodes N1 =4024 N2 =4023 N3 =4081 destination nodes N1'=5024 N2'=5023 N3'=5081 After Initialization phase impactor occurs in right position and size N3 N1 N2 Rigid body plate compress cube element (loaded via CONLO) File: TRSFM_NPOS_MOVE_ROTA.pc Impactor transformation = Rotation and scaling 40

41 MADYMO Coupling 41

42 MADYMO Coupling MADYMO Coupling Coupling with MADYMO 6.3 version PAMCRASH and MADYMO may have different UNIT system, automatic unit conversion done Supported platforms : Linux ia32 (GlibC 2.3) Linux ia64 (GlibC 2.3) HP-UX ia64 42

43 MADYMO Coupling Example AMS Frontal Impact UNIT MM KG MS KELVIN Powered by VCP 2.5 MADYMO 50th percentile Hybrid III UNIT S.I. 43

44 Material 22 and 5 - DDM Dilatation Damage Measure 44

45 DDM card : SECFO / 2VOLFRAC NAME ddm 1 DDM PART 1 END DDM Dilatation Damage Measure Example-1 DDM test on material 22 Volumetric tension test VEL3D DDM evaluated on the complete PART. DDM Critical value Vcrit Result of the test - element stress pressure DDM value is 1 When element pressure falls below Vcrit the DDM records current volume fraction form the initialy specified set DDM value (DDM for our case takes the value 0 or 1 because only one element is considered) Element stress Files: DDM_mat22.pc 45

46 DDM Dilatation Damage Measure Example-2 DDM test on material 5 DDM card : SECFO / 2VOLFRAC NAME ddm 1 DDM -0.5 PART 1 END DDM evaluated on the complete PART. Compression tension test of 12 elements model fixed VEL3D Result of the test - element stress pressure The DDM records the current volume fraction that falls bellow the Vcrit value of the pressure DDM value 0.25 DDM value in case 1 = 4/12 = DDM value in case 2 = 3/12 = 0.25 Where: Volume of one element = 1 Volume of complete part = 12 DDM Critical value Vcrit Element stress Files: DDM_mat5_12elem.pc Case 1 : Pressure of 4 elements falls below Vcrit Pressure of 3 elements Case 2 : falls below Vcrit 46

Presentation of PAM-CRASH v2004. Part 1: Solver News

Presentation of PAM-CRASH v2004 Part 1: Solver News. 1 Overview New Options Elements Materials Others Quality Numerical precision and robustness 2 CRASH/SAFE 2G Evolution V2002: Basic reengineering Common