diff options
12 files changed, 7848 insertions, 0 deletions
diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap10_Curves.ipynb b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap10_Curves.ipynb new file mode 100644 index 00000000..b87fe16a --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap10_Curves.ipynb @@ -0,0 +1,1023 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "# Chapter 10: Curves" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 379 pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('Tangent length =', 58.45305445925609)\n", + "('Length long of cord=', 114.35142994976763)\n", + "('Length of curve =', 115.19173063162576)\n", + "('chainage of commencement =', 1262.046945540744)\n", + "('chainage of tangency =', 1377.2386761723697)\n", + "('apex distance =', 6.143663587883047)\n", + "('versed sine of curve is', 6.009409798203436)\n" + ] + } + ], + "source": [ + "#ch-10 page 379 pb-1\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "r=275;\n", + "t=24;\n", + "l=1320.5;\n", + "\n", + "tl=r*math.tan((t/2)*(math.pi/180));\n", + "print('Tangent length =',tl);\n", + "llc=2*r*math.sin((t/2)*(math.pi/180));\n", + "print('Length long of cord=',llc);\n", + "loc=(math.pi*r*t/180);\n", + "print('Length of curve =',loc)\n", + "coc=l-tl;\n", + "ct=coc+loc;\n", + "print('chainage of commencement =',coc);\n", + "print('chainage of tangency =',ct);\n", + "k=math.cos((t/2)*math.pi/180);\n", + "ad=r*((1/k)-1);\n", + "print('apex distance =',ad)\n", + "k1=math.cos((t/2)*(math.pi/180))\n", + "vsc=r*(1-k1);\n", + "print('versed sine of curve is',vsc);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 379,380 pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('radius of curve ', 573.0)\n", + "('Tangent length =', 153.5348872630333)\n", + "('Length of curve =', 300.02209841782525)\n", + "('Length long of cord=', 296.60662568748876)\n", + "('chainage of commencement =', 2606.4651127369666)\n", + "('chainage of tangency =', 2906.487211154792)\n", + "('length of each half =', 148.30331284374438)\n", + "('O30=', 18.738622298863106, 'O60=', 16.374481794326243, 'O90=', 12.412299602376265, 'O120=', 6.81817453294525, 'O148.3=', 0.0)\n" + ] + } + ], + "source": [ + "#ch-10 page 379,380 pb-2\n", + "from __future__ import division\n", + "\n", + "import math\n", + "ac=45.5;cb=75.5;\n", + "#a\n", + "\n", + "t=cb-ac;\n", + "l1=1719;\n", + "l=2760;\n", + "\n", + "#b\n", + "r=l1/3;\n", + "print('radius of curve ',r);\n", + "\n", + "#c\n", + "tl=r*math.tan((t/2)*(math.pi/180));\n", + "print('Tangent length =',tl);\n", + "#d\n", + "loc=(math.pi*r*t/180);\n", + "print('Length of curve =',loc)\n", + "#e\n", + "llc=2*r*math.sin((t/2)*(math.pi/180));\n", + "print('Length long of cord=',llc);\n", + "\n", + "#f,g\n", + "coc=l-tl;\n", + "ct=coc+loc;\n", + "print('chainage of commencement =',coc);\n", + "print('chainage of tangency =',ct);\n", + "\n", + "#h\n", + "\n", + "half=0.5*llc;\n", + "print('length of each half =',half);\n", + "\n", + "ini=30;\n", + "\n", + "k=math.sqrt(r*r-(half*half));\n", + "o=r-k\n", + "k1=r-o;\n", + "O30=(math.sqrt(r*r-(ini*ini)))-k1;\n", + "O60=(math.sqrt(r*r-(2*ini*2*ini)))-k1;\n", + " \n", + "O90=(math.sqrt(r*r-(3*ini*3*ini)))-k1;\n", + "O120=(math.sqrt(r*r-(4*ini*4*ini)))-k1;\n", + "Oh=(math.sqrt(r*r-(half*half)))-k1;\n", + "\n", + "print('O30=',O30,'O60=',O60,'O90=',O90,'O120=',O120,'O148.3=',Oh);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### ch-10 page 381 pb-3" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('Tangent length =', 150.2288093231531)\n", + "('Length of curve =', 278.55454861829503)\n", + "('chainage of T1=', 360.00119067684693)\n", + "('chainage of T2=', 638.555739295142)\n", + "('chainage covered=', 630.0011906768469)\n", + "('Length of final sub cord=', 8.55454861829503)\n", + "('first ofset=', 1.5)\n", + "('second ofset=', 3.0)\n", + "('tenth ofset=', 0.5496946010193738)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "a=126.8;\n", + "t=180-a;\n", + "r=300;\n", + "#b\n", + "tl=r*math.tan((t/2)*(math.pi/180));\n", + "print('Tangent length =',tl);\n", + "\n", + "#c\n", + "loc=(math.pi*r*t/180);\n", + "print('Length of curve =',loc)\n", + "\n", + "#d\n", + "l=510.23;\n", + "ct1=l-tl;\n", + "ct2=ct1+loc;\n", + "\n", + "print('chainage of T1=',ct1);\n", + "print('chainage of T2=',ct2);\n", + "\n", + "#f\n", + "n=9;\n", + "b=30;\n", + "cc=ct1+270;\n", + "lfsc=ct2-cc;\n", + "print('chainage covered=',cc);\n", + "print('Length of final sub cord=',lfsc);\n", + "\n", + "O1=(b*b)/(2*r);\n", + "O2=(b*b)/r;\n", + "\n", + "O10=(lfsc*(b+lfsc))/(2*r);\n", + "\n", + "print('first ofset=',O1);\n", + "print('second ofset=',O2);\n", + "print('tenth ofset=',O10);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 382 pb-4" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "(165.95962740164575, 158.58327092428829)\n", + "('Radius R=', 143.72525333242524)\n", + "('Tangent length BT1=', 82.97981370082289)\n", + "('Tangent length CT1=', 67.02018629917711)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "ab=30;bc=90;cd=140;\n", + "l1=250;l2=150;l3=325;\n", + "\n", + "abc=210-bc;\n", + "t1=0.5*abc;\n", + "bcd=270-cd;\n", + "t2=0.5*bcd;\n", + "t3=180-(t1+t2);\n", + "\n", + "\n", + "k=(math.sin(t2*(math.pi/180)))/(math.sin(t3*(math.pi/180)));\n", + "OB=l2*k;\n", + "k1=(math.sin(t1*(math.pi/180)))/(math.sin(t3*(math.pi/180)));\n", + "OC=l2*k1;\n", + "print(OB,OC);\n", + "R=OB*(math.sin(t1*(math.pi/180)));\n", + "print('Radius R=',R);\n", + "\n", + "BT1=OB*(math.cos(t1*(math.pi/180)));\n", + "CT1=OC*(math.cos(t2*(math.pi/180)));\n", + "\n", + "print('Tangent length BT1=',BT1);\n", + "print('Tangent length CT1=',CT1);\n", + "\n", + "\n", + " \n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 383 pb-5" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('R1=', 368.61561236693893)\n", + "('length of arc T1T2=', 209.43951023931953)\n", + "('length of arc T2T3=', 386.0133666034928)\n", + "('chainage of T1=', 792.8203230275509)\n", + "('chainage of T3=', 1388.2731998703632)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "r=400;\n", + "t1=15;t2=30;t3=60;\n", + "ct=900;\n", + "l=320;\n", + "BT2=r*(math.tan((t1)*math.pi/180));\n", + "CT2=l-BT2;\n", + "\n", + "r1=(CT2)/(math.tan((t2)*math.pi/180));\n", + "\n", + "print('R1=',r1);\n", + "t1t2=(math.pi*r*t2)/(180);\n", + "\n", + "t2t3=(math.pi*r1*t3)/(180);\n", + "\n", + "print('length of arc T1T2=',t1t2);\n", + "print('length of arc T2T3=',t2t3);\n", + "\n", + "\n", + "ct1=ct-BT2;\n", + "ct3=ct1+t1t2+t2t3;\n", + "\n", + "print('chainage of T1=',ct1);\n", + "print('chainage of T3=',ct3);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 384 pb-6" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('R2=', 1092.8203230275515)\n", + "('length of arc T1T2=', 209.43951023931953)\n", + "('length of arc T2T3=', 572.1993830861634)\n", + "('chainage of point of reverse curvature =', 1709.4395102393196)\n", + "('chainage of finishing point T3=', 2281.638893325483)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "r1=400;\n", + "t1=30;d=200;\n", + "ct1=1500;\n", + "k=1-(math.cos(t1*(math.pi/180)))\n", + "T1G=r1*(k);\n", + "\n", + "r2=(d-T1G)/k;\n", + "print('R2=',r2);\n", + "\n", + "t1t2=(math.pi*r1*t1)/180;\n", + "t2t3=(math.pi*r2*t1)/180;\n", + "print('length of arc T1T2=',t1t2);\n", + "print('length of arc T2T3=',t2t3);\n", + "\n", + "ct2=ct1+t1t2;\n", + "ct3=ct2+t2t3;\n", + "\n", + "print('chainage of point of reverse curvature =',ct2);\n", + "print('chainage of finishing point T3=',ct3);\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 385 pb-7" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('EF=', 228.7451827250347)\n", + "('chainage of T1=', 701.0877129159065)\n", + "('chainage of D=', 1015.2469782748858)\n", + "('chainage of T2', 1137.4200259144889)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "a1=135;a2=145;\n", + "t1=180-a1;\n", + "t2=180-a2;\n", + "t3=180-(t1+t2);\n", + "r1=400;r2=200;\n", + "ct=1000;\n", + "\n", + "ED=r1*(math.tan((t1/2)*(math.pi/180)));\n", + "\n", + "FD=r2*(math.tan((t2/2)*(math.pi/180)));\n", + "\n", + "EF=ED+FD;\n", + "\n", + "print('EF=',EF);\n", + "\n", + "BE=EF*(math.sin(t2*(math.pi/180)))/(math.sin(t3*(math.pi/180)));\n", + "\n", + "BF=EF*(math.sin(t1*(math.pi/180)))/(math.sin(t3*(math.pi/180)))\n", + "\n", + "\n", + "ct1=ct-(BE+ED);\n", + "\n", + "cd=ct1+((math.pi*r1*t1)/(180));\n", + "\n", + "ct2=cd+((math.pi*r2*t2)/(180));\n", + "\n", + "print('chainage of T1=',ct1);\n", + "print('chainage of D=',cd);\n", + "print('chainage of T2',ct2);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 386 pb-8" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('Radius R=', 272.7765415715475)\n", + "('angle Theta=', 109.0)\n", + "('curve length T1D=', 145.2058875651141)\n", + "('curve length DT2=', 192.814375291381)\n", + "('chainage of T1=', 1305.430383812096)\n", + "('chainage of T2=', 1643.4506466685912)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "t1=30.5;\n", + "t2=40.5;\n", + "EF=175;\n", + "cb=1500;\n", + "\n", + "k1=math.tan((t1/2)*(math.pi/180));\n", + "k2=math.tan((t2/2)*(math.pi/180));\n", + "\n", + "r=EF/(k1+k2);\n", + "print('Radius R=',r);\n", + "\n", + "et1=r*k1;\n", + "ft2=r*k2;\n", + "\n", + "t3=180-(t1+t2);\n", + "print('angle Theta=',t3);\n", + "k3=(math.sin(t2*(math.pi/180)))/(math.sin(t3*(math.pi/180)));\n", + "k4=(math.sin(t1*(math.pi/180)))/(math.sin(t3*(math.pi/180)));\n", + "\n", + "be=EF*k3;\n", + "bf=EF*k4;\n", + "\n", + "t1d=(math.pi*r*t1)/180;\n", + "dt2=(math.pi*r*t2)/180;\n", + "\n", + "print('curve length T1D=',t1d);\n", + "print('curve length DT2=',dt2);\n", + "\n", + "ct1=cb-(be+et1);\n", + "\n", + "ct2=ct1+t1d+dt2;\n", + "print('chainage of T1=',ct1);\n", + "print('chainage of T2=',ct2)\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 387 pb-9" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('theta 3=', 10.649036741314365)\n", + "('Radius R=', 829.124128893828)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "t1=80-70;\n", + "l=50;\n", + "k=1/(math.cos(20*(math.pi/180)));\n", + "\n", + "k1=k*(math.sin(t1*(math.pi/180)));\n", + "t3=math.asin(k1);\n", + "t3=t3*(180/(math.pi));\n", + "print('theta 3=',t3);\n", + "\n", + "t3=180-t3;\n", + "t2=180-(t3+t1);\n", + "\n", + "r=l*(math.sin(t1*(math.pi/180)))/(math.sin(0.6*(math.pi/180)));\n", + "print('Radius R=',r);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 388 pb-10" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('radius of circular curvature=', 402.713280728911)\n", + "('length of transistion curve =', 90.83333333333334)\n", + "('spiral angle=', 6.461627773592511)\n", + "('central angle=', 47.07674445281498)\n", + "('length of circular curve =', 330.88702808033025)\n", + "('shift of curve =', 0.8536568115234379)\n", + "('tangent length =', 278.4161466916694)\n", + "('chainage of 1st tangent point =', 871.5838533083306)\n", + "('chainage of 2nd tangent point =', 1384.1375480553274)\n", + "('chainage of 1st junction point =', 962.417186641664)\n", + "('chainage of 2nd junction point =', 1293.3042147219942)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "sp=80;\n", + "v=(sp*1000)/(60*60);\n", + "cr=(1/8);\n", + "g=9.81;\n", + "a=60;\n", + "\n", + "#a\n", + "\n", + "r=(v*v)/(g*cr);\n", + "print('radius of circular curvature=',r);\n", + "\n", + "#b\n", + "k=0.3;\n", + "l=(v*v*v)/(k*r);\n", + "print('length of transistion curve =',l);\n", + "\n", + "sa=l/(2*r);\n", + "sa=sa*(180/(math.pi));\n", + "print('spiral angle=',sa);\n", + "ca=a-(2*sa);\n", + "print('central angle=',ca);\n", + "\n", + "lcc=(math.pi*r*ca)/180;\n", + "print('length of circular curve =',lcc);\n", + "\n", + "s=(l*l)/(24*r);\n", + "print('shift of curve =',s);\n", + "ag=a/2;\n", + "t=(r+s)*(math.tan(ag*(math.pi/180)))+(l/2);\n", + "print('tangent length =',t);\n", + "#c\n", + "cip=1150;\n", + "c1t=cip-t;\n", + "c1j=c1t+l;\n", + "c2j=c1j+lcc;\n", + "c2t=c2j+l;\n", + "\n", + "print('chainage of 1st tangent point =',c1t);\n", + "print('chainage of 2nd tangent point =',c2t);\n", + "\n", + "print('chainage of 1st junction point =',c1j);\n", + "print('chainage of 2nd junction point =',c2j);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 389 pb-11" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('radius =', 343.8)\n", + "('tangent length =', 108.39972361659453)\n", + "('curve length =', 210.01546889247766)\n", + "('chainage of 1st point =', 1471.6002763834056)\n", + "('chainage of 2nd point =', 1681.6157452758832)\n", + "('length of final sub chord =', 8.399723616594429)\n", + "('chainage covered=', 1660)\n", + "('length of final sub chord', 21.615745275883228)\n", + "('deflection angle for initial sub chord =', 43.52844633883164, 'min')\n", + "('deflection angle for full chord', 2.5910642019719075, 'min')\n", + "('deflection angle for final sub chord', 1.8669261261094798, 'min')\n", + "('total deflection angle=', 17.5)\n", + "('apex distance =', 16.6843132234107)\n", + "('versed sine of curve =', 15.91211233275958)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "a=145;\n", + "cpi=1580;\n", + "de=5;\n", + "pi=30;\n", + "lct=0.00555;\n", + "\n", + "da=180-a;\n", + "\n", + "r=(1719)/5;\n", + "\n", + "print('radius =',r);\n", + "\n", + "#a\n", + "\n", + "tl=r*(math.tan((da/2)*(math.pi/180)));\n", + "print('tangent length =',tl);\n", + "\n", + "#b\n", + "\n", + "cl=(math.pi*r*da)/180;\n", + "print('curve length =',cl);\n", + "\n", + "#c\n", + "\n", + "c1t=cpi-tl;\n", + "print('chainage of 1st point =',c1t);\n", + "\n", + "#d\n", + "c2t=c1t+cl;\n", + "print('chainage of 2nd point =',c2t);\n", + "\n", + "#e\n", + "lisc=1480-c1t;\n", + "print('length of final sub chord =',lisc);\n", + "#f\n", + "n=6;\n", + "ini=30;\n", + "cc=1480+(n*30);\n", + "print('chainage covered=',cc);\n", + "#g\n", + "lfsc=c2t-cc;\n", + "print('length of final sub chord',lfsc);\n", + "#h\n", + "dasc=((c2t+100)*lisc)/(r);\n", + "print('deflection angle for initial sub chord =',dasc,'min');\n", + "#i\n", + "dafc=((c2t+100)*pi)/r;\n", + "print('deflection angle for full chord',dafc/60,'min');\n", + "#j\n", + "dafsc=((c2t+100)*lfsc)/r;\n", + "print('deflection angle for final sub chord',dafsc/60,'min');\n", + "\n", + "#k\n", + "\n", + "tda=da/2;\n", + "print('total deflection angle=',tda);\n", + "\n", + "\n", + "#l\n", + "k=1/(math.cos(tda*(math.pi/180)));\n", + "ad=r*(k-1);\n", + "print('apex distance =',ad);\n", + "\n", + "vs=r*(1-(math.cos(tda*(math.pi/180))));\n", + "print('versed sine of curve =',vs);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 391 pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('length of vertical curve =', 240.0)\n", + "('chainage of A', 430.0)\n", + "('chainage of C', 670.0)\n", + "('Rl of A', 374.9)\n", + "('Rl of C', 374.66)\n", + "('Rl of E', 374.78)\n", + "('Rl of F', 375.14)\n", + "('tangent correction at the apex =', 0.36000000000001364)\n", + "('tangent correction at 1st,2nd,3rd,4th,5th,6th, points', 0.01, 0.04, 0.09, 0.16, 0.25, 0.36)\n", + "RL of the points on grade\n", + "(375.0, 375.1, 375.20000000000005, 375.30000000000007, 375.4000000000001, 375.5000000000001)\n", + "RL of the points on curve\n", + "(374.99, 375.06, 375.11000000000007, 375.14000000000004, 375.1500000000001, 375.1400000000001)\n", + "Rls of points on the grade right side\n", + "(375.36, 375.22, 375.08000000000004, 374.94000000000005, 374.80000000000007)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "pi=20;\n", + "cb=550;\n", + "rlb=375.5;\n", + "g1=0.5;\n", + "g2=-0.7;\n", + "#a\n", + "vc=((g1-g2)*20)/0.1;\n", + "\n", + "print('length of vertical curve =',vc);\n", + "\n", + "#b,c\n", + "\n", + "ca=cb-(vc/2);\n", + "cc=ca+vc;\n", + "print('chainage of A',ca);\n", + "print('chainage of C',cc);\n", + "\n", + "#d,e,f,g\n", + "\n", + "rla=rlb-((g1*0.5*vc)/100);\n", + "rlc=rlb-((-g2*0.5*vc)/100);\n", + "rle=0.5*(rla+rlc);\n", + "rlf=0.5*(rlb+rle);\n", + "\n", + "print('Rl of A',rla);\n", + "print('Rl of C',rlc);\n", + "print('Rl of E',rle);\n", + "print('Rl of F',rlf);\n", + "#h\n", + "tc=rlb-rlf;\n", + "print('tangent correction at the apex =',tc);\n", + "\n", + "#i\n", + "tc1=((g1-g2)*(pi*pi))/(400*0.5*vc);\n", + "tc2=((g1-g2)*(2*pi*2*pi))/(400*0.5*vc);\n", + "tc3=((g1-g2)*(3*pi*3*pi))/(400*0.5*vc);\n", + "tc4=((g1-g2)*(4*pi*4*pi))/(400*0.5*vc);\n", + "tc5=((g1-g2)*(5*pi*5*pi))/(400*0.5*vc);\n", + "tc6=((g1-g2)*(6*pi*6*pi))/(400*0.5*vc);\n", + "print('tangent correction at 1st,2nd,3rd,4th,5th,6th, points',tc1,tc2,tc3,tc4,tc5,tc6);\n", + "\n", + "#j\n", + "rp=(g1*pi)/100;\n", + "\n", + "rl1=rla+rp;\n", + "rl2=rl1+rp;\n", + "rl3=rl2+rp;\n", + "rl4=rl3+rp;\n", + "rl5=rl4+rp;\n", + "rl6=rl5+rp;\n", + "print('RL of the points on grade');\n", + "print(rl1,rl2,rl3,rl4,rl5,rl6)\n", + "\n", + "#k\n", + "rlc1=rl1-tc1;\n", + "rlc2=rl2-(tc2);\n", + "rlc3=rl3-(tc3);\n", + "rlc4=rl4-(tc4);\n", + "rlc5=rl5-(tc5);\n", + "rlc6=rl6-(tc6);\n", + "\n", + "print('RL of the points on curve');\n", + "print(rlc1,rlc2,rlc3,rlc4,rlc5,rlc6);\n", + "\n", + "#l\n", + "\n", + "fp=0.14;\n", + "\n", + "rlg5=rlb-fp;\n", + "rlg4=rlg5-fp;\n", + "rlg3=rlg4-fp;\n", + "rlg2=rlg3-fp;\n", + "rlg1=rlg2-fp;\n", + "\n", + "print('Rls of points on the grade right side');\n", + "print(rlg5,rlg4,rlg3,rlg2,rlg1);\n", + "\n", + "\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-10 page 393 pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 13, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('length of vertical curve', 300.0)\n", + "('RL of A=', 252.0)\n", + "('RL of C=', 251.25)\n", + "('RL of E=', 251.625)\n", + "('RL of F=', 251.0625)\n", + "RL on the grade on the side AB \n", + "(251.7, 251.39999999999998, 251.09999999999997, 250.79999999999995)\n", + "RL on grade on side BC\n", + "(250.65, 250.8, 250.95000000000002, 251.10000000000002)\n", + "tangent correction from expression \n", + "(-0.0225, -0.09, -0.2025, -0.36)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "cb=400;\n", + "rlb=250.5;\n", + "pi=30;\n", + "g1=-1.0;\n", + "g2=0.5;\n", + "g=0.1;\n", + "ga=20;\n", + "#a\n", + "vc=(g1-g2)/g;\n", + "vc=-vc*ga;\n", + "print('length of vertical curve',vc);\n", + "\n", + "#b,c\n", + "ca=cb-(0.5*vc);\n", + "cc=ca+vc;\n", + "\n", + "#d,e,f,g\n", + "\n", + "rla=rlb+((0.5*vc)/100);\n", + "\n", + "rlc=rlb+((0.5*0.5*vc)/100);\n", + "\n", + "rle=0.5*(rla+rlc);\n", + "\n", + "rlf=0.5*(rle+rlb);\n", + "\n", + "print('RL of A=',rla);\n", + "print('RL of C=',rlc);\n", + "print('RL of E=',rle);\n", + "print('RL of F=',rlf);\n", + "\n", + "#h\n", + "fp=pi/100;\n", + "\n", + "rl1=rla-fp;\n", + "rl2=rl1-fp;\n", + "rl3=rl2-fp;\n", + "rl4=rl3-fp;\n", + "print('RL on the grade on the side AB ');\n", + "print(rl1,rl2,rl3,rl4);\n", + "\n", + "#i\n", + "\n", + "rp=(0.5*pi)/100;\n", + "\n", + "rls4=rlb+rp\n", + "rls3=rls4+rp\n", + "rls2=rls3+rp\n", + "rls1=rls2+rp\n", + "\n", + "print('RL on grade on side BC');\n", + "print(rls4,rls3,rls2,rls1);\n", + "\n", + "#j\n", + "\n", + "y1=((g1-g2)*(pi*pi))/(cb*0.5*vc);\n", + "y2=((g1-g2)*(2*pi*2*pi))/(cb*0.5*vc);\n", + "y3=((g1-g2)*(3*pi*3*pi))/(cb*0.5*vc);\n", + "y4=((g1-g2)*(4*pi*4*pi))/(cb*0.5*vc);\n", + "\n", + "print('tangent correction from expression ');\n", + "print(y1,y2,y3,y4);\n", + "\n" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 2", + "language": "python", + "name": "python2" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 2 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython2", + "version": "2.7.11" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap1_Introduction.ipynb b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap1_Introduction.ipynb new file mode 100644 index 00000000..a5f2c912 --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap1_Introduction.ipynb @@ -0,0 +1,1467 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "# Chapter 1: Introduction\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem1, pg 25" + ] + }, + { + "cell_type": "code", + "execution_count": 24, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "true length= 327.4905\n" + ] + } + ], + "source": [ + "\n", + "\n", + "l=20; #chain length\n", + "e=0.03; #error\n", + "l1=l+e; #L'\n", + "ml=327; #measured length\n", + "truel=(l1/l)*(ml) #true length\n", + "print (\"true length=\",truel)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem2, pg 25" + ] + }, + { + "cell_type": "code", + "execution_count": 25, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "amount of error= 0.20083682008368697\n" + ] + } + ], + "source": [ + "\n", + "\n", + "l1=20; #chain 1 length\n", + "e=0.05; #error\n", + "l11=l1+e; \n", + "ml1=1200; #measured lenght\n", + "tl=(l11/l1)*ml1; #true lenght of line\n", + "\n", + "l2=30; #chain 2 length\n", + "ml2=1195; #measured length\n", + "\n", + "l21=(tl/ml2)*l2; \n", + "ae=l21-l2; #amount of error\n", + "print('amount of error=',ae)\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem3, pg 25" + ] + }, + { + "cell_type": "code", + "execution_count": 26, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "true length1= 901.35\n", + "true length 2= 678.3750000000001\n", + "true distance= 1579.7250000000001\n" + ] + } + ], + "source": [ + "\n", + "\n", + "l1=20\n", + "e=(0.06/2) #consider mean elongation\n", + "l11=l1+e;\n", + "ml=900;\n", + "tl=(l11/l1)*ml;\n", + "print('true length1=',tl)\n", + "l2=20;\n", + "e2=(0.06+0.14)/2;\n", + "l21=20+e2;\n", + "ml2=1575-ml;\n", + "\n", + "tl2=(l21/l2)*ml2;\n", + "print('true length 2=',tl2)\n", + "td=tl+tl2;\n", + "print('true distance=',td)\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem4, pg26" + ] + }, + { + "cell_type": "code", + "execution_count": 27, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "distance between stations on map= 35.0 centimeters\n", + "true distance on ground = 1750.0 meters\n" + ] + } + ], + "source": [ + "\n", + "\n", + "s=100;\n", + "dsm=3500;\n", + "adsm=dsm/s;\n", + "\n", + "print('distance between stations on map=',adsm,'centimeters')\n", + "\n", + "actuals=50;\n", + "td=adsm*actuals;\n", + "\n", + "print('true distance on ground =',td,'meters')\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 5, pg 26" + ] + }, + { + "cell_type": "code", + "execution_count": 28, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "actual area present= 132.01840894148586 square cm\n", + "true area= 212286.9217619987 square meters\n" + ] + } + ], + "source": [ + "\n", + "\n", + "present=19.5\n", + "actual=20;\n", + "cm1=actual/present;\n", + "cm12=(actual*actual)/(present*present);\n", + "pm=125.5;\n", + "apm=pm*cm12;\n", + "print('actual area present=',apm,'square cm');\n", + "\n", + "cm=40;\n", + "cm2=cm*cm;\n", + "\n", + "area=cm2*apm;\n", + "scale=(20.05*20.05)/(20*20);\n", + "ta=scale*area;\n", + "print('true area=',ta,'square meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 6, pg 27" + ] + }, + { + "cell_type": "code", + "execution_count": 29, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + " for n=1\n", + "the temperature correction is 0.00396 meters\n", + "the pull corretion is 0.002380952380952381 meters\n", + "the sag correction is -0.0026680499999999995 meters\n", + "the total correction is 0.0036729023809523816 meters\n", + "the true length is 780.0954954619046\n", + " for n=2\n", + "the temperature correction is 0.00396 meters\n", + "the pull corretion is 0.002380952380952381 meters\n", + "the sag correction is -0.0006670124999999999 meters\n", + "the total correction is 0.005673939880952382 meters\n", + "the true length is 780.1475224369049\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "\n", + "L=30;\n", + "t0=20;\n", + "p0=10;\n", + "pm=15;\n", + "tm=32;\n", + "a=0.03;\n", + "al=11/(1000000);\n", + "E=2.1*(1000000);\n", + "w=0.693;\n", + "ml=780;\n", + "n=1;\n", + "print(' for n=1')\n", + "ct=al*L*(tm-t0);\n", + "print('the temperature correction is',ct,'meters');\n", + "\n", + "cp=(pm-p0)*L/(a*E);\n", + "print('the pull corretion is ',cp,' meters');\n", + "\n", + "cs=-L*w*w/(24*pm*pm*n*n);\n", + "print('the sag correction is ',cs,'meters');\n", + "\n", + "e=ct+cp+cs;\n", + "print('the total correction is ',e,'meters');\n", + "\n", + "l1=L+e;\n", + "\n", + "truelength=(l1/L)*ml;\n", + "print('the true length is ',truelength);\n", + "\n", + "n=2;\n", + "\n", + "print(' for n=2')\n", + "ct=al*L*(tm-t0);\n", + "print('the temperature correction is',ct,'meters');\n", + "\n", + "cp=(pm-p0)*L/(a*E);\n", + "print('the pull corretion is ',cp,' meters');\n", + "\n", + "cs=-L*w*w/(24*pm*pm*n*n);\n", + "print('the sag correction is ',cs,'meters');\n", + "\n", + "e=ct+cp+cs;\n", + "print('the total correction is ',e,'meters');\n", + "\n", + "l1=L+e;\n", + "\n", + "truelength=(l1/L)*ml;\n", + "print('the true length is ',truelength);\n", + "\n", + "\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 7, pg 28" + ] + }, + { + "cell_type": "code", + "execution_count": 30, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "the temperature correction is 0.0021999999999999997 meters\n", + "the pull corretion is -0.002380952380952381 meters\n", + "the sag correction is -0.0013333333333333335 meters\n", + "the total correction is -0.001514285714285715 meters\n", + "the horizontal distance is 19.998485714285714\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "L=20;\n", + "t0=20;\n", + "p0=15;\n", + "p=10;\n", + "tm=30;\n", + "a=0.02;\n", + "al=11/(1000000);\n", + "E=2.1*(1000000);\n", + "w=0.4;\n", + "\n", + "n=1;\n", + "ct=al*L*(tm-t0);\n", + "print('the temperature correction is',ct,'meters');\n", + "\n", + "cp=(p-p0)*L/(a*E);\n", + "print('the pull corretion is ',cp,' meters');\n", + "\n", + "cs=-L*w*w/(24*p*p*n*n);\n", + "print('the sag correction is ',cs,'meters');\n", + "\n", + "e=ct+cp+cs;\n", + "print('the total correction is ',e,'meters');\n", + "\n", + "hd=L+e;\n", + "\n", + "print('the horizontal distance is ',hd);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 8, pg 29" + ] + }, + { + "cell_type": "code", + "execution_count": 31, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "for p=5 case\n", + "the temperature correction is 0.00165 meters\n", + "the pull corretion is 0.0 meters\n", + "the sag correction is -0.02178 meters\n", + "the total correction is -0.020130000000000002 meters\n", + "the horizontal distance is 29.97987\n", + "for p=11 case\n", + "the temperature correction is 0.00165 meters\n", + "the pull corretion is 0.004285714285714286 meters\n", + "the sag correction is -0.0045000000000000005 meters\n", + "the total correction is 0.001435714285714285 meters\n", + "the horizontal distance is 30.001435714285716\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "L=30;\n", + "t0=20;\n", + "p0=5;\n", + "tm=25;\n", + "a=0.02;\n", + "al=11/(1000000);\n", + "E=2.1*(1000000);\n", + "float(E);\n", + "float(al);\n", + "w1=22;\n", + "w=0.66;\n", + "n=1;\n", + "\n", + "p=5;\n", + "print('for p=5 case');\n", + "\n", + "ct=al*L*(tm-t0);\n", + "float(ct);\n", + "print('the temperature correction is',ct,'meters');\n", + "\n", + "cp=(p-p0)*L/(a*E);\n", + "print('the pull corretion is ',cp,' meters');\n", + "\n", + "cs=-L*w*w/(24*p*p*n*n);\n", + "print('the sag correction is ',cs,'meters');\n", + "\n", + "e=ct+cp+cs;\n", + "print('the total correction is ',e,'meters');\n", + "\n", + "hd=L+e;\n", + "\n", + "print('the horizontal distance is ',hd);\n", + "\n", + "p=11;\n", + "print('for p=11 case');\n", + "\n", + "ct=al*L*(tm-t0);\n", + "print('the temperature correction is',ct,'meters');\n", + "\n", + "cp=(p-p0)*L/(a*E);\n", + "print('the pull corretion is ',cp,' meters');\n", + "\n", + "cs=-L*w*w/(24*p*p*n*n);\n", + "print('the sag correction is ',cs,'meters');\n", + "\n", + "e=ct+cp+cs;\n", + "print('the total correction is ',e,'meters');\n", + "\n", + "hd=L+e;\n", + "\n", + "print('the horizontal distance is ',hd);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 9, pg 30" + ] + }, + { + "cell_type": "code", + "execution_count": 32, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "the temperature correction is 0.00264 meters\n", + "the pull corretion is 0.003492063492063492 meters\n", + "the sag correction is -0.001171875 meters\n", + "the total correction is 0.004960188492063492 meters\n", + "the true length is 680.1686464087301\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "L=20;\n", + "t0=20;\n", + "p0=5;\n", + "pm=16;\n", + "tm=32;\n", + "a=0.03;\n", + "al=11/(1000000);\n", + "E=2.1*(1000000);\n", + "w=0.6;\n", + "ml=680;\n", + "n=1;\n", + "\n", + "\n", + "ct=al*L*(tm-t0);\n", + "print('the temperature correction is',ct,'meters');\n", + "\n", + "cp=(pm-p0)*L/(a*E);\n", + "print('the pull corretion is ',cp,' meters');\n", + "\n", + "cs=-L*w*w/(24*pm*pm*n*n);\n", + "print('the sag correction is ',cs,'meters');\n", + "\n", + "e=ct+cp+cs;\n", + "print('the total correction is ',e,'meters');\n", + "\n", + "l1=L+e;\n", + "\n", + "truelength=(l1/L)*ml;\n", + "print('the true length is ',truelength);" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 10, pg 31" + ] + }, + { + "cell_type": "code", + "execution_count": 33, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "the temperature correction is 0.0061600000000000005 meters\n", + "the pull corretion is -0.0033333333333333335 meters\n", + "the sag correction is -0.008979994074074075 meters\n", + "the total correction is -0.006153327407407408 meters\n", + "the correctt distance is 679.8505620486773\n" + ] + } + ], + "source": [ + "\n", + "\n", + "\n", + "L=28;\n", + "t0=20;\n", + "p0=10;\n", + "pm=5;\n", + "tm=40;\n", + "a=0.02;\n", + "al=11/(1000000);\n", + "E=2.1*(1000000);\n", + "w1=470;\n", + "ml=680;\n", + "n=1;\n", + "\n", + "w=(470*28)/30;\n", + "w=w/1000;\n", + "\n", + "ct=al*L*(tm-t0);\n", + "print('the temperature correction is',ct,'meters');\n", + "\n", + "cp=(pm-p0)*L/(a*E);\n", + "print('the pull corretion is ',cp,' meters');\n", + "\n", + "cs=-L*w*w/(24*pm*pm*n*n);\n", + "print('the sag correction is ',cs,'meters');\n", + "\n", + "e=ct+cp+cs;\n", + "print('the total correction is ',e,'meters');\n", + "\n", + "l1=L+e;\n", + "\n", + "dis=(l1/L)*ml;\n", + "print('the correctt distance is ',dis);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 11, pg32" + ] + }, + { + "cell_type": "code", + "execution_count": 34, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "from fig p.1.1\n", + "87\n", + " the value of EF is 135.34797293685585 meters\n", + " the value of DF is 103.68256255569626 meters\n", + " the value of EG is 205.85953773426738 meters\n" + ] + } + ], + "source": [ + "#ch-1, problems on obstacles in chaining, page-32,pb-1\n", + "\n", + "from __future__ import division\n", + "\n", + "import math;\n", + "\n", + "print('from fig p.1.1')\n", + "DE=87;\n", + "print(DE);\n", + "EF=float(87/(math.cos(50*(math.pi/180))))\n", + "\n", + "DF=87*(math.tan(50*(math.pi/180)))\n", + "\n", + "EG=87/(math.cos(65*(math.pi/180)))\n", + "\n", + "\n", + "print(' the value of EF is ',EF,'meters');\n", + "\n", + "print(' the value of DF is ',DF,'meters');\n", + "\n", + "print(' the value of EG is ',EG,'meters'); \n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 12, pg 33" + ] + }, + { + "cell_type": "code", + "execution_count": 35, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "width of river is 227.23577649516116 meters\n" + ] + } + ], + "source": [ + "#ch-1 page-33, pb-2\n", + "import math\n", + "\n", + "\n", + "x=(380.0285/2.5754);\n", + "\n", + "PA=x;\n", + "AQ=367-x;\n", + "al=180-(36.45+86.55);\n", + "bt=86.35-40-35;\n", + "\n", + "TA=AQ*math.tan(46*(math.pi/180));\n", + "\n", + "print('width of river is ',TA,'meters');\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### Problem 13, pg 34" + ] + }, + { + "cell_type": "code", + "execution_count": 36, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "width of river is 316.63370603933663 meters\n" + ] + } + ], + "source": [ + "# cha-1 page-34 pb-3\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "x=(849.224)/2.6196\n", + "\n", + "\n", + "\n", + "PA=x;\n", + "AQ=517-x;\n", + "al=78-33.67;\n", + "bt=180-(43.333+78);\n", + "\n", + "TA=AQ*math.tan(58.66*(math.pi/180));\n", + "\n", + "print('width of river is ',TA,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 14, pg35" + ] + }, + { + "cell_type": "code", + "execution_count": 37, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "chainage of C is 277.08203230275507 meters\n" + ] + } + ], + "source": [ + "# cha-1 page-34,35 pb-4\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "al=288.5-(48.5+180);\n", + "bt=90-48.5;\n", + "BAC=360-41.5;\n", + "\n", + "AC=40*(math.tan(60*(math.pi/180)));\n", + "\n", + "A=207.8;\n", + "\n", + "C=A+AC;\n", + "\n", + "print('chainage of C is',C,'meters');\n", + " \n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### Problem 15, pg36" + ] + }, + { + "cell_type": "code", + "execution_count": 38, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "width of the river is 74.99999999999999 meters\n" + ] + } + ], + "source": [ + "\n", + "\n", + "import math\n", + "BB=287.25;\n", + "MC=62.25;\n", + "al=(BB-180)-MC;\n", + "BM=75;\n", + "BC=BM*(math.tan(45*(math.pi/180)))\n", + "\n", + "print('width of the river is ',BC,'meters')\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### Problem 16, pg 36" + ] + }, + { + "cell_type": "code", + "execution_count": 39, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "AB= 241.8677324489565\n" + ] + } + ], + "source": [ + "#CH-1 PAGE-36 PB-6;\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "\n", + "AC=250;\n", + "AD=300;\n", + "DB=150;\n", + "BC=100;\n", + "DC=DB+BC;\n", + "\n", + "cosal=(AD*AD+DC*DC-(AC*AC))/(2*AD*DC);\n", + "\n", + "AB=math.sqrt((AD*AD+DB*DB)-2*(AD*DB*cosal));\n", + "\n", + "print('AB=',AB);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 17, pg37" + ] + }, + { + "cell_type": "code", + "execution_count": 40, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "al 63.43494882292201\n", + "bt= 26.56505117707799\n", + "k= 0.5\n", + "chinage of c is 375.5 meters\n" + ] + } + ], + "source": [ + "# ch-1 page-36,37 pb-7\n", + "\n", + "from __future__ import division\n", + "\n", + "\n", + "import math\n", + "\n", + "BE=50;\n", + "AB=25;\n", + "AEC=157.5-67.5;\n", + "\n", + "al=math.atan2(BE,AB);\n", + "al=al*(180/math.pi);\n", + "\n", + "print('al',al)\n", + "\n", + "bt=90-al;\n", + "print('bt=',bt);\n", + "k=(math.tan(bt*math.pi/180))\n", + "\n", + "print('k=',k)\n", + "BC=BE/k;\n", + "C=275.5+BC;\n", + "print('chinage of c is',C,'meters')\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 18, pg38" + ] + }, + { + "cell_type": "code", + "execution_count": 41, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "measured length is 79.71623152917896 meters\n", + "true length is 79.61658623976749 meters\n" + ] + } + ], + "source": [ + "#ch-1 page -37,38 pb-1\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "\n", + "a=17.5;\n", + "b=19.3;\n", + "c=17.8;\n", + "d=13.6;\n", + "e=12.9;\n", + "\n", + "da=2.35;\n", + "db=4.20;\n", + "dc=2.95;\n", + "dd=1.65;\n", + "de=3.25;\n", + "\n", + "AB=math.sqrt((a*a)-(da*da));\n", + "BC=math.sqrt((b*b)-(db*db));\n", + "CD=math.sqrt((c*c)-(dc*dc));\n", + "DE=math.sqrt((d*d)-(dd*dd));\n", + "EF=math.sqrt((e*e)-(de*de));\n", + "\n", + "total=AB+BC+CD+DE+EF;\n", + "print('measured length is ',total,'meters');\n", + "\n", + "e=0.025;\n", + "l=20;\n", + "l1=l-e;\n", + "ml=total;\n", + "\n", + "tl=(l1/l)*ml;\n", + "\n", + "print('true length is ',tl,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 19, pg 38" + ] + }, + { + "cell_type": "code", + "execution_count": 42, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "measured length is 531.2592044589876 meters\n", + "true length is 532.587352470135 meters\n" + ] + } + ], + "source": [ + "#ch-1 page -38 pb-2\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "ab=550;\n", + "AB=ab*(math.cos(15*(math.pi/180)));\n", + "\n", + "l=20;\n", + "e=0.05;\n", + "l1=l+e;\n", + "ml=AB;\n", + "print('measured length is ',ml,'meters');\n", + "\n", + "tl=(l1/l)*ml;\n", + "\n", + "print('true length is ',tl,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 20, pg39" + ] + }, + { + "cell_type": "code", + "execution_count": 43, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "horizontal distance 1 is 275.74617084341827 meters\n", + "horizontal distance 2 is 278.61041325879694 meters\n", + "horizontal distance 3 is 279.8856909525744 meters\n" + ] + } + ], + "source": [ + "#ch-1 page -38,39 pb-3\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "ab=280;\n", + "\n", + "AB1=ab*(math.cos(10*(math.pi/180)));\n", + "\n", + "print('horizontal distance 1 is ',AB1,'meters');\n", + "\n", + "cosal=(10/(math.sqrt(101)));\n", + "\n", + "AB2=ab*cosal;\n", + "\n", + "print('horizontal distance 2 is ',AB2,'meters');\n", + "\n", + "bb=8;\n", + "AB3=math.sqrt(ab*ab-(bb*bb));\n", + "\n", + "print('horizontal distance 3 is ',AB3,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 21, pg40" + ] + }, + { + "cell_type": "code", + "execution_count": 44, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "measured length is 101.35201880331583 meters\n", + "true horizontal distance is 101.26755878764641 meters\n" + ] + } + ], + "source": [ + "#ch-1 page -39,40 pb-4\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "a=28.7;\n", + "b=23.4;\n", + "c=20.9;\n", + "d=29.6;\n", + "\n", + "ag=5;\n", + "bg=7;\n", + "cg=10;\n", + "dg=12;\n", + "\n", + "AB=a*(math.cos(ag*(math.pi/180)));\n", + "\n", + "BC=b*(math.cos(bg*(math.pi/180)));\n", + "\n", + "CD=c*(math.cos(cg*(math.pi/180)));\n", + "\n", + "DE=d*(math.cos(dg*(math.pi/180)));\n", + "\n", + "total=AB+BC+CD+DE;\n", + "\n", + "ml=total;\n", + "\n", + "print('measured length is ',ml,'meters');\n", + "\n", + "l=30;\n", + "e=0.025;\n", + "l1=l-e;\n", + "\n", + "tl=(l1/l)*ml;\n", + "\n", + "print('true horizontal distance is ',tl,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 22, pg 40" + ] + }, + { + "cell_type": "code", + "execution_count": 45, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "theta1= 30.009552668941378\n", + "theta2= 106 degrees 32.534711618974654 minutes\n" + ] + } + ], + "source": [ + "#ch-1 page -40 pb-1\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "\n", + "a=23;\n", + "b=16.5;\n", + "c=12;\n", + "\n", + "\n", + "t1=math.acos((a*a+b*b-(c*c))/(2*a*b));\n", + "t1=t1*(180/math.pi);\n", + "\n", + "print('theta1=',t1);\n", + "\n", + "t2=math.acos((c*c+b*b-(a*a))/(2*c*b));\n", + "t2=t2*(180/math.pi);\n", + "dg=int(t2)\n", + "mi=t2-int(t2)\n", + "mi=(mi*60);\n", + "print('theta2=',dg,'degrees',mi,'minutes');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 23, pg 41" + ] + }, + { + "cell_type": "code", + "execution_count": 46, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "theta1= 5 degrees 46.94403663966165 minutes\n", + "theta2= 165 degrees 26.421472313304548 minutes\n" + ] + } + ], + "source": [ + "#ch-1 page -40,41 pb-2\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "\n", + "a=257;\n", + "b=156;\n", + "c=103;\n", + "\n", + "\n", + "t1=math.acos((a*a+b*b-(c*c))/(2*a*b));\n", + "t1=t1*(180/math.pi);\n", + "\n", + "dg1=int(t1)\n", + "mi1=t1-int(t1)\n", + "mi1=(mi1*60);\n", + "print('theta1=',dg1,'degrees',mi1,'minutes');\n", + "\n", + "\n", + "t2=math.acos((c*c+b*b-(a*a))/(2*c*b));\n", + "t2=t2*(180/math.pi);\n", + "dg=int(t2)\n", + "mi=t2-int(t2)\n", + "mi=(mi*60);\n", + "print('theta2=',dg,'degrees',mi,'minutes');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 24, pg 42" + ] + }, + { + "cell_type": "code", + "execution_count": 47, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "RF is 0.025\n", + "length of scale is 15.000000000000002 meters\n" + ] + } + ], + "source": [ + "#CH-1 PAGE-42 PB-1;\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "\n", + "sc=100;\n", + "a=2.5;\n", + "m=6;\n", + "\n", + "RF=(a/sc);\n", + "\n", + "print('RF is ',RF);\n", + "\n", + "length=RF*m*sc;\n", + "\n", + "print('length of scale is ',length,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 25, pg 42" + ] + }, + { + "cell_type": "code", + "execution_count": 48, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "RF= 0.0002\n", + "length of final scale is 700.0\n" + ] + } + ], + "source": [ + "#CH-1 PAGE-42,43 PB-2;\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "\n", + "sc=100;\n", + "area=93750;\n", + "l=6.0;\n", + "b=6.25;\n", + "\n", + "cm2=(area)/(l*b);\n", + "\n", + "cm=math.sqrt(cm2);\n", + "RF=1/(sc*cm);\n", + "\n", + "print('RF=',RF);\n", + "\n", + "leng=14;\n", + "leng=leng*cm;\n", + "\n", + "print('length of final scale is ',leng);\n", + "\n", + "\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### Problem 26, pg 43" + ] + }, + { + "cell_type": "code", + "execution_count": 49, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "RF= 0.00025\n", + "length of scale is 600.0 meters\n" + ] + } + ], + "source": [ + "#CH-1 PAGE-43 PB-3;\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "l=1.2;\n", + "al=30;\n", + "al=al/100;\n", + "sc=1000;\n", + "\n", + "\n", + "RF=(al)/(sc*l);\n", + "print('RF=',RF);\n", + "\n", + "\n", + "cm1=(1/RF)/(100);\n", + "\n", + "lsc=15;\n", + "cm15=lsc*cm1;\n", + "\n", + "print('length of scale is ',cm15,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### Problem 27, pg44" + ] + }, + { + "cell_type": "code", + "execution_count": 50, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "1cm= 30.0\n", + "RF= 0.03333333333333333\n", + "length of scale is 13.333333333333334 CENTIMETERS\n" + ] + } + ], + "source": [ + "#CH-1 PAGE-44 PB-4;\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "\n", + "sc=100;\n", + "hect=10000;\n", + "area=0.45*hect;\n", + "\n", + "cm1=(area)/5;\n", + "cm=math.sqrt(cm1);\n", + "\n", + "print('1cm=',cm);\n", + "RF=1/(cm);\n", + "print('RF=',RF);\n", + "\n", + "\n", + "maxl=400;\n", + "\n", + "los=(RF*maxl);\n", + "\n", + "print('length of scale is',los,'CENTIMETERS');\n" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 2", + "language": "python", + "name": "python2" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 2 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython2", + "version": "2.7.11" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap2_Chain-Surveying.ipynb b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap2_Chain-Surveying.ipynb new file mode 100644 index 00000000..e3e0c595 --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap2_Chain-Surveying.ipynb @@ -0,0 +1,127 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "# Chapter 2: Chain Surveying" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### pg-56, pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('max length of offset should be', 6.8842279474019135, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "ag=5;\n", + "giv=0.03;\n", + "\n", + "L=20;\n", + "l=(giv*L/(math.sin(ag*math.pi/180)));\n", + "\n", + "\n", + "AB=l;\n", + "\n", + "BC=AB*(math.sin(ag*(math.pi/180)));\n", + "BC=BC/20;\n", + "\n", + "print('max length of offset should be',l,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### ch-2 page-56, pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('required displacement perpendicular to chain is', 0.0020556978681392835, 'meters')\n", + "('displacement parallel ot chain is', 0.07850393436441575, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "\n", + "AD=AB=15;\n", + "ag=3;\n", + "AC=15*(math.cos(ag*(math.pi/180)))\n", + "\n", + "CD=AB-AC\n", + "sc=10;\n", + "\n", + "CD=CD/sc;\n", + "\n", + "print('required displacement perpendicular to chain is',CD,'meters');\n", + "\n", + "\n", + "BC=AB*(math.sin(ag*(math.pi/180)));\n", + "\n", + "BC=BC/sc;\n", + "print('displacement parallel ot chain is',BC,'meters');\n", + "\n" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 2", + "language": "python", + "name": "python2" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 2 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython2", + "version": "2.7.11" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap3_Compass-Traversing.ipynb b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap3_Compass-Traversing.ipynb new file mode 100644 index 00000000..af93c556 --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap3_Compass-Traversing.ipynb @@ -0,0 +1,1695 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "# Chapter 3: Compass Traversing" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3, section 3.10, pg85, problem 1" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "N 45 degrees 30.0 minutes E\n", + "S 54 degrees 15.0 minutes E\n", + "S 42 degrees 15.0 minutes W\n", + "N 39 degrees 30.0 minutes W\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "import math\n", + "\n", + "#(a)\n", + "WCB_AB=45+(30/60)\n", + "QB_AB=WCB_AB\n", + "mins=(QB_AB-int(QB_AB))*60\n", + "deg=int(QB_AB)\n", + "print \"N\",deg,\"degrees\",mins,\"minutes E\"\n", + " \n", + "#(b)\n", + "WCB_BC=125+(45/60)\n", + "QB_BC=180-WCB_BC\n", + "mins=(QB_BC-int(QB_BC))*60\n", + "deg=int(QB_BC)\n", + "print \"S\",deg,\"degrees\",mins,\"minutes E\"\n", + "\n", + "#(c)\n", + "WCB_CD=222+(15/60)\n", + "QB_CD=WCB_CD-180\n", + "deg=int(QB_CD)\n", + "mins=(QB_CD-deg)*60\n", + "print \"S\",deg,\"degrees\",mins,\"minutes W\"\n", + "\n", + "#(d)\n", + "WCB_DE=320+(30/60)\n", + "QB_DE=360-WCB_DE\n", + "deg=int(QB_DE)\n", + "mins=(QB_DE-deg)*60\n", + "print \"N\",deg,\"degrees\",mins,\"minutes W\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3,section 3.10,problem 2,pg 85" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "216 degrees 30.0 minutes\n", + "136 degrees 30.0 minutes\n", + "26 degrees 45.0 minutes\n", + "319 degrees 45.0 minutes\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "import math\n", + "\n", + "#(a)\n", + "QB_AB=36+(30/60)\n", + "WCB_AB=180+QB_AB\n", + "mins=(WCB_AB-int(WCB_AB))*60\n", + "deg=int(WCB_AB)\n", + "print deg,\"degrees\",mins,\"minutes\"\n", + " \n", + "#(b)\n", + "QB_BC=43+(30/60)\n", + "WCB_BC=180-QB_BC\n", + "mins=(WCB_BC-int(WCB_BC))*60\n", + "deg=int(WCB_BC)\n", + "print deg,\"degrees\",mins,\"minutes\"\n", + "\n", + "\n", + "#(c)\n", + "QB_CD=26+(45/60)\n", + "WCB_CD=QB_CD\n", + "mins=(WCB_CD-int(WCB_CD))*60\n", + "deg=int(WCB_CD)\n", + "print deg,\"degrees\",mins,\"minutes\"\n", + "\n", + "#(d)\n", + "QB_DE=40+(15/60)\n", + "WCB_DE=360-QB_DE\n", + "mins=(WCB_DE-int(WCB_DE))*60\n", + "deg=int(WCB_DE)\n", + "print deg,\"degrees\",mins,\"minutes\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3,section 3.11,problem 1,pg 85" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "130 degrees 30.0 minutes\n", + "325 degrees 15.0 minutes\n", + "30 degrees 30.0 minutes\n", + "240 degrees 45.0 minutes\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "import math\n", + "\n", + "#(a)\n", + "FB_AB=310+(30/60)\n", + "BB_AB=FB_AB-180\n", + "mins=(BB_AB-int(BB_AB))*60\n", + "deg=int(BB_AB)\n", + "print deg,\"degrees\",mins,\"minutes\"\n", + " \n", + "#(b)\n", + "FB_BC=145+(15/60)\n", + "BB_BC=FB_BC+180\n", + "mins=(BB_BC-int(BB_BC))*60\n", + "deg=int(BB_BC)\n", + "print deg,\"degrees\",mins,\"minutes\"\n", + "\n", + "#(c)\n", + "FB_CD=210+(30/60)\n", + "BB_CD=FB_CD-180\n", + "mins=(BB_CD-int(BB_CD))*60\n", + "deg=int(BB_CD)\n", + "print deg,\"degrees\",mins,\"minutes\"\n", + "#(d)\n", + "FB_DE=60+(45/60)\n", + "BB_DE=FB_DE+180\n", + "mins=(BB_DE-int(BB_DE))*60\n", + "deg=int(BB_DE)\n", + "print deg,\"degrees\",mins,\"minutes\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3,section 3.11,problem 2,pg 86" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "N 30 degrees 30.0 minutes W\n", + "S 40 degrees 30.0 minutes E\n", + "N 60 degrees 15.0 minutes E\n", + "S 45 degrees 30.0 minutes W\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "import math\n", + "\n", + "#(a)\n", + "FB_AB=30+(30/60)\n", + "BB_AB=FB_AB\n", + "mins=(BB_AB-int(BB_AB))*60\n", + "deg=int(BB_AB)\n", + "print \"N\",deg,\"degrees\",mins,\"minutes W\"\n", + "\n", + "#(b)\n", + "FB_AB=40+(30/60)\n", + "BB_AB=FB_AB\n", + "mins=(BB_AB-int(BB_AB))*60\n", + "deg=int(BB_AB)\n", + "print \"S\",deg,\"degrees\",mins,\"minutes E\"\n", + "\n", + "#(c)\n", + "FB_AB=60+(15/60)\n", + "BB_AB=FB_AB\n", + "mins=(BB_AB-int(BB_AB))*60\n", + "deg=int(BB_AB)\n", + "print \"N\",deg,\"degrees\",mins,\"minutes E\"\n", + "\n", + "#(d)\n", + "FB_AB=45+(30/60)\n", + "BB_AB=FB_AB\n", + "mins=(BB_AB-int(BB_AB))*60\n", + "deg=int(BB_AB)\n", + "print \"S\",deg,\"degrees\",mins,\"minutes W\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### chapter 3,section 3.11,problem 3,pg 86" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "220 degrees 30.0 minutes\n", + "130 degrees 45.0 minutes\n", + "325 degrees 45.0 minutes\n", + "35 degrees 30.0 minutes\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "#(a)\n", + "BB_AB=40+(30/60)\n", + "FB_AB=BB_AB+180\n", + "mins=(FB_AB-int(FB_AB))*60\n", + "deg=int(FB_AB)\n", + "print deg,\"degrees\",mins,\"minutes\"\n", + " \n", + "#(b)\n", + "BB_BC=310+(45/60)\n", + "FB_BC=BB_BC-180\n", + "mins=(FB_BC-int(FB_BC))*60\n", + "deg=int(FB_BC)\n", + "print deg,\"degrees\",mins,\"minutes\"\n", + " \n", + "\n", + "#(c)\n", + "BB_CD=145+(45/60)\n", + "FB_CD=BB_CD+180\n", + "mins=(FB_CD-int(FB_CD))*60\n", + "deg=int(FB_CD)\n", + "print deg,\"degrees\",mins,\"minutes\"\n", + " \n", + "\n", + "#(d)\n", + "BB_DE=215+(30/60)\n", + "FB_DE=BB_DE-180\n", + "mins=(FB_DE-int(FB_DE))*60\n", + "deg=int(FB_DE)\n", + "print deg,\"degrees\",mins,\"minutes\"\n", + " \n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3,section 3.11,problem 4,pg 86" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "S 30 degrees 30.0 minutes E\n", + "N 40 degrees 15.0 minutes W\n", + "S 60 degrees 45.0 minutes W\n", + "N 45 degrees 30.0 minutes E\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "#(a)\n", + "BB_AB=30+(30/60)\n", + "FB_AB=BB_AB\n", + "mins=(FB_AB-int(FB_AB))*60\n", + "deg=int(FB_AB)\n", + "print \"S\",deg,\"degrees\",mins,\"minutes E\"\n", + "\n", + "#(b)\n", + "BB_BC=40+(15/60)\n", + "FB_BC=BB_BC\n", + "mins=(FB_BC-int(FB_BC))*60\n", + "deg=int(FB_BC)\n", + "print \"N\",deg,\"degrees\",mins,\"minutes W\"\n", + "\n", + "#(c)\n", + "BB_CD=60+(45/60)\n", + "FB_CD=BB_CD\n", + "mins=(FB_CD-int(FB_CD))*60\n", + "deg=int(FB_CD)\n", + "print \"S\",deg,\"degrees\",mins,\"minutes W\"\n", + "\n", + "#(d)\n", + "BB_DE=45+(30/60)\n", + "FB_DE=BB_DE\n", + "mins=(FB_DE-int(FB_DE))*60\n", + "deg=int(FB_DE)\n", + "print \"N\",deg,\"degrees\",mins,\"minutes E\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3, section 3.12, pg87, problem 1" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "truebearing of AB= 130 degrees 15.0 minutes\n", + "magnetic bearing of AB= 219 degrees 0.0 minutes\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "magneticbearing=135+0.5\n", + "declination=5+0.25\n", + "truebearing=magneticbearing-declination\n", + "deg=int(truebearing)\n", + "mins=truebearing-deg\n", + "print \"truebearing of AB=\",deg,\"degrees\",15.0,\"minutes\"\n", + "\n", + "truebearing=210+(45/60)\n", + "declination=8+(15/60)\n", + "magnetic_bearing=truebearing+declination\n", + "deg=int(magnetic_bearing)\n", + "mins=magnetic_bearing-deg\n", + "print \"magnetic bearing of AB=\",deg,\"degrees\",mins,\"minutes\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3, section 3.12, pg87, problem 2" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Required true bearing= S 40 degrees 30.0 minutes W\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "RB_CD=30+(15/60)\n", + "WCB_CD=180+RB_CD\n", + "declination=10+(15/60)\n", + "TB=WCB_CD+declination\n", + "truebearing=TB-180\n", + "deg=int(truebearing)\n", + "mins=(truebearing-deg)*60\n", + "print \"Required true bearing=\",\"S\",deg,\"degrees\",mins,\"minutes\",\"W\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3, section 3.12, pg88, problem 3" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Magnetic bearing= 312 degrees 45.0 minutes\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "magneticbearing=320+(30/60)\n", + "declination=3+(30/60)\n", + "truebearing=magneticbearing-declination\n", + "declination2=4+(15/60)\n", + "MB=truebearing-declination2\n", + "deg=int(MB)\n", + "mins=(MB-deg)*60\n", + "print \"Magnetic bearing=\",deg,\"degrees\",mins,\"minutes\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3, section 3.12, pg88, problem 4" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "4 degrees 30.0 minutes E\n", + "5 degrees 45.0 minutes W\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "#(a)\n", + "magneticbearing=175+(30/60)\n", + "magneticdeclination=180-magneticbearing\n", + "deg=int(magneticdeclination)\n", + "mins=(magneticdeclination-deg)*60\n", + "print deg,\"degrees\",mins,\"minutes E\"\n", + "\n", + "#(b)\n", + "\n", + "magneticdeclination=5+(45/60)\n", + "deg=int(magneticdeclination)\n", + "mins=(magneticdeclination-deg)*60\n", + "print deg,\"degrees\",mins,\"minutes W\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3, section 3.13, pg88, problem 1" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "AngleAOB= 109 degrees 45.0 minutes\n", + "AngleBOC= 80 degrees 30.0 minutes\n", + "AngleCOD= 89 degrees 45.0 minutes\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "bearingOB=140+(15/60)\n", + "bearingOA=30+(30/60)\n", + "angleAOB=bearingOB-bearingOA\n", + "deg=int(angleAOB)\n", + "mins=(angleAOB-deg)*60\n", + "print \"AngleAOB=\",deg,\"degrees\",mins,\"minutes\"\n", + "\n", + "bearingOC=220+(45/60)\n", + "angleBOC=bearingOC-bearingOB\n", + "deg=int(angleBOC)\n", + "mins=(angleBOC-deg)*60\n", + "print \"AngleBOC=\",deg,\"degrees\",mins,\"minutes\"\n", + "\n", + "bearingOD=310+(30/60)\n", + "angleCOD=bearingOD-bearingOC\n", + "deg=int(angleCOD)\n", + "mins=(angleCOD-deg)*60\n", + "print \"AngleCOD=\",deg,\"degrees\",mins,\"minutes\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### chapter 3, section 3.13, pg89, problem 2" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Interior angle B= 105 degrees 15.0 minutes\n", + "Interior angle C= 99 degrees 45.0 minutes\n", + "Exterior angle D= 260 degrees 15.0 minutes\n", + "Interior angle D= 99 degrees 45.0 minutes\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "interiorB=(45+(30/60))+180-(120+(15/60))\n", + "deg=int(interiorB)\n", + "mins=(interiorB-deg)*60\n", + "print \"Interior angle B=\",deg,\"degrees\",mins,\"minutes\"\n", + "\n", + "interiorC=(120+(15/60))+180-(200+(30/60))\n", + "deg=int(interiorC)\n", + "mins=(interiorC-deg)*60\n", + "print \"Interior angle C=\",deg,\"degrees\",mins,\"minutes\"\n", + "\n", + "exteriorD=(280+(45/60))+180-(200+(30/60))\n", + "deg=int(exteriorD)\n", + "mins=(exteriorD-deg)*60\n", + "print \"Exterior angle D=\",deg,\"degrees\",mins,\"minutes\"\n", + "\n", + "interiorD=360-(260+(15/60))\n", + "deg=int(interiorD)\n", + "mins=(interiorD-deg)*60\n", + "print \"Interior angle D=\",deg,\"degrees\",mins,\"minutes\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### section 3.13, problem 3" + ] + }, + { + "cell_type": "code", + "execution_count": 13, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "80 degrees 30.0 minutes\n", + "200 degrees 30.0 minutes\n", + "320 degrees 30.0 minutes\n", + "80 degrees 30.0 minutes\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "import math\n", + "\n", + "\n", + "\n", + "FB_AB=80+(30/60)\n", + "FB_BC=FB_AB+180-60\n", + "FB_CA=FB_BC-180+300\n", + "\n", + "\n", + "deg1=int(FB_AB)\n", + "mins1=(FB_AB-deg1)*60\n", + "deg2=int(FB_BC)\n", + "mins2=(FB_BC-deg2)*60\n", + "deg3=int(FB_CA)\n", + "mins3=(FB_CA-deg3)*60\n", + "\n", + "\n", + "print deg1,\"degrees\",mins1,\"minutes\";\n", + "print deg2,\"degrees\",mins2,\"minutes\";\n", + "print deg3,\"degrees\",mins3,\"minutes\";\n", + "print deg1,\"degrees\",mins1,\"minutes\";\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### section 3.13, problem 3" + ] + }, + { + "cell_type": "code", + "execution_count": 14, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "120 degrees 30.0 minutes\n", + "210 degrees 30.0 minutes\n", + "300 degrees 30.0 minutes\n", + "30 degrees 30.0 minutes\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "import math\n", + "\n", + "\n", + "\n", + "FB_AB=120+(30/60)\n", + "FB_BC=FB_AB+180-90\n", + "FB_CD=FB_BC-180+270\n", + "FB_DA=FB_CD-180-90\n", + "\n", + "deg1=int(FB_AB)\n", + "mins1=(FB_AB-deg1)*60\n", + "deg2=int(FB_BC)\n", + "mins2=(FB_BC-deg2)*60\n", + "deg3=int(FB_CD)\n", + "mins3=(FB_CD-deg3)*60\n", + "deg4=int(FB_DA)\n", + "mins4=(FB_DA-deg4)*60\n", + "\n", + "print deg1,\"degrees\",mins1,\"minutes\";\n", + "print deg2,\"degrees\",mins2,\"minutes\";\n", + "print deg3,\"degrees\",mins3,\"minutes\";\n", + "print deg4,\"degrees\",mins4,\"minutes\";\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### Chapter 3, section 3.13, pg 91, problem 5" + ] + }, + { + "cell_type": "code", + "execution_count": 15, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "30 degrees 0.0 minutes\n", + "318 degrees 0.0 minutes\n", + "246 degrees 0.0 minutes\n", + "174 degrees 0.0 minutes\n", + "102 degrees 0 minutes\n", + "30 degrees 0 minutes\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "\n", + "interiorB=540/5\n", + "FB_AB=30+(0/60)\n", + "FB_BC=FB_AB+180+interiorB\n", + "FB_CD=FB_BC-180+interiorB\n", + "FB_DE=FB_CD-180+interiorB\n", + "FB_EA=FB_DE+180-360+interiorB\n", + "FB_AB=FB_EA+180-360+interiorB\n", + "deg1=int(FB_AB)\n", + "mins1=(FB_AB-deg1)*60\n", + "deg2=int(FB_BC)\n", + "mins2=(FB_BC-deg2)*60\n", + "deg3=int(FB_CD)\n", + "mins3=(FB_CD-deg3)*60\n", + "deg4=int(FB_DE)\n", + "mins4=(FB_DE-deg4)*60\n", + "deg5=int(FB_EA)\n", + "mins5=0\n", + "deg6=int(FB_AB)\n", + "mins6=0\n", + "print deg1,\"degrees\",mins1,\"minutes\"\n", + "print deg2,\"degrees\",mins2,\"minutes\"\n", + "print deg3,\"degrees\",mins3,\"minutes\"\n", + "print deg4,\"degrees\",mins4,\"minutes\"\n", + "print deg5,\"degrees\",mins5,\"minutes\"\n", + "print deg6,\"degrees\",mins6,\"minutes\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### pg 92, prob6" + ] + }, + { + "cell_type": "code", + "execution_count": 16, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "exterior angle A 150 degrees 15 minutes\n", + "interior angle A 209 degrees 45 minutes\n", + "Exterior angle B 309 degrees 45 minutes\n", + "interior angle B 50 degrees 15 minutes\n", + "interior angle C 95 degrees 15 minutes\n", + "interior angle D 102 degrees 15 minutes\n", + "interior angle E 82 degrees 30 minutes\n", + "540 degrees 540.0 degrees\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "BB_AB=330+(15/60)\n", + "BB_BC=200+(30/60)\n", + "BB_CD=115+(45/60)\n", + "BB_DE=38+(0/60)\n", + "BB_EA=300+(30/60)\n", + "\n", + "exteriorA=BB_EA-(150+(15/60))\n", + "interiorA=360-exteriorA\n", + "exteriorB=BB_AB-(20+(30/60))\n", + "interiorB=360-exteriorB\n", + "interiorC=(295+(45/60))-BB_BC\n", + "interiorD=218-BB_CD\n", + "interiorE=(120.5)-BB_DE\n", + "\n", + "deg1=int(exteriorA)\n", + "mins1=int((exteriorA-deg1)*60)\n", + "deg2=int(interiorA)\n", + "mins2=int((interiorA-deg2)*60)\n", + "deg3=int(exteriorB)\n", + "mins3=int((exteriorB-deg3)*60)\n", + "deg4=int(interiorB)\n", + "mins4=int((interiorB-deg4)*60)\n", + "deg5=int(interiorC)\n", + "mins5=int((interiorC-deg5)*60)\n", + "deg6=int(interiorD)\n", + "mins6=int((interiorD-deg6)*60)\n", + "deg7=int(interiorE)\n", + "mins7=int((interiorE-deg7)*60)\n", + "\n", + "n=5\n", + "check=(2*n-4)*90\n", + "summ=interiorA+interiorB+interiorC+interiorD+interiorE\n", + "\n", + "print \"exterior angle A\",deg1,\"degrees\",mins1,\"minutes\"\n", + "print \"interior angle A\",deg2,\"degrees\",mins2,\"minutes\"\n", + "print \"Exterior angle B\",deg3,\"degrees\",mins3,\"minutes\"\n", + "print \"interior angle B\",deg4,\"degrees\",mins4,\"minutes\"\n", + "print \"interior angle C\",deg5,\"degrees\",mins5,\"minutes\"\n", + "print \"interior angle D\",deg6,\"degrees\",mins6,\"minutes\"\n", + "print \"interior angle E\",deg7,\"degrees\",mins7,\"minutes\"\n", + "print check,\"degrees\",summ,\"degrees\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### pg 93, prob7" + ] + }, + { + "cell_type": "code", + "execution_count": 17, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Interior angle A 58 degrees 45 minutes\n", + "Interior angle B 105 degrees 30 minutes\n", + "Interior angle C 109 degrees 30 minutes\n", + "Interior angle D 86 degrees 15 minutes\n", + "360 degrees 360.0 degrees\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "FB_AB=45+(30/60)\n", + "FB_BC=60+(0/60)\n", + "FB_CD=10+(30/60)\n", + "FB_DA=75+(45/60)\n", + "\n", + "\n", + "BB_AB=45+(30/60)\n", + "BB_BC=60+(0/60)\n", + "BB_CD=10+(30/60)\n", + "BB_DA=75+(45/60)\n", + "\n", + "\n", + "interiorA=180-(FB_AB+BB_DA)\n", + "interiorB=(FB_BC+BB_AB)\n", + "interiorC=180-(BB_BC+FB_CD)\n", + "interiorD=(FB_DA+BB_CD)\n", + "\n", + "\n", + "deg1=int(interiorA)\n", + "mins1=int((interiorA-deg1)*60)\n", + "deg2=int(interiorB)\n", + "mins2=int((interiorB-deg2)*60)\n", + "deg3=int(interiorC)\n", + "mins3=int((interiorC-deg3)*60)\n", + "deg4=int(interiorD)\n", + "mins4=int((interiorD-deg4)*60)\n", + "\n", + "\n", + "n=4\n", + "check=(2*n-4)*90\n", + "summ=interiorA+interiorB+interiorC+interiorD\n", + "\n", + "print \"Interior angle A\",deg1,\"degrees\",mins1,\"minutes\"\n", + "print \"Interior angle B\",deg2,\"degrees\",mins2,\"minutes\"\n", + "print \"Interior angle C\",deg3,\"degrees\",mins3,\"minutes\"\n", + "print \"Interior angle D\",deg4,\"degrees\",mins4,\"minutes\"\n", + "print check,\"degrees\",summ,\"degrees\"\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### pg 93, prob8" + ] + }, + { + "cell_type": "code", + "execution_count": 18, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Interior angle A= 79 degrees 30 minutes\n", + "Interior angle B= 99 degrees 30 minutes\n", + "exterior angle C= 258 degrees 15 minutes\n", + "Interior angle C= 101 degrees 45 minutes\n", + "exterior angle D= 170 degrees 15 minutes\n", + "Interior angle D= 189 degrees 45 minutes\n", + "Interior angle E= 70 degrees 30 minutes\n", + "540 degrees 541.0 degrees\n", + "error= 1 degrees\n", + "-12.0 minutes\n", + "corrected values are:\n", + "Interior angle A= 79 degrees 18.0 minutes\n", + "Interior angle B= 99 degrees 18.0 minutes\n", + "Interior angle C= 101 degrees 33.0 minutes\n", + "Interior angle D= 189 degrees 33.0 minutes\n", + "Interior angle E= 70 degrees 18.0 minutes\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "BB_AB=330+(0/60)\n", + "BB_BC=48+(0/60)\n", + "BB_CD=127+(45/60)\n", + "BB_DE=120+(0/60)\n", + "BB_EA=229+(30/60)\n", + "\n", + "FB_AB=150+(0/60)\n", + "FB_BC=230+(30/60)\n", + "FB_CD=306+(15/60)\n", + "FB_DE=298+(0/60)\n", + "FB_EA=49+(30/60)\n", + "\n", + "\n", + "interiorA=BB_EA-FB_AB\n", + "interiorB=BB_AB-FB_BC\n", + "exteriorC=FB_CD-BB_BC\n", + "interiorC=360-(258+(15/60))\n", + "exteriorD=FB_DE-BB_CD\n", + "interiorD=360-exteriorD\n", + "interiorE=BB_DE-FB_EA\n", + "\n", + "deg1=int(interiorA)\n", + "mins1=int((interiorA-deg1)*60)\n", + "deg2=int(interiorB)\n", + "mins2=int((interiorB-deg2)*60)\n", + "deg3=int(exteriorC)\n", + "mins3=int((exteriorC-deg3)*60)\n", + "deg4=int(interiorC)\n", + "mins4=int((interiorC-deg4)*60)\n", + "deg5=int(exteriorD)\n", + "mins5=int((exteriorD-deg5)*60)\n", + "deg6=int(interiorD)\n", + "mins6=int((interiorD-deg6)*60)\n", + "deg7=int(interiorE)\n", + "mins7=int((interiorE-deg7)*60)\n", + "\n", + "n=5\n", + "check=(2*n-4)*90\n", + "summ=interiorA+interiorB+interiorC+interiorD+interiorE\n", + "\n", + "print \"Interior angle A=\",deg1,\"degrees\",mins1,\"minutes\"\n", + "print \"Interior angle B=\",deg2,\"degrees\",mins2,\"minutes\"\n", + "print \"exterior angle C=\",deg3,\"degrees\",mins3,\"minutes\"\n", + "print \"Interior angle C=\",deg4,\"degrees\",mins4,\"minutes\"\n", + "print \"exterior angle D=\",deg5,\"degrees\",mins5,\"minutes\"\n", + "print \"Interior angle D=\",deg6,\"degrees\",mins6,\"minutes\"\n", + "print \"Interior angle E=\",deg7,\"degrees\",mins7,\"minutes\"\n", + "print check,\"degrees\",summ,\"degrees\"\n", + "\n", + "error=541-540\n", + "correction=(-60/5)\n", + "print \"error=\",error,\"degrees\"\n", + "print correction,\"minutes\"\n", + "\n", + "correctedvalue1=mins1+correction\n", + "correctedvalue2=mins2+correction\n", + "correctedvalue4=mins4+correction\n", + "correctedvalue6=mins6+correction\n", + "correctedvalue7=mins7+correction\n", + "\n", + "print \"corrected values are:\"\n", + "print \"Interior angle A=\",deg1,\"degrees\",correctedvalue1,\"minutes\"\n", + "print \"Interior angle B=\",deg2,\"degrees\",correctedvalue2,\"minutes\"\n", + "print \"Interior angle C=\",deg4,\"degrees\",correctedvalue4,\"minutes\"\n", + "print \"Interior angle D=\",deg6,\"degrees\",correctedvalue6,\"minutes\"\n", + "print \"Interior angle E=\",deg7,\"degrees\",correctedvalue7,\"minutes\"\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### pg 95, prob1" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Interior angle A= 44 degrees 0 minutes\n", + "Interior angle B= 26 degrees 30 minutes\n", + "exterior angle C= 200 degrees 15 minutes\n", + "Interior angle C= 159 degrees 45 minutes\n", + "Interior angle D= 42 degrees 15 minutes\n", + "Interior angle E= 267 degrees 30 minutes\n", + "540 degrees 540.0 degrees\n", + "242.75 correct 330.25 correct\n", + "corrected values are:\n", + "FB_AB= 194 degrees 15 minutes\n", + "FB_BC= 40 degrees 45 minutes\n", + "FB_CD= 20 degrees 30 minutes\n", + "FB_DE= 242 degrees 45 minutes\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "BB_AB=13+(0/60)\n", + "BB_BC=222+(30/60)\n", + "BB_CD=200+(30/60)\n", + "BB_DE=62+(45/60)\n", + "BB_EA=147+(45/60)\n", + "\n", + "FB_AB=191+(45/60)\n", + "FB_BC=39+(30/60)\n", + "FB_CD=22+(15/60)\n", + "FB_DE=242+(45/60)\n", + "FB_EA=330+(15/60)\n", + "\n", + "#(a)\n", + "interiorA=FB_AB-BB_EA\n", + "interiorB=FB_BC-BB_AB\n", + "exteriorC=BB_BC-FB_CD\n", + "interiorC=360-(200+(15/60))\n", + "interiorD=FB_DE-BB_CD\n", + "interiorE=FB_EA-BB_DE\n", + "\n", + "deg1=int(interiorA)\n", + "mins1=int((interiorA-deg1)*60)\n", + "deg2=int(interiorB)\n", + "mins2=int((interiorB-deg2)*60)\n", + "deg3=int(exteriorC)\n", + "mins3=int((exteriorC-deg3)*60)\n", + "deg4=int(interiorC)\n", + "mins4=int((interiorC-deg4)*60)\n", + "deg6=int(interiorD)\n", + "mins6=int((interiorD-deg6)*60)\n", + "deg7=int(interiorE)\n", + "mins7=int((interiorE-deg7)*60)\n", + "\n", + "n=5\n", + "check=(2*n-4)*90\n", + "summ=interiorA+interiorB+interiorC+interiorD+interiorE\n", + "\n", + "print \"Interior angle A=\",deg1,\"degrees\",mins1,\"minutes\"\n", + "print \"Interior angle B=\",deg2,\"degrees\",mins2,\"minutes\"\n", + "print \"exterior angle C=\",deg3,\"degrees\",mins3,\"minutes\"\n", + "print \"Interior angle C=\",deg4,\"degrees\",mins4,\"minutes\"\n", + "print \"Interior angle D=\",deg6,\"degrees\",mins6,\"minutes\"\n", + "print \"Interior angle E=\",deg7,\"degrees\",mins7,\"minutes\"\n", + "print check,\"degrees\",summ,\"degrees\"\n", + "\n", + "#(b)\n", + "\n", + "print FB_DE,\"correct\",FB_EA,\"correct\"\n", + "\n", + "AB=FB_EA-180+interiorA\n", + "BC=(194+(15/60))-180+interiorB\n", + "CD=(40+(45/60))+180-exteriorC\n", + "DE=(20+(30/60))+180+interiorD\n", + "\n", + "deg1=int(AB)\n", + "mins1=int((AB-deg1)*60)\n", + "deg2=int(BC)\n", + "mins2=int((BC-deg2)*60)\n", + "deg3=int(CD)\n", + "mins3=int((CD-deg3)*60)\n", + "deg4=int(DE)\n", + "mins4=int((DE-deg4)*60)\n", + "\n", + "print \"corrected values are:\"\n", + "print \"FB_AB=\",deg1,\"degrees\",mins1,\"minutes\"\n", + "print \"FB_BC=\",deg2,\"degrees\",mins2,\"minutes\"\n", + "print \"FB_CD=\",deg3,\"degrees\",mins3,\"minutes\"\n", + "print \"FB_DE=\",deg4,\"degrees\",mins4,\"minutes\"\n", + "\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### pg 95, prob1" + ] + }, + { + "cell_type": "code", + "execution_count": 20, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "included angle A= 79 degrees 30 minutes\n", + "included angle B= 99 degrees 30 minutes\n", + "included angle C= 101 degrees 45 minutes\n", + "exterior angle D= 171 degrees 15 minutes\n", + "included angle D= 188 degrees 45 minutes\n", + "exterior angle D= 289 degrees 30 minutes\n", + "included angle E= 70 degrees 30 minutes\n", + "540 degrees 540.0 degrees\n", + "68.25 correct 148.75 correct 248.25 correct\n", + "correction= 1\n", + "corrected values are:\n", + "FB_AB= 68 degrees 15 minutes\n", + "FB_CD= 227 degrees 0 minutes\n", + "BB_CD= 47 degrees 0 minutes\n", + "FB_DE= 218 degrees 15 minutes\n", + "BB_DE= 38 degrees 15 minutes\n", + "AB=100m, BC=100m,CD=50m, scale=20m for plot\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "BB_AB=248+(15/60)\n", + "BB_BC=326+(15/60)\n", + "BB_CD=46+(0/60)\n", + "BB_DE=38+(15/60)\n", + "BB_EA=147+(45/60)\n", + "\n", + "FB_AB=68+(15/60)\n", + "FB_BC=148+(45/60)\n", + "FB_CD=224+(30/60)\n", + "FB_DE=217+(15/60)\n", + "FB_EA=327+(45/60)\n", + "\n", + "#(a)\n", + "includedA=-FB_AB+BB_EA\n", + "includedB=-FB_BC+BB_AB\n", + "includedC=BB_BC-FB_CD\n", + "includedD=360-(171+(15/60))\n", + "exteriorD=FB_DE-BB_CD\n", + "exteriorE=FB_EA-BB_DE\n", + "includedE=360-(289+(30/60))\n", + "\n", + "deg1=int(includedA)\n", + "mins1=int((includedA-deg1)*60)\n", + "deg2=int(includedB)\n", + "mins2=int((includedB-deg2)*60)\n", + "deg3=int(includedC)\n", + "mins3=int((includedC-deg3)*60)\n", + "deg4=int(exteriorD)\n", + "mins4=int((exteriorD-deg4)*60)\n", + "deg5=int(includedD)\n", + "mins5=int((includedD-deg5)*60)\n", + "deg6=int(exteriorE)\n", + "mins6=int((exteriorE-deg6)*60)\n", + "deg7=int(includedE)\n", + "mins7=int((includedE-deg7)*60)\n", + "\n", + "n=5\n", + "check=(2*n-4)*90\n", + "summ=includedA+includedB+includedC+includedD+includedE\n", + "\n", + "print \"included angle A=\",deg1,\"degrees\",mins1,\"minutes\"\n", + "print \"included angle B=\",deg2,\"degrees\",mins2,\"minutes\"\n", + "print \"included angle C=\",deg3,\"degrees\",mins3,\"minutes\"\n", + "print \"exterior angle D=\",deg4,\"degrees\",mins4,\"minutes\"\n", + "print \"included angle D=\",deg5,\"degrees\",mins5,\"minutes\"\n", + "print \"exterior angle D=\",deg6,\"degrees\",mins6,\"minutes\"\n", + "print \"included angle E=\",deg7,\"degrees\",mins7,\"minutes\"\n", + "print check,\"degrees\",summ,\"degrees\"\n", + "\n", + "#(b)\n", + "\n", + "print FB_AB,\"correct\",FB_BC,\"correct\",BB_AB,\"correct\"\n", + "\n", + "\n", + "FB_BC=(328+(45/60))-(326+(15/60))\n", + "FB_CD=(224+(30/60))+FB_BC\n", + "BB_CD=227-180\n", + "correctionatD=1 \n", + "FB_DE=(217+(15/60))+1\n", + "BB_DE=FB_DE-180\n", + "\n", + "deg1=int(FB_AB)\n", + "mins1=int((FB_AB-deg1)*60)\n", + "deg2=int(FB_CD)\n", + "mins2=int((FB_CD-deg2)*60)\n", + "deg3=int(BB_CD)\n", + "mins3=int((BB_CD-deg3)*60)\n", + "deg4=int(FB_DE)\n", + "mins4=int((FB_DE-deg4)*60)\n", + "deg5=int(BB_DE)\n", + "mins5=int((BB_DE-deg5)*60) \n", + "\n", + "print \"correction=\",correctionatD\n", + "print \"corrected values are:\"\n", + "print \"FB_AB=\",deg1,\"degrees\",mins1,\"minutes\"\n", + "print \"FB_CD=\",deg2,\"degrees\",mins2,\"minutes\"\n", + "print \"BB_CD=\",deg3,\"degrees\",mins3,\"minutes\"\n", + "print \"FB_DE=\",deg4,\"degrees\",mins4,\"minutes\"\n", + "print \"BB_DE=\",deg5,\"degrees\",mins5,\"minutes\"\n", + "print \"AB=100m, BC=100m,CD=50m, scale=20m for plot\"\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### pg 100, prob3" + ] + }, + { + "cell_type": "code", + "execution_count": 21, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "59.0 correct 139.5 correct 239.0 correct\n", + "correctionatC= 2.5\n", + "correctionatD= 1.25\n", + "correctionatE= 0.25\n", + "corrected values are:\n", + "BB_CD= 217.75 BB_DE= 209.25 BB_EA= 138.75\n", + "FB_CD= 217 degrees 45 minutes\n", + "FB_DE= 209 degrees 15 minutes\n", + "FB_EA= 318 degrees 45 minutes\n", + "declination= -10 degrees W\n", + "true bearing values:\n", + "BB_AB= 229.0\n", + "BB_BC= 309.5\n", + "BB_CD= 27.75\n", + "BB_DE= 19.0\n", + "BB_EA= 128.75\n", + "FB_AB= 49\n", + "FB_BC= 129.5\n", + "FB_CD= 207.75\n", + "FB_DE= 199.25\n", + "FB_EA= 308.75\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "BB_AB=239+(00/60)\n", + "BB_BC=317+(00/60)\n", + "BB_CD=36+(30/60)\n", + "BB_DE=29+(00/60)\n", + "BB_EA=138+(45/60)\n", + "\n", + "FB_AB=59+(00/60)\n", + "FB_BC=139+(30/60)\n", + "FB_CD=215+(15/60)\n", + "FB_DE=208+(0/60)\n", + "FB_EA=318+(30/60)\n", + "\n", + "print FB_AB,\"correct\",FB_BC,\"correct\",BB_AB,\"correct\"\n", + "\n", + "correctionatC=2+(30/60)\n", + "FB_CD=(215+(15/60))+correctionatC\n", + "correctionatD=1+(15/60)\n", + "FB_DE=208+correctionatD\n", + "correctionatE=(15/60)\n", + "FB_EA=(318+(30/60))+correctionatE\n", + "\n", + "\n", + "deg2=int(FB_CD)\n", + "mins2=int((FB_CD-deg2)*60)\n", + "deg4=int(FB_DE)\n", + "mins4=int((FB_DE-deg4)*60)\n", + "deg5=int(FB_EA)\n", + "mins5=int((FB_EA-deg5)*60) \n", + "\n", + "print \"correctionatC=\",correctionatC\n", + "print \"correctionatD=\",correctionatD\n", + "print \"correctionatE=\",correctionatE\n", + "print \"corrected values are:\"\n", + "print \"BB_CD=\",217.75,\" BB_DE=\",209.25,\" BB_EA=\",138.75\n", + "print \"FB_CD=\",deg2,\"degrees\",mins2,\"minutes\"\n", + "print \"FB_DE=\",deg4,\"degrees\",mins4,\"minutes\"\n", + "print \"FB_EA=\",deg5,\"degrees\",mins5,\"minutes\"\n", + "print \"declination=\",-10,\"degrees W\"\n", + "\n", + "BB_AB=239+(00/60)-10\n", + "BB_BC=317+(00/60)-10+correctionatC\n", + "BB_CD=36+(30/60)-10+correctionatD\n", + "BB_DE=29+(00/60)-10\n", + "BB_EA=138+(45/60)-10\n", + "\n", + "FB_AB=59-10\n", + "FB_BC=(139+(30/60))-10\n", + "FB_CD=(215+(15/60))-10+correctionatC\n", + "FB_DE=(208+(0/60))-10+correctionatD\n", + "FB_EA=(318+(30/60))-10+correctionatE\n", + "\n", + "print \"true bearing values:\"\n", + "print \"BB_AB=\",BB_AB \n", + "print \"BB_BC=\",BB_BC\n", + "print \"BB_CD=\",BB_CD\n", + "print \"BB_DE=\",BB_DE\n", + "print \"BB_EA=\",BB_EA\n", + "\n", + "print \"FB_AB=\",FB_AB\n", + "print \"FB_BC=\",FB_BC\n", + "print \"FB_CD=\",FB_CD\n", + "print \"FB_DE=\",FB_DE\n", + "print \"FB_EA=\",FB_EA\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### pg 102, prob4" + ] + }, + { + "cell_type": "code", + "execution_count": 22, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "45.5 correct 60.0 correct 45.5 correct\n", + "correctionatC= 0.666666666667\n", + "correctionatD= 1.5\n", + "corrected values are:\n", + "BB_CD=N 4.83 W BB_BC=N 60 degrees W\n", + "FB_CD=N 4.83 W FB_DA=N 85 degrees W\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "BB_AB=45+(30/60)\n", + "BB_BC=60+(40/60)\n", + "BB_CD=3+(20/60)\n", + "BB_DA=85+(00/60)\n", + "\n", + "\n", + "FB_AB=45+(30/60)\n", + "FB_BC=60+(0/60)\n", + "FB_CD=5+(30/60)\n", + "FB_DA=83+(30/60)\n", + "\n", + "\n", + "print FB_AB,\"correct\",FB_BC,\"correct\",BB_AB,\"correct\"\n", + "\n", + "correctionatC=-0+(40/60)\n", + "FB_CD=(5+(30/60))+correctionatC\n", + "correctionatD=1+(30/60)\n", + "FB_DA=83+(30/60)+correctionatD\n", + "\n", + "\n", + "\n", + "deg2=int(FB_CD)\n", + "mins2=int((FB_CD-deg2)*60)\n", + "deg4=int(FB_DA)\n", + "mins4=int((FB_DA-deg4)*60)\n", + " \n", + "\n", + "print \"correctionatC=\",correctionatC;\n", + "print \"correctionatD=\",correctionatD\n", + "\n", + "print \"corrected values are:\";\n", + "print \"BB_CD=N\",4.83,\"W\",\" BB_BC=N\",60,\"degrees W\";\n", + "print \"FB_CD=N\",4.83,\"W\",\" FB_DA=N\",85,\"degrees W\"\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### pg 102, prob4" + ] + }, + { + "cell_type": "code", + "execution_count": 23, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "80.0 correct 40.5 correct 80.0 correct\n", + "correctionatB= 0.75\n", + "correctionatC= 0.5\n", + "corrected values are:\n", + "BB_AB=N 40.5 E BB_BC=N 80 degrees E\n", + "FB_CD=N 20 E FB_DA=S 80 degrees E\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "BB_AB=41+(15/60)\n", + "BB_BC=79+(30/60)\n", + "BB_CD=20+(0/60)\n", + "BB_DA=80+(00/60)\n", + "\n", + "\n", + "FB_AB=40+(30/60)\n", + "FB_BC=80+(45/60)\n", + "FB_CD=19+(30/60)\n", + "FB_DA=80+(00/60)\n", + "\n", + "\n", + "print FB_DA,\"correct\",FB_AB,\"correct\",BB_DA,\"correct\";\n", + "\n", + "correctionatB=-0+(45/60)\n", + "FB_BC=(80+(45/60))+correctionatB\n", + "correctionatC=0+(30/60)\n", + "FB_CD=19+(30/60)+correctionatC\n", + " \n", + "\n", + "print \"correctionatB=\",correctionatB;\n", + "print \"correctionatC=\",correctionatC;\n", + "\n", + "print \"corrected values are:\";\n", + "print \"BB_AB=N\",40.5,\"E\",\" BB_BC=N\",80,\"degrees E\";\n", + "print \"FB_CD=N\",20,\"E\",\" FB_DA=S\",80,\"degrees E\"\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### pg 104, prob6" + ] + }, + { + "cell_type": "code", + "execution_count": 24, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "included angle A= 79 degrees 45 minutes\n", + "included angle B= 99 degrees 30 minutes\n", + "included angle C= 101 degrees 45 minutes\n", + "exterior angle D= 171 degrees 30 minutes\n", + "included angle D= 188 degrees 30 minutes\n", + "exterior angle D= 289 degrees 30 minutes\n", + "included angle E= 70 degrees 30 minutes\n", + "540 degrees 540.0 degrees\n", + "59.0 correct 139.5 correct 239.0 correct\n", + "correction= 1.25\n", + "corrected values are:\n", + "BB_BC= 319.5 BB_CD= 73.75 degrees BB_DE= 29.25 degrees\n", + "FB_CD= 217.75 FB_DE= 209.25 degrees FB_EA= 318.75 degrees\n" + ] + } + ], + "source": [ + "\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "BB_AB=239+(00/60)\n", + "BB_BC=317+(0/60)\n", + "BB_CD=36+(30/60)\n", + "BB_DE=29+(00/60)\n", + "BB_EA=138+(45/60)\n", + "\n", + "FB_AB=59+(0/60)\n", + "FB_BC=139+(30/60)\n", + "FB_CD=215+(15/60)\n", + "FB_DE=208+(0/60)\n", + "FB_EA=318+(30/60)\n", + "\n", + "#(a)\n", + "includedA=-FB_AB+BB_EA\n", + "includedB=-FB_BC+BB_AB\n", + "includedC=BB_BC-FB_CD\n", + "includedD=360-(171+(30/60))\n", + "exteriorD=FB_DE-BB_CD\n", + "exteriorE=FB_EA-BB_DE\n", + "includedE=360-(289+(30/60))\n", + "\n", + "deg1=int(includedA)\n", + "mins1=int((includedA-deg1)*60)\n", + "deg2=int(includedB)\n", + "mins2=int((includedB-deg2)*60)\n", + "deg3=int(includedC)\n", + "mins3=int((includedC-deg3)*60)\n", + "deg4=int(exteriorD)\n", + "mins4=int((exteriorD-deg4)*60)\n", + "deg5=int(includedD)\n", + "mins5=int((includedD-deg5)*60)\n", + "deg6=int(exteriorE)\n", + "mins6=int((exteriorE-deg6)*60)\n", + "deg7=int(includedE)\n", + "mins7=int((includedE-deg7)*60)\n", + "\n", + "n=5\n", + "check=(2*n-4)*90\n", + "summ=includedA+includedB+includedC+includedD+includedE\n", + "\n", + "print \"included angle A=\",deg1,\"degrees\",mins1,\"minutes\"\n", + "print \"included angle B=\",deg2,\"degrees\",mins2,\"minutes\"\n", + "print \"included angle C=\",deg3,\"degrees\",mins3,\"minutes\"\n", + "print \"exterior angle D=\",deg4,\"degrees\",mins4,\"minutes\"\n", + "print \"included angle D=\",deg5,\"degrees\",mins5,\"minutes\"\n", + "print \"exterior angle D=\",deg6,\"degrees\",mins6,\"minutes\"\n", + "print \"included angle E=\",deg7,\"degrees\",mins7,\"minutes\"\n", + "print check,\"degrees\",summ,\"degrees\"\n", + "\n", + "#(b)\n", + "\n", + "print FB_AB,\"correct\",FB_BC,\"correct\",BB_AB,\"correct\"\n", + "\n", + "\n", + "\n", + "FB_CD=(215+(15/60))+(2+(30/60))\n", + "BB_CD=(37+(45/60))\n", + "correctionatD=(1+(15/60)) \n", + "FB_DE=(208+(0/60))+correctionatD\n", + "FB_EA=(318+(30/60))+(0+(15/60))\n", + "\n", + " \n", + "print \"correction=\",correctionatD;\n", + "print \"corrected values are:\";\n", + "print \"BB_BC=\",319.5, \"BB_CD=\",73.75,\"degrees\", \"BB_DE=\",29.25,\"degrees\";\n", + "print \"FB_CD=\",217.75, \"FB_DE=\",209.25,\"degrees\", \"FB_EA=\",318.75,\"degrees\";\n", + "\n", + "\n" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 2", + "language": "python", + "name": "python2" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 2 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython2", + "version": "2.7.11" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap5_Levelling.ipynb b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap5_Levelling.ipynb new file mode 100644 index 00000000..a73f9fc8 --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap5_Levelling.ipynb @@ -0,0 +1,903 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "# Chapter 5: Levelling" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### ch-5 page 151, pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('correct reading on A =', 2.524965929375, 'meters')\n", + "('correct reading of B =', 1.75499327, 'meters')\n", + "('true difference is', 0.769972659375, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "a=150;\n", + "b=100;\n", + "ar=2.525;\n", + "br=1.755;\n", + "\n", + "sc=1000;\n", + "d=(a*a)/(sc*sc);\n", + "\n", + "A=0.0673*d*d;\n", + "\n", + "fa=ar-A;\n", + "\n", + "print('correct reading on A =',fa,'meters');\n", + "\n", + "\n", + "d=(b*b)/(sc*sc);\n", + "\n", + "B=0.0673*d*d;\n", + "fb=br-B;\n", + "\n", + "print('correct reading of B =',fb,'meters');\n", + "\n", + "AB=fa-fb;\n", + "print('true difference is',AB,'meters');\n", + "\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page 152, pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('heigght of lighthouse is', 60.57000000000001, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "\n", + "d=30;\n", + "sc=1000;\n", + "\n", + "h=0.0673*d*d;\n", + "\n", + "print('heigght of lighthouse is',h,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page 152, pb-3" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('D=', 27.25696334003587)\n", + "('dimp of horizon', 0.0042789581381531975, 'degrees')\n", + "('dimp of horizon', 14.709974521760092, 'minutes')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "h=50;\n", + "\n", + "d=math.sqrt(h/0.0673);\n", + "print('D=',d);\n", + "\n", + "r=6370;\n", + "dip=d/r;\n", + "print('dimp of horizon',dip,'degrees');\n", + "\n", + "dip1=dip*((180*60)/math.pi)\n", + "print('dimp of horizon',dip1,'minutes');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page 152,153, pb-4" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('distance between man and object is', 39.44664791774385, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "h1=50;\n", + "h2=10;\n", + "c=0.0673;\n", + "\n", + "d1=math.sqrt(h1/c);\n", + "\n", + "d2=math.sqrt(h2/c);\n", + "\n", + "dis=d1+d2;\n", + "\n", + "print('distance between man and object is',dis,'meters');\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### ch-5 page-153, pb-5" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('height of the hill is ', 309.46147646724046, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "h1=10;\n", + "c=0.0673\n", + "d1=math.sqrt(h1/c);\n", + "\n", + "d2=d1-80; #since d1+d2=80;\n", + "h2=c*d2*d2;\n", + "\n", + "print('height of the hill is ',h2,'meters');\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page-153,154 pb-6" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('distance AB =', 86.24055457549457, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "h1=100;\n", + "h2=150;\n", + "\n", + "r2=12880;\n", + "c=(6/7)*(1000/r2);\n", + "d1=math.sqrt(h1/c)\n", + "d2=math.sqrt(h2/c)\n", + "\n", + "d=d1+d2;\n", + "print('distance AB =',d,'meters');\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page-154 pb-7" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('R=', 20.00000000000007)\n", + "('sensitiveness of bubble is ', 20.626499999999925, 'seconds')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "r1=2.550;\n", + "r2=2.500;\n", + "\n", + "s=r1-r2;\n", + "d=0.002;\n", + "D=100;\n", + "n=5;\n", + "r=(n*d*D/s);\n", + "\n", + "print('R=',r);\n", + "\n", + "alp=(s/(n*D))*206265;\n", + "\n", + "print('sensitiveness of bubble is ',alp,'seconds');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page-154,155 pb-8" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('error is ', 0.01939252902819189, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "n=2;\n", + "D=100;\n", + "alp=20;\n", + "\n", + "\n", + "s=(alp*n*D)/206265;\n", + "\n", + "print('error is ',s,'meters');\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page-156, pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('true level of difference is', 1.115, 'meters')\n", + "('RL of B =', 124.435, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "a=2.245;\n", + "b=3.375;\n", + "AB=b-a;\n", + "\n", + "ap=1.955;\n", + "bp=3.055;\n", + "\n", + "dAB=bp-ap;\n", + "\n", + "tl=(AB+dAB)/2;\n", + "print('true level of difference is',tl,'meters')\n", + "rla=125.55;\n", + "rlb=rla-tl;\n", + "\n", + "\n", + "print('RL of B =',rlb,'meters');\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page 157, pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('true RL of B', 524.065, 'meters')\n", + "('combined corrction for 500m=', 0.016825, 'meters')\n", + "('collimation error per 100m=', -0.0023599999999999997, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "aa=1.155;\n", + "ab=2.595;\n", + "ba=0.985;\n", + "bb=2.415;\n", + "\n", + "td=((ab-aa)+(bb-ba))/2\n", + "\n", + "rla=525.5;\n", + "rlb=rla-td;\n", + "dab=500;\n", + "print('true RL of B',rlb,'meters');\n", + "\n", + "dab1=dab/1000;\n", + "\n", + "correct=0.0673*dab1*dab1;\n", + "print('combined corrction for 500m=',correct,'meters');\n", + "\n", + "sc=100;\n", + "a=1.155;\n", + "e=-(0.0118*sc)/(dab);\n", + "\n", + "\n", + "print('collimation error per 100m=',e,'meters')\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page 157,158, pb-3" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('true difference between A and B is ', 0.33999999999999997, 'meters')\n", + "('amount of collimation error =', -0.015000000000000124, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "aa=1.725;\n", + "ab=1.370;\n", + "ba=1.560;\n", + "bb=1.235;\n", + "\n", + "A=aa-ab;\n", + "B=ba-bb;\n", + "\n", + "AB=(A+B)/2;\n", + "\n", + "print('true difference between A and B is ',AB,'meters');\n", + "\n", + "CB=bb;\n", + "CA=CB+AB;\n", + "\n", + "OCA=1.560;\n", + "e=OCA-CA;\n", + "\n", + "print('amount of collimation error =',e,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### ch-5 page 158,159, pb-4" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('apparent difference of level between A and B is', 0.52, 'meters')\n", + "('apparent difference of level at B', 0.8999999999999999, 'meters')\n", + "('true differece of level=', 0.71)\n", + "('correction to be applied at A is =', -0.18999999999999995)\n", + "('RL of B=', 449.29, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "aa=1.725;\n", + "ab=2.245;\n", + "ba=2.145;\n", + "bb=3.045;\n", + "\n", + "AB=200;\n", + "rla=450;\n", + "\n", + "\n", + "aAB=ab-aa;\n", + "\n", + "print('apparent difference of level between A and B is',aAB,'meters');\n", + "\n", + "dB=bb-ba\n", + "\n", + "print('apparent difference of level at B',dB,'meters')\n", + "\n", + "td=(aAB+dB)/2;\n", + "\n", + "print('true differece of level=',td);\n", + "\n", + "CB=bb;\n", + "\n", + "CA=CB-td;\n", + "\n", + "e=ba-CA;\n", + "\n", + "print('correction to be applied at A is =',e)\n", + "\n", + "rlb=rla-td;\n", + "\n", + "print('RL of B=',rlb,'meters')\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page 185,186 pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 13, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "0.52\n", + "('apparent difference of level =', 0.8999999999999999, 'meters')\n", + "('true difference of level=', 0.71, 'meters')\n", + "('true reading on A=', 2.335, 'meters')\n", + "('collimation error =', -0.18999999999999995, 'meters')\n", + "('RL of B=', 449.29, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "aa=1.725;\n", + "ab=2.245;\n", + "ba=2.145;\n", + "bb=3.045;\n", + "dAB=200;\n", + "rla=450.0;\n", + "AB=ab-aa;\n", + "print(AB)\n", + "adif=bb-ba\n", + "\n", + "print('apparent difference of level =',adif,'meters');\n", + "\n", + "#a\n", + "td=(AB+adif)/2;\n", + "print('true difference of level=',td,'meters')\n", + "#b\n", + "\n", + "tb=bb;\n", + "ta=bb-td;\n", + "\n", + "print('true reading on A=',ta,'meters');\n", + "\n", + "#c\n", + "\n", + "e=ba-ta;\n", + "\n", + "print('collimation error =',e,'meters');\n", + "\n", + "#d\n", + "\n", + "rlb=rla-td;\n", + "print('RL of B=',rlb,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page 186,187 pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 14, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('correct staff reading on B should be =', 1.0650000000000002, 'meters')\n", + "('collimation error is ', 0.08499999999999974, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "import math\n", + "\n", + "ma=1.585;\n", + "mb=1.225;\n", + "aa=1.425;\n", + "ab=1.150;\n", + "\n", + "dAB=100;\n", + "\n", + "#a\n", + "td=ma-mb;\n", + "B=aa-td;\n", + "\n", + "print('correct staff reading on B should be =',B,'meters');\n", + "\n", + "#c\n", + "\n", + "\n", + "e=ab-B;\n", + "print('collimation error is ',e,'meters')\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page 187 pb-3" + ] + }, + { + "cell_type": "code", + "execution_count": 15, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "first setting\n", + "('true difference is', 0.08499999999999996, 'meters')\n", + "('apparent difference of level =', 0.06999999999999984, 'meters')\n", + "second setting\n", + "('collimation error is', 0.015000000000000124, 'meters')\n", + "('correction at A=', 0.0015000000000000126, 'meters')\n", + "('correction at B=', 0.01650000000000014, 'meters')\n" + ] + } + ], + "source": [ + "#ch-5 page 187 pb-3\n", + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "dAB=100;\n", + "\n", + "aa=1.875;\n", + "ab=1.790;\n", + "\n", + "le=10;\n", + "\n", + "ba=1.630;\n", + "bb=1.560;\n", + "\n", + "\n", + "td=aa-ab;\n", + "\n", + "apd=ba-bb;\n", + "print('first setting')\n", + "print('true difference is',td,'meters');\n", + "print('apparent difference of level =',apd,'meters');\n", + "\n", + "print('second setting');\n", + "\n", + "A=ba-td;\n", + "\n", + "e1=bb-A\n", + "\n", + "cA=(le/dAB)*e1\n", + "cB=((le+dAB)/dAB)*e1\n", + "print('collimation error is',e1,'meters')\n", + "print('correction at A=',cA,'meters')\n", + "print('correction at B=',cB,'meters')\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page 163 pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 16, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "(0.645, 1.115, 0.7650000000000001, 0.23499999999999988, 0.85, 3.6100000000000003)\n", + "(2.835, 1.1949999999999998, 0.625, 1.375, 6.029999999999999)\n", + "('k=', -2.4200000000000017)\n", + "('k1=', -2.419999999999999)\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "import math\n", + "\n", + "bs1=2.375;bs2=2.835;bs3=0.435;\n", + "is1=1.730;is2=0.615;is3=2.070;is4=1.835;is5=1.630;\n", + "is6=2.255;\n", + "fs1=3.450;fs2=0.985;fs3=3.630;\n", + "\n", + "sbs=bs1+bs2+bs3;\n", + "sis=is1+is2+is3+is4+is5+is6;\n", + "sfs=fs1+fs2+fs3;\n", + "\n", + "r1=bs1-is1;\n", + "r2=is1-is2;\n", + "r3=bs2-is3;\n", + "r4=is3-is4;\n", + "r5=is4-fs2;\n", + "sr=r1+r2+r3+r4+r5;\n", + "print(r1,r2,r3,r4,r5,sr);\n", + "\n", + "\n", + "f1=bs2;\n", + "f2=is5-bs3;\n", + "f3=fs3-is6;\n", + "f4=is6-is5\n", + "sf=f1+f2+f3+f4;\n", + "print(f1,f2,f4,f3,sf);\n", + "\n", + "k=sbs-sfs\n", + "print('k=',k);\n", + "k1=sr-sf\n", + "print('k1=',k1);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-5 page 163,164 pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 17, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "(0.9049999999999998, 1.12, 1.7349999999999999, 1.365, 5.125)\n", + "(1.4749999999999999, 1.465, 0.665, 1.29, 4.8950000000000005)\n", + "('k=', 0.22999999999999954)\n", + "('k1=', 0.22999999999999954)\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "import math\n", + "\n", + "bs1=3.150;bs2=3.860;bs3=0.470;\n", + "is1=2.245;is2=2.125;is3=0.760;is4=1.935;is5=3.225;\n", + "fs1=1.125;fs2=2.235;fs3=3.890;\n", + "\n", + "sbs=bs1+bs2+bs3;\n", + "sis=is1+is2+is3+is4+is5;\n", + "sfs=fs1+fs2+fs3;\n", + "\n", + "r1=bs1-is1;\n", + "r2=is1-fs1;\n", + "r3=bs2-is2;\n", + "r4=is2-is3;\n", + "\n", + "sr=r1+r2+r3+r4;\n", + "print(r1,r2,r3,r4,sr);\n", + "\n", + "\n", + "f1=fs2-is3;\n", + "f2=is4-bs3;\n", + "f3=is5-is4;\n", + "f4=fs3-is5;\n", + "sf=f1+f2+f3+f4;\n", + "print(f1,f2,f4,f3,sf);\n", + "\n", + "k=sbs-sfs\n", + "print('k=',k);\n", + "k1=sr-sf\n", + "print('k1=',k1);\n" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 2", + "language": "python", + "name": "python2" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 2 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython2", + "version": "2.7.11" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap7_Computation-of-Area.ipynb b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap7_Computation-of-Area.ipynb new file mode 100644 index 00000000..25fe5def --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap7_Computation-of-Area.ipynb @@ -0,0 +1,802 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "# Chapter 7: Computation of Area" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### ch-7 page 207 pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "chainage 0 and 20\n", + "('area=', 420.0)\n", + "chainage 20 and 65\n", + "('area=', 2250.0)\n", + "chainage 65 and 110\n", + "('area=', 1305.0)\n", + "chainage 90 and 110\n", + "('area=', 600.0)\n", + "chainage 40 and 90\n", + "('area=', 2000.0)\n", + "chainage 0 and 40\n", + "('area=', 400.0)\n", + "('area of field =', 6975.0)\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "\n", + "print('chainage 0 and 20')\n", + "a1=0;b1=20;\n", + "\n", + "base=b1-a1;\n", + "o1=0;o2=42;\n", + "mo1=(o2+o1)/2;\n", + "\n", + "ae1=base*mo1;\n", + "print('area=',ae1);\n", + "\n", + "print('chainage 20 and 65')\n", + "a1=20;b1=65;\n", + "\n", + "base=b1-a1;\n", + "o1=58;o2=42;\n", + "mo2=(o2+o1)/2;\n", + "\n", + "ae2=base*mo2;\n", + "print('area=',ae2);\n", + "\n", + "\n", + "print('chainage 65 and 110')\n", + "a1=65;b1=110;\n", + "\n", + "base=b1-a1;\n", + "o1=0;o2=58;\n", + "mo3=(o2+o1)/2;\n", + "\n", + "ae3=base*mo3;\n", + "print('area=',ae3);\n", + "\n", + "\n", + "print('chainage 90 and 110')\n", + "a1=90;b1=110;\n", + "\n", + "base=b1-a1;\n", + "o1=0;o2=60;\n", + "mo4=(o2+o1)/2;\n", + "\n", + "ae4=base*mo4;\n", + "print('area=',ae4);\n", + "\n", + "print('chainage 40 and 90')\n", + "\n", + "a1=40;b1=90;\n", + "\n", + "base=b1-a1;\n", + "o1=60;o2=20;\n", + "mo5=(o2+o1)/2;\n", + "\n", + "ae5=base*mo5;\n", + "print('area=',ae5);\n", + "\n", + "print('chainage 0 and 40')\n", + "a1=0;b1=40;\n", + "\n", + "base=b1-a1;\n", + "o1=20;o2=0;\n", + "mo6=(o2+o1)/2;\n", + "\n", + "ae6=base*mo6\n", + "print('area=',ae6);\n", + "\n", + "\n", + "area=ae1+ae2+ae3+ae4+ae5+ae6;\n", + "\n", + "print('area of field =',area);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### ch-7 page 209,210 pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "chainage 15.5 and 27.5\n", + "('area=', 135.0)\n", + "chainage 15.5 and 50\n", + "('area=', 905.625)\n", + "chainage 50 and 75.5\n", + "('area=', 835.125)\n", + "chainage 75.5 and 86.7\n", + "('area=', 198.80000000000004)\n", + "chainage 86.7 and 90\n", + "('area=', 17.324999999999985)\n", + "chainage 60 and 90\n", + "('area=', 532.5)\n", + "chainage 35.5 and 60\n", + "('area=', 490.0)\n", + "chainage 27.5 and 35.5\n", + "('area=', 60.0)\n", + "('ap,ae=', 3022.05, 152.325)\n", + "('total area of field =', 2869.7250000000004)\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "\n", + "print('chainage 15.5 and 27.5')\n", + "a1=15.5;b1=27.5;\n", + "\n", + "base=b1-a1;\n", + "o1=0;o2=22.5;\n", + "mo1=(o2+o1)/2;\n", + "\n", + "ae1=base*mo1;\n", + "ap1=0;\n", + "an1=ae1;\n", + "print('area=',ae1);\n", + "\n", + "print('chainage 15.5 and 50')\n", + "a1=15.5;b1=50;\n", + "\n", + "base=b1-a1;\n", + "o1=22.5;o2=30;\n", + "mo2=(o2+o1)/2;\n", + "\n", + "ae2=base*mo2;\n", + "ap2=ae2;\n", + "an2=0;\n", + "print('area=',ae2);\n", + "\n", + "\n", + "print('chainage 50 and 75.5')\n", + "a1=50;b1=75.5;\n", + "\n", + "base=b1-a1;\n", + "o1=30;o2=35.5;\n", + "mo3=(o2+o1)/2;\n", + "\n", + "ae3=base*mo3;\n", + "ap3=ae3;\n", + "an3=0;\n", + "print('area=',ae3);\n", + "\n", + "\n", + "print('chainage 75.5 and 86.7')\n", + "a1=75.5;b1=86.7;\n", + "\n", + "base=b1-a1;\n", + "o1=35.5;o2=0;\n", + "mo4=(o2+o1)/2;\n", + "\n", + "ae4=base*mo4;\n", + "ap4=ae4;\n", + "an4=0;\n", + "print('area=',ae4);\n", + "\n", + "print('chainage 86.7 and 90')\n", + "\n", + "a1=86.7;b1=90;\n", + "\n", + "base=b1-a1;\n", + "o1=0;o2=10.5;\n", + "mo5=(o2+o1)/2;\n", + "\n", + "ae5=base*mo5;\n", + "ap5=0;\n", + "an5=ae5;\n", + "print('area=',ae5);\n", + "\n", + "print('chainage 60 and 90')\n", + "a1=60;b1=90;\n", + "\n", + "base=b1-a1;\n", + "o1=10.5;o2=25.0;\n", + "mo6=(o2+o1)/2;\n", + "\n", + "ae6=base*mo6\n", + "ap6=ae6;\n", + "an6=0;\n", + "print('area=',ae6);\n", + "\n", + "print('chainage 35.5 and 60')\n", + "a1=35.5;b1=60;\n", + "\n", + "base=b1-a1;\n", + "o1=25;o2=15;\n", + "mo7=(o2+o1)/2;\n", + "\n", + "ae7=base*mo7\n", + "ap7=ae7;\n", + "an7=0;\n", + "print('area=',ae7);\n", + "\n", + "print('chainage 27.5 and 35.5')\n", + "a1=27.5;b1=35.5;\n", + "\n", + "base=b1-a1;\n", + "o1=15;o2=0;\n", + "mo8=(o2+o1)/2;\n", + "\n", + "ae8=base*mo8\n", + "ap8=ae8;\n", + "an8=0\n", + "print('area=',ae8);\n", + "\n", + "an=an1+an2+an3+an4+an5+an6+an7+an8;\n", + "ap=ap1+ap2+ap3+ap4+ap5+ap6+ap7+ap8;\n", + "\n", + "area=ap-an;\n", + "print('ap,ae=',ap,an)\n", + "print('total area of field =',area);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-7 page 214 pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Mid ordinate rule\n", + "('required area is', 188.0, 'square meters')\n", + "average ordinate rule\n", + "('required area is', 161.14285714285714, 'sqare meters')\n", + "trapezoidal rule\n", + "('required area is ', 188.0, 'square meters')\n", + "simpsons rule\n", + "('required area is ', 196.66666666666669, 'square meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "\n", + "dis=10;\n", + "a=0;g=0;\n", + "b=2.5;c=3.5;d=5;e=4.6;f=3.2;\n", + "\n", + "print('Mid ordinate rule');\n", + "\n", + "h1=(a+b)/2;\n", + "h2=(b+c)/2;\n", + "h3=(c+d)/2;\n", + "h4=(d+e)/2;\n", + "h5=(e+f)/2;\n", + "h6=(f+g)/2;\n", + "area=dis*(h1+h2+h3+h4+h5+h6);\n", + "\n", + "print('required area is',area,'square meters');\n", + "\n", + "print('average ordinate rule');\n", + "dis=10;\n", + "p=6;\n", + "bl=dis*p;\n", + "no=7;\n", + "\n", + "\n", + "area2=bl*(a+b+c+d+e+f+g)/no;\n", + "\n", + "print('required area is',area2,'sqare meters');\n", + "\n", + "print('trapezoidal rule');\n", + "\n", + "\n", + "area3=(dis/2)*(2*(a+b+c+d+e+f+g));\n", + "\n", + "print('required area is ',area3,'square meters');\n", + "print('simpsons rule');\n", + "\n", + "area4=(dis/3)*(4*(b+d+f)+2*(c+e));\n", + "print('required area is ',area4,'square meters');\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-7 page 216 pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "trapezoidal rule\n", + "('required area is ', 820.125, 'square meters')\n", + "simpsons rule\n", + "(756.0, 57.375)\n", + "('required area is ', 813.375, 'square meters')\n" + ] + } + ], + "source": [ + "#ch-7 page 216 pb-2\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "print('trapezoidal rule');\n", + "\n", + "o1=3.5;o2=4.3;o3=6.75;o4=5.25;o5=7.5;o6=8.8;o7=7.9;\n", + "o8=6.4;o9=4.4;o10=3.25;\n", + "\n", + "dis=15;\n", + "\n", + "area1=(dis/2)*(o1+o10+(2*(o2+o3+o4+o5+o6+o7+o8+o9)));\n", + "\n", + "print('required area is ',area1,'square meters');\n", + "\n", + "print('simpsons rule')\n", + "\n", + "A1=dis/3*(o1+o9+4*(o2+o4+o6+o8)+2*(o3+o5+o7));\n", + "\n", + "A2=dis/2*(o10+o9);\n", + "\n", + "area2=A1+A2;\n", + "print(A1,A2)\n", + "\n", + "print('required area is ',area2,'square meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### cha 7 page -216 pb-3" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "trapezoidal rule\n", + "(89.5, 106.49999999999999, 158.0)\n", + "('total area=', 354.0, 'meters')\n", + "simpsons rule\n", + "(89.66666666666667, 102.33333333333333, 157.33333333333334)\n", + "('total area is ', 349.33333333333337, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "o1=2.5;o2=3.8;o3=4.6;o4=5.2;o5=6.1;o6=4.7;o7=5.8;o8=3.9;o9=2.20;\n", + "\n", + "d1=5;\n", + "d2=10;\n", + "d3=20;\n", + "\n", + "\n", + "print('trapezoidal rule')\n", + "\n", + "del1=(d1/2)*(o1+o5+2*(o2+o3+o4));\n", + "del2=(d2/2)*(o5+o7+2*(o6));\n", + "del3=(d3/2)*(o7+o9+2*(o8));\n", + "\n", + "total1=del1+del2+del3;\n", + "print(del1,del2,del3)\n", + "\n", + "print('total area=',total1,'meters');\n", + "\n", + "print('simpsons rule')\n", + "\n", + "de1=(d1/3)*(o1+o5+4*(o2+o4)+2*(o3));\n", + "de2=(d2/3)*(o5+o7+4*(o6));\n", + "de3=(d3/3)*(o7+o9+4*(o8));\n", + "\n", + "\n", + "total2=de1+de2+de3;\n", + "print(de1,de2,de3)\n", + "\n", + "print('total area is ',total2,'meters')\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha 7 page -225 pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('A=', 748.0)\n", + "('required area is', 748.0, 'meters')\n" + ] + } + ], + "source": [ + "#cha 7 page -225 ;pb-1\n", + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "ir=9.377;\n", + "fr=3.336;\n", + "m=100;\n", + "c=23.521;\n", + "\n", + "n=1;\n", + "\n", + "a1=m*(fr-ir+10*(n)+c);\n", + "\n", + "a2=m*(fr-ir-10*(n)+c);\n", + "\n", + "print('A=',a2);\n", + "print('required area is',a2,'meters');\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha 7 page -225,226 pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('A=', 81460.00000000001)\n", + "('required area is', 81460.00000000001, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "ir=8.652;\n", + "fr=6.798;\n", + "c=20;\n", + "m=100;\n", + "n=1;\n", + "\n", + "sc=100;\n", + "\n", + "a2=m*(fr-ir-10*(n)+c);\n", + "\n", + "a2=a2*sc;\n", + "\n", + "print('A=',a2);\n", + "print('required area is',a2,'meters');\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha 7 page -226 pb-3" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('required area is', 9747.499999999998, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "ir=4.855;\n", + "fr=8.754;\n", + "m=100;\n", + "\n", + "n=0;\n", + "c=0;\n", + "sc=25\n", + "a=m*(fr-ir)\n", + "a=a*sc;\n", + "print('required area is',a,'meters');\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-7 page-226 pb-4" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "case 1\n", + "('A=', 100.0)\n", + "('M=', 100.0)\n", + "case 2\n", + "('required area is', 1347.0)\n", + "('area of zero circle is', 2122.0)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "print('case 1')\n", + "\n", + "ir=3.415;\n", + "fr=4.415;\n", + "n=0;\n", + "c=0;\n", + "sc=16; #1cm^2=16m^2;\n", + "h=10000;\n", + "ag=0.16*h;\n", + "\n", + "am=ag/sc;\n", + "print('A=',am);\n", + "\n", + "m=am/(fr-ir);\n", + "\n", + "print('M=',m);\n", + "\n", + "print('case 2');\n", + "\n", + "fr_ir=2.25;\n", + "c=21.22;\n", + "n=1\n", + "\n", + "a1=m*(fr_ir-10+c);\n", + "print('required area is',a1);\n", + "\n", + "area=m*c;\n", + "\n", + "print('area of zero circle is',area);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "\n", + "### cha-7 page-227 pb-5" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('area of zero circle is', 1221.0, 'square centimeters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "l=10;b=15;\n", + "a1=l*b;\n", + "\n", + "ir=0.686;\n", + "fr=9.976;\n", + "n=2;\n", + "m=100;\n", + "\n", + "marea=150;\n", + "\n", + "c=(marea/100)+10.710;\n", + "\n", + "area=m*c;\n", + "print('area of zero circle is',area,'square centimeters');\n", + "\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-7 page-228 pb-6" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "case 1\n", + "('area of figure is', 705.0, 'square cm')\n", + "case 2\n", + "('area of figure is', -1357.0, 'sq cm')\n", + "('C=', 20.62)\n", + "('area of zero circle is', 2062.0, 'square cm')\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "print('case 1')\n", + "n=1;\n", + "c=0;\n", + "m=100;\n", + "fr=4.825;\n", + "ir=7.775;\n", + "area1=m*(fr-ir+10*n)\n", + "\n", + "print('area of figure is',area1,'square cm');\n", + "\n", + "print('case 2')\n", + "fr=8.755;\n", + "ir=2.325;\n", + "m=100;\n", + "n=2;\n", + "area2=m*(fr-ir-10*n+c)\n", + "\n", + "print('area of figure is',area2,'sq cm')\n", + "c=(area1/m)+13.570;\n", + "print('C=',c)\n", + "\n", + "areac=m*c;\n", + "print('area of zero circle is',areac,'square cm');\n", + "\n" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 2", + "language": "python", + "name": "python2" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 2 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython2", + "version": "2.7.11" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap8_Computation-of-Volume.ipynb b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap8_Computation-of-Volume.ipynb new file mode 100644 index 00000000..197ff573 --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap8_Computation-of-Volume.ipynb @@ -0,0 +1,713 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "# Chapter 8: Computation of volume" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-8 page-241 pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "(10.215, 14.84375, 28.433749999999996, 34.375, 23.633750000000003, 16.23375, 9.58375)\n", + "by trapezoidal rule\n", + "('V=', 5096.775, 'meter cube')\n", + "by prismoidal rule\n", + "('V=', 5143.25, 'meter cube')\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "h1=0.90;h2=1.25;h3=2.15;h4=2.50;h5=1.85;h6=1.35;h7=0.85;\n", + "\n", + "b=10;\n", + "sh=1.5;\n", + "\n", + "h=40;\n", + "\n", + "d1=(b+(sh*h1))*h1;\n", + "d2=(b+(sh*h2))*h2;\n", + "d3=(b+(sh*h3))*h3;\n", + "d4=(b+(sh*h4))*h4;\n", + "d5=(b+(sh*h5))*h5;\n", + "d6=(b+(sh*h6))*h6;\n", + "d7=(b+(sh*h7))*h7;\n", + "\n", + "print(d1,d2,d3,d4,d5,d6,d7)\n", + "print('by trapezoidal rule');\n", + "v=(h/2)*(d1+d7+2*(d2+d3+d4+d5+d6));\n", + "print('V=',v,'meter cube');\n", + "\n", + "print('by prismoidal rule');\n", + "\n", + "v1=(h/3)*(d1+d7+4*(d2+d4+d6)+2*(d3+d5));\n", + "\n", + "print('V=',v1,'meter cube');\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-8 page-241,242 pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "(19.736842105263158, 29.487179487179485, 19.047619047619047, 23.214285714285715)\n", + "(7.125, 11.845, 7.68, 13.78, 16.5, 7.125)\n", + "from chainage 0 to 50\n", + "(70.3125, 179.3165955128205)\n", + "from chainage 50 to 100\n", + "(78.74925714285715, 174.63782051282053)\n", + "from chainage 100 to 150\n", + "(73.14285714285714, 213.2455)\n", + "from chainage 150 to 200\n", + "(221.13535714285715, 159.94642857142858)\n", + "from chainage 200 to 250\n", + "590.625\n", + "('total cutting =', 1033.9649714285715)\n", + "('total fitting=', 727.1463445970695)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "h1=0.75;h2=1.15;h3=0.80;h4=1.30;h5=1.5;h6=0.75\n", + "b=8;sh=2;\n", + "\n", + "x1=(50*h1)/(h1+h2);\n", + "x2=(50*h2)/(h2+h3);\n", + "x3=(50*h3)/(h4+h3);\n", + "x4=(50*h4)/(h4+h5);\n", + "print(x1,x2,x3,x4);\n", + "\n", + "a1=(b+(sh*h1))*h1;\n", + "a2=(b+(sh*h2))*h2;\n", + "a3=(b+(sh*h3))*h3;\n", + "a4=(b+(sh*h4))*h4;\n", + "a5=(b+(sh*h5))*h5;\n", + "a6=(b+(sh*h6))*h6;\n", + "print(a1,a2,a3,a4,a5,a6)\n", + "\n", + "print('from chainage 0 to 50');\n", + "c1=(a1/2)*(x1);\n", + "f1=(a2/2)*(x2+0.79);\n", + "print(c1,f1);\n", + "\n", + "\n", + "\n", + "print('from chainage 50 to 100');\n", + "f2=(a2/2)*(x2);\n", + "c2=(a3/2)*(x3+1.46);\n", + "print(c2,f2);\n", + "\n", + "print('from chainage 100 to 150');\n", + "c3=(a3/2)*(x3);\n", + "f3=(a4/2)*(30.95);\n", + "print(c3,f3);\n", + "\n", + "print('from chainage 150 to 200');\n", + "f4=(a4/2)*(x4);\n", + "c4=(a5/2)*(x4+3.59);\n", + "print(c4,f4);\n", + "\n", + "print('from chainage 200 to 250');\n", + "c5=((a1+a5)/2)*50;\n", + "\n", + "print(c5);\n", + "\n", + "tc=c1+c2+c3+c4+c5;\n", + "tf=f1+f2+f3+f4;\n", + "\n", + "print('total cutting =',tc);\n", + "print('total fitting=',tf);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-8 page-244 pb-3" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "(26.25, 11.5625, 8.5025, 15.9225, 22.44, 30.9225, 6.5625)\n", + "('volume=', 5472.125)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "h=50;\n", + "h1=2.50;h2=1.25;h3=0.95;h4=1.65;h5=2.20;h6=2.85;h7=0.75;\n", + "b=8;sh=1;\n", + "\n", + "\n", + "a1=(b+(sh*h1))*h1;\n", + "a2=(b+(sh*h2))*h2;\n", + "a3=(b+(sh*h3))*h3;\n", + "a4=(b+(sh*h4))*h4;\n", + "a5=(b+(sh*h5))*h5;\n", + "a6=(b+(sh*h6))*h6;\n", + "a7=(b+(sh*h7))*h7;\n", + "\n", + "print(a1,a2,a3,a4,a5,a6,a7);\n", + "\n", + "v=(h/3)*(a1+a7+4*(a2+a4+a6)+2*(a3+a5));\n", + "\n", + "\n", + "print('volume=',v);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-8 page-245 pb-4" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "according to trapezoidal rule\n", + "('volume =', 330375.0)\n", + "according to prismoidal rule\n", + "('volume =', 330250.0)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "a1=2050;a2=8400;a3=16300;a4=24600;a5=31500;\n", + "\n", + "h=5;\n", + "\n", + "print('according to trapezoidal rule')\n", + "\n", + "v1=(h/2)*(a1+a5+2*(a2+a3+a4));\n", + "\n", + "print('volume =',v1);\n", + "\n", + "print('according to prismoidal rule')\n", + "\n", + "v2=(h/3)*(a1+a5+4*(a2+a4)+2*(a3));\n", + "\n", + "print('volume =',v2)\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-8 page-245,246 pb-5" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "bottom section\n", + "('area A1=', 1200)\n", + "mid section\n", + "('area A2=', 2000.0)\n", + "top section\n", + "('area A3=', 30000)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "print('bottom section')\n", + "L=40;\n", + "B=30;\n", + "a1=L*B;\n", + "print('area A1=',a1)\n", + "\n", + "print('mid section')\n", + "b=40;\n", + "sh=2.5;\n", + "\n", + "l=L+2*2*sh;\n", + "b=B+2*2*sh;\n", + "a2=l*b;\n", + "print('area A2=',a2);\n", + "\n", + "print('top section')\n", + "sh=5;\n", + "\n", + "l1=L+2*sh;\n", + "b1=B*2*2*sh;\n", + "a3=l1*b1;\n", + "print('area A3=',a3)\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-8 page-246,247 pb-6" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "first section\n", + "(6.666666666666666, 5.454545454545455)\n", + "('area A1=', 20.36363636363636)\n", + "second section\n", + "(7.777777777777778, 6.363636363636363)\n", + "('area A2=', 33.494949494949495)\n", + "third section\n", + "(8.88888888888889, 7.2727272727272725)\n", + "('area A3=', 48.64646464646465)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "b=8;\n", + "h=2;\n", + "n=10;\n", + "s=1;\n", + "\n", + "print('first section');\n", + "b1=(b/2)+((n*s)/(n-s))*(h+(b/(2*n)));\n", + "b2=(b/2)+((n*s)/(n+s))*(h-(b/(2*n)));\n", + "\n", + "a1=0.5*((((b/(2*s))+h))*(b1+b2)-((b*b)/(2*s)));\n", + "print(b1,b2)\n", + "\n", + "print('area A1=',a1);\n", + "\n", + "print('second section');\n", + "b=8;h=3;n=10;s=1;\n", + "\n", + "\n", + "b1=(b/2)+((n*s)/(n-s))*(h+(b/(2*n)));\n", + "b2=(b/2)+((n*s)/(n+s))*(h-(b/(2*n)));\n", + "\n", + "a2=0.5*((((b/(2*s))+h))*(b1+b2)-((b*b)/(2*s)));\n", + "print(b1,b2)\n", + "\n", + "print('area A2=',a2);\n", + "\n", + "print('third section');\n", + "b=8;h=4;n=10;s=1;\n", + "\n", + "\n", + "b1=(b/2)+((n*s)/(n-s))*(h+(b/(2*n)));\n", + "b2=(b/2)+((n*s)/(n+s))*(h-(b/(2*n)));\n", + "\n", + "a3=0.5*((((b/(2*s))+h))*(b1+b2)-((b*b)/(2*s)));\n", + "print(b1,b2)\n", + "\n", + "print('area A3=',a3);\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-8 page-247,248 pb-7" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "first section\n", + "(6.666666666666667, 5.454545454545454)\n", + "('area A1=', 11.36363636363636)\n", + "second section\n", + "(8.75, 5.833333333333333)\n", + "('area A2=', 26.041666666666664)\n", + "third section\n", + "(7.428571428571429, 5.777777777777778)\n", + "('area A3=', 17.920634920634917)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "b=10;\n", + "h=1;\n", + "n=10;\n", + "s=1;\n", + "\n", + "print('first section');\n", + "b1=(b/2)+((n*s)/(n-s))*(h+(b/(2*n)));\n", + "b2=(b/2)+((n*s)/(n+s))*(h-(b/(2*n)));\n", + "\n", + "a1=0.5*((((b/(2*s))+h))*(b1+b2)-((b*b)/(2*s)));\n", + "print(b1,b2)\n", + "\n", + "print('area A1=',a1);\n", + "\n", + "print('second section');\n", + "b=10;h=2;n=5;s=1;\n", + "\n", + "\n", + "b1=(b/2)+((n*s)/(n-s))*(h+(b/(2*n)));\n", + "b2=(b/2)+((n*s)/(n+s))*(h-(b/(2*n)));\n", + "\n", + "a2=0.5*((((b/(2*s))+h))*(b1+b2)-((b*b)/(2*s)));\n", + "print(b1,b2)\n", + "\n", + "print('area A2=',a2);\n", + "\n", + "print('third section');\n", + "b=10;h=1.5;n=8;s=1;\n", + "\n", + "\n", + "b1=(b/2)+((n*s)/(n-s))*(h+(b/(2*n)));\n", + "b2=(b/2)+((n*s)/(n+s))*(h-(b/(2*n)));\n", + "\n", + "a3=0.5*((((b/(2*s))+h))*(b1+b2)-((b*b)/(2*s)));\n", + "print(b1,b2)\n", + "\n", + "print('area A3=',a3);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-8 page-248 pb-8" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "at station 1\n", + "('area=', 14.25)\n", + "at station 2\n", + "('area=', 18.975)\n", + "('v=', 830.625, 'cp=', 0.20833333333333262)\n", + "('correct volume =', 830.4166666666666)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "print('at station 1');\n", + "h=1;h1=2.55;h2=0.95;b=9;b1=7.5;b2=5.25;\n", + "w1=b1+b2\n", + "a=(((h/2)*(b1+b2))+((b/4)*(h1+h2)));\n", + "print('area=',a)\n", + "\n", + "print('at station 2');\n", + "h=1.5;h1=2.8;h2=1.35;b=9;b1=8.1;b2=4.75;\n", + "\n", + "a1=(((h/2)*(b1+b2))+((b/4)*(h1+h2)));\n", + "d=50;\n", + "k=10.01;\n", + "v=(d/2)*(a+a1);\n", + "w2=b1+b2\n", + "print('area=',a1)\n", + "h2=1;\n", + "h1=1.5;\n", + "cp=(d/12)*(h1-h2)*(w2-w1);\n", + "\n", + "\n", + "cv=v-cp;\n", + "print('v=',v,'cp=',cp)\n", + "print('correct volume =',cv);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-8 page-249 pb-9" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "section 1\n", + "(7.03125, 1.0416666666666667)\n", + "section 2\n", + "(9.03125, 0.375)\n", + "('vc=', 401.5625, 'vf=', 35.41666666666667)\n", + "('corrected volume (in cutting)=', 400.5208333333333)\n", + "('corrected volume(in filling)', 34.027777777777786)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "print('section 1')\n", + "b=10;n=5;s=1;s1=2;\n", + "d=50;h1=0.5;h2=0.7;\n", + "\n", + "ac=0.5*(((0.5*b+n*h1)*(0.5*b+n*h1))/(n-s));\n", + "\n", + "af=0.5*(((0.5*b-n*h1)*(0.5*b-n*h1))/(n-s1));\n", + "\n", + "print(ac,af)\n", + "\n", + "\n", + "print('section 2')\n", + "\n", + "\n", + "ac1=0.5*(((0.5*b+n*h2)*(0.5*b+n*h2))/(n-s));\n", + "\n", + "af1=0.5*(((0.5*b-n*h2)*(0.5*b-n*h2))/(n-s1));\n", + "D=50;\n", + "print(ac1,af1)\n", + "vc=((ac+ac1)/2)*D;\n", + "vf=((af+af1)/2)*D;\n", + "\n", + "print('vc=',vc,'vf=',vf);\n", + "\n", + "D=50;\n", + "pcc=(D/(12*(n-s)))*(n*n*(h1-h2)*(h1-h2));\n", + "\n", + "\n", + "pcf=(D/(12*(n-s1)))*(n*n*(h1-h2)*(h1-h2));\n", + "\n", + "\n", + "cvc=vc-pcc;\n", + "cvf=vf-pcf;\n", + "\n", + "print('corrected volume (in cutting)=',cvc);\n", + "\n", + "print('corrected volume(in filling)',cvf)\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### cha-8 page-251 pb-10" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "at station 1\n", + "(40.2, 82.475)\n", + "('area =', 21.137499999999996)\n", + "at station 2\n", + "(53.7, 105.675)\n", + "('area =', 25.987499999999997)\n", + "volume by average end area rule\n", + "('volume=', 1178.1249999999998)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "a1=0;a2=3.0;\n", + "b1=2.20;b2=5.50;\n", + "c1=1.75;c2=3.0;\n", + "d1=1.5;d2=0;\n", + "e1=4.75;e2=5.25;\n", + "f1=6.40;f2=7.30;\n", + "g1=0;g2=3.0;\n", + "\n", + "print('at station 1')\n", + "sp=(e1*d2)+(f1*e2)+(d2*f2)+(c1*d2)+(b1*c2)+(a1*b2);\n", + "\n", + "sq=(e2*d1)+(e1*f2)+(f1*g2)+(d1*c2)+(c1*b2)+(b1*a2);\n", + "\n", + "area1=0.5*(sp-sq)\n", + "area1=abs(area1);\n", + "print(sp,sq)\n", + "print('area =',area1)\n", + "\n", + "a1=0;a2=3.0;\n", + "b1=3.1;b2=5.25;\n", + "c1=2.20;c2=3.0;\n", + "d1=2;d2=0;\n", + "e1=5.25;e2=6;\n", + "f1=7.40;f2=8.5;\n", + "g1=0;g2=3.0;\n", + "print('at station 2')\n", + "sp1=(e1*d2)+(f1*e2)+(d2*f2)+(c1*d2)+(b1*c2)+(a1*b2);\n", + "\n", + "sq1=(d1*e2)+(e1*f2)+(f1*g2)+(d1*c2)+(c1*b2)+(b1*a2);\n", + "print(sp1,sq1)\n", + "\n", + "\n", + "area2=0.5*(sp1-sq1)\n", + "area2=abs(area2);\n", + "print('area =',area2)\n", + "\n", + "print('volume by average end area rule')\n", + "v=50*((area1+area2)/2);\n", + "print('volume=',v)\n" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 2", + "language": "python", + "name": "python2" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 2 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython2", + "version": "2.7.11" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap9_Theodolite-Traversing.ipynb b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap9_Theodolite-Traversing.ipynb new file mode 100644 index 00000000..3153bd76 --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chap9_Theodolite-Traversing.ipynb @@ -0,0 +1,502 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "# Chapter 9: Theodolite Traversing" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-9 page 302 pb-1" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "(65.11978202514794, -63.50741753704864, -51.91315691660193)\n", + "(38.20554919114786, 168.9587462008847, -30.579186368382416)\n", + "(50.300792428502625, -176.58510902365015)\n", + "('distance DA=', 183.60955979422673)\n", + "('bearing of DA=', 74.10023981818601)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "l1=75.5;\n", + "l2=180.5\n", + "l3=60.25\n", + "\n", + "t1=30.4;t2=69.4;t3=30.5;\n", + "t2=180-t2;\n", + "t3=180-t3;\n", + "\n", + "Lc1=l1*math.cos(t1*(math.pi/180))\n", + "Lc2=l2*math.cos(t2*(math.pi/180))\n", + "Lc3=l3*math.cos(t3*(math.pi/180))\n", + "\n", + "Ls1=l1*math.sin(t1*(math.pi/180))\n", + "Ls2=l2*math.sin(t2*(math.pi/180))\n", + "Ls3=-l3*math.sin(t3*(math.pi/180))\n", + "\n", + "print(Lc1,Lc2,Lc3);\n", + "print(Ls1,Ls2,Ls3);\n", + "Lc4=-Lc1-Lc2-Lc3;\n", + "Ls4=-Ls1-Ls2-Ls3;\n", + "\n", + "print(Lc4,Ls4);\n", + "\n", + "t4=-math.atan(Ls4/Lc4);\n", + "t4=t4*(180/math.pi);\n", + "\n", + "l4=math.sqrt(Lc4*Lc4+Ls4*Ls4);\n", + "\n", + "print('distance DA=',l4);\n", + "print('bearing of DA=',t4);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-9 page 304 pb-2" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('length of DA is', 145.80634036039953, 'or', 12.053659639600497)\n", + "when length of DA ,L=145.8\n", + "('bearing at AB is=N', 82.44640641462031)\n", + "when length of DA ,L=12.04\n", + "0.999661660714\n", + "('bearing at AB is=N', 1.4904797844587976)\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "l1=100;\n", + "l2=80;\n", + "l3=60;\n", + "\n", + "t2=39.5;t3=40.5;t4=49.75;\n", + "\n", + "L2=l2*math.cos(t2*(math.pi/180));\n", + "L3=l3*math.cos(t3*(math.pi/180));\n", + "\n", + "D2=l2*math.sin(t2*(math.pi/180));\n", + "D3=l3*math.sin(t3*(math.pi/180));\n", + "\n", + "l41=(157.86+math.sqrt(157.86*157.86-4*1757.5))/2;\n", + "l42=(157.86-math.sqrt(157.86*157.86-4*1757.5))/2;\n", + "\n", + "print('length of DA is',l41,'or',l42);\n", + "\n", + "print('when length of DA ,L=145.8')\n", + "\n", + "k=math.cos(t4*(math.pi/180))\n", + "k1=(L2+L3-(k*l41))/100;\n", + "t1=math.acos(k1);\n", + "t1=t1*(180/(math.pi))\n", + "print('bearing at AB is=N',t1)\n", + "\n", + "\n", + "print('when length of DA ,L=12.04')\n", + "\n", + "k=math.cos(t4*(math.pi/180))\n", + "k1=(L2+L3-(k*l42))/100;\n", + "k1=k1+0.004;\n", + "t11=math.acos(k1);\n", + "t11=t11*(180/(math.pi))\n", + "print(k1)\n", + "print('bearing at AB is=N',t11)\n", + "\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-9 page 305 pb-3" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('latitude of AB,CD,DE are', 86.59373062437335, -57.03044742000232, -25.250000000000007)\n", + "('Depature of AB,CD,DE are', 51.00760547755076, -48.708603624763775, -43.73428289111415)\n", + "(-4.313283204371025, 41.43528103832716)\n", + "('length of BC=', 348.51410778926174)\n", + "('length of EA=', 317.28203276885586)\n" + ] + } + ], + "source": [ + "#ch-9 page 305 pb-3\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "l1=100.5;l3=75;l4=50.5;\n", + "t1=30.5;t2=45;t3=40.5;t4=60;t5=40.25;\n", + "\n", + "\n", + "L1=l1*math.cos(t1*(math.pi/180))\n", + "L3=-l3*math.cos(t3*(math.pi/180))\n", + "L4=-l4*math.cos(t4*(math.pi/180))\n", + "\n", + "print('latitude of AB,CD,DE are',L1,L3,L4);\n", + "D1=l1*math.sin(t1*(math.pi/180))\n", + "D3=-l3*math.sin(t3*(math.pi/180))\n", + "D4=-l4*math.sin(t4*(math.pi/180))\n", + "print('Depature of AB,CD,DE are',D1,D3,D4);\n", + "\n", + "L2_L5=-(L1+L3+L4);\n", + "D2_D5=-(D1+D3+D4);\n", + "print(L2_L5,D2_D5)\n", + "\n", + "k=0.117;\n", + "l5=(L2_L5+D2_D5)/(k);\n", + "\n", + "k1=0.763;\n", + "\n", + "l2=(k1*l5)-L2_L5;\n", + "l2=l2/0.707;\n", + "\n", + "print('length of BC=',l2);\n", + "print('length of EA=',l5);\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-9 page 307 pb-4" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('latitudes of AQ,QR,RB are', 65.21958064293187, -61.02257873940248, -36.93176700776003)\n", + "('Depature of AQ,QR,RB are', 38.03493526693723, 52.118205878497236, -65.27667719549248)\n", + "('length of AB=', 41.114514530539196, 'meters')\n" + ] + } + ], + "source": [ + "#ch-9 page 307 pb-4\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "l1=75.5;l2=80.25;l3=75;\n", + "t1=30.25;t2=40.5;t3=60.5;\n", + "\n", + "\n", + "L1=l1*math.cos(t1*(math.pi/180))\n", + "L2=-l2*math.cos(t2*(math.pi/180))\n", + "L3=-l3*math.cos(t3*(math.pi/180))\n", + "print('latitudes of AQ,QR,RB are',L1,L2,L3);\n", + "\n", + "\n", + "D1=l1*math.sin(t1*(math.pi/180))\n", + "D2=l2*math.sin(t2*(math.pi/180))\n", + "D3=-l3*math.sin(t3*(math.pi/180))\n", + "print('Depature of AQ,QR,RB are',D1,D2,D3);\n", + "\n", + "L4=-(L1+L2+L3);\n", + "D4=-(D1+D2+D3);\n", + "\n", + "l4=math.sqrt(L4*L4+(D4*D4));\n", + "\n", + "print('length of AB=',l4,'meters');\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-9 page 308 pb-5" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('latitudes of BQ,QP,PA are', -96.23556979807799, -172.32583208830516, -61.552945012933385)\n", + "('Depature of BQ,QP,PA are', 115.71069572705566, -101.50767259214082, -108.79446199248746)\n", + "('length of AB=', 343.3992171422471, 'meters')\n", + "('bearing of AB=', 15.989201746570728)\n", + "('PAB=', 44.51079825342927, 'QBA=', 66.23920174657073)\n" + ] + } + ], + "source": [ + "#ch-9 page 308 pb-5\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "l1=150.5;l2=200;l3=125;\n", + "t1=50.25;t2=30.5;t3=60.5;\n", + "\n", + "\n", + "L1=-l1*math.cos(t1*(math.pi/180))\n", + "L2=-l2*math.cos(t2*(math.pi/180))\n", + "L3=-l3*math.cos(t3*(math.pi/180))\n", + "print('latitudes of BQ,QP,PA are',L1,L2,L3);\n", + "\n", + "\n", + "D1=l1*math.sin(t1*(math.pi/180))\n", + "D2=-l2*math.sin(t2*(math.pi/180))\n", + "D3=-l3*math.sin(t3*(math.pi/180))\n", + "print('Depature of BQ,QP,PA are',D1,D2,D3);\n", + "\n", + "L4=-(L1+L2+L3);\n", + "D4=-(D1+D2+D3);\n", + "\n", + "l4=math.sqrt(L4*L4+(D4*D4));\n", + "\n", + "print('length of AB=',l4,'meters');\n", + "\n", + "t4=math.atan(D4/L4);\n", + "t4=t4*(180/math.pi);\n", + "print('bearing of AB=',t4);\n", + "\n", + "PAB=t3-t4;\n", + "QBA=t1+t4;\n", + "\n", + "print('PAB=',PAB,'QBA=',QBA);\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-9 page 308 pb-6" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('latitudes of AB,BC,CD,DE are', 121.76738460229168, 106.68654829016975, -133.98333444865727, 19.805712703281312)\n", + "('Depature of AB,BC,CD,DE are', 45.52695956373076, 186.6627451151656, 78.9222154403895, 118.35427218446777)\n", + "('length of EA=', 444.4100422146986, 'meters')\n", + "('bearing of EA=', 75.09947760257306)\n", + "('bearing from F to C is =', 5.818201574554788)\n", + "('distance from F to C is =', 172.2028708809785)\n" + ] + } + ], + "source": [ + "#ch-9 page 308 pb-6\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "l1=130;l2=215;l3=155.5;l4=120;\n", + "t1=20.5;t2=60.25;t3=30.5;t4=80.5;\n", + "\n", + "\n", + "L1=l1*math.cos(t1*(math.pi/180))\n", + "L2=l2*math.cos(t2*(math.pi/180))\n", + "L3=-l3*math.cos(t3*(math.pi/180))\n", + "L4=l4*math.cos(t4*(math.pi/180))\n", + "print('latitudes of AB,BC,CD,DE are',L1,L2,L3,L4);\n", + "\n", + "\n", + "D1=l1*math.sin(t1*(math.pi/180))\n", + "D2=l2*math.sin(t2*(math.pi/180))\n", + "D3=l3*math.sin(t3*(math.pi/180))\n", + "D4=l4*math.sin(t4*(math.pi/180))\n", + "print('Depature of AB,BC,CD,DE are',D1,D2,D3,D4);\n", + "\n", + "L5=-(L1+L2+L3+L4);\n", + "D5=-(D1+D2+D3+D4);\n", + "\n", + "l5=math.sqrt(L5*L5+(D5*D5));\n", + "\n", + "print('length of EA=',l5,'meters');\n", + "\n", + "t5=math.atan(D5/L5);\n", + "t5=t5*(180/math.pi);\n", + "print('bearing of EA=',t5);\n", + "\n", + "FA=l5/2;\n", + "l6=FA;\n", + "t6=t5;\n", + "L6=-l6*math.cos(t6*(math.pi/180))\n", + "D6=-l6*math.sin(t6*(math.pi/180))\n", + "\n", + "L7=-(L1+L2+L6)\n", + "D7=-(D1+D2+D6)\n", + "\n", + "t7=math.atan(D7/L7);\n", + "t7=t7*(180/math.pi);\n", + "print('bearing from F to C is =',t7);\n", + "\n", + "l7=math.sqrt(L7*L7+(D7*D7));\n", + "\n", + "print('distance from F to C is =',l7);\n", + "\n", + "\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### ch-9 page 308 pb-7" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('latitudes of AB,DE,EA are', 362.50000000000006, -538.0859250785912, -574.3378222288467)\n", + "('Depature of AB,DE,EA are', 627.8684177437179, -782.9198791909149, -27.587423899772766)\n", + "('t2-t3=', 100.07865810778766)\n", + "('Bearing of BC is', 63.0)\n", + "('Bearing of CD is', 37.12098009569709)\n" + ] + } + ], + "source": [ + "#ch-9 page 308 pb-7\n", + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "l1=725;l2=1050;l3=1250;l4=950;l5=575;\n", + "t1=60;t4=55.5;t5=2.75;\n", + "\n", + "\n", + "L1=l1*math.cos(t1*(math.pi/180))\n", + "L4=-l4*math.cos(t4*(math.pi/180))\n", + "L5=-l5*math.cos(t5*(math.pi/180))\n", + "print('latitudes of AB,DE,EA are',L1,L4,L5);\n", + "\n", + "\n", + "D1=l1*math.sin(t1*(math.pi/180))\n", + "D4=-l4*math.sin(t4*(math.pi/180))\n", + "D5=-l5*math.sin(t5*(math.pi/180))\n", + "print('Depature of AB,DE,EA are',D1,D4,D5);\n", + "\n", + "t2_t3=math.acos(0.1750);\n", + "t2_t3=180-(t2_t3*(180/math.pi));\n", + "\n", + "print('t2-t3=',t2_t3);\n", + "\n", + "t3=math.asin(0.6035);\n", + "t3=t3*(180/math.pi);\n", + "t2=t2_t3-t3;\n", + "t2=math.ceil(t2);\n", + "\n", + "print('Bearing of BC is',t2);\n", + "print('Bearing of CD is',t3);\n" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 2", + "language": "python", + "name": "python2" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 2 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython2", + "version": "2.7.11" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chapter11_Tacheometric-Traversing.ipynb b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chapter11_Tacheometric-Traversing.ipynb new file mode 100644 index 00000000..b203f9a0 --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/chapter11_Tacheometric-Traversing.ipynb @@ -0,0 +1,616 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "# Chapter 11: Tacheometric Surveying" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "\n", + "### section 11.7 , pg 413, problem 1" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "RL of instrument axis= 764.345 m\n", + "RL of D= 784.042 m\n", + "Distance of CD=147.097m\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "import math\n", + "\n", + "retiftoi=100\n", + "fplusd=0.15\n", + "s1=2.450-1.150\n", + "thetha1=5+(20/60)\n", + "v1=(100*1300*math.sin(10+(40/60))/2)+(0.15*math.sin(5+(20/60)));\n", + "s2=1.5\n", + "thetha2=8+(12/60)\n", + "V2=21.197\n", + "d2=147.097\n", + "RL=750.500+1.8+12.045 \n", + "RLD=RL+V2-1.5\n", + "print \"RL of instrument axis=\",RL,\"m\"\n", + "print \"RL of D=\", RLD,\"m\"\n", + "print \"Distance of CD=147.097m\" \n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### section 11.7, pg 415, problem 2" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "RL of axis when isnt. at P= 265.109\n", + "RL of A= 280.38\n", + "RL at B= 298.021\n", + "RL of B= 296.571\n", + "Distance between A and B= 118.009\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "v1=7.534\n", + "v2=16.871\n", + "v3=15.326\n", + "RLatp=255.750+v1+1.825\n", + "RLofA=265.109+v2-1.6\n", + "RLatB=280.380+v3+2.315\n", + "RLofB=298.021-1.450\n", + "D3=118.009\n", + "print \"RL of axis when isnt. at P=\", RLatp\n", + "print \"RL of A=\", RLofA\n", + "print \"RL at B=\", RLatB\n", + "print \"RL of B=\", RLofB\n", + "print \"Distance between A and B=\", D3\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### section 11.7 , pg 413, problem 1" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "10.494\n", + "RL of axis when isnt. at A= 462.449\n", + "RL of A= 461.104\n", + "RL at B= 487.151\n", + "RL of B= 485.601\n", + "RL of C 510.533\n" + ] + } + ], + "source": [ + "\n", + "from __future__ import division\n", + "import math\n", + "\n", + "\n", + "v1=10.494\n", + "d1=108.989\n", + "V2=24.807\n", + "d2=176.514\n", + "v3=25.652\n", + "d3=145.477\n", + "RL=450.500+1.455+v1 \n", + "RLofA=462.449-1.345\n", + "RLofB=462.449+24.807-1.655\n", + "RLatB=487.151\n", + "RLofC=RLofB+v3-2.250+1.53\n", + "print v1\n", + "print \"RL of axis when isnt. at A=\", RL\n", + "print \"RL of A=\", RLofA\n", + "print \"RL at B=\", RLatB\n", + "print \"RL of B=\", RLofB\n", + "print \"RL of C\", RLofC\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### ch-11 page 416 pb-4" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "in 1st observation\n", + "('v1,d1=', 9.386067902413853, 119.26130043570826)\n", + "in 2nd observation\n", + "('v2,d2=', 26.26555359446006, 145.25041419362984)\n", + "('RL of A=', 159.18106790241387)\n", + "('RL of B=', 175.81555359446008)\n", + "('difference of level AB=', 104.0330138511747, 'meters')\n", + "('gradient of AB is 1 in', 6.254056529136824)\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "c=100;\n", + "h=1.55;\n", + "rlo=150;\n", + "ra1=1.155;ra2=1.755;ra3=2.355;\n", + "rb1=1.250;rb2=2;rb3=2.750;\n", + "t1=30.5;t2=75.5;\n", + "a1=4.5;a2=10.25;\n", + "\n", + "print('in 1st observation')\n", + "v1=c*(ra3-ra1)*(math.sin(9*(math.pi/180)));\n", + "v1=v1/2;\n", + "d1=c*(ra3-ra1)*(math.cos(a1*(math.pi/180)))*(math.cos(a1*(math.pi/180)));\n", + "print('v1,d1=',v1,d1);\n", + "\n", + "print('in 2nd observation');\n", + "\n", + "v2=c*(rb3-rb1)*(math.sin(20.5*(math.pi/180)));\n", + "v2=v2/2;\n", + "d2=c*(rb3-rb1)*(math.cos(a2*(math.pi/180)))*(math.cos(a2*(math.pi/180)));\n", + "print('v2,d2=',v2,d2);\n", + "\n", + "rl=rlo+h;\n", + "rla=rl+v1-ra2;\n", + "rlb=rl+v2-rb2;\n", + "\n", + "print('RL of A=',rla);\n", + "print('RL of B=',rlb);\n", + "\n", + "t=t2-t1;\n", + "AB=math.sqrt((d1*d1+d2*d2)-2*(d1*d2*(math.cos(t*(math.pi/180)))));\n", + "print('difference of level AB=',AB,'meters');\n", + "\n", + "dab=rlb-rla;\n", + "gab=AB/dab;\n", + "print('gradient of AB is 1 in',gab);\n", + "\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### ch-11 page 418 pb-5" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('v1,v2=', 31.256671980047464, 31.1867536226639)\n", + "('h1,h2=', 2.0188558936750263, 1.5699268991777582)\n", + "('RL of A=', 418.7244721262775)\n", + "('RL of B=', 419.24331947815836)\n", + "('distance between A an B is', 323.2978586242886)\n", + "('gradient of PA is 1 in ', 5.567473732648181)\n", + "('gradient of PB is 1 in ', 4.68342893110529)\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "h=1.5;\n", + "a1=10;a2=12;\n", + "c=100;\n", + "ra1=1.150;ra2=2.050;ra3=2.950;\n", + "rb1=0.855;rb2=1.605;rb3=2.355;\n", + "rlp=450.5;\n", + "\n", + "\n", + "\n", + "v1=c*(ra3-ra1)*(math.sin(a1*(math.pi/180)));\n", + "\n", + "v2=c*(rb3-rb1)*(math.sin(a2*(math.pi/180)));\n", + "\n", + "h1=ra2*(math.cos(a1*(math.pi/180)));\n", + "h2=rb2*(math.cos(a2*(math.pi/180)));\n", + "\n", + "print('v1,v2=',v1,v2);\n", + "print('h1,h2=',h1,h2);\n", + "\n", + "rlai=rlp+h;\n", + "\n", + "rla=rlai-v1-h1;\n", + "rlb=rlai-v2-h2;\n", + "\n", + "print('RL of A=',rla);\n", + "print('RL of B=',rlb);\n", + "\n", + "d1=c*(ra3-ra1)*(math.cos(a1*(math.pi/180)))-ra2*(math.sin(a1*(math.pi/180)));\n", + "d2=c*(rb3-rb1)*(math.cos(a2*(math.pi/180)))-rb2*(math.sin(a2*(math.pi/180)));\n", + "\n", + "dab=d1+d2;\n", + "print('distance between A an B is',dab);\n", + "gpa=d1/(rlp-rla);\n", + "gpb=d2/(rlp-rlb);\n", + "\n", + "print('gradient of PA is 1 in ',gpa);\n", + "print('gradient of PB is 1 in ',gpb);\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### ch-11 page 419 pb-6\n" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "(96.98463103929541, 158.78462024097664, 117.67570175629913)\n", + "('latitudes of AB,BC,CD=', 83.56478621811925, -121.63607598835735, -83.20928669276485)\n", + "('depatures of AB,BC,CD ', 49.22342087003188, 102.06478649968226, -83.20928669276483)\n", + "(121.28057646300294, -68.07892067694931)\n", + "('Bearing of DA=', 29.30698225670086)\n", + "('length DA=', 139.08169422226874)\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "c=100;\n", + "ra1=1.25;ra2=1.75;ra3=2.25;\n", + "rb1=0.95;rb2=1.75;rb3=2.55;\n", + "rc1=1.55;rc2=2.15;rc3=2.75;\n", + "a1=10;a2=5;a3=8;\n", + "\n", + "ab=c*(ra3-ra1)*(math.cos(a1*(math.pi/180)))*(math.cos(a1*(math.pi/180)));\n", + "bc=c*(rb3-rb1)*(math.cos(a2*(math.pi/180)))*(math.cos(a2*(math.pi/180)));\n", + "cd=c*(rc3-rc1)*(math.cos(a3*(math.pi/180)))*(math.cos(a3*(math.pi/180)));\n", + "\n", + "print(ab,bc,cd);\n", + "\n", + "lab=ab*(math.cos(30.5*(math.pi/180)));\n", + "lbc=-bc*(math.cos(40*(math.pi/180)));\n", + "lcd=-cd*(math.cos(45*(math.pi/180)));\n", + "print('latitudes of AB,BC,CD=',lab,lbc,lcd);\n", + "\n", + "dab=ab*(math.sin(30.5*(math.pi/180)));\n", + "dbc=bc*(math.sin(40*(math.pi/180)));\n", + "dcd=-cd*(math.sin(45*(math.pi/180)));\n", + "print('depatures of AB,BC,CD ',dab,dbc,dcd);\n", + "\n", + "lc=-(lab+lbc+lcd);\n", + "ls=-(dab+dbc+dcd);\n", + "\n", + "print(lc,ls)\n", + "k=-ls/lc;\n", + "t=math.atan(k);\n", + "t=t*(180/(math.pi));\n", + "\n", + "print('Bearing of DA=',t);\n", + "DA=math.sqrt(lc*lc+ls*ls);\n", + "print('length DA=',DA);\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### ch-11 page 419 pb-7" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "('Distance AC=', 158.62738402665204)\n", + "('Distance BD=', 189.49088179672577)\n", + "('total latitude of C=', 18.46481737819161)\n", + "('total depature of C=', 21.113710931586226)\n", + "('total latitude of D=', 9.659924163502069)\n", + "-15.6914002615\n", + "('total depature of D=', 308.2914002614939)\n", + "('length CD=', 329.52276617048415, 'meters')\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "h1=1.48;h2=1.42;c=100;\n", + "ra1=0.77;ra2=1.60;ra3=2.43;\n", + "rb1=0.86;rb2=1.84;rb3=2.82;\n", + "a1=12.166;a2=10.5;\n", + "la=112.82;da=106.4;\n", + "lb=198.5;db=292.6;\n", + "ac=c*(ra3-ra1)*(math.cos(a1*(math.pi/180)))*(math.cos(a1*(math.pi/180)));\n", + "bd=c*(rb3-rb1)*(math.cos(a2*(math.pi/180)))*(math.cos(a2*(math.pi/180)));\n", + "\n", + "print('Distance AC=',ac);\n", + "print('Distance BD=',bd);\n", + "lac=-ac*(math.cos(53.5*(math.pi/180)));\n", + "tlc=la+lac;\n", + "print('total latitude of C=',tlc);\n", + "\n", + "dac=ac*(math.sin(53.5*(math.pi/180)));\n", + "da=-da;\n", + "tdc=da+dac;\n", + "print('total depature of C=',tdc);\n", + "\n", + "lbd=-bd*(math.cos(4.75*(math.pi/180)));\n", + "tld=lb+lbd;\n", + "print('total latitude of D=',tld);\n", + "\n", + "db=-db;\n", + "ddb=-bd*(math.sin(4.75*(math.pi/180)));\n", + "tdd=-(db+ddb);\n", + "print(ddb)\n", + "print('total depature of D=',tdd);\n", + "\n", + "dx=tdc+tdd;\n", + "cx=tlc-tld;\n", + "\n", + "CD=math.sqrt(dx*dx+cx*cx);\n", + "print('length CD=',CD,'meters');\n", + "\n", + "\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### chapter 11, section 11.8, pg 422, example 1" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "distance = 262.890670554\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "c=600\n", + "fplusd=0.5\n", + "s=3\n", + "n=6.860\n", + "distance= (c*s/n)+ fplusd\n", + "print \"distance =\",distance\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### chapter 11, section 11.8, pg423, eg2" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "5.142\n", + "RL of A= 259.692\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "d=65.340\n", + "x=4.5\n", + "y= math.tan(x)\n", + "v=5.142\n", + "RLofA=255.500+v-0.950\n", + "print v\n", + "print \"RL of A=\", RLofA\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### chapter 11, section 11.8, pg423, eg2" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "distance between B and BM= 49.706\n", + "RL of B= 515.398\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "s1=2\n", + "h1=0.655\n", + "v1=6.578\n", + "RL=v1+h1+510.5\n", + "v2=1.085\n", + "d2=12.396\n", + "h2=1.25\n", + "RLofB=RL-v2-h2\n", + "d=37.31+12.396\n", + "print \"distance between B and BM=\",d\n", + "print \"RL of B=\", RLofB\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### chapter 11, section 11.8, pg423, eg2" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": { + "collapsed": false + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "n= 15.9100040177\n" + ] + } + ], + "source": [ + "from __future__ import division\n", + "\n", + "import math\n", + "\n", + "d=124.45\n", + "c=1000\n", + "s=2\n", + "fplusd=0.3\n", + "thetha=(5+(6/30))\n", + "n=1980/d\n", + "print \"n=\",n\n" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 2", + "language": "python", + "name": "python2" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 2 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython2", + "version": "2.7.11" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/screenshots/13509393_889732434488692_1124174787_o.png b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/screenshots/13509393_889732434488692_1124174787_o.png Binary files differnew file mode 100644 index 00000000..6a4574cb --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/screenshots/13509393_889732434488692_1124174787_o.png diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/screenshots/13517827_889732431155359_775192974_o.png b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/screenshots/13517827_889732431155359_775192974_o.png Binary files differnew file mode 100644 index 00000000..10a4c1c7 --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/screenshots/13517827_889732431155359_775192974_o.png diff --git a/SURVYNG_AND_LEVELLING__by_N.N.BASAK/screenshots/13549112_889732614488674_584774303_o.png b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/screenshots/13549112_889732614488674_584774303_o.png Binary files differnew file mode 100644 index 00000000..7ecc2c63 --- /dev/null +++ b/SURVYNG_AND_LEVELLING__by_N.N.BASAK/screenshots/13549112_889732614488674_584774303_o.png |
