Parsing NASA thermodynamic data

AIM: To Simulate the Birdstrike on Aero-engine using Explicit Analysis in LS-Dyna. KEYWORDS: Include, Elastic, Non-Linear Plasticity, Boundary conditions, Elform, etc., SOFTWARE: LS-Dyna OBJECTIVES: 1. The given file contains all the 4 parts which as to separate and save part wise. 2. The suitable material…

OBJECTIVE: To calculate the Mooney Rivlin and Ogden material constants and compare the both using stress-strain data from a dogbone specimen tensile test with 100 percent strain. The given material data is the engineering stress-strain in MPa/(mm/mm). The comparison should be shown from the d3hsp file and using simulation.…

Link to the report

PFAhttps://docs.google.com/document/d/1dokHwqnAJ2L53b0WeqtpwUWlgnhLG4XQUnU1xKUxwt8/edit?usp=sharing    

PFA   https://docs.google.com/document/d/1Z7EfgNIC3iXig0Hu92_wsRWGU6cJ3jAJsmY203fPaG4/edit?usp=sharing  

PFA https://docs.google.com/document/d/1YBD6ltYcw3Oe_LIcaDxwy-eHrFkUXIG4noaXD8ounKQ/edit?usp=sharing

PFA, https://docs.google.com/document/d/1xn778ASSeILLstC-IfoqDC0gICGJj9HI-hWqRYh5Scs/edit?usp=sharing

PFA, https://docs.google.com/document/d/1PKGREE4PWnmuHHvGXD4TOqd6Rg10lj98OQTANRGKTuQ/edit?usp=sharing

https://docs.google.com/document/d/15P6gLgD_lQc2IoyIm9ytG7ubYsVds4cBUNfzDhUnSIM/edit?usp=sharing   Aim: To conduct the Drop test on Solid block in ls Dyna Procedure: Material Definition          *MAT_020 - Which is used to make the part as rigid by checking the CMO as 1…

 Aim: solving the given equation by using both explicit and Implicit method   Solution: import matplotlib.pyplot as plt import numpy as np def diff(u): #a=((u**3)/3 + u) #b=(u**2)+1 a=((u**3)+(9*u**2)+(4*u)) b=((3*u**2)+(18*u)+4) return [a,b] u=0 f1=0.01 error=1 u_1=[] de_u=[] f=0 f_int=[] f_int.append(0) u_1.append(0)…

https://docs.google.com/document/d/1oTxftwfRjC6zDAdVhaPS9Xa8Dc_MXVT2-NitJTwr9nA/edit?usp=sharing  

Aim:      To write a code and simulation of forward kinematics of 2R robotic arm   Introduction:             l1=length of link 1 l2=length of link 2 θ1=inclination angle of link 1 θ2=inclination angle of link 2 x0,y0 origin points (x0,x1),(y0,y1) coordinate…

# # set a [tk_getOpenFile -title "select the file"] # set file [open "test.xml" r] # for {set i 0} {$i<18} {incr i} { # gets $file # } # cd "C:/Users/Skill-Lync/Downloads/cONNECTOR PROJECT" # set g [open "test_file.txt" w] # while {[eof $file] == 0} { # set val {} # set a [gets $file] # set val [maximum $a] # puts $g…

i

By using tcl we will generate an industry-standard report       set f [open "test3.csv" w] set systemTime [clock seconds] puts $f "Username,$env(username)" puts $f "Time, [clock format $systemTime -format %H:%M:%S]" puts $f "Date, [clock format $systemTime -format {%A the %d of %B %Y}]" puts $f "Solver profile,[hm_getsolver]"…

Aim:         To write a program to simulate the transient behaviour of a simple pendulum and to create an animation of it's motion. In the above equation, g = gravity in m/s2, L = length of the pendulum in m, m = mass of the ball in kg, b=damping coefficient.  (d^2 θ/dt^2) +(dθ/dt)*b/m+(g/l)*sin⁡θ=0…

 Aim:         To write a code for data analysis of engine file and to find the spcific fuel consumption,power output of the engine   Code for data analysis   import math import numpy as n1 import matplotlib.pyplot as plt import sys try: open('engine_ .out','r') print('hi') except : print('File…

Aim:    To write a Matlab program to calculate drag force against a cyclist.   Code for drag force with velocity clear all close all clc rho = 1.2; a = 0.1; cd = 0.8; v = 1:50; f = rho * a * v.^2 * cd/2; plot(v,f,'LineWidth',2) xlabel('VELOCITY') ylabel('DRAG FORCE')            …

hello

  Aim:         To Write a Python program to calculate drag force against a cyclist.    Case:1 Drag force based on varying velocity   import matplotlib.pyplot as plt import numpy as np a=2 vs=[5,6,7,10,12,15,18] cd=0.8 p=78 f=[] for v in vs: f.append(cd*0.5*p*v*v*a) plt.plot(vs,f)…

AIM:       To simulate Rankine cycle using the MATLAB       1-2   The fluid is pumped to the boiler at high pressure and it is isentropic process 2-3   Heat is supplied on the fluid through the boiler 3-4   The fluid is expanded through isentropically through…

Aim:     To write a code in Matlab to parse the NASA thermodynamic data file and then calculate thermodynamic properties of particular gas species in the file.   Code for file parsing:   clear all close all clc f=fopen('THERMO.dat') for i=1:5 a=fgetl(f); end k=1; c=0; i=0 while(~feof(f)) while(i<5)…

Aim:    To write a code for the PV diagram of the otto cycle and to calculate thermal efficiency of the cycle using MATLAB   Introduction of OTTO cycle     Process (0-1) Air intake piston moves TDC to BDC Process (1-2)Isentopic compression Process (2-3)Constant volume heat addition Process (3-4)Isentopic…

Aim:      to write a code for minimum pressure for breaking ice with air cushion vehicle and to find optimum relaxation factor               p=Cushion pressure               h=Thickness of the ice field          …

Aim:    To write a program in Matlab to simulate the forward kinematics of a 2R Robotic Arm     l1=length of link 1 l2=length of link 2 θ1=inclination angle of link 1 θ2=inclination angle of link 2 x0,y0 origin points (x0,x1),(y0,y1) coordinate of link 1 (x1,x2),(y1,y2) coordinate of link…

Aim:       To write code in MATLAB to optimize the stalagmite function and find the global maxima of the function Introduction of ga Algorithm:       In a genetic algorithm, a population of the solution to an optimization problem is evolved toward better solutions. Each…

Aim:         To write a program to simulate the transient behaviour of a simple pendulum and to create an animation of its motion.   In the above equation, g = gravity in m/s2, L = length of the pendulum in m, m = mass of the ball in kg, b=damping coefficient.  (d^2 θ/dt^2) +(dθ/dt)*b/m+(g/l)*sin⁡θ=0…

Objective:      To compare the bending performance for three types of beams i.e Rectangular beam, I section beam and C section and also compare the results of stress, displacement and strain for all the three case and to find the best of three   Case 1 (Rectangular beam)   Boundary condition :…

Aim:    To write Matlab code for curve fitting   Code for the curve fitting   clear all close all clc cpdata= load('data'); temp = cpdata(:,1); cp = cpdata(:,2); plot(temp,cp,'linewidth',4) hold on abc = polyfit(temp,cp,3); dupcp = polyval(abc,temp); abc = polyfit(temp,cp,1); ducp = polyval(abc,temp);…

Aim:        To Minimize the following function        5-(x-2)^2 -2(y-1)^2         subject to the following constraint         x + 4y = 3     u(x,y)=5-(x-2)^2-2(y-1)^2 const=x+4y-3 f(x,y)=u(x,y)+ʎ(const) f(x,y)=5-(x-2)^2-2(y-1)^2+ʎ(x+4y-3)…

Aim:         To write a program to simulate the transient behaviour of a simple pendulum and to create an animation of its motion.   In the above equation, g = gravity in m/s2, L = length of the pendulum in m, m = mass of the ball in kg, b=damping coefficient.  (d^2 θ/dt^2) +(dθ/dt)*b/m+(g/l)*sin⁡θ=0…

 Aim:         To write a code for PV diagram of OTTO cycle and to calculate thermal efficiency of the cycle using PYTHON   Introduction of OTTO cycle     Process (0-1) Air intake piston moves TDC to BDC Process (1-2)Isentopic compression Process (2-3)Constant volume heat addition…