Objective: To write a program to simulate the transient behavior of a simple pendulum.   Program: #ODE import math import matplotlib.pyplot as plt import numpy as np from scipy.integrate import odeint b = 0.05 #damping co-eff g = 9.81 #gravity l = 1 #length of pendulum m = 1 #mass of the pendulum ct = 1 theta_0…

Objective: To clean up and mesh IP substrate models for structural analysis.   Procedure: 1. Model 1: The model was checked for any geometric errors and all topological errors were fixed. Quality criterions were set as stated in the problem using the "quality criteria" window. To extract the meshed midplane of…

Objective: To clean up and mesh a hood geometry with an element size of 5mm.    Procedure: 1. Inner hood panel: The inner hood panel was isolated and checked for any geometric errors. All the topological errors were fixed by using various commands such as \"topo\", \"delete\", \"coons\", \"paste\", \"release…

Objective: To compare the simulation results of the base setup and improved shell element properties.   Procedure: The Radioss starter file was imported and termination time (T_stop) which is located in eng run card was set as 55ms.  In order to get 20 animation steps (/ANIM/DT) the frequency value found under…

Objective: To mesh a cylinder model with element size as 1 unit Procedure: 1. Initial Set-up: The entity mode was changed to PID. \"Draw CONS\" option was ticked in the quality criterion window. Geometry was checked was any irregularities using check command.   2. PID Set-up: The cylinder model was split into different…

Objective: To clean up and mesh the BMW M6 model.    Procedure: 1. Fixing topological errors:    1.1. Styling lines were made visible using \"macros\" command and then deleted.    1.2. Logo: Overlapping surfaces were deleted and the logo was projected to the hood using \"project\" command.…

Objective:  To mesh a bumper assembly with a given mesh size. To compare the simulation results of a crash tube by making changes in contact properties and notches in the model. Procedure: Bumper assembly mesh: \'Toggle edge\' command under \'quick edit\' was used to connect edges. \'Extend\' command under \'surface…

Objective: To simulate forward kinematics of a 2R Robotic Arm.   Program: # 2R robotic arm import math import matplotlib.pyplot as plt L1 = 1 #Length of first arm L2 = 1 #Lemgth of second arm Theta_start = 0 #Start angle of both rotating arm Theta_end = math.pi/2 #End angle of both rotating arm n=10 #No. of divisions…

Objective: To create 2D and 3D mesh for the given models.   Procedure: 1. Housing: Hub: The hub model was 2D meshed using \'automesh\' command by giving the appropriate mesh criteria. Tria elements with an element size of 5mm were used for the mesh.         All 2D mesh errors were fixed by…

Objective: To compare the effect of using different material cards on a model for simulations. The model being tested is a flat plate subjected to impact from a sphere.   Model 1: (Input parameters) Law 2 (Johnson-Cook model) which is used for elastoplastic materials was used for model 1.  The input parameters…

Objective: To write a program to calculate cushion pressure needed to break a given thickness of ice with the Newton-Raphson method using Python.   Program: #Newton-Raphson import math import matplotlib.pyplot as plt import numpy as np r = 40 #size of air cushion b = 0.5 #width of ice wedge s = 150*144 #tensile strength…

Objective: To clean up and mesh a turbocharger model. Procedure:  1. Clean up:    1.1. Delete overlapping part:           The glitched section shown in the left figure occurs when 2 surfaces overlap each other. This can be solved by deleting one of the surfaces.    …

Objective: To create a data visualizer tool using python.   Program: #data analysis import math import matplotlib.pyplot as plt import numpy as np filename = input(\'Enter file name \') if filename.endswith(\'.out\'): crank=[] pressure=[] max_pressure=[] min_pressure=[] mean_temp=[] min_temp=[] max_temp=[] volume=[]…

Objective: To write a program to plot the PV diagram of an otto cycle.    Program: #otto cycle import math import matplotlib.pyplot as plt import numpy as np P1 = 101325 #presure at stage 1 T1 = 500 #temperature at stage 1 T3 = 2500 #temperature after combustion (stage 3) G = 1.4 #heat capacity ratio (gamma)…

Objective: To clean up the geometry and mesh three models.   Procedure: 1.Model 1 (Practice): Geometry clean up: The model was switched from \'auto\' mode to \'topo\' mode from the visualization toolbar for viewing the errors in geometry. \'Filler surf\' command under \'quick edit\' was used to patch missing surfaces.…

Objective: To perform crash analysis on Chrysler Neon side model.   Procedure:  The unit system used for this analysis was [Kg mm ms] A cylindrical rigid wall was created along the z-axis at the side of the car. The input parameters given are as follows: Slide=2 Dsearch=1000 Fric=0.1 Diameter=254 A self contact…

Objective: To perform crash analysis on Chrysler Neon frontal model.   Procedure: The unit system used for this analysis was [Kg mm ms] A planar rigid wall was created in front of the car along the x-axis with a friction value of 0.1. A self contact interface was created by setting the entire model as both master…

Objective: To write a program to perform curve fitting using a Cp vs Temperature data file.   Program: #curve_fit import math import matplotlib.pyplot as plt import numpy as np from scipy.optimize import curve_fit def fn_linear(t,a,b): #linear return a*t+b def fn_cubic(t,a,b,c,d): #cubic return a*pow(t,3)+b*pow(t,2)+c*t+d…