FLOW IS TURBULENT  AS Reynolds no. changes velocity changes Re = rho*U_inlet*D/muD = cylinder diameter - 0.02mmu = dynamic viscosity - 1.849*10^-5 kg/m-srho = density - 1.184 kg/m^3data from (https://www.me.psu.edu/cimbala/me433/Links/Table_A_9_CC_Properties_of_Air.pdf)@ temperature 25 degrees Initail run with…

Forward Kinamatics Of 2R Robotic Arm Robotic arm is made of two link, link1 rotate 0 to 180 degree and link2 rotate 30 to 150 degree. %forward kinematics of 2R robotics clear all close all clc %inputs l1=0.7;%length of arm1 l2=0.3;%length of arm2 theta1=linspace(0,180,10);%arm1 angle of rotion with interval theta2=linspace(30,150,40);%arm2…

The NACA four-digit wing sections define the profile by: First digit describing maximum camber as percentage of the chord. Second digit describing the distance of maximum camber from the airfoil leading edge in tens of percents of the chord. Last two digits describing maximum thickness of the airfoil as percent of the…

import numpy as np #import numpy module for linspace command from scipy.integrate import odeint #import ode intrigate from scipy.integrate import matplotlib.pyplot as plt #import matplotlib.pyplot for graph import math #import math module #function to make ode function to be sloved def model(theta,t,b,g,l,m): theta1=theta[0]…

close all clear all clc R=8.314; f1=fopen(\'THERMO.dat\',\'r\'); %open thermo.dat from same directory of matlab code f2=fopen(\' Molecular mass.txt\',\'w\'); %open notepad save molecular mass and write line1=fgetl(f1);%fgetl is use to read first line of file for i=1:4 fgetl(f1) end T=linspace(300,1000,5000); for i=1:52…

Simulation shows Mass of air flow inside the cyclinder. Rate of flow inside the cyclinder varies by valve lift. Simulation calculates mass flow inlet and flow pattern inside cyclinder at various valve lift.   BOUNDARY CONDITIONS @inlet total pressure inlet @exit  static pressure    inlet pressure >…

import numpy as np #import module for array import matplotlib.pyplot as plt #import module for graph from scipy.optimize import curve_fit #import curvefit from scipy.optimize #function for linear , square and cubic fit def func1(t,a,b): return a*t+b def func2(t,a,b,c): return a*pow(t,2)+b*t+c def func3(t,a,b,c,d): return…

import math import matplotlib.pyplot as plt def engine_kinematics(bore,stroke,conn_rod,cr,start_crank,end_crank): import math bore = 0.1 stroke=0.1 conn_rod=0.15 cr=12 a=stroke/2 R=conn_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 = 50…

Genetic Algorithm is an optimization method based on the mechanics of natural genetics and natural selection.Genetic Algorithm mimics the principle of natural genetics and natural selection to constitute search and optimization procedures.GA is used for scheduling to find the near to optimum solution in a short time. In…

clear all close all clc %inputs b=0.1; %damping co-efficient g=9.81; %acceleration due to gravity l=1; %length of pendulum m=0.5; %mass of ball %intial condition theta_0 = [0;5]; %time points t_span = linspace(0,10,500); %time 0 to 10 sec with 500 intervals [t,results]=ode45(@(t,theta) odefunc(t, theta, b,g,l,m),t_span,theta_0);%ode45…

import math import matplotlib.pyplot as plt import numpy as np try: #open file and read if not able to read file not found values = open(\'engine_data.out\',\'r\') except FileNotFoundError: print( \"File not recognized. Please provide a valid CONVERGE output file\") data=[] # splitting the line and removing the line start…

The NACA four-digit wing sections define the profile by: First digit describing maximum camber as percentage of the chord. Second digit describing the distance of maximum camber from the airfoil leading edge in tens of percents of the chord. Last two digits describing maximum thickness of the airfoil as percent of the…

import matplotlib.pyplot as plt import numpy as np import plotly.graph_objects as go def function(p,h,beta,r,sigma): #equation term1 = pow(p,3)*(1-pow(beta,2)) term2 = (0.4*h*pow(beta,2) - (sigma*pow(h,2))/pow(r,2))*pow(p,2) term3 = (pow(sigma,2)*pow(h,4)*p)/(3*pow(r,4)) term4 = pow((sigma*pow(h,2)/(3*pow(r,2))),3) return…

A centrifugal pump is a mechanical device designed to move a fluid by means of the transfer of rotational energy from one or more driven rotors, called impellers.  Fluid enters the rapidly rotating impeller along its axis and is cast out by centrifugal force along its circumference through the impeller’s…

DESIGN OF PLANETARY GEAR Motion study of a Planetary Gear.A planetary gear set is made up of three types of gears; a sun gear, planet gears, and a ring gear. The sun gear is located at the center, and transmits torque to the planet gears which are typically mounted on a moveable carrier. The planet gears orbit around the…

MBD Simulation of Piston Assembly          Wrist Pin Offset ` Crank Speed (Rpm) 1              "0" mm 2000 2   10 mm Positive  2000 3   10 mm Negative  2000             Angular Velocity…

Aim : To simulates IC Engine Valve Train with two different Cam.(3mm & 6mm lift) Valve Lift.    The contact force between  Cam and Push Rod Pushrod and Rocker Arm Rocker Arm and Valve   The time for which the valves remained open was optimised only for single engine speed. However, as the engine…

Aim : Internal Geneva Gear Mechanism          Perform motion analysis by rotating the driver wheel at 10rpm and 20rpm at get results with Frames as per characterstics of simulation.  Optain contact forces and angular velocity of gears.   The Geneva Mechanism or Maltese cross…

 Aim:  Frequency analysis of rotating shaft. Introduction :Frequency analysis helps determine critical speed based on resonant frequency and control vibrations in system. Theory:The critical speed is the theoretical angular velocity that excites the natural frequency of a rotating object, such…