import numpy as np import matplotlib.pyplot as plt from scipy.optimize import curve_fit def funcl(t, a, b): return a*t +b def func3(t, a, b, c,d): return a*pow(t,3) + b*pow(t,2) + c*t + d def read_file(): temp = [] cp = [] for line in open(\'data\',\'r\'): values = line.split(\',\') temp.append(float(values[0])) cp.append(float(values[1]))…

  import matplotlib.pyplot as plt import math import numpy as np from numpy import trapz as trp try: file = open(\'engine_data.out\',\'r\') except: print(\' \"File not recognized. Please provide a valid CONVERGE output file\"\') file = \'engine_data.out\' row_count = 0 column_count = 0 count_header =0 column_title…

Aim: To perform static analysis on the plates with holes of different geometries and compare the results of stress , strain and displacement. Geometric Model of the plates: Case 1: Dimensions of the Plate : Length=300mmHeight=120mmThickness=30mmCircular Hole at the center:Diameter=60mm 3D Model of the Plate: Case 2: Dimensions…

Case 1 : For relaxation factor = 1   import matplotlib.pyplot as plt def f(p,h): x1 = pow(p,3)*(1-pow(b,2)) x2 = (0.4*h*pow(b,2)) - ((s*pow(h,2))/pow(r,2)) x3 = (pow(s,2))*(pow(h,4))/(3*pow(r,4)) x4 = s*pow(h,2)/(3*pow(r,2)) x5 = pow(x4,3) function = x1 + pow(p,2)*x2 + p*x3 - x5 return function def fprime(p,h): x11…

Aim: To compare the bending performance and suggest best out of the following three beams. Case 1 : Case 2: Case 3: Introduction: A beam is a structural element that primarily resists loads applied laterally to the beam's axis. Its mode of deflection is primarily by bending. The loads applied to the beam result in reaction…

Aim:  To create a  model of a rotating flywheel on a shaft and identify 5 critical frequencies by generating 5 mode shapes. Boundary Conditions: Material Used = Alloy Steel Fixtures = bearing supports on either side of the shaft at stepped part of the shaft A Shaft with rotating flywheel is model using solidworks…

clear all close all clc % Input Variables b=0.05; %damping coefficient g=9.81; %acceleration due to gravity l=1; %length of the string m=1; %mass of the pendulum bob %Initial_position theta_0 = [0; 3]; %Time t_span = linspace(0 ,10, 500); %Solve ODE [t,results] = ode45(@(t,theta) ode_func(t, theta, b, g, l ,m), t_span,…

Aim: To write a code in Matlab to parse the NASA thermodynamic data file and then calculate thermodynamic properties of various gas species. Theory: NASA came up with polynomials that can be used to evaluate thermodynamic properties such as Cp, H and S. They have also documented the coefficients that are required to evaluate…

Cylinder diameter = 0.02m Length = 0.05m Base Velocity= 10 m/s density of air= 1.225 kg/m3 dynamic viscosity= 1.81*10-5\">10−510-5 (at 15° c) Calculation:  Re = rho*V*D/mu D = cylinder diameter mu = dynamic viscosity rho = density   at 10m/s Reynold number=13535     Case 1: Increased…

import math import matplotlib.pyplot as plt def engine_kinematics( bore, stroke , con_rod , cr , start_crank , end_crank) : a = stroke/2 R = con_rod/a V_s = (math.pi/4)*pow(bore,2)*stroke V_c = V_s/(cr - 1) sc = math.radians(start_crank) ec = math.radians(end_crank) num_values = 100 dtheta = (ec-sc)/(num_values -1) V =[]…

Aim :- To perform a buckling analysis on a cyclonic separator stand.   To create a stiffener to the legs of the stand and determine the increase in buckling factor of safety. To perform parametric design study to run simulations at different positions of the stiffener and obtain the optimum position for the stiffener…

Design Model of IC Engine Valve Train : Cam : Push Rod: Rocker Arm:   Valve: Valve Mount: Assembly Model of IC Engine Valve Train: Simulation: Simulation is done for following two cases: Sl.No CAM Lift (mm) Speed (RPM)  Material  1 3.5 1500 Cast Carbon Steel 2 6 1500 Cast Carbon Steel During  the motion…