Assignment 3-2D Element Formulation Challenge

AIM: To run the simulation for the given component with and without changing the recommended element parameters using a raioss solver.

OBJECTIVE:

• To change the run time to 55ms and to change the number of animation steps to 60.
• To run the basic simulation with default element parameters.
• To check the energy error and mass error for the base simulation.
• To plot the graphs for kinetic energy,internal energy,hourglass energy,total energy,contact energy and rigid wall forces.
• To run new simulation with recommended element parameters and to plot the graphs again.

Recommended 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.

PROCEDURE:

Phase 1- Import the starter file into the solver deck radioss

• While operating hypermesh software. A user profile window will pop up. Switch to the radioss user profile as shown in below figure and start importing the model into the solver deck.

• Now after switch the user profile to radioss.import the model (starter file) into the radioss solver deck.
• To import the model, go to standard panel>>Import>>Import solver Deck.

• Import solver deck option should be selected to import the radioss starter file

• Here the import browser will appear as shown above figure. select the appropriate starter file to import into GUI.
• Switch the file type to RADIOSS block and import the model into GUI.

• Here the model is in wireframe mode,switch to the shaded mode as shown in below.
• To switch the shaded mode,go to visualization tab>shaded elements and mesh lines.

shaded element and mesh lines mode

Phase 2- Edit ENG_RUN File [Tstop] and ANIM_DT [Tfrequency]:

• As per our requirement,we need to change the Tstop value and Tfrequency value.
• To change the values,go to model browser>control cards>edit the values by selecting the appropriate control card.

1) Run time [Tstop]

• Run time is the total time taken for the sound wave to travel through the entire length of the model.
• Time for sound to travel length of rail formula is deltat=L/C

Where delta t= time step

L= length of the rail

C= speed of the sound

Here the Tstop=55ms value has been already given in the question,so need to calculate.

Enter the value 55 in ENG_RUN File which is shown in the below fig.

2) ANIM_DT[Tfrequency]:

As per our requirement,engine file have to generate between 25-60 animation files.

Here we are taking value 60,the solver will create 60 animation files within the termination time.

Tfrequency= Runtime/ Number of animation steps

Tfrequency= 55/60

Tfrequency= 0.9166

Now go the model browser>select appropriate control card >ANIM_DT> enter value as 0.9166 as shown in the figure.

Phase 3- Run the base simulation with and without changing the recommended shell parameters

case 1- without changing the recommended shell parameters.

• Run the base simulation with default shell parameters as shown below,don't assign recommended shell parameters.

• Now run the simulation,to run simulation,go to analysis panel as shown in fig.
• Go to analysis panel>radioss>select the input file>save it in different folder and rename it as case1>run.
• Check the include connectors,if there are any connectors in the model.The connectors will also be taken ino account.
• Type -nt 4 in options tab,this will make the simulation faster.
• Where nt indicates no of threads,4 indicates assigning the task to 4 cores in the system.

• Here save the input file in a different folder,and rename it as case-1 which is shown in below.
• Now run the solver by hitting on radioss.

• After completing the simulation,the radioss will pop up a solver window stating Radioss job completed which indicates the simulation has been completed which is shown below fig.

• Now verify the Tstop and Tfrequency values are properly assigned in engine.rad file [00001.rad]
• Open engine .rad file [00001.rad] with notepad as shown in below fig.

• Next verify the required number of animation steps obtained in the file location or not.
• In the begining the specified number of anmation steos was assigned as 60.
• Here the number of animation steps obtained are 61 which is shown in below fig.

Phase 4- Do energy error and mass energy error checks

• Now engine output file [00001.out] has been generated after analysis.
• The engine output file [00001.out] contains the data like

1. Time step
2. Energy Error
3. Mass Error
4. Internal energy error
5. Kinetic energy error
6. contact energy error

• Here check the mass error and energy error.
• The obtained values for energy error,mass error,internal energy error,kinetic energy error and contact energy error has been shown in below fig.
• The acceptable energy error range is -15% to +5%.
• The acceptable mass error rnage is 0 to 2%.

Scroll down to see the maximum value obtained for the energy error,kinetic energy error,internal energy error and mass error as shown in below figure.

• The maximum energy error is = -10.3%
• The maximum internal energy error = 0.1932 x
• The maximum kinetic error = 0.3563 x
• The maximum mass error = 0.1653 x

Note:

1. The acceptable energy error range is -15% to +5%.
2. The acceptable mass error rnage is 0 to 2%.

Now these two error values are acceptable,they are not going beyond the standard values.

Phase 5- Postprocessing

1) Review the simulation using- Hyperview.

2) Plot the graphs using- Hypergraph 2D

Review the simulation using Hyperview:

• Hyperview allows for loading and viewing result obtained from several sources.
• Based on the solver type of the files and the results you would like to visualize and analyze,there are different ways to load the input deck and their corresponding result into hyperview.
• First to begin in the postprocessing in the hypermesh split the screen as shown below.
• Import the animation file.h3d into hyperview.

• And then activate the client Hyperview

• To access the load model panel.
• Select load model button from the hyperview panel and open .h3d file as shown in below fig

Opening .h3d File.

After loading the model into hyperview,it will be represented as shown in fig.

• After importing the .h3d file into the GUI. enable the contour.
• The contour tool create contour plots of a model graphically visualize the analysis results.
• To enable contour. Go to results toolbar> contour.

• Now switch to the von misses -stress in result type and select the component,select the averaging method as simple and then click apply as shown in below.

• After applying,Run the simulation,the simulation animation is shown in below.

2) Plot the grapth using- Hypergraph 2D

• Now plot the graph using hypergraph 2D,we are plotting the graphs to see what is happening in the rail component.
• Hypergraph 2D is a powerful data analysis and plotting tool with interfaces to many popular file formats.
• It is sophisicated math engine capable of processing even the most complex mathematical expressions.
• Hypergraph 2D combines these features with hogh quality presentation output and customization capablities to create a complete data analysis system for any organization.

• Here switch to the hypergraph 2D to plot the graphs.
• To switch,go to client selector>>choose the hypergraph 2D as shown above.

• The first graph is plotted for the rigid wall forces.
• To plot the graph,go to hypermesh 2D>> data file>> element_formulation-shell-3_assignmentT01>> Apply.

• Here the graph has been plotted for the rigid wall forces as shown in below fig.
• There are various parameters in rigid wall like normal force.Tangent force.FX-X Total force,FY-Y total force,FZ-Z Total force,total resultant force.
• We are selecting only total resultant force because we dont know the exact axis of the component placed,so that's why we are selecting total resultant force and then plotting the graph.

• While plotting the graph,units profile window will pop up,there select solver units and then plot the graph because the graph will be plotted according to the solver units.

Case 1- Without assigning the recommended parameters in the property.

Rigid Wall Forces:

• Total resultant force=
• Rigid wall forces are the forces which are acting on the rigid wall.
• Here the rigid wall forces is at higher value in the begining.
• In begining there is no impact in the simulation,but in the graph the rigid wall forces are at the higher values.
• This is because,the solver didn't record in that time,so the rigid wall force is at higher value in the begining.
• When the rail component hits the rigid wall,the rigid wall forces increases from the begining.
• It goes on decreasing,because the deformation is happening to the rail component when it hits to the rigid wall.

[Note: The solver starts to record when the rail component hits on the rigid wall,this is the reason why in the begining,the rigid wall forces is higher value or increases and goes on decreasing.]

Next plot the graphs for the kinetic energy,internal energy,total energy,contact energy and hourglass energy.Those graphs are plotted and shown in the below fig.

Kinetic Energy:

Kinetic energy is at peak in the starting,why because we have huge mass which is shown in below.

To check mass,go to tool>>mass details>>mass,mass,COG,Inertia.

To check velocity,go to solver brower>>INIVEL Card>>Initial velocity.

• Kinetic energy is lower and decreases,why because there is a velocity applied to the component,so its getting deformed and the kinetic energy decreases.
• The formula for kinetic energy= 1/2 MV*2
• Given data
• Mass= 503881.28179g
• Velocity= -15.6ms
• KE= 1/2 X (503881.28179)(-15.6)*2
• KE = 61312274.37
• The value calculated here is shown in the fig.
• With the calculation,we came to know why kinetic energy is increasing in the begining itself,why because the initial velocity is given for the component.
• The kinetic energy decreases,why because the velocity decreases( the rail component starts to deform,after complete deformation,it will be in rest) when the component comes to a rest,so velocity becomes zero and the kinetic energy decreaes.

Internal Energy:

• The formula for internal energy is E=Q+W
• Here the heat is neglected,cause there is no heat transfer,we will be having only workdone.
• W=FXD
• Here the internal energy is in zero,because there is no deformation initially,so the internal energy is in zero and starts from zero and goes on increasing.
• The internal energy increases because the dformation is happening there is a displacement,so the internal energy increases,when the displacement or deformation occurs.

Total Energy:

• Here total energy starts from the higher value,in starting itself its increasing why because the total energy is sum of kinetic energy+contact energy+hourglass+internal energy.
• All energies are in initial condition except kinetic energy,so kinetic energy is only there.
• Total energy= kinetic energy+0+0+0.
• Kinetic energy is directly proportional to the total energy.
• SO the toatal energy is increasing in the begining.
• Due to 10.3% of energy error,the total energy is decreasing.

• Why it is decreasing.let's see mathematically
• energy error=10.3%
• 10.3/100 x61314900.
• =0.103 x 61314900
• = 6315434.7
• Now subtract this value with maximum total energy error.
• Total energy error maximum value = 61314900
• =61314900-631543.7
• =54999465.3
• We got an minimum value which is shown in above fig.

Contact Energy:

• Here the contact energy starts from origin,cause the rail component is in intial condition.
• Intially there is no contact between the rail component and rigid wall.
• There is also no deformation intially but when the component goes and hits on the rigid wall,the contact energy starts increasing.
• When the deformation or dispalcement happens to the component,the contact energy increases cause the crash tube or rail component comes within the contact to the rigid wall,so the contact energy goes on increasing.

Hourglass Energy:

• delta t> t - hourglass exsist.
• delta
• Where deltat- shockwave
• t- time taken for the information.
• Here the hourglass energy is starting from origin because the recommended parameters has not been specified.
• The hourglass energy goes on increasing because to balance the rest of energies.An artificial energy is required that is called hourglass energy and thats's the reason why hourglass energy goes on increasing from origin.
• A lot of hourglass energy is needed,so its goes on increasing.

Case 2- With assigning the recommended parameters in the property.

• Now enter the recommended parameters in properties called Hat-section,close-out section which is given in the question.
• The recommended parameters have been entered in the properties as shown below fig.

• Similrly follow the same process and do the same what we have done for the case 1.
• Now run the simulation for the case 2.

• Now go open the engine output file to check whether there is improvement or not.
• There is improvement when compared to the values obtained from case 1.
• Here the energy and mass error values are acceptable and start to run the simulation.

• Now run the simulation,to run simulation we need a additional window.In the main page menu,go and split the window into 8 as shown in fig.
• In any new window switch to hyprview and then load the .h3d file.The loaded .h3d file is shown in fig.

• Now run the simulation and in the result type switch to the von misses stress as shown in below fig.

• Now plot the graphs for the energies, to see the changes in the results with respect to the case 1 plot the graphs for all the energies.
• Here the graph has been plotted for the rigid wall forces which is shown in below fig.
• Similarly plot the graphs for all the energies.

• Here there is little bit improvement after assigning the recommended shell parameters when compared to case 1.
• Here, there is no hourglass effect due to recommended parameters,so the rigid wall forces curve is again increasing,if you compare the case 1 and case 2 simulation animation properly,you will get an idea.
• The graph has been plotted for all the energies like kinetic energy,internal energy,contact energy,total energy anf hourglass energy which is shown in below fig.

• Here there is a little bit of improvement in all the enrgies when comapred to the case-1
• Now plot the graph for the hourglass energy to see the result,whether there is any improvement in hourglass energy or not.

• Here the hourglass energy is constant,why because the recommended parameter Ishell= 24[QEPH 4 noded shell]
• It will be having 4 nodes in the element,so hourglass won't happen,it will be constant.
• In previous case,it was under integrated node,so there was hourglass deformation.

Final Image-

Here in the case 2 graphs which is shown below,there is no much change.

Comparison of case 1 and case 2 simulation.

S.No.	Characteristic	case-1	case-2
1.	Maximum energy error	-10.3%	-0.1%
2.	Maximum kinetic energy	0.3563E+08	0.3941E+08
3.	Maximum internal energy	0.1932E+08	0.2186E+08
4.	Total energy	Toatal energy increases and slighly decreases due to energy error	Total energy remains constant
5.	Hourglass energy	Hourglass energy increases due to balance all the energies	Hourglass energy remain constant due to recommended shell parameters.

Result:

• Hence the simulation has been runned without assigning recommended shell parameters.
• Hence the simulation has been runned with assigned recommended shell parameters.
• The graphs have been plotted for case-1 and case-2 successfully.
• While comparing case-1 and case-2,case-2 is effective,cause there is no hourglass deformation.

Conclusion and learning outcome:

• I conclude that case case-2 with assigning recommended shell parameters is effective.
• In this challenge,i came to know.
• How to run a simulation with and without assigning shell parameters.
• How to do postprocessing (plotting the graphs for the various parameters).
• How to understand the curves that are increasing and decreasing in the graphs.
• How to compare the obtained values with the various parameters.