CODE : import matplotlib.pyplot as plt import math import numpy as np beeta=0.5 r=40 sigma=150 def func(p,h): term1=pow(p,3)*(1-pow(beeta,2)) term2=(0.4*h*pow(beeta,2)-sigma*pow(h,2)/pow(r,2))*pow(p,2) term3=pow(sigma,2)*pow(h,4)/3/pow(r,4)*p term4=pow(sigma*pow(h,2)/3/pow(r,2),3) eq=term1+term2+term3+term4 return eq def…

1. What does popt and pcov mean?  Popt = An array that holds the value of the coefficients of variables in a function. pcov = An array that stores the covariance (cov) of the coefficients of the variables in the given function. Covariance is the variance between the two variables required of tge curve fit. 2.What…

CODE : \"\"\" OTTO CYCLE SIMULATOR \"\"\" import math import matplotlib.pyplot as plt import numpy as np def engine_kinematics(bore, stroke, con_rod, cr , start_crank, end_crank): import matplotlib.pyplot as plt # Geometric parameters a = stroke / 2 R = con_rod / a # Volume parameters V_s = math.pi*(1/4)*pow(bore,2)*stroke…

Frequency Analysis of a rotating shaft  Model :   Frequency Analysis : Steps   Create new study and select frequency study. Select Stainless steel material for shaft. Select Aluminium material for flywheel. Fix both ends using the fixture geometry option. Mesh to a gridsize of 5mm. Run the simulation.…

CODE: R=8.314 f=fopen(\'THERMO.dat\',\'r\'); Header=fgetl(f); global_temp_range=fgetl(f); global_temp_range=str2num(global_temp_range); line1=fgetl(f); line2=fgetl(f); line3=fgetl(f); for i=1:53 line0=fgetl(f); chemical_element=strtok(line0); a=strfind(line0,\'G \'); local_temp=line0(a+1:end-1); local_temp=str2num(local_temp);…

CODE : \"\"\" ENGINE DATA ANALYSIS \"\"\" \"\"\" 0-Crank 4-Mean_Temp 8-Mass 12-C_p 16-Gamma 16.Dyn_Visc 1-Pressue 5-Max_temp 9- Desnsity 13-C_v 2-Max_Pres 6-Min_temp 10-Integrated_HR 14-Gamma 3-Min_Pres 7-Volume 11-HR_Rate 15-Kin_Visc \"\"\" Enter_x_parameter=8 Enter_y_parameter=9 \"\"\" Compatability Check \"\"\" if Enter_x_parameter…

CODE: import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt import math # function that returns dz/dt def model(theta,t,b,g,l,m): theta1=theta[0] theta2=theta[1] dtheta1_dt=theta2 dtheta2_dt=-(b/m)*theta2 - (g/l) * math.sin(theta1) dtheta_dt = [dtheta1_dt,dtheta2_dt] return dtheta_dt b=0.02…

CODE : b=0.5; l=1; m=1; g=9.81; theta_0=[0;4]; t_span=linspace(0,20,400); [t,result]=ode45(@(t,theta)ode_function(b,g,l,m,t,theta),t_span,theta_0); figure(1) hold on plot(t,result(:,1),\'color\',\'r\'); plot(t,result(:,2),\'color\',\'b\'); xlabel(\'Time\') ylabel(\'Displacement / velocity\') legend(\'Displacement\',\'velocity\')…