PERFORMING SIDE CRASH ANALYSIS OF A CAR ON A POLE
OBJECTIVES
- To understand the different Kinematic Conditions.
- To understand the different Loading Conditions.
- To assign recommended contact interfaces parameters for crash analysis.
- To understand the Time step Control for crash analysis.
- To learn about the Output Requests for post processing.
- To learn about the Check Debug for post processing
- To plot the different contours formed with use of various contact interfaces and parameter change.
Check with the unit system, either use [Mg mm s] or [kg mm ms]
- This can be done using the starter _le neon_front_0000.rad in notepad as shown below.

- Import the file neon_side_reduced_0000.rad
- It will appear in wireframe mode.
- Click on shaded mesh to view the surfaces in shaded form.

Create an appropriate interface, friction 0.2 and recommended parameters.
- There are different contact interfaces are present for each part. But it is good to have lower number of interfaces so that we will get faster simulation.
- Delete the available interfaces in the original model, and create a new one with the Type 7 interface.
- Recommended properties as follows

- To create the interface of Type 7 you can create in Hypermesh or Hypercrash
- To create in Hypermesh right-click > create >contact>Type 7, apply appropriate properties and friction as 0.2


- Igap (Determines how the size of the gap is calculated): 3 Variable gap + gap scale correction of the computed gap + size of the mesh is taken into account to
- Avoid initial penetrations during self-contact.
- Fscale_gap(Gap scale factor) : 0.8
- Gapmin(Minimum gap for activation of the interface ): 0.5
- Inacti(Action to take if initial penetrations exist): 6 Gap is variable with time but initial penetration is computed
- Istf(Affects how the stiffness of the interface is calculated): 4 Kmin=(Km, Ks)
- Iform(Friction formulation): 2 Stiffness
- Stmin(Minimum stiffness to use in the interface): 1
- Idel(What to do with slave nodes and master segments if an element fails (deleted) that they are attached to): 2 when an element is deleted, the corresponding segment is removed from the master side of the interface
- Fric(coefficient of friction ): 0.2


Make sure of no penetrations and intersections
- To run the analysis, we should make sure that there is no intersection and penetrations between two bodies of parts
- To check for penetrations and intersections in a hypermesh.



- As there are no penetrations and intersections available in a model. We can run the analysis. If the quality check shows any intersection or penetration we can fix it by using an auto fix tool.
- If the auto fix is failed to clear the penetrations or intersections then we have to correct as per the manual method.
Check rigid bodies if any issues
As checked in the model, there are no issues with the rigid bodies in the given model.
Create a rigid wall with friction 0.1
- We can create in both Hypermesh and Hypercrash
- To create a CYLINRICAL rigid wall in Hypermesh go to Create> Rigid Wall
- Select the type of the rigid wall you want as a cylinder, infinite plane, or sphere.
- For this application, we need an infinite plane type of rigid wall.
- Select the base node from which you want to create the rigid wall.
- Enter the value for D search as 1500 and friction as 0.1.


Comparing the weight with a full-scale model for weight balancing
- The weight of the full car model is 1216.55 kg as shown in below image.

- The weight of the reduced car side model is 166 kg.
- The center of gravity of the model should be in line with the C.G of the full car, because in order to compare with the full model we should have the C.G in the same place. The location of the C.G is at the position when the line from the bottom of the driver seat meets with the central line. So in order to get the C.G we have to arrange the mass in certain way.

- The total mass of the Full car model is 700kg, so we need to match these mass in our model and attain the C.G at Correct place by Distributing those mass in front and rear.
- The current model as the total mass of 166kg still 534kg has to be added to the model.
- We can add the remaining mass on the floor panel 240kg is added to the RH Side of the car and 144kg is added to the rear part of car side and 100Kg is added to frontal side for attempting 700Kg challenge, so that C.G lies inline from the bottom of the driver seat meets with the central line.

Apply initial velocity if 35mph
- As we are following the system [kg mm ms] we need to input the value of velocity in mm/ms.
- So after conversion the velocity 35mph= 15.65mm/ms.
- In Hypercrash go to Load Case > Initial Velocity > Select the whole model and apply velocity in X-Direction as 15.65mm/ms as shown in the below image.

- In Hypermesh go to Boundary Condiotions > Initial Velocity > Select the whole model and apply velocity in X-Direction as 15.65mm/ms as shown in the below image.

Use model checker to ensure good quality.
- Before going to the actual run we should ensure the quality of the model.
- In Hypermesh go to Tool> Model checker> Radioss Block > Then the system will show the error and warning. Try to find out the error and correct then or use the autocorrect tool.

- In Hypercrash go to Quality> Model checker> Run. Then the system will show the error and warning. Try to find out the error and correct then or use the autocorrect tool.

- Timestep changes, Runtime, Animation Frequency changes.

- A constant nodal time step
- T_scal:0.67
- Time step:0.5
- Termination Time: 80ms.
- Before running the simulation is better to do the model checker and make sure there is no error in the simulation
Output request
- Sectional force in the rails at the location of indicated node.
- Time History Output of the following nodes are required to measure Intrusion at B pillar, hinge pillar and fuel tank region
- 123561 in Fuel tank region
- 123475 B pillar,
- 1240 Hinge
- And time history output of node 337773 node in the side Door is also taken to measure the peak velocity of the door.

Section Creation
In Hypercrash go to Data History> section> select component> select 3 nodes to create a section.

Defined Sections
After Analysis Block

Debugging
- For debugging we have to see neon_side_reduced_0001.rad file

The Energy Error computed by RADIOSS is a percentage.
- If the error is negative, it means that some energy has been dissipated.
- Negative Energy Error since it is not counted in the energy balance. The normal amount of Hourglass energy is about 10% to15%.
- If the error is positive, there is an energy creation. In case of using QEPH shell formulation or fully integrated elements, the Energy
- Error can be slightly positive since there is no Hourglass energy and the computation is much more accurate. An error of 1% or 2% will be acceptable.
The Mass error is 0. 0.7647E-02 Kg which is acceptable as no change is happened to mass.
The Simulation time for frontal crash is much more.
- Action to take if initial penetration. Removal of initial penetration where possible. Elsewhere, reduce to less than 30% of defined gap value and adjust the gap by using in Inacti.
- When we are not using no failure criteria simulation time required is quite less as compared to others.
As per the above image, the energy error is -1.6% and mass error is 0.76E-02. So we can go for results as the values are in an acceptable range.




All Energies



- Section Force at B Pillar RH = 3.5 KN

- Displacement at node at B Pillar = 325mm
- Velocity at node at B Pillar = 3mm/ms

- Displacement at node at Hinge Pillar = 50mm
- Velocity at node at Hinge Pillar = 2 mm/ms

- Displacement at node at Fuel tank = 700 mm
- Velocity at node at Fuel tank = 5 mm/ms
- Peak Velocity at Inner Node on Door
- Peak Velocity at node at Door = 17 mm/ms

CONCLUSION:
- Since the Centre of gravity of the car and the mass addition plays a very critical role we matched the CG of a full car model to Get Proper results.
- The intrusion of the B-pillar Fuel Tank region and the hinge is been measured
- B pillar plays a major role in reducing the side-impact and protecting the passengers from its effect, it redirects the force to roof and the cross beams so that the shock is reduced and the passengers are protected.
- Side Impact bars in the door helps to reduce the intrusion in B-pillar and Fuel tank region
- Peak Velocity of Door Starts at 17mm/ms but we can't take that as Peak velocity Because that's is applied velocity of car at Starting , so the Peak Velocity can be considered only after the impact
- Peak Velocity of the Door is 17mm/ms
- The Mass error is 7647E-02 i.e 0.76% which is acceptable as no change is happened to mass.
- To plot the different contours formed with use of various material laws and parameter change.
- The sectional force in the bumper to rail for both RH and LH side plays an important role in absorbing the max impact.
- The role of shotgun design also plays an important role to keep the impact less transmitting to the interior parts of the vehicle.
- The max impact is taken by the frontal cross-sections hence the forces in the A-pillar section for both RH and LH side are very less.
- The position of the Center of Gravity of a vehicle is vital in the crashworthiness, so maintaining the weight balancing of a vehicle is the skill of design engineer.
- As we are using the half model and with added masses, if we can go for full vehicle body analysis then more realistic results can be obtained.
NOTE- .h3d file is having big size its not uploading. Kindly find it on following link
https://drive.google.com/file/d/17x3WKXui3OpYwBNmhqv0JQW0M4cUT7BX/view?usp=sharing