initial commit / add all books

22 files changed, 154 insertions, 0 deletions

diff --git a/548/CH5/EX5.01/5_01.sce b/548/CH5/EX5.01/5_01.sce
new file mode 100755
index 000000000..e7e1be089
--- /dev/null
+++ b/548/CH5/EX5.01/5_01.sce
@@ -0,0 +1,10 @@
+pathname=get_absolute_file_path('5_01.sce')

+filename=pathname+filesep()+'5_01data.sci'

+exec(filename)

+L=q*c*Cl;disp(L,"L=","L=q*c*Cl","Lift per unit span L:")

+D=q*c*Cd;disp(D,"D=","D=q*c*Cd","Drag per unit span D:")

+M=q*c*Cm*c;disp(M,"M=","M=q*c*Cm*c","Moment per unit span M:")

+printf("\Answer:\n")

+printf("\n\Lift about the quarter chord,per unit span: %f N\n\n",L)

+printf("\n\Drag about the quarter chord,per unit span: %f N\n\n",D)

+printf("\n\moment about the quarter chord,per unit span: %f N.m\n\n",M)

diff --git a/548/CH5/EX5.02/5_02.sce b/548/CH5/EX5.02/5_02.sce
new file mode 100755
index 000000000..3209722ac
--- /dev/null
+++ b/548/CH5/EX5.02/5_02.sce
@@ -0,0 +1,6 @@
+pathname=get_absolute_file_path('5_02.sce')

+filename=pathname+filesep()+'5_02data.sci'

+exec(filename)

+printf("\Answer:\n")

+printf("\n\angle of attack for 700 N lift: %f degree\n\n",a)

+printf("\n\angle of attack for zero lift:: %f degree\n\n",a1)

diff --git a/548/CH5/EX5.03/5_03.sce b/548/CH5/EX5.03/5_03.sce
new file mode 100755
index 000000000..c72fba469
--- /dev/null
+++ b/548/CH5/EX5.03/5_03.sce
@@ -0,0 +1,6 @@
+pathname=get_absolute_file_path('5_03.sce')

+filename=pathname+filesep()+'5_03data.sci'

+exec(filename)

+Cp=(P1-P)/q;disp(Cp,"Cp=","Cp=(P1-P)/q","pressure coefficient Cp :")

+printf("\Answer:\n")

+printf("\n\pressure coefficient at this point of wing: %f \n\n",Cp)
\ No newline at end of file
diff --git a/548/CH5/EX5.04/5_04.sce b/548/CH5/EX5.04/5_04.sce
new file mode 100755
index 000000000..29758bbf1
--- /dev/null
+++ b/548/CH5/EX5.04/5_04.sce
@@ -0,0 +1,6 @@
+pathname=get_absolute_file_path('5_04.sce')

+filename=pathname+filesep()+'5_04data.sci'

+exec(filename)

+Cp=(P1-P)/q;disp(Cp,"Cp=","Cp=(P1-P)/q","pressure coefficient Cp :")

+printf("\Answer:\n")

+printf("\n\pressure coefficient : %f \n\n",Cp)
\ No newline at end of file
diff --git a/548/CH5/EX5.05/5_05.sce b/548/CH5/EX5.05/5_05.sce
new file mode 100755
index 000000000..fda942924
--- /dev/null
+++ b/548/CH5/EX5.05/5_05.sce
@@ -0,0 +1,6 @@
+pathname=get_absolute_file_path('5_05.sce')

+filename=pathname+filesep()+'5_05data.sci'

+exec(filename)

+Cp=Cpo/(sqrt(1-M^2));disp(Cp,"Cp=","Cp=Cpo/(sqrt(1-M^2))","pressure coefficient Cp :")

+printf("\Answer:\n")

+printf("\n\pressure coefficient : %f \n\n",Cp)
\ No newline at end of file
diff --git a/548/CH5/EX5.06/5_06.sce b/548/CH5/EX5.06/5_06.sce
new file mode 100755
index 000000000..628c60be9
--- /dev/null
+++ b/548/CH5/EX5.06/5_06.sce
@@ -0,0 +1,6 @@
+pathname=get_absolute_file_path('5_06.sce')

+filename=pathname+filesep()+'5_06data.sci'

+exec(filename)

+P1=(q*Cp)+P;disp(P1,"P1=","P1=q*Cp+p","pressure at this point P1:")

+printf("\Answer:\n")

+printf("\n\pressure at this point : %f N/m^2\n\n",P1)
\ No newline at end of file
diff --git a/548/CH5/EX5.07/5_07.sce b/548/CH5/EX5.07/5_07.sce
new file mode 100755
index 000000000..aa4217a34
--- /dev/null
+++ b/548/CH5/EX5.07/5_07.sce
@@ -0,0 +1,7 @@
+pathname=get_absolute_file_path('5_07.sce')

+filename=pathname+filesep()+'5_07data.sci'

+exec(filename)

+V2=V*((Cp1-Cp2)+(V1/V)^2)^0.5;

+disp(V2,"V2=","V2=V*((Cp1-Cp2)+(V1/V)^2)^0.5","velocity at point 2 V2:")

+printf("\Answer:\n")

+printf("\n\Velocity at point 2: %f m/s\n\n",V2)
\ No newline at end of file
diff --git a/548/CH5/EX5.08/5_08.sce b/548/CH5/EX5.08/5_08.sce
new file mode 100755
index 000000000..d67f76e2d
--- /dev/null
+++ b/548/CH5/EX5.08/5_08.sce
@@ -0,0 +1,6 @@
+pathname=get_absolute_file_path('5_08.sce')

+filename=pathname+filesep()+'5_08data.sci'

+exec(filename)

+Cn=integrate('1-0.95*y','y',0,1.0)-integrate('1-300*y^2 ','y',0,0.1)-integrate('-2.2277+2.2277*y','y',0.1,1.0)

+printf("\Answer:\n")

+printf("\n\Normal force coefficient : %f \n\n",Cn)

diff --git a/548/CH5/EX5.09/5_09.sce b/548/CH5/EX5.09/5_09.sce
new file mode 100755
index 000000000..e2ed1cc69
--- /dev/null
+++ b/548/CH5/EX5.09/5_09.sce
@@ -0,0 +1,6 @@
+pathname=get_absolute_file_path('5_09.sce')

+filename=pathname+filesep()+'5_09data.sci'

+exec(filename)

+Cl=Co/(sqrt(1-M^2));disp(Cl,"Cl=","Cl=Co/(sqrt(1-M^2))","Lift coefficient Cl :")

+printf("\Answer:\n")

+printf("\n\Lift coefficient at Mach 0.7: %f \n\n",Cl)
\ No newline at end of file
diff --git a/548/CH5/EX5.10/5_10.JPG b/548/CH5/EX5.10/5_10.JPG
new file mode 100755
index 000000000..c7b55ca9c
--- /dev/null
+++ b/548/CH5/EX5.10/5_10.JPG
diff --git a/548/CH5/EX5.10/5_10.sce b/548/CH5/EX5.10/5_10.sce
new file mode 100755
index 000000000..285a6a71c
--- /dev/null
+++ b/548/CH5/EX5.10/5_10.sce
@@ -0,0 +1,15 @@
+pathname=get_absolute_file_path('5_10.sce')

+filename=pathname+filesep()+'5_10_data.sci'

+exec(filename)

+clf();

+i = 1;

+while(i<=length(M)) 

+    Cpcr(i)=(2/(y*M(i)^2))*[[(2+(y-1)*M(i)^2)/(y+1)]^(y/(y-1))-1]

+    Cpmin(i)=Cpomin/sqrt(1-M(i)^2);

+   i = i+1; 

+end 

+xlabel("Mach Number");

+ylabel("Cp");

+plot2d(M,Cpcr,2);

+plot2d(M,Cpmin);

+disp("The intersection point of both the graphs i.e approx 0.74 is the critical Mach no of the NACA-0012 airfoil.")

diff --git a/548/CH5/EX5.11/5_11.sce b/548/CH5/EX5.11/5_11.sce
new file mode 100755
index 000000000..1bea26bfe
--- /dev/null
+++ b/548/CH5/EX5.11/5_11.sce
@@ -0,0 +1,8 @@
+pathname=get_absolute_file_path('5_11.sce')

+filename=pathname+filesep()+'5_11_data.sci'

+exec(filename)

+L=q*c*Cl;disp(L,"L=","L=q*c*Cl","lift per unit span for mach 3 :")

+Dw=q*c*Cd;disp(Dw,"Dw=","Dw=q*c*Cd","Wave drag per unit span for mach 3 :")

+L1=q1*c*Cl1;disp(L1,"L1=","L1=q1*c*Cl1","lift per unit span for mach 2:")

+Dw1=q1*c*Cd1;disp(Dw1,"Dw1=","Dw1=q1*c*Cd1","Wave drag per unit span for mach 2 :")

+

diff --git a/548/CH5/EX5.12/5_12.sce b/548/CH5/EX5.12/5_12.sce
new file mode 100755
index 000000000..7cefdef1e
--- /dev/null
+++ b/548/CH5/EX5.12/5_12.sce
@@ -0,0 +1,10 @@
+pathname=get_absolute_file_path('5_12.sce')

+filename=pathname+filesep()+'5_12_data.sci'

+exec(filename)

+a=L*(M^2-1)^0.5/(4*q*S);

+disp(a,"a=","a=L*(M^2-1)^0.5/(4*q*S)","angle of attack at sea level:")

+a1=L*(M^2-1)^0.5/(4*q1*S);

+disp(a1,"a1=","a1=L*(M^2-1)^0.5/(4*q1*S)","angle of attack at 10 Km:")

+printf("\Answer:\n")

+printf("\n\angle of attack at sea level: %f degree\n\n",a*180/%pi)

+printf("\n\angle of attack at 10 Km: %f degree\n\n",a1*180/%pi)
\ No newline at end of file
diff --git a/548/CH5/EX5.13/5_13.sce b/548/CH5/EX5.13/5_13.sce
new file mode 100755
index 000000000..1786eae93
--- /dev/null
+++ b/548/CH5/EX5.13/5_13.sce
@@ -0,0 +1,7 @@
+pathname=get_absolute_file_path('5_13.sce')

+filename=pathname+filesep()+'5_13data.sci'

+exec(filename)

+L=q*S*4*a/sqrt(M^2-1);

+disp(L,"L=","L=q*S*4*a/sqrt(M^2-1)","Lift exerted on airplane L:")

+printf("\Answer:\n")

+printf("\n\Lift exerted on airplane: %f N\n\n",L)
\ No newline at end of file
diff --git a/548/CH5/EX5.14/5_14.sce b/548/CH5/EX5.14/5_14.sce
new file mode 100755
index 000000000..46c7ea514
--- /dev/null
+++ b/548/CH5/EX5.14/5_14.sce
@@ -0,0 +1,7 @@
+pathname=get_absolute_file_path('5_14.sce')

+filename=pathname+filesep()+'5_14data.sci'

+exec(filename)

+Cl=L/(q*S);

+disp(Cl,"Cl=","Cl=L/(q*S)","Lift coefficient Cl:")

+printf("\Answer:\n")

+printf("\n\Lift coefficient: %f \n\n",Cl)
\ No newline at end of file
diff --git a/548/CH5/EX5.15/5_15.sce b/548/CH5/EX5.15/5_15.sce
new file mode 100755
index 000000000..210c2b446
--- /dev/null
+++ b/548/CH5/EX5.15/5_15.sce
@@ -0,0 +1,8 @@
+pathname=get_absolute_file_path('5_15.sce')

+filename=pathname+filesep()+'5_15data.sci'

+exec(filename)

+Cdi=Cl^2/(%pi*e*AR);disp(Cdi,"Cdi=","Cdi=Cl^2/(%pi*e*AR)","induced drag coefficient Cdi:")

+Di=q*S*Cdi;disp(Di,"Di=","Di=q*S*Cdi","induced drag Di:")

+printf("\Answer:\n")

+printf("\n\induced drag coefficient: %f \n\n",Cdi)

+printf("\n\induced drag: %f N\n\n",Di)
\ No newline at end of file
diff --git a/548/CH5/EX5.16/5_16.sce b/548/CH5/EX5.16/5_16.sce
new file mode 100755
index 000000000..23062367e
--- /dev/null
+++ b/548/CH5/EX5.16/5_16.sce
@@ -0,0 +1,6 @@
+pathname=get_absolute_file_path('5_16.sce')

+filename=pathname+filesep()+'5_16data.sci'

+exec(filename)

+Dt=(Cd+Cdi)*S*(D*V^2/2);disp(Dt,"Dt=","Dt=(Cd+Cdi)*S*q","total drag Di:")

+printf("\Answer:\n")

+printf("\n\Totl drag: %f N\n\n",Dt)
\ No newline at end of file
diff --git a/548/CH5/EX5.17/5_17.sce b/548/CH5/EX5.17/5_17.sce
new file mode 100755
index 000000000..aebe84a61
--- /dev/null
+++ b/548/CH5/EX5.17/5_17.sce
@@ -0,0 +1,8 @@
+pathname=get_absolute_file_path('5_17.sce')

+filename=pathname+filesep()+'5_17data.sci'

+exec(filename)

+Cl=a1*(a-a2);disp(Cl,"Cl=","Cl=a1(a-a2)","lift coefficient Cl:")

+Cd=cd+Cl^2/(%pi*e*AR);disp(Cd,"Cd=","Cd=cd+Cl^2/(%pi*e*AR)","total drag coefficient Cd=:")

+printf("\Answer:\n")

+printf("\n\Lift coefficient: %f \n\n",Cl)

+printf("\n\Totl drag coefficient: %f \n\n",Cd)
\ No newline at end of file
diff --git a/548/CH5/EX5.18/5_18.sce b/548/CH5/EX5.18/5_18.sce
new file mode 100755
index 000000000..e8f065a60
--- /dev/null
+++ b/548/CH5/EX5.18/5_18.sce
@@ -0,0 +1,6 @@
+pathname=get_absolute_file_path('5_18.sce')

+filename=pathname+filesep()+'5_18data.sci'

+exec(filename)

+Di=q*S*Cdi;disp(Di,"Di=","Di=q*S*Cdi","induced drag on one wing Di:")

+printf("\Answer:\n")

+printf("\n\Induced drag exerted on both the wings: %f N\n\n",2*Di)

diff --git a/548/CH5/EX5.19/5_19.sce b/548/CH5/EX5.19/5_19.sce
new file mode 100755
index 000000000..e07be4064
--- /dev/null
+++ b/548/CH5/EX5.19/5_19.sce
@@ -0,0 +1,4 @@
+pathname=get_absolute_file_path('5_19.sce')

+filename=pathname+filesep()+'5_19data.sci'

+exec(filename)

+disp("comparing the results of part a and b we can see the high-aspect ratio wing experiences a 26% higher increase in Cl than the low-aspect ratio wing.")
\ No newline at end of file
diff --git a/548/CH5/EX5.20/5_20.sce b/548/CH5/EX5.20/5_20.sce
new file mode 100755
index 000000000..5707adcfd
--- /dev/null
+++ b/548/CH5/EX5.20/5_20.sce
@@ -0,0 +1,9 @@
+pathname=get_absolute_file_path('5_20.sce')

+filename=pathname+filesep()+'5_20data.sci'

+exec(filename)

+V1=sqrt(2*Wt/(D*S*Clm));disp(V1,"V1=","V1=sqrt(2*Wt/(D*S*Clm))","stalling speed for full fuel tank V1:")

+V2=sqrt(2*Wf/(D*S*Clm));disp(V2,"V2=","V2=sqrt(2*Wf/(D*S*Clm))","stalling speed for empty fuel tank V1:")

+printf("\Answer:\n")

+printf("\n\stalling speed for full fuel tank : %f m/s\n\n",V1)

+printf("\n\stalling speed for empty fuel tank : %f m/s\n\n",V2)

+

diff --git a/548/CH5/EX5.21/5_21.sce b/548/CH5/EX5.21/5_21.sce
new file mode 100755
index 000000000..c68868168
--- /dev/null
+++ b/548/CH5/EX5.21/5_21.sce
@@ -0,0 +1,7 @@
+pathname=get_absolute_file_path('5_21.sce')

+filename=pathname+filesep()+'5_21data.sci'

+exec(filename)

+V=sqrt(2*Wt/(D*S*Clm));disp(V,"V=","V=sqrt(2*Wt/(D*S*Clm))","stalling speed for Boeing 727 V:")

+printf("\Answer:\n")

+printf("\n\stalling speed for full fuel tank : %f m/s\n\n",V)

+disp("stalling speed for lockhead F-104 is a much higher value than the Boeing 727.","comparison with stalling speed for full fuel tank of example 5.20:")
\ No newline at end of file