summaryrefslogtreecommitdiff
path: root/49/CH6
diff options
context:
space:
mode:
Diffstat (limited to '49/CH6')
-rwxr-xr-x49/CH6/EX6.1/ex1.sce16
-rwxr-xr-x49/CH6/EX6.10/ex10.sce11
-rwxr-xr-x49/CH6/EX6.2/ex2.sce21
-rwxr-xr-x49/CH6/EX6.3/ex3.sce31
-rwxr-xr-x49/CH6/EX6.4/ex4.sce40
-rwxr-xr-x49/CH6/EX6.5/ex5.sce27
-rwxr-xr-x49/CH6/EX6.6/ex6.sce28
-rwxr-xr-x49/CH6/EX6.7/ex7.sce15
-rwxr-xr-x49/CH6/EX6.8/ex8.sce14
-rwxr-xr-x49/CH6/EX6.9/ex9.sce10
10 files changed, 213 insertions, 0 deletions
diff --git a/49/CH6/EX6.1/ex1.sce b/49/CH6/EX6.1/ex1.sce
new file mode 100755
index 000000000..08e15fd36
--- /dev/null
+++ b/49/CH6/EX6.1/ex1.sce
@@ -0,0 +1,16 @@
+//CHAPTER 6 _ PRESSURE AND SOUND MEASUREMENT
+//Caption : MANOMETERS
+// Example 1 // Page 329
+D1=0.1 //('Enter the diameter of well =:')
+D2=0.01 //('Enter the diameter of the tube =:')
+g=9.81;
+pho_air=1.23 //('Enter the density of air in kg/m^3 =:')
+pho_liquid=1200 //('Enter the density of liquid in manometer =:')
+h=1 //('Enter the height by which liquid decreases in smaller area arm when exposed to the nominal pressure of p2 =:')
+// Let the pressure difference is represented by P=p1-p2
+disp("The pressure difference is given by:")
+disp("P=h*(1+((D2/D1)^2)*g*(pho_liquid-pho_air))")
+P=h*(1+((D2/D1)^2)*g*(pho_liquid-pho_air))*10^-3;
+printf('So the pressure difference is given by %1.2f kPa \n',P)
+
+
diff --git a/49/CH6/EX6.10/ex10.sce b/49/CH6/EX6.10/ex10.sce
new file mode 100755
index 000000000..66ed239ec
--- /dev/null
+++ b/49/CH6/EX6.10/ex10.sce
@@ -0,0 +1,11 @@
+//CHAPTER 6 _ PRESSURE AND SOUND MEASUREMENT
+//Caption : Sound Measurement
+// Example 10// Page 370
+Lp1=75 //('enter the sound level first machine=:')
+Lp2=77 //('enter the sound level second machine=:')
+Lp3=79 //('enter the sound level third machine=:')
+disp("Since the noise levels are incoherent,the total sound pressure is the sum of the mean square value of the individual sound pressures")
+disp("Lp_total=10*log10(10^(Lp1/10)+10^(Lp2/10)+10^(Lp3/10))")
+Lp_total=10*log10(10^(Lp1/10)+10^(Lp2/10)+10^(Lp3/10));
+printf('The total sound pressure is %ddB',Lp_total)
+//decibles are normally rounded off to the nearest integers \ No newline at end of file
diff --git a/49/CH6/EX6.2/ex2.sce b/49/CH6/EX6.2/ex2.sce
new file mode 100755
index 000000000..d99616966
--- /dev/null
+++ b/49/CH6/EX6.2/ex2.sce
@@ -0,0 +1,21 @@
+//CHAPTER 6 _ PRESSURE AND SOUND MEASUREMENT
+//Caption : MANOMETERS
+// Example 2 // Page 329
+pho_l=900
+disp("pho_l=900 ") //('Enter the density of the fluid =:')
+Pa= 500000
+disp("Pa= 500000 ") //('Enter the air pressure =:')
+t=298
+disp("t=298 ") //('Air is at what temperature(in deg cent) =:')
+R=287;
+disp("R=287;")
+g=9.81;
+T=t+273;
+disp("pho_a=Pa/(R*T);")
+pho_a=Pa/(R*T);
+printf('The density of air is %fd kg/m^3 \n',pho_a)
+h=.2 //('Enter the difference in the height of the fluid in the manometer=:')
+disp("Pres_diff=(g*h)*(pho_l-pho_a)")
+Pres_diff=(g*h)*(pho_l-pho_a)*10^-3
+printf('The differential pressure is %1.2f kPa\n',Pres_diff)
+
diff --git a/49/CH6/EX6.3/ex3.sce b/49/CH6/EX6.3/ex3.sce
new file mode 100755
index 000000000..1da1ada8b
--- /dev/null
+++ b/49/CH6/EX6.3/ex3.sce
@@ -0,0 +1,31 @@
+//CHAPTER 6 _ PRESSURE AND SOUND MEASUREMENT
+//Caption : Elastic Transducers
+// Example 3 // Page 337
+Sa=1000
+disp("Sa=1000") //('Enter the sensitivity of LVDT =:')
+//Properties of diaphragm
+E=200*10^9 //('Enter the value of modulus of elasticity=:')
+disp("E=200*10^9 ")
+v=0.3 //('Enter the Poissons ratio=:')
+disp("v=0.3 ")
+d=0.2 //('Enter the diameter of diaphragm=:')
+disp("d=0.2 ")
+R=d*(1/2);
+P_max=2*10^6 //('What is the maximum pressure?')
+disp("P_max=2*10^6 ")
+p=7800 //('What is the density of steel?')
+disp("Thickness is given by:")
+disp("t=(3*P_max*R^4*(1-v^4)/(4*E))^(1/4);")
+t=(3*P_max*R^4*(1-v^4)/(4*E))^(1/4)
+T=t*1000;
+printf('Thickness is %1.1f mm\n',T)
+//To calculate the lowest pressure in kPa which may be sensed by this instrument , resolution and the natural frequency of the diaphragm
+y=.001 //('Enter the l)east value of measurement=:')
+p_min=(y*16*E*t^3)/(3*R^4*(1-v^2)*Sa)
+printf('So the minimum pressure and resolution is %d Pa \n',p_min)
+f=(10.21/R^2)*((E*t^2)/(12*(1-v^2)*p))^(1/2)
+printf('The natural frequency of diaphragm is %fd Hz',f)
+
+
+
+
diff --git a/49/CH6/EX6.4/ex4.sce b/49/CH6/EX6.4/ex4.sce
new file mode 100755
index 000000000..6e4b6f240
--- /dev/null
+++ b/49/CH6/EX6.4/ex4.sce
@@ -0,0 +1,40 @@
+//CHAPTER 6 _ PRESSURE AND SOUND MEASUREMENT
+//Caption : Design of Pressure Transducers
+// Example 4 // Page 338
+p_max=10*10^6 //('Enter the capacity of the transducer=:')
+D=.05 //('Enter the diameter of diaphragm=:')
+R=D/2;
+v=0.3; // poissons ratio
+E=200*10^9;
+// We know that
+// y=3pR^4(1-v^2)/16t^3E
+// if y<t/4, the non linearity is restricted to 0.3%
+//So t is given by
+t=(3*p_max*R^4*(1-v^2)/(4*E))^(1/4)
+disp(t)
+printf('thickness comes out to be %fd m\n',t);
+Sr_max=(3*p_max*R^2)/(4*t^2)
+printf('So the max radial stress is %fd Pa\n',Sr_max)
+printf('The given fatigue strength is 500MPa\n' )
+if Sr_max > 500*10^6 then
+ disp("The diaphragm must be redesigned");
+ t1=((3*p_max*R^2)/(4*500*10^6))^(1/2);
+printf('The required thickness is %fd m\n',t1)
+
+else
+ disp("The design is OK");
+end
+// Let the voltage ratio be represented by Err
+Err=(820*p_max*R^2*(1-v^2))/(E*(t1^2))
+printf('The voltage ratio is %fd\n', Err)
+// For maximum power dissipation
+PT=1
+RT=120
+Ei=2*(PT*RT)^(1/2);
+disp("Let the sensitivity of the transducer be represented by ss")
+ss=(820*R^2*(1-v^2)*Ei)/(E*t1^2)
+printf('sensitivity is %fd\n', ss)
+// Part c
+S_LVDT=(ss*16*t^3*E)/(3*R^4*(1-v^2)*Ei)
+printf('SENSITIVITY OF LVDT IS %fd \n',S_LVDT)
+
diff --git a/49/CH6/EX6.5/ex5.sce b/49/CH6/EX6.5/ex5.sce
new file mode 100755
index 000000000..4d0ab154e
--- /dev/null
+++ b/49/CH6/EX6.5/ex5.sce
@@ -0,0 +1,27 @@
+//CHAPTER 6 _ PRESSURE AND SOUND MEASUREMENT
+//Caption : Pressure Gage
+// Example 5 // Page 347
+p_max=10*10^6 //('Enter the maximum differential pressure')
+fn=20000 //(' Enter the frequency')
+E=200*10^9; // modulus of elasticity
+v=0.3; // poissons ratio
+p=7800 // density of steel
+disp("Let t/R be represented by TR ")
+TR=((3*p_max*(1-v^2))/(4*E))^(1/4)
+// we know R^2/t = r2t=10.21(Et^2/12(1-v^2)p)^0.5/R^2 using it , we have
+r2t=(10.21*sqrt(E/(12*(1-v^2)*p)))/fn
+R=TR*r2t;
+printf('value of R is %fd m\n', R)
+
+t=R*TR;
+printf(' value of t is %fd m \n',t)
+
+eo=8.85*10^-12
+er=1.0006;
+d=.001 //('Enter the distance between the plates of capacitor=:')
+S=-(eo*er*%pi*R^2)/d^2;
+// variation of capacitor distance with respect to pressure is given by
+q=(3*R^4*(1-v^2))/(16*E*t^3)
+// total sensitivity of the pressure transducer is given by
+sensitivity=S*q*10^18;
+printf(' So the total sensitivity of the pressure transducer is given by %1.2f pF/MPa\n',sensitivity) \ No newline at end of file
diff --git a/49/CH6/EX6.6/ex6.sce b/49/CH6/EX6.6/ex6.sce
new file mode 100755
index 000000000..2c601e1d9
--- /dev/null
+++ b/49/CH6/EX6.6/ex6.sce
@@ -0,0 +1,28 @@
+//CHAPTER 6 _ PRESSURE AND SOUND MEASUREMENT
+//Caption : High Pressure Measurement
+// Example 6 // Page 357
+R1=100 //(' Enter the resistance of Mangnin wire=:')
+disp("R1=100")
+b=25*10^-12; // standard for mangnin
+disp("b=25*10^-12;")
+disp("u=0.5")
+u=0.5 //(' enter the uncertainty in measuring pressure for gage=:')
+// to calculate maximum uncertainty in differential pressure
+udp=u*(10-0.1)*10^6/100;
+uR=R1*b*udp;
+printf('So the maximum uncertainty in measuring resistance is %fd ohm \n',uR)
+//to calculate the output bridge voltage for 10 MPa
+Ei=5 //('enter the input voltage=:')
+disp("p1=0.1*10^6")
+disp("R2=R1*(1+b*p1)")
+disp("p2=10*10^6 ")
+p1=0.1*10^6 //('enter the pressure at which bridge is assumed to be balanced=:')
+R2=R1*(1+b*p1)
+p2=10*10^6 //('enter the pressure at which output voltage is to be calculated=:')
+R3=R1*(1+b*p2);
+dR=R3-R2;
+r=1;
+Eo=(r*dR*Ei)/((1+r)^2*R2)
+printf(' The output bridge voltage is %fd volt\n',Eo)
+
+
diff --git a/49/CH6/EX6.7/ex7.sce b/49/CH6/EX6.7/ex7.sce
new file mode 100755
index 000000000..3e32f4d67
--- /dev/null
+++ b/49/CH6/EX6.7/ex7.sce
@@ -0,0 +1,15 @@
+//CHAPTER 6 _ PRESSURE AND SOUND MEASUREMENT
+//Caption : McLeod Gage
+// Example 7 // Page 362
+disp("Vb=150*10^-6")
+disp("d=1.5*10^-3")
+disp("a=%pi*d^2/4;")
+Vb=150*10^-6 //('enter the volume of the Mc Leod gage=:')
+d=1.5*10^-3 //('enter the diameter of capillary=:')
+a=%pi*d^2/4;
+p=40*10^-6 //('enter the pressure for which the gage reading is to be noted=:')
+//y=(-p*area_cap+sqrt((p*area_cap)^2-4*p*area_cap*Vb))/(2*area_cap);
+l=p*a;
+
+y=(sqrt(l^2+(4*l*Vb))-l)/(2*a)
+printf('The gage reading comes out to be %fd mof Hg\n',y)
diff --git a/49/CH6/EX6.8/ex8.sce b/49/CH6/EX6.8/ex8.sce
new file mode 100755
index 000000000..2d7eabefd
--- /dev/null
+++ b/49/CH6/EX6.8/ex8.sce
@@ -0,0 +1,14 @@
+//CHAPTER 6 _ PRESSURE AND SOUND MEASUREMENT
+//Caption : Knudsen Gage
+// Example 8 // Page 363
+disp("Td=40")
+disp("Tv=300")
+disp("p=2*10^-6")
+Td=40 //('enter the temperature difference=:')
+Tv=300 //('enter the gas temperature at which the force has to be calculated=:')
+p=2*10^-6 //('enter the pressure(in m of Hg)=:')
+pa=p*13600*9.81;
+k=4*10^-4; // knudsen constant
+F=(pa*Td)/(k*Tv);
+printf('So the required force is %1.1f N',F)
+
diff --git a/49/CH6/EX6.9/ex9.sce b/49/CH6/EX6.9/ex9.sce
new file mode 100755
index 000000000..1a846c6de
--- /dev/null
+++ b/49/CH6/EX6.9/ex9.sce
@@ -0,0 +1,10 @@
+//CHAPTER 6 _ PRESSURE AND SOUND MEASUREMENT
+//Caption : Sound Measurement
+// Example 9// Page 369
+disp("Lp=104")
+Lp=104 //('enter the sound pressure level in decibles=:')
+disp("pa=20*10^-6;")
+disp("p=sqrt(10^(Lp/10)*pa^2);")
+pa=20*10^-6; // rms pressure threshold of hearing
+p=sqrt(10^(Lp/10)*pa^2);
+printf('root mean square sound pressure is %1.3fPa\n',p)
generated by cgit (git 2.34.1) at 2024-12-20 14:21:15 +0000