Assignment 7-Side Pole Crash Simulation Challenge

AIM :- 

To create a deck setup for the Neon side crash -BIW and requst for TH (Time History) plots

OBJECTIVES:-

• Check unit system and either follow[Mg mm s] or [Kg mm ms].
• Create appropriate interface ,friction 0.2 and recommended parameters.
• Make sure of no penetrations and intersection,Correct rigid bodies if any issues.
• Create rigid wall with friction 0.1.
• Added masses to reach target weight 700kg while getting CG about the required range.
• Initial velocity 35 mph.
• Time-step :0.5 to 0.1 microseconds.
• Run 80 ms.

Output requests:

• Sectional force in the cross member.
• Intrusion at B pillar,hinge pillar and fuel tank region.Provide recommendation on what can help to reduce Fuel tank intrusion.
• Peak velocity of inner node of the door.

THEORY:-

 what is crash ?

A crash test is a form of destructive testing usually performed in order to ensure safe design standards in crashworthiness and crashcompatability for various modes of transportation (see automobile safety) or related systems and components.

Crash is the explicit analysis. Explicit dynamics analysis is used to determine the dynamic response of a structure due to stress wave propagation, impact or rapidly changing time- dependent loads. Momentum exchange between moving bodies and inertial effects are usually important aspects of the type of analysis being conducted. This type of analysis can also be used to model mechanical phenomena that are highly nonlinear. Nonlinearities may stem from the materials, (hyperelasticity, plastic flows, failure), from contact (high speed collisions and impact) and from the geometric deformation (buckling and collapse). Events with time scales of less than 1 second (usually of order milliseconds) are efficiently simulated with this type of analysis. Based on all these conditions the user will select an explicit solver to solve the analysis. Typical examples of explicit analysis are crash analysis, drop test, bird strike, fan blade out…etc. Obtaining an answer from an explicit simulation can be simple — it can be more challenging to get the correct answer. The predictions of a simulation are usually very difficult to confirm using traditional techniques, making it crucial that the analysis is carried out carefully and thoroughly. Nonlinear analysis requires a lot of experience and a careful understanding of what you want to achieve, as well as a careful understanding of the software's analytical capabilities. When choosing a method for a dynamic problem, considerations must be made regarding the length of time the response can take place compared to the explicit method's stability limit. 

Types of crashes tested

There are three types of crash frontal , side , and roof crash

                                FRONTAL CRASH

         

                                  SIDE POLE CRASH                                                                                          SIDE CRASH

                                     ROOF CRASH

Types of SIDE crashes tested

There are two types of side crash

Side-impact tests: these forms of accidents have a very significant likelihood of fatality, as cars do not have a significant crumple zone to absorb the impact forces before an occupant is injured.

Pole-impact tests: A difficult test which places a large amount of force on a small proportion on the side of the vehicle

Here we are doing the side crash with pole (pole impact test) forces are going to act as shown in the below fig

METHODOLODGY:-

1. IMPORT THE FILE 
2. SET UP FOR THE FRONTAL CRASH OF NEON FRONT
3. OUTPUT REQUESTS
4. RESULTS and simulation
5. plotting graphs

 

PROCEDURE:-

1:-IMPORT THE FILE 

• Open the hyper mesh software and select the profile as radioss 
• Go to file and click on the import and select the import solver deck 
• click on the file option and select the file _0000.rad
  and click on the import as shown in fig 1.1

                        Fig 1.1

• After importing the model the model is as shown fig 1.2

                                                                   Fig 1.2

2:-SET UP FOR THE FRONTAL CRASH OF NEON FRONT

1 :- Check unit system and either follow[Mg mm s] or [Kg mm ms].

• In engine cards we can see the BEGINE_CARD there we can check the units followed in the model as shown in fig 2.1

                          Fig 2.1

• Units used in the model is [Kg mm ms]

2 :- Create appropriate interface ,friction 0.2 and recommended parameters.

• Here we have to create 2 contact interfaces one is for self interface and other one is to get forces near bumper area . Both interfaces are having Recommendations interfaces with gap mim 0.5 and friction 0.2
• And the recommended parameters are shown in the below fig 2.2

                                                                 Fig 2.2

• Here we are using only TYPE 7 interface for the type 7 interface these are parameters used
• To create the type 7 interface we can use hypermesh or hypercrash in hypermesh as shown in fig 2.3 

                                 Fig 2.3

• Parameters are set as recommonded shown in fig 2.4

• For the self contact interface we have to give master and slave as car components as shown in the fig 2.5 and Fig 2.6 
• Where the Blue colour is the master and red colour is the slave

                                                                             Fig 2.5

 

 

                                                                          Fig 2.6

3 :- Make sure of no penetrations and intersection,Correct rigid bodies if any issues.

• To check the penetrations and intersection we have to open the hypercrash 
• To open the hypercrash click on the applications and select the hypercrash and set the units to [Kg mm ms] an run
• And import the model by dragging the file into the hypercrash 
• Cllick on the quality and select the check all intersections the window shows as below fig 2.7

 

                           Fig 2.7

• There are no penetrations and intersection in the model
• And also there are no issues in the rigid bodies 

 

4 :- Create rigid wall with friction 0.1.

• Create Rigid wall by Solver > Create > Rwall > cylinder .

• Give sliding with friction of 0.1.

• Search distance 2000 (Distance for slave search).

 

• XM YM ZM are the coordinates the point where the rigidwall will be created 
• And the wall is created as shown in the fig 2.8

                                                             Fig 2.8

5 :- Mass Balancing: In order to maintain CG of car we have to add the mass for missing components. In given model there are so many parts are missing because of node limit in Studentversion.

• Check CG by Menu > Mass > Balancing > Show CG point. Before adding mass the CG point shown in the figure below

                                                                                                               CG point Before Mass adding fig 2.9

• Add the mass for missing components such as seat, doors, Back seats etc and change the CG point location, and ensure the total mass of car should be 700Kg  (Load case > Added mass)

                                                                                                               CG point After Mass adding fig 2.10

6 :- Initial velocity: In order move car towards  rigid wall. Given velocity 35mph convert to m/s so we have 15.646 m/s.

• Right click in tab area > creatw > I-L > Initial velocity,
• apply velocity to all components in Y direction shown in the fig. below.
• In order move car towards  rigid wall. Given velocity 35mph convert to m/s so we have 15.646 m/s
• Apply velocity to all components in X direction shown in the fig.2.12 below

                                                                                                               Initial velocity fig 2.11

7 :- Time-step :0.5 to 0.1 microseconds ,Run 80 ms.

• Input the Run time 80ms that is Final time for the run. in engine cards ENG_RUN
• Time step: Time step defaults for all elements which is 0.003ms (Solver > Create > Engine keywords > DT >DT) . If the time step of a simulation reaches ΔTmin">ΔTmin,then the following actions will happen depending on which element type is controlling the time step. For /BRICK or /QUAD elements, the simulation stops and a restart file is created. For /SHELL or /SH_3N elements, the element controlling the time step is deleted from the simulation. For all other element types, ΔTmin">ΔTminis not used.

8 :- Set up model checker 

• Before creating animation it is always best practice to eliminate error to get accurate results or some times due to these errors animation files dosen't gets created.(Hypercrach>>Quality>>Model check>>Run) As we see in  fig 2.12  there were o errrors & some warning  found while on running model checker. all warnings are not harmful but errors must be eleminated before running simulation. 

                                                 Fig 2.12

3:-OUTPUT REQUESTS

1. Sectional force in the cross member.
2. Intrusion at B pillar,hinge pillar and fuel tank region.Provide recommendation on what can help to reduce Fuel tank intrusion.
3. Peak velocity of inner node of the door.

1. Sectional force in the cross member.

• We have to create the sectional forces on cross members
• First we have to create moving frames on the indicated nodes 

FRAME :- Frame is a type of reference coordinate system where the quantities defined follow its own origin.Here, unlike the moving skew, the moving frame is attached with the node which is selected as origin.Calculates the motion of the body relative to the frame defined.Frames are used only in notable options like imposed velocity, imposed displacement, TH/node etc.,

• To create a frame Solver > frame > mov 

• We have to select the origin node that which node is indicated in the output requests and X-AXIS or Z-AXIS and XY PLAIN we have to give that which the frame should move fame is created as shown in the fig 3.1

                                                               Fig 3.1

• And then we have to create the sections on the indicated node

• A section is a set of nodes and a set of elements. A section is a cut in the structure, where forces, moments or energy will be computed and stored in output files (using /TH/SECTIO). In Radioss /SECT, /SECT/PARAL and /SECT/CIRCLE can be used to define a section.

• To create the sections path shown in fig 3.2

                                     Fig 3.2

• We have to select the N1 ,N2 ,and N3 as which the direction of the car is moving as shown in fig 3.3

]

                                                                              Fig 3.3

• In that we have to select the elements around to create a section over there as shown in the fig 3.4

                                                                                 Fig 3.4

• And same for the all indicated nodes where we have to create the sections as shown in the fig 3.5

                                                                         Fig 3.5

2. Intrusion at B pillar,hinge pillar and fuel tank region.Provide recommendation on what can help to reduce Fuel tank intrusion.

• We are goping to create the spring elements on the B pillar,hinge pillar and fuel tank region as shown in the fig 3.6

                                                                    Fig 3.6

• By crerating the springs on the B pillar,hinge pillar and fuel tank region we can know the intrustions on the perticular region 
• The stiffness and mass of the spring should be very less which is not considerable 

   3. Peak velocity of inner node of the door.

• We have to find the peak velocity of inner node of the door 
• So that we are going to calculate the peak velocity of the node 337773 as shown in the fig 3.7

                                                                        fig 3.7

• To get the results outputs we have to create the TH for all the sections ,spring elements , and node 

 4:-RESULTS and simulation

• Before running the simulation make sure to delete the unused and run the model checker TOOLS>MODEL CHECKER>RADIOSS BLOCK
• If there is any errors or incompatable kinamatic contitions fix the error and run the simulation  ANALYSIS>RADIOSS>RADIOSS  simulation shown in GIF 4.1

                                                                                                                  GIF 4.1

• After running the simulation check for the energy error and mass error as shown in fig 4.1

                                                                      Fig 4.1

• energy error is -1.2% and mass error is 0.7606E-02 the model is good
• And the time taken for the simulation is 3778.97 s and number of cycles is  80001

• We can see these results in the 0001.out file
• Displacement is shown in GIF 4.2 and vonmisses in GIF 4.3 of the model

                                                                                         GIF 4.2

                                                                                     GIF 4.3

5:-plotting graphs

• To plot graph change it to the HYPER GRAPH 2D and select the T01 file

5.1 :- Graph 5.1 shows the  IE ,KE ,CE ,HE ,TE to time

                                                                       Graph 5.1

• Initially IE is zero and its increasingw,r,t to time until the the complete crash happens. because  at initial stage we have maximum velocity,Kinetic energy is directly proportional to the Velocity. so kinetic energy is max at initial stage after hitting the KE is observed by plate so internal energy increase until the component reaches complete crash after that there is no energy observed (internal energy ΔU = q+w)
• we have  Total Energy = Internal Energy + Kinetic energy + External Work, Here TE is not depending on external work because there is no any force acting on car BIW only the car is moving with velocity.as IE increasing and KE decreasing in same path so the TE is remains Constance. TE is slightly decreasing because of Contact energy i.e Contact energy is the energy acting opposite to the system to prevent the penetration acting on the system due to applied velocity.
• Hourglass energy is constant because both are controlled by elements formulation. here hourglass is controlledby Physical stabilization i.e changing material level in ISHELL is 24 QEPH Shell formulation.

5.2 :- SECTIONAL FORCES GRAPHS

• Force which is coming from the Doors and ots frame is received by Cross members. to find the this force we have sections defined. From the section forces graph, 

• We have created the 2 sections on 2 cross members 
• Graph 5.2 shows the forces on the sections on the cross members 

 

                                                                                Graph 5.2

• Max force on Cross section 1 is 5.9871 KN at 16ms
• Max force on Cross section 2 is 4.4289 KN at 19ms

5.3 :- Intrusions at B pillar,hinge pillar and fuel tank region

• Intrusion is change in space or decreasing length horizontally of car driver compartment. so due to Intrusions we can say safety of driver. intrusion is measured by creating spring element with less mass and stiffness so the one can check spring compression, from compression will give change in length of spring. in the fig. 11.6 shows the graph for spring compression as Intrusion.

                                                                                Graph 5.3

• 9900037 spring elements isHinge piller max intrustion is 504.821 mm
• 9900038 spring elements is under B piller max intrutions is 802.983 mm
• 9900039 spring elements is under fuel tank region max intrustion is 776.721 mm

5.4 :- Peak velocity of inner node of door

• To know the velocity of particular node so that we can study effect of sudden impact on velocity. here because impact the node 337773 is accelerated and increased velocity to 16.1601 m/s at 5.5 ms. here node is away from the pole center and at the time of impact its acts like centrifugal force of mass m, Radious R (from center of rigid pole) Due to this angular velocity will added so that increased velocity can see. shown in the Graph 5.4 Peak velocity is 16.1601 m/s.

                                                                                Graph 5.4

CONCLUSION :- 

• Intrusion of frame near A pillar is maximum ,Intrusion frame ner  driver seat is minimum & Angle of collision  plays important role as in this case it was ner perpendicular so intrusion of frame near a pillar is observed maximum. 
• We can see the forces on members are more 
• Centre of Gravity plays important role in calculation of forces which can be adjusted by adding mass