Assignment 3-2D Element Formulation Challenge

2D Element Formulation Challenge

AIM: To compare the results with base simulation and improved shell element properties. By using a crash beam file we have to simulate the model with base simulation and compare it with improved shell element properties.

PROCEDURE:

• Import the given crash beam file to the Hypermesh window. And change the view to shaded elements to see the meshed component clearly.

• Now, to change the run time to 55ms, go to Cards Collector >> ENG_RUN >> Tstop >> change to desired value ( here the given default value was 60.0 ). It has been changed to 55ms from 60ms

• Now, the number of animation steps during has to be set between min.25 to max.60. We select the value of 50 animation steps. Now, to set this 50 animation steps, the Tfreq value has to be changed.
  
  Tfreq = ( Total run time / no. of. animation steps required )
  
  So, Tfreq = (55/50) = 1.1
  
  Now, to set the Tfreq value, go to ENG_ANIM_DT >> Tfreq >> Set the desired value (Given default was 1.0, It was changed to 1.1)

• With the Tstop, Tfreq and animation steps set, we can now export the corresponding file to the radioss solver and proceed with the analysis to check the energy error and mass error. Make sure to enter the ( -nt 4 ) In the options tab (to access number of cores in the hardware which is hardware specific).

• Now, while checking the results, we can observe that the energy error starts from -0.0% and goes upto -10.3% and the Mass error remains the same 1.659E throughout the simulation which is within the range of acceptance.

• Now, we can proceed with Ploting rigid wall forces, Internal energy, Hourglass energy, Contact energy and total energy of the simulation.
  For this, we need to export the output file and view them in the Hyperview and Hypergraph windows to see the animation and the energy plots.
  
  
  To Open the animation file, we need to search for .h3d file and open it in the Hyperview tab to view the animation. Once it is opened we can see the Contour Options to View the various phenomenon, Here we need to check the Displacment and the Hourglass energy. Firstly the Hour glass energy is observed.

• we can see that the maximum displacement was about 761mm.

• Now, to see the Global Variables, we can split the window and open up Hypergraph in the other window. Here we need to import the T01 ( Time History File ) to view all the global variables like Hourglass energy, internal energy etc.
  
  Once the File is imported, Fix time along X axis and set the global variables in the Y axis and observe the results :

• We can see that the Kinetic Energy is maximum at the starting since a huge force is applied to the crash tube which possess certain mass and then it gradually decreases since the component deforms. Similarly, we can see that the Internal energy is minimum while starting and it gradually increases as the component starts to deform and absorb the energy.
• Since the elements used were of reduced integration, there are hour glass energy which increases from the starting. The Total energy also slightly decreases from the beginning and this is due to the energy error and the hour glass energy produced.

Rigid wall force with time

• We can see that the Total Resultant Force is in a non linear manner. There are certain spikes in between which indicates certain spike in force which is hitting the rigidwall.

Now, Let's see what the results are when we change the Properties of the elements as recommended and check the results.

Here are the changes made to the Properties:

Parameters		Comment
Ishell =24		QEPH 4 nodes shells are best combination of cost and accuracy.
Ismstr=2		Full geometric nonlinearities (default)with possible small strain formulation activation in RADIOSS Engine.
Ish3n=2		Standard 3 noded shell (C0)with modification for large rotation(default)
N=5		Number of integration points set to 5 for accuracy bending.
Ithick=1		Thickness change is taken into account for accuracy.
Iplas=1		Iterative plasticity for good accuracy.

• Once all these were changed, and the Tstop set to 55ms and the Tfreq set to 1.1 ( for 50 animation steps ) now let's run the solver and check the mass and energy error.

• Changing the values for the Hat section and close out panel properties, only for, Ishell, Ismstr, Ish3n, N, Ithick, Iplas with improved shell element properties.

• Here we can see that the Energy error is -0.1% which is a great sign and the mass is almost the same 1.659E. So, this is a well accepted result. Further we can proceed with the contour plots and the global variables plot and cross check the results with the previous condition.

• To check the Contour plots and animation, Open the .h3d file in the Hyperview and switch to contour option. Let's see the displacement created in the simulation.

• We can see that the maximum displacement is found to be 781mm. Which is 20mm more diplacement than base simulation.

• Now, Open the T01 ( time history file ) of this phenomena in the Hypergraph where the element properties were changed to check the Global variables and the Resultant forces.
  
  After Opening the T01 file, Fix Time along x-axis and Fix the global variables Internal energy, Kinetic energy, Hourglass energy and Total energy along y-axis and plot the graph.

• We can clearly see that there is no Hour glass energy, Kinetic energy falls and the Internal Energy gains throughout the time period. The Total Energy seems to be consistent throughout the Crashing.

Rigid wall force with time

• We can see that the Total Resultant force is slightly above the base simulation Resultant force and it is behaving in non linear manner.

1. TOTAL ENERGY : It is evident that the Total energy is consistent in the second case whereas in the default property the Total energy was decreasing. This is due to the fact that the Energy error changed from -10.3% to -0.1% when the element properties were changed.

2. HOUR GLASS ENERGY : The Hourglass energy doesn't exist in the second case whereas the default element has some hour glass energy. This is due to the type of integration wherein the default follows reduced integration method while in the other case it is Improved Reduced Integration wherein it has the capacity of adding artificial stiffness directly to the material to avoid the Hour glass energy.

3. CONTACT ENERGY : It is almost the same in both the cases.

4. INTERNAL ENERGY : The Internal energy also increases throughout the simulation in both the cases.

5. RIGID WALL FORCES : The Resultant rgidwall force was dropped in the second case than the resultant force in the default element property.

6. DISPLACEMENT : There is a slight increase in maximum displacement of elements in the second case which is roughly about 20mm when compared to the default property elements. This might be due to the better integration method followed in the recommended element property.

• RESULT :
  
  Thus, the simulation has been carried out with both the default element property and with recommended element property and results like global variables, rigid wall force, hourglass energy, energy errors, total energy, displacement were compared.