Comparing default shell properties with recommended properties through simulation in HyperView and plotting in HyperGraph

OBJECTIVE: To compare the results with base simulation and improved shell elements properties by performing energy, mass error checks, plotting graph in HyperGraph and animating in HyperView.

SIMULATION 1: With default Shell Properties.

Imported Element_Formulation-Shell-3_assignment_0000.rad file. The property was not changed.

Changed Runtime to 55milli seconds and animation frequency to 2.2 to get 25 animation files.

SIMULATION 2: With recommended shell properties, same Runtime and animation frequency as above.

After assigning recommended properties:

ENERGY ENERGY AND MASS  ERRORS:

Energy error must be negative and decreasing. The error must be less than 15% at the end of the simulation (5% for a good model). If error increases a little and then decreases it may be acceptable.If the error increases, this means that a problem has occurred. Then the error may grow to 99.9% with a computation failure message.

So ideally the error should be atmost -15% and better to have -5%. The closer to 0% error the better. For QEPH shell formulation or fully integrated elements, the error can be slightly positive since there is no hourglass energy and computation is much more accurate. An error of 1-2% is acceptable. Otherwise, any positive error means which means energy has been created and it is a model issue.

Opened Element_Formulation-Shell-3_assignment_0001.out file to check the energy error and mass errors.

At the end of simulation1 error is shown below:

ENERGY ERROR AND MASS ERROR: Simulation 1

-10.3% AND MASS ERROR 0.1659E-03

Since energy error is negative it indicates energy has been dissipated from the system. 

Also since energy error is less than -15%, so the model is ACCEPTABLE for further calculations.

ENERGY ERROR AND MASS ERROR: Simulation 2

-0.1% and 0.1659E-03

ELEMENT QEPH(4node) Ishell 24 has inbuilt Hour Glass stabilization due to which there is a reduce in energy errors in the second simulation.

PLOTTING OF VARIOUS GRAPH:

SIMULATION 1 AND 2: HOURGLASS

From the above plot we can see HG in simulation 2 is constant and zero. This is mainly due to the selection of QEPH element.

HG energy is controlled by two methods.

HG PERTURBATION METHOD: For element suspected to HG modes, internal HG forces are applied to resist the HG mode deformation.

Stiffness Method: Here Extra stiffness is applied on the material itself even before HG error can appear so there is HG energy generation in this method. QEPH elements use this method.

SIMULATION 2 AND 1:INTERNAL ENERGY

            SIMULATION 2                                              SIMULATION 1

Internal energy is the sum of kinetic and potential energy. In simulation 1 change in internal energy is almost linear and in 2 we see a curve. At the beginning KE was high and IE was zero due to zero deformation and at the end, KE goes to zero and IE becomes high and constant but still less than KE mainly due to energy error. we can see zero KE only if we increase run time more than 55ms.

If we increase runtime to 85ms, KE becomes zero and IE become maximum and remain constant.

SIMULATION 2 AND 1:CONTACT ENERGY

SIMULATION 1                                SIMULATION 2

In both cases, we see around 17ms, a sudden increase in contact energy. This is mainly due to buckling. In simulation2, this was captured better.

SIMULATION 1 AND 2:TOTAL ENERGY

  SIMULATION 1                             SIMULATION 2

 

  SIMULATION 1 AND 2:RIGID WALL FORCES

    SIMULATION 1                                                     SIMULATION 2

Initially, there was huge force, but due to internal deformation this force gradually decreased and friction force also got generated. The total resultant force is a sum of normal force and tangential force.

  SIMULATION 1 AND 2:VON MISES STRESS

       SIMULATION 1                                                            SIMULATION 2

Both simulations show the same value of VON MISES stress but node at which max force is acting got changed. The reason could be due to the recommended properties applied for 2nd simulation. At 716 node there is no stress in both simulations.

CONCLUSION:

With default properties, element deformation does not look natural. To overcome this, shell property was changed according to recommendations. QEPH ISHELL 24, with inbuilt HG stabilization, was used to capture deformations and transfer of stresses more effectively.

 

Radioss file link:

https://drive.google.com/file/d/1gmZ4m4_4lIipgBnZP38uyTw71gZlxmKQ/view

 

 

 


Projects by LUV KUMAR

OBJECTIVE: To simulate the side crash analysis of the given neon model. STEPS: UNIT SYSTEM: During whole analysis, Kg KN mm ms was followed. INTERFACES: Deleted already defined interfaces and created Type 7 global interface. The properties for the Type 7 interface a Read more

FRONTAL CRASH ANALYSIS
LUV KUMAR · 2019-10-31 04:05:45

OBJECTIVE: To perform a frontal crash simulation using Radioss as a solver. STEPS: View components: In the NEON FRONT model there are 75 components. ASSIGNING MATERIAL ID, PROPERTY ID, AND CARD IMAGE: For rigid bodies, card image is none and prop id, mat id is uns Read more

OBJECTIVE: Studying Interfaces and boundary condition in Hypercrash, Hypermesh and meshing bumper assembly MESHING OF BUMPER ASSEMBLY: Topological error like faces not connected was removed using extend option under surface edit option. Auto mesh under 2D option wa Read more

OBJECTIVE: Through this challenge, we will study different material laws in Radioss and will select the best material law for the given scenario by simulating and plotting graph. STEPS CASE1: Changing the model name from FAILURE_JOHNSON_0000 to Law2_epsmax_failure.R Read more

OBJECTIVE: To clean the geometry, extract mid surface, mesh the extracted surface and assign thickness to given models. LINK TO HYPERMESH FILES https://drive.google.com/file/d/1Op9RLc8y85U2qtlLyBNh5GYCglVHuNqN/view   PROCEDURE: Model was imported and to see Read more

OBJECTIVE: Surface mesh and volume mesh generation on given housing and arm bracket model. LINK TO HYPER MESH RADIOSS FILES: https://drive.google.com/file/d/1h0tOn4b9e9EDoVFqaGuKCjnWRgrHxUvL/view 3D MESHING: 3D Elements should be used when all dimensions are compara Read more

NEWTON RAPHSON METHOD
LUV KUMAR · 2019-09-13 07:11:16

OBJECTIVE: To find minimum pressure using Newton-Raphson method PROGRAMMING EXPLAINED: Two functions are defined. f and fprime. These function are defined to calculate next iteration value using NR method. To calculate pressure values for different thickness, for a Read more

File Parsing in Python
LUV KUMAR · 2019-09-13 06:25:32

OBJECTIVE: To demonstate file parsing using Python. when we run the program, first thing appears is  to enter the filename,if we speicify correct name then we can see the rest output otherwise a message will show file not found. To store different properties li Read more

SciPy Curve Fitting
LUV KUMAR · 2019-09-12 14:32:22

OBJECTIVE: To curve fit given data using linear and cubic polynomial. For curve fitting, we imported curve_fit from scipy library. Defined two function, func and func1 which are used for fitting curve. Defined another function read_file to read temperature and cp v Read more


Loading...

The End