Assignment 7 - Side Crash Neon BIW
Side Crash Neon BIW
Aim: To create deck setup, analyse the results with mass balance for Side Crash Neon BIW and request for time history plots.
Objective:
1. To check unit system and either follow[Mg mm s] or [Kg mm ms]
2. To create an appropriate interface ,friction 0.2 and recommended parameters.
3. To make sure of no penetrations and intersections.
4. To correct rigid bodies if any issues.
5. To create rigid wall with friction 0.1.
6. To compare the model weight with the full scale 300k nodes model and use added masses to reach target weight 700kg while getting CG about the required range.
7. To create initial velocity 35 mph.
8. To input Time-step :0.5 to 0.1 microseconds and Run 80 ms
Introduction:
Crash worthiness is one of the core topics in the passive safety of vehicles and plays an important role in the condition that the impact cannot be avoided. Generally, the analysis of crashworthiness is based on the related crash responses, i.e. the displacement, velocity and acceleration, of critical parts of a vehicle in full car crash tests. However, these tests are required appropriate facilities, one or more cars with measuring devices, experienced staff and a long time to prepare. It means they are complicated, expensive, long-lasting and therefore not easy to realize. This is especially true in the early stage of vehicle design. Therefore, vehicle designers and researchers made a lot of effort to build numerical models to describe the crash processes. Up to now, various technologies are used to model the vehicle crash. Typical crash models may be classified into three broad categories
1) Detailed nonlinear finite element models.
2) Multi body models and multi body based lumped parameter models.
3) Functional approximation or response surface models.
In the proceeding of the study, a piece wise linear model is proposed to represent the vehicle-rigid wall Side pole crash. Compared to existing models, this model is developed based on the analysis of crash responses and therefore can reflect the performance of vehicle structures in crashes.the vehicle structure would transmit the impact load in a proper way, i.e. load-carrying path. In addition, there are some weak components arranged as crumple zone. During the crash, the crumple zone will deform and absorb energy and ensure the rest parts of the vehicle are safe. In the local design, the material and shape of beams, shells and connectors are also considered to achieve better energy-absorbing performance [1].
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.
Setup Procedure
1. Import the file:
In the data we have two files
1.1. neon_front_0000.rad (starter input file): Which has all data required for the solver such as model information, elements, nodes, loads, boundary condition etc.
1.2. neon_front_0001.rad (engine input file): Which has all data required for the radioss engine i.e speed , time step, run time, what all data we want as output.
To import file,
a. change user profile to radioss Because we are using solver as Radioss
b. File>Import>Solver deck. we can see the new tab where options are available for browsing files. Browse the file FIRST_RUN_0000.rad (starter input file),
c. Click on import so we can see imported file with required for the solver such as model information, elements, nodes, loads, boundary condition etc. Show in the fig. 1.1 (both hypermesh and Crash
Fig.1.1: Import file
2. Setup for crash: Here some setup are done by Crash and remaining setups are done by Hypermesh.
2.1. Check the unit system: To check for the unit system used in the model by clicking on model>cards>BEGIN_CARD it's in Kg-mm-ms.
3. Create interface: Here required 2 contact interfaces one is for self interface and other one is to get forces near the bumper area where it is difficult to create sections. Both interfaces have Recommendations interfaces with gap mim 0.5 and friction 0.2. 0.2.
3.1. Self interface: Self interface is created in HyperCrash. To create
a. Right click on tab area>create new > A-c> Contact Interface >Type 7
b. Edit the and apply all Recommendations interfaces. and select the self impact option..
3.2. Interface for Bumper: Similarly create Contact interface. In fig 3.1 shows the both interfaces.
Fig.3.1: Contact Interface
4. Check Initial penetrations and intersection: To check Click on Quality>check Intersection on tree > select parts > click on check penetration. In the figure 4.1 can see the procedure and results
Fig.4.1: Initial penetrations and intersection
5. Mass Balancing: In order to maintain CG of a car we have to add the mass for missing components. In the given model there are so many parts missing because of the node limit in the Student Version.
a. 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
b. 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
6. Initial velocity: In order to move the car towards a rigid wall. Given velocity 35mph converts to m/s so we have 15.646 m/s.
a. Right click in tab area > create > I-L > Initial velocity,
b. apply velocity to all components in Y direction shown in the fig. below.
Initial velocity
7. Run Time And Time step:
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">ΔTmin is not used.
8. Create rigid wall: Create Rigid wall by Solver > Create > Rwall > PCYL.
Give sliding with friction of 0.1.
Search distance 1000 (Distance for slave search).
XM, YM, ZM is the location where we need Rigid wall, the axis is in Z direction.
Rigid Pole
8. Create Sections: Sections are created in components to get the exact forces at a particular section. Total 2 sections are required in cross members.
To create a section, we need to define a moving frame.
Create a moving frame: Use the Systems panel to create rectangular, cylindrical and spherical coordinate systems. Use this function when you want to define nodes, loads and constraints in a different coordinate system. the following steps to create moving frame by Create by Node Reference Sub panel.
a. Solver > frame > mov
b. Select the coordinate system type you want to create.
c. Select the node that defines the origin by selecting it in the modeling window, or by entering the node ID.
d. x-axis node: Select the node that defines the x direction by selecting it in the modeling window, or by entering the node ID.
e. xy-plane node: The following Output requests are asked in the given problem.Select the node that defines the xy plane by selecting it in the modeling window, or by entering the node ID.
f. relative size: Specify how large, in model units, systems are displayed relative to the model size.
In fig. 8.1 shows the moving frame
Fig.8.1 Moving Frame
Create Section: 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.
Generally, to define the section, the following points are required:
a. A cutting plane (which defined either by 3 points or by parallelogram or by a circular disk)
b. A reference point to compute forces
c. A direction of the section
d. In /SECT, the cutting plane is defined by a group of elements and its orientation by a group of nodes.
Review the section by pressing 'Q'
Similarly define all the sections shown in the figure 8.2
Fig.8.2 Section
In the figure 9.1 Shows the complete Setup model and ready for simulation.
Fig.9.1: Complete model
10. Output requests: In output requests we are going to request all required output which we need. To request Solver > Create > TH > select things to record in result.
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.
Run simulation Analysis > Radioss
Simulation Results:
In fig 10.1 shows the simulation for Side pole crash of car BIW.
Fig. 10.1: Side pole crash
From the simulation we can see that. Initially B Pillar hit the pole with a velocity of 15.646 m/s. The B pillar is going to start hitting then the forces are transferred to the left and right and forward direction in roof rails and cross members at bottom. In each stage the forces are going to reduce because energy absorbed in each stage. .
Load-carrying Path and Crash Process: Generally, vehicles have a high relative velocity in the crashes. For this reason,
11. Graphs:
11.1 Energy: Energy, in physics, the capacity for doing work. It may exist in potential, kinetic, thermal, electrical, chemical, nuclear, or other various forms.
In FEA we have Total Energy, Internal Energy, Kinetic energy, Contact energy and hourglass energy
Total Energy (TE): Total Energy is the sum of Internal Energy, Kinetic energy & External Work
Internal Energy (IE): The internal energy of a thermodynamic system is the energy contained within it. It is the energy necessary to create or prepare the system in any given internal state. It does not include the kinetic energy of motion of the system as a whole, nor the potential energy of the system as a whole due to external force fields, including the energy of displacement of the surroundings of the system. It keeps account of the gains and losses of energy of the system that are due to changes in its internal state
Kinetic energy (KE): In physics, the kinetic energy (KE) of an object is the energy that it possesses due to its motion. It is defined as the work needed to accelerate a body of a given mass from rest to its stated velocity.
Contact energy (CE): Contact energy is the energy acting opposite to the system to prevent the penetration acting on the system due to applied load case.
Hourglass energy (HE): This is usually done by adding extra stiffness or viscosity across the diagonal which adds extra strain energy into the system and we call this as Hourglass Energy.
Energy Errors:
Energy error is 0 to -2.4 % and mass error is 0% this is due to energy being dissipated from the system. In case of under integrated elements, the Hourglass Energy can explain a negative Energy Error since it is not counted in the energy balance.
Results:
11.2. In Side pole crash we have following energy graphs.
1. Initially Car is at rest so there is no any energies, once car moving with velocity of 15.646 m/s towards Side pole for Side pole crash with maximum KE and IE of 86244.3 `(Kg-mm^2)/(ms^2)`
2. When a car is hitting the rigid wall the velocity of every part is going to drop w.r.t to time with different magnitude which is dependent on the distance from the rigid pole. shown in the fig. below.
3. Initially IE is zero and it's increasing w,r,t to time until the complete crash happens. because at the 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 increases 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 the car BIW only the car is moving with velocity.
as IE increases and KE decreases in the same path so the TE remains Constant. 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.
4. Hourglass energy is constant because both are controlled by element formulation. Here hourglass is controlled by Physical stabilization i.e changing material level in ISHELL is 24 QEPH Shell formulation.
In the fig. 11.2 shown energy graph
Fig. 11.2: Energy Graph
From the graph we have
a. KE max= 86244.3 Kg-mm^2/ms^2 at time 0 ms.
KE min= 63103.9Kg-mm^2/ms^2 at time 79.5 ms.
b. Internal energy max= 21760 Kg-mm^2/ms^2 at time 79.5 ms^2
Internal energy min= 0 Kg-mm^2/ms^2 at time 0ms^2
c. Contact energy max= 1015 Kg-mm^2/ms^2 at time 79.5 ms^2
d. Hourglass energy max = 67.8 Kg-mm^2/ms^2 at time 79.5 ms^2
11.3 Rigid Pole forces: Rigid pole forces are the reaction forces for the applied forces.
1. Initially The rigid pole forces, once the car BIW with velocity 15.646 m/s reaches rigid pole it starts producing the reaction forces, in graph we can see there is sudden increasing of rigid pole forces due to B pillar touching to rigid pole and it decreasing slightly.
2. Rwall forces are increasing to max (64.5KN at 3.5ms) value when all the forces are transferring to RH and LH roof rails and Roof Cross members, once rails start deforming, Rwall forces are decreasing because Rails are absorbing energy. also velocity decreasing due to decreasing velocity the rigid pole reactions forces are decreasing.
3. Remaining energy absorbed by Roof and floor panel cross members. we can see less Rwall forces are dropping down because velocity dropping. from 60ms the Rwall forces are increasing because all maximum deformation are already happen and and remaining forces are distributed to Rwall shown in the fig 11.3.
Fig. 11.3: Rigid Wall forces.
11.4. Forces At Sections: Here we are going to review the Load-carrying Path and forces at particular section, we have total 2 sections.
11.4.1. Axial force received Cross members: Force which is coming from the Doors and its frame is received by Cross members. to find this force we have sections defined. From the section forces graph,
a. Force received Cross member 1 is 6.5KN @ 15 ms
b. Force received at Cross member 2 is 4.7KN @ 17 ms
After that the section forces are produced because all parts touching the Rigid pole are deformed max amount, remaining forces are distributed again to the Cross members.
Fig. 11.4.1: Axial force received Cross Members
11.6. 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 the driver. intrusion is measured by creating spring elements 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.
Initially all springs have a length of 1401.307 mm in fig. below we have shown elements ID
from graph we have compression of spring with IDs
9900037 is compressed 891mm (63%) fuel tank region
9900038 is compressed 914mm (65%) Max at below B pillar
9900039 is compressed 615mm (43%) Door below
9900040 is compressed 400mm (28%) hinge pillar
Here we have more Intrusions. It is not acceptable.
To reduce Intrusions at fuel tank region we can provide cross member,
Fig. 11.6: Intrusions Graph.
11.7. Peak velocity of inner node of the door: To know the velocity of particular node so that we can study the effect of sudden impact on velocity. here because impact the node 337773 is accelerated and increased velocity to 16.2 m/s at 5.5 ms. Here the node is away from the pole center and at the time of impact it acts like centrifugal force of mass m, Radius R (from center of rigid pole) . Due to this angular velocity will be added so that increased velocity can be seen. shown in the fig. 11.7 Peak velocity is 16.2 m/s.
Fig. 11.6: Peak velocity of inner node of the door
Conclusion: The overall objective of the work was to simulate a Side Crash and validate the results of the simulations obtained from the crash-test. Simulation was performed using the Radioss software package.
The graphical results obtained all showed that the test models behavior were shown
Due to the limited number of nodes in the student version some of the biw, engine blocks, suspension system etc are removed, so results are not as in real condition.
It has been observed that there is Maximum deformation of the cabin and also there was Maximum intrusion at given nodes in the driver compartment. it's not acceptable. to prevent intrusion by providing cross members at top and bottom.
References
[1] A MODEL OF VEHICLE-FIXED BARRIER FRONTAL CRASH AND ITS APPLICATION IN THE ESTIMATION OF CRASH KINEMATICS, Zuolong, Wei Hamid Reza, Krimi Kjell Gunnar, Robbersmyr University of Agder Norway Paper Number 15-0161