{"id":2139,"date":"2022-03-13T12:19:51","date_gmt":"2022-03-13T12:19:51","guid":{"rendered":"https:\/\/blog.metu.edu.tr\/eresmech\/?page_id=2139"},"modified":"2022-12-06T01:50:00","modified_gmt":"2022-12-06T01:50:00","slug":"3-8","status":"publish","type":"page","link":"https:\/\/blog.metu.edu.tr\/eresmech\/mechanisms\/ch3\/3-8\/","title":{"rendered":"3-8"},"content":{"rendered":"<div id=\"pl-gb2139-69d5b79774fa7\"  class=\"panel-layout\" ><div id=\"pg-gb2139-69d5b79774fa7-0\"  class=\"panel-grid panel-no-style\" ><div id=\"pgc-gb2139-69d5b79774fa7-0-0\"  class=\"panel-grid-cell\" ><div id=\"panel-gb2139-69d5b79774fa7-0-0-0\" class=\"so-panel widget widget_sow-editor panel-first-child panel-last-child widgetopts-SO\" data-index=\"0\" ><div\n\t\t\t\n\t\t\tclass=\"so-widget-sow-editor so-widget-sow-editor-base\"\n\t\t\t\n\t\t>\n<div class=\"siteorigin-widget-tinymce textwidget\">\n\t<h1><b>3.8<\/b> Numerical Solution of the Loop Closure Equations<\/h1>\n<p>The equations derived in the previous sections for position analysis can easily be used to obtain numerical results when the values of the link lengths are known. One can use a scientific calculator to compute the values of the dependant variables for a given value of the input parameter.\u00a0 If the analysis for the whole range of the input variable is to be performed, some means of programming or a package programme may be useful. In this section examples for different ways of numerical solution will be shown using.<\/p>\n<ul>\n<li>Scientific calculator,<\/li>\n<li>Spread-sheet package program such as EXCEL\u00ae<\/li>\n<li>A programming language such as BASIC, PASCAL, FORTRAN or C,<\/li>\n<li>A package program for solving equations such as MathCAD\u00ae or\u00a0 MATLAB\u00ae<\/li>\n<\/ul>\n<p>It must be noted that the numerical solution of mechanisms is not restricted by the examples shown. Advances in the programming techniques and hardware has a great influence on the numerical methods and new software and hardware capabilities may change the approach while keeping the basic mathematics the same.<\/p>\n<p>Teaching a programming language or a package program is\u00a0<span style=\"font-family: Arial, Helvetica, sans-serif\">our aim<\/span>. Therefore it is strongly recommended to the reader to refer to any one of the well written texts for a programming language or package program for a complete understanding of commands and features of these programs.<\/p>\n<p><span style=\"color: #ff0000\"><strong>Example 1:<\/strong><\/span><\/p>\n<p align=\"center\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-981 aligncenter\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-1.gif\" alt=\"\" width=\"476\" height=\"249\" \/><\/p>\n<p style=\"text-align: left\" align=\"center\">For the mechanism shown: a<sub>2<\/sub>\u00a0= 50 mm; a<sub>3<\/sub>\u00a0= 250 mm; b<sub>3<\/sub>\u00a0= 120 mm; a<sub>1<\/sub> = 20 mm and \u03b3<sub>3<\/sub>\u00a0= 30\u00b0<span style=\"font-size: 13.3333px\">. <\/span>Determine the coordinates of point C when \u03b8<sub>12<\/sub> = 60\u00b0.<\/p>\n<p>Since:<\/p>\n<p style=\"text-align: center\">sin\u03b8<sub>13<\/sub> = (a<sub>2<\/sub>sin\u03b8<sub>12<\/sub> \u2212 a<sub>1<\/sub>)\/a<sub>3<\/sub><\/p>\n<p style=\"text-align: center\">sin\u03b8<sub>13<\/sub> = (50 sin60\u00b0 \u2212 20)\/250 = 0.0932<\/p>\n<p>&nbsp;<\/p>\n<p>which yields \u03b8<sub>13<\/sub> = 5.3\u00b0 or 174.7\u00b0<\/p>\n<p>From the figure, we note that we must select \u03b8<sub>13<\/sub>\u00a0= 174.7\u00b0. Then:<\/p>\n<p style=\"text-align: center\">s<sub>14<\/sub>\u00a0= a<sub>2<\/sub> cos\u03b8<sub>12<\/sub> \u2212 a<sub>3<\/sub> cos\u03b8<sub>13<\/sub>\u00a0= 50 cos60\u00b0 \u2212 250 cos 174.6\u00b0\u00a0= 273.9 mm<\/p>\n<p style=\"text-align: center\">x<sub>C<\/sub>\u00a0= x + b<sub>3<\/sub>cos(\u03b8<sub>13<\/sub> \u2212 \u03b3<sub>3<\/sub>) = 273.9 + 120 cos(174.7\u00b0 \u2212 30\u00b0) =\u00a0176.0 mm<\/p>\n<p style=\"text-align: center\">y<sub>C<\/sub>\u00a0= a<sub>1<\/sub>\u00a0+ b<sub>3<\/sub>sin(\u03b8<sub>13<\/sub> \u2212 \u03b3<sub>3<\/sub>)\u00a0= 20 + 120 sin(174.7\u00b0 \u2212 30\u00b0) =\u00a089.4mm<\/p>\n<p>If the calculator is programmable, the above equations can all be programmed as successive operations and the co-ordinates of point C can be obtained for different crank angles. Otherwise, the same calculation must be typed several times. The values of the fixed parameters and the input parameter must all be stored and recalled when necessary. The computation will very much depend on the calculator that is available.<\/p>\n<p><span style=\"color: #ff0000\"><strong>Example 2:<\/strong><\/span><\/p>\n<p align=\"center\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-982\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-2.gif\" alt=\"\" width=\"460\" height=\"553\" \/><\/p>\n<p>Consider the four-bar mechanism shown which has the link lengths:<\/p>\n<p align=\"left\">a<sub>1<\/sub>\u00a0= 70 mm; \u00a0a<sub>2<\/sub>\u00a0= 35 mm; a<sub>3<\/sub>\u00a0= 62.3 mm; a<sub>4<\/sub>\u00a0= 56 mm; |BC| = c<sub>3<\/sub>\u00a0= 84.1mm; |AC| = b<sub>3<\/sub>\u00a0= 126.6mm<\/p>\n<p>You are to determine\u00a0<span style=\"font-family: Symbol\">q<\/span><sub>13<\/sub>,\u00a0<span style=\"font-family: Symbol\">q<\/span><sub>14<\/sub>\u00a0and the co-ordinates of point C for 0\u00b0 &lt;\u00a0<span style=\"font-family: Symbol\">q<\/span><sub>12<\/sub>&lt; 360\u00b0 with 5\u00b0\u00a0increments. The necessary equations are:<\/p>\n<p style=\"text-align: center\">x<sub>A<\/sub> = s cos\u03d5 = a<sub>2<\/sub>cos\u03b8<sub>12<\/sub> \u2212 a<sub>1<\/sub> \u00a0\u00a0\u00a0\u00a0\u00a0(Horizontal component of\u00a0<strong>AB<sub>0<\/sub><\/strong>)<\/p>\n<p style=\"text-align: center\">y<sub>A<\/sub> = s sin\u03d5 = a<sub>2<\/sub>sin\u03b8<sub>12<\/sub> \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0(Vertical component of\u00a0<strong>AB<sub>0<\/sub><\/strong>)<\/p>\n<p style=\"text-align: center\">s = <span class=\"katex-eq\" data-katex-display=\"false\"> \\displaystyle \\sqrt{{{{\\text{x}}_{\\text{A}}}^{2}+{{\\text{y}}_{\\text{A}}}^{2}}} <\/span>\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0(\u00a0magnitude of\u00a0<strong>AB<sub>0<\/sub><\/strong>)<\/p>\n<p style=\"text-align: center\">\u03d5 = atan2(x<sub>A<\/sub>, y<sub>A<\/sub>)\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0(angular position of\u00a0<strong>AB<sub>0<\/sub>\u00a0<\/strong>with respect to positive x axis)<\/p>\n<p style=\"text-align: center\">\u03bc = \u2220ABB<sub>0<\/sub> = cos<sup>-1<\/sup><span class=\"katex-eq\" data-katex-display=\"false\"> \\displaystyle \\left[ {\\frac{{{{\\text{a}}_{3}}^{2}+{{\\text{a}}_{4}}^{2}-{{\\text{s}}^{2}}}}{{2{{\\text{a}}_{3}}{{\\text{a}}_{4}}}}} \\right] <\/span><\/p>\n<p style=\"text-align: center\">\u03c8 = \u2220AB<sub>0<\/sub>B = cos<sup>-1<\/sup><span class=\"katex-eq\" data-katex-display=\"false\"> \\displaystyle \\left[ {\\frac{{{{\\text{a}}_{4}}^{2}+{{\\text{s}}^{2}}-{{\\text{a}}_{3}}^{2}}}{{2{{\\text{a}}_{4}}{{\\text{s}}}}}} \\right] <\/span><\/p>\n<p style=\"text-align: center\">\u03b8<sub>14<\/sub> = \u03d5 \u2212 \u03c8<\/p>\n<p style=\"text-align: center\">\u03b8<sub>13<\/sub>\u00a0= \u03b8<sub>14<\/sub>\u00a0\u2212 \u03bc<\/p>\n<p style=\"text-align: center\">x<sub>C<\/sub> = x<sub>A<\/sub>\u00a0+ a<sub>1<\/sub>\u00a0+ b<sub>3<\/sub>cos(\u03b3 + \u03b8<sub>14<\/sub>) \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0(Horizontal component of the distance <strong>A<sub>0<\/sub>C<\/strong>)<\/p>\n<p style=\"text-align: center\">y<sub>C<\/sub> = y<sub>A<\/sub>\u00a0+ b<sub>3<\/sub>sin(\u03b3 + \u03b8<sub>13<\/sub>) \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0(Vertical component of the distance <strong>A<sub>0<\/sub>C<\/strong>)<\/p>\n<p style=\"text-align: left\">where:<\/p>\n<p style=\"text-align: center\">\u03b3 = \u2220ABC = cos<sup>-1<\/sup><span class=\"katex-eq\" data-katex-display=\"false\"> \\displaystyle \\left[ {\\frac{{{{\\text{a}}_{3}}^{2}+{{\\text{b}}_{3}}^{2}-{{\\text{c}}_{3}}^{2}}}{{2{{\\text{a}}_{3}}{{\\text{b}}_{3}}}}} \\right] <\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-2141\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2022\/03\/3-8_clip_image031.gif\" alt=\"\" width=\"1029\" height=\"750\" \/><\/p>\n<p>For this example we will use of\u00a0<strong>EXCEL<\/strong>, a spread sheet program. In a spread sheet program you have a number of cells which can be considered as small boxes. Each cell is in a certain column (designated by letters) and a certain row (designated by a number). Thus a cell has a name such as A1 (first row first column) or C25 (25th row third column), etc. For each cell you can place some number or you can compute and place a new number generated from the numbers which are already placed\u00a0 into other cells. During this computation besides the simple mathematical operations such as addition and multiplication, the spreadsheet has several built-in functions. For example you can square the content of the cell A1 and multiply it will the cosine of the content of the cell A2 and place the new computed number into cell A3. Now, this future can be used very effectively for the analysis of mechanisms.<\/p>\n<p>Place the values of the fixed link lengths a<sub>1<\/sub>, a<sub>2<\/sub>, a<sub>3<\/sub>, a<sub>4<\/sub>\u00a0into cells A1 to A4 respectively. Into the cell A10, place the value of the angle \u03b8<sub>12<\/sub>\u00a0(0). Now multiply the value in cell A10 by \u03c0\/180 to obtain \u03b8<sub>12<\/sub>\u00a0in radians and place it into the cell B10 by simply typing\u00a0<strong>=A10 * PI ()\/180<\/strong>\u00a0into the cell B10. (all the angles must be in radians). You can also convert the angle in degrees into radians by using the build-in function in Excel \u201c<strong>RADIANS()<\/strong>\u201d by typing \u201c=RADIANS(A10)\u201d into cell B10. \u00a0Eq (a) is typed into the cell C10 by typing \u201c=$A$2*cos(B10)-$A$1\u201d. The $ sign must be used for fixed lengths ($ sign in front of column or row designation fixes that row or column. If $ sign is placed on both column and row that cell is fixed. i.e. the cell address becomes an absolute address instead of relative cell address. This is important when the cells are to be copied into other cells Instead of typing the dollar sign, after the cell name is typed you can press the \u201cF4\u201d key). In a similar fashion, you can type in all the formulas from B10 to M10 as shown in row 8 of the spread-sheet , remembering that the formulas are to be solved from top to bottom and from left to right. Normally you will see the resulting numerical values in the cells. Once these equations are typed, if you change the value of \u03b8<sub>12<\/sub>\u00a0in cell A10, the remaining dependent variables will all change automatically. Thus, you can determine the value of \u03b8<sub>13<\/sub>, \u03b8<sub>14<\/sub>\u00a0, and the coordinates of point C\u00a0 for any value of the input crank angle.<\/p>\n<div>\n<p>Now, let us place a different value of \u03b8<sub>12<\/sub>\u00a0into the cell A9. Next use\u00a0<strong>COPY<\/strong>\u00a0and\u00a0<strong>PASTE<\/strong>\u00a0option in the\u00a0<strong>EDIT <\/strong>menu to copy the contents of the cells from B10 to M10 into the cells B11 to M11. Note that the formulas will be copied and if you have not fixed the cell with the $ sign, what is A10, B10&#8230;F10 in the equations will all be A11, B11,&#8230;F11 in this row while $A$1 is kept the same.\u00a0 Next, let us place values of for every 15\u00b0\u00a0increments from 0\u00b0\u00a0to 360\u00b0\u00a0in column A starting with A10 (you can use 1\u00b0\u00a0interval as well. The only reason why 15\u00b0\u00a0interval is used is to show the result\u00a0<span style=\"font-family: Arial, Helvetica, sans-serif\">on<\/span> one screen only.) and copy the formulas in row 7 into all the rows . You will thus obtain<\/p>\n<\/div>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-984 aligncenter\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-4.gif\" alt=\"\" width=\"700\" height=\"479\" \/><\/p>\n<p>the value of all the dependent variables and the co-ordinates of point C for every crank angle you have selected!<\/p>\n<p>In order to visualise the motion, one can use the chart option of the spread-sheet program and plot \u03b8<sub>13<\/sub>, \u03b8<sub>14<\/sub>\u00a0and \u03bc as function of \u03b8<sub>12<\/sub>, or plot y<sub>c<\/sub>\u00a0versus x<sub>c<\/sub>\u00a0to determine the coupler path C as shown in the following figures<\/p>\n<p>Note that one can use the closed form method instead of step-wise method as well.<\/p>\n<p align=\"center\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-985\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-5.gif\" alt=\"\" width=\"600\" height=\"260\" \/><\/p>\n<p align=\"center\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-986\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-6.gif\" alt=\"\" width=\"642\" height=\"474\" \/><\/p>\n<p style=\"text-align: center\" align=\"center\">path of coupler point C<\/p>\n<p>To simulate the mechanism, let us also copy cells A10 to M10 into cells A6 to M6. We shall use the Control Toolbox feature of Excel. Usually the toolbox is not switched on. To switch this toolbox on, click \u201ctools-customize\u201d menu, place a check mark on \u201cControl Toolbox\u201d (You can also place a check mark on \u201cVisual Basic\u201d since we shall utilize this toolbox as well). First click onto the blue triangle-yellow ruler icon to switch to the design mode. There are several different buttons that can be used. We shall use the \u201cSpin Button\u201d which has two arrows up and down. If you are in the design mode, when you click onto the spin button the arrow cursor will change to plus sign. When you are near the cell A7 Right click on the mouse button and form a rectangular box . You will see that the spin button will be formed when you release the right button of the mouse. Right click onto this button shape you have formed and from the pop-up menu select \u201cproperties\u201d. You can change the colour, size of the button. One important feature is that you can link this button to a cell. In the \u201cLinked Cell\u201d box type A6, the cell which contains the value of the independent variable (crank angle in degrees). Into the \u201cMax\u201d cell type \u201c360\u201d and into the \u201cSmallChange\u201d cell type an integer step size (say 2). Close the menu and on the Control ToolBox click onto the Blue Triangle icon to close the design mode. \u0130f you click the left arrow on the spin button the value in cell A6 will decrease by 2 and if you click the right arrow, it will increase by 2 and the cells B6-M6 will also change accordingly.<\/p>\n<p>Let us now go to the cell Q2 and R2 and type 0 for the coordinates of A<sub>0<\/sub>. Next into cell Q3 and R3 type the coordinates of point A which is \u201c=A2*Cos(B6)\u201d and \u201c=A2*Sin(B6)\u201d. Into cells Q9, R9 type the coordinates of point B<sub>0<\/sub>\u00a0(\u201c=A1\u201d and \u201c0\u201d). Type the coordinates of point B into cells Q4 and R4 as \u201c=Q9+A4*Cos(J6)\u201d and \u201cR9+A4*Sin(J6)\u201d. If you type \u201c=L6\u201d and \u201c=M6\u201d into the cells Q5, R5 you have copied the coordinates of point C into these cells. Into the cells Q7, R7 and Q8, R8 copy the coordinates of points A and B Thus obtaining a table as shown in the following figure (actually in the cells you will be seeing the result as numbers). These are the coordinates of the relevant points on the mechanism.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-2142\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2022\/03\/3-8_clip_image040.gif\" alt=\"\" width=\"501\" height=\"217\" \/><\/p>\n<p>Next click on the chart icon. Select \u201cXY Scatter\u201d as the Chart type. As sub-type select \u201cscatter with data points connected by lines\u201d and click next. Select series and add series. For x values, select cells Q2 to Q8 which are the x values of points A<sub>0<\/sub>ABB<sub>0<\/sub>\u00a0and C. For y values select cells R2 to R8 which are the y values. You will obtain a chart similar to the one as shown\u00a0<span style=\"font-family: Arial, Helvetica, sans-serif\">below.<\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-2143 size-full aligncenter\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2022\/03\/3-8_clip_image042.gif\" alt=\"\" width=\"362\" height=\"241\" \/><\/p>\n<p>Although the four-bar mechanism is drawn, figure is not exactly satisfactory. First click on the chart and select Chart options from the pop-up menu. Select gridlines for both x and y axes. From the pop-up menu also select\u00a0<strong>&#8220;format plot area&#8221;<\/strong>to change the colour of the plot area from grey to some other colour (say no colour). Next click on the x and y axes and on the scales menu change the minimum and maximum -40&lt;x&lt;100 and -40&lt;y&lt;140 and change major unit to 20 for both x and y axis. On the screen try to make the grid lines look as a square. In such a case you have made the x and y axis scales the same and the mechanism proportion is more real. Left Click on the mechanism link and select \u201cFormat Data Point\u201d. In \u201cPatterns\u201d menu select the colour and the weight of the line that represents the links. For the marker shape select the circle and for the foreground select \u201cNo Colour\u201d. These circles will represent the revolute joints between the links. You can also plot the coupler curve shown in Figure as another series on this chart. Now the mechanism will be seen as shown. When you click on the spin button, the input crank angle value will change and so will the positions of all the links of the mechanism drawn on the chart. You can also use different colours for each link by defining A<sub>0<\/sub>A, ABCA and BB<sub>0<\/sub>\u00a0as different series.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-2144 size-full aligncenter\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2022\/03\/3-8_clip_image044.gif\" alt=\"\" width=\"344\" height=\"366\" \/><\/p>\n<p><b><a href=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/four-bar.xls\">Click here to download the Excel file for this example<\/a><\/b>.<\/p>\n<p>The analysis and simulation of mechanisms using Excel is explained in more detail in\u00a0<strong><a href=\"https:\/\/blog.metu.edu.tr\/eresmech\/mechanisms\/ek2\/\">Appendix II<\/a>.<\/strong><\/p>\n<\/div>\n<\/div><\/div><\/div><\/div><div id=\"pg-gb2139-69d5b79774fa7-1\"  class=\"panel-grid panel-no-style\" ><div id=\"pgc-gb2139-69d5b79774fa7-1-0\"  class=\"panel-grid-cell panel-grid-cell-empty\" ><\/div><div id=\"pgc-gb2139-69d5b79774fa7-1-1\"  class=\"panel-grid-cell panel-grid-cell-mobile-last\" >        <div id=\"panel-gb2139-69d5b79774fa7-1-1-0\" class=\"so-panel widget widget_wylwidget panel-first-child panel-last-child widgetopts-SO\" data-index=\"1\" ><div class=\"panel-widget-style panel-widget-style-for-gb2139-69d5b79774fa7-1-1-0\" >        <h3 class=\"widget-title\">Excel for Mechanical Engineers<\/h3>        <div class=\"lyte-wrapper lidget\" style=\"width:711px; height:400px; min-width:200px; max-width:100%;\"><div class=\"lyMe\" id=\"YLW_rckKx2deqRA\"><div id=\"lyte_rckKx2deqRA\" data-src=\"https:\/\/img.youtube.com\/vi\/rckKx2deqRA\/hqdefault.jpg\" class=\"pL\"><div class=\"play\"><\/div><div class=\"ctrl\"><div class=\"Lctrl\"><\/div><\/div><\/div><\/div><noscript><a href=\"https:\/\/youtu.be\/rckKx2deqRA\"><img decoding=\"async\" src=\"https:\/\/img.youtube.com\/vi\/rckKx2deqRA\/hqdefault.jpg\" alt=\"\" \/><\/a><\/noscript><\/div>\n        <div><\/div>\n        <\/div><\/div>        <\/div><div id=\"pgc-gb2139-69d5b79774fa7-1-2\"  class=\"panel-grid-cell panel-grid-cell-empty\" ><\/div><\/div><div id=\"pg-gb2139-69d5b79774fa7-2\"  class=\"panel-grid panel-no-style\" ><div id=\"pgc-gb2139-69d5b79774fa7-2-0\"  class=\"panel-grid-cell\" ><div id=\"panel-gb2139-69d5b79774fa7-2-0-0\" class=\"so-panel widget widget_sow-editor panel-first-child panel-last-child widgetopts-SO\" data-index=\"2\" ><div\n\t\t\t\n\t\t\tclass=\"so-widget-sow-editor so-widget-sow-editor-base\"\n\t\t\t\n\t\t>\n<div class=\"siteorigin-widget-tinymce textwidget\">\n\t<p><span style=\"color: #ff0000\"><strong>Example 3:<\/strong><\/span><\/p>\n<p>In this example we shall use MathCAD\u00ae which is a program for solving equations. There are other programs such as Mathematica\u00ae, Eureka\u00ae , MATLAB\u00ae, etc. which can essentially perform similar equation solving. In spread-sheet approach on the screen you see the cells and the numbers stored in the cells. In case of MathCAD the same basic idea is utilised. However the cells or boxes in which you store the numbers are not shown directly but\u00a0 they exist as a label name that the user assigns. For example, when you type: a:=120 , b:=50, \u03b8:=20, etc. The program generates boxes with labels a, b, attached to the cells in the computer which stores the values you have typed. (This assignment is done by typing the label and then typing \u201c:\u201d, from which you see \u201c:=\u201d on the screen and the computer expects an assignment to be made to the label typed).\u00a0 There is also a cell with label <span style=\"font-family: Symbol\">p<\/span>\u00a0which already has the value 3.141592654 that you can use. Similar to the spread-sheet, besides the simple mathematical operations there are several built-in functions as well. Now, if you type: x<sub>a<\/sub>:=a*cos(\u03b8<span style=\"font-family: Symbol\">*\u03c0<\/span>\/180)+b on the screen you will see (to type x<sub>a<\/sub>\u00a0type &#8220;<strong>x.a<\/strong>&#8221; on the screen. To type the Greek characters you can use\u00a0<strong>&#8220;Greek toolbar<\/strong>&#8221; or use [Ctrl] G to change any latin character to its Greek equivalent i.e typing\u00a0<strong>&#8220;q [Ctrl] G &#8220;<\/strong>will result in &#8220;<span style=\"font-family: Symbol\"><strong>q<\/strong><\/span>&#8220;<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-1003\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/eqn218-7.gif\" alt=\"\" width=\"201\" height=\"84\" srcset=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/eqn218-7.gif 134w, https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/eqn218-7-100x42.gif 100w\" sizes=\"auto, (max-width: 201px) 100vw, 201px\" \/><\/p>\n<p>A new cell with label x<sub>a<\/sub>\u00a0is made available and into this cell the numerical value is placed using the value of the cells a, b and\u00a0<span style=\"font-family: Symbol\">q<\/span>\u00a0and the operations indicated. If you want to see the contents of the cell simply type x<sub>a<\/sub>\u00a0=. If you had a:=120, b:=50, \u03b8:=20 assigned beforehand, on the screen you will see: x<sub>a<\/sub>\u00a0= 143.969.\u00a0 Hence you can see the operations that you have performed as you would write it as a text, and you see the values in certain cells whenever necessary. The formulas that you type are shown in a more realistic form.<\/p>\n<p>In figure (a) below, an adjustable pump is shown. Crank A<sub>0<\/sub>A (which is constructed as an eccentric) is driven by an electric motor through the worm gear. The schematic drawing of the mechanism is shown in\u00a0<span style=\"font-family: Arial, Helvetica, sans-serif\">f<\/span>igure (b). The<span style=\"font-family: Arial, Helvetica, sans-serif\">\u00a0<\/span>stroke adjustment is obtained by moving the location of the pivot point B<sub>0<\/sub>\u00a0by means of an adjustment screw.<\/p>\n<p>Note that A<sub>0<\/sub>ABB<sub>0<\/sub> forms a four-bar mechanism. Unlike the previous examples the fixed link A<sub>0<\/sub>B<sub>0<\/sub>\u00a0is not horizontal and B<sub>0<\/sub>\u00a0is to the left of A<sub>0<\/sub>.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-989\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-9.gif\" alt=\"\" width=\"792\" height=\"552\" \/><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-990\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-10.gif\" alt=\"\" width=\"875\" height=\"458\" \/><\/p>\n<p align=\"center\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-991\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-11.gif\" alt=\"\" width=\"486\" height=\"382\" \/><\/p>\n<p>We can write the vector equation:<\/p>\n<p style=\"text-align: center\"><strong>B<sub>0<\/sub>A<\/strong> = <strong>B<sub>0<\/sub>A<sub>0<\/sub><\/strong>\u00a0+ <strong>A<sub>0<\/sub>A<\/strong><\/p>\n<p style=\"text-align: center\">se<sup>i\u03d5<\/sup>\u00a0 = b<sub>1 <\/sub>\u2212 is<sub>1\u00a0<\/sub>+ a<sub>2<\/sub>e<sup>i<\/sup><sup>\u03b8<sub>12<\/sub><\/sup><\/p>\n<p>or, equating the real and imaginary parts:<\/p>\n<p style=\"text-align: center\">s cos\u03d5\u00a0= b<sub>1<\/sub>\u00a0+ a<sub>2<\/sub>cos\u03b8<sub>12<\/sub><br \/>\ns sin\u03d5 = s<sub>1<\/sub>\u00a0+ a<sub>2<\/sub>sin\u03b8<sub>12<\/sub><\/p>\n<p align=\"left\">From which we can solve for s and \u03d5. Then using the triangle B<sub>0<\/sub>A B, the angle \u03c8\u00a0can be determined using the cosine theorem:<\/p>\n<p style=\"text-align: center\">\u03c8 = cos<sup>-1<\/sup><span class=\"katex-eq\" data-katex-display=\"false\"> \\displaystyle \\left[ {\\frac{{{{\\text{a}}_{4}}^{2}+{{\\text{s}}^{2}}-{{\\text{a}}_{3}}^{2}}}{{2{{\\text{a}}_{4}}{{\\text{s}}}}}} \\right] <\/span><\/p>\n<p>Now the angle \u03b8<sub>14<\/sub>\u00a0is:<\/p>\n<p style=\"text-align: center\">\u03b8<sub>14<\/sub> = \u03d5 \u2212 \u03c8 \u2212 \u03c0\/2<\/p>\n<p align=\"center\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-992\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-12.gif\" alt=\"\" width=\"621\" height=\"390\" \/><\/p>\n<p>For the double slide (links 4 and 5), the loop equation is:<\/p>\n<p style=\"text-align: center\">c<sub>4<\/sub>e<sup>i<\/sup><sup>\u03b8<sub>14<\/sub><\/sup>\u00a0+ is<sub>4<\/sub>e<sup>i<\/sup><sup>\u03b8<sub>14<\/sub><\/sup>\u00a0= \u2212is<sub>1<\/sub> \u2212 s<sub>15<\/sub>\u00a0+ b<sub>1<\/sub><\/p>\n<p>\u00a0which yields:<\/p>\n<p style=\"text-align: center\">s<sub>4<\/sub>\u00a0= \u2212(s<sub>1<\/sub> + c<sub>4<\/sub>sin\u03b8<sub>14<\/sub>)\/cos\u03b8<sub>14<\/sub><\/p>\n<p style=\"text-align: center\">s<sub>15<\/sub> = \u2212c<sub>4<\/sub>cos\u03b8<sub>14<\/sub> + s<sub>4<\/sub>sin\u03b8<sub>14<\/sub> + b<sub>1<\/sub><\/p>\n<p align=\"left\">In MathCAD you will write these equations exactly the same way. But before performing any computation, you must define the numerical values of the variables you are going to use either as input or as a result of some computation. The MathCAD sheet for the input crank angle \u03b8<sub>12<\/sub>\u00a0= 120\u00b0\u00a0will look as follows:<\/p>\n<hr size=\"2\" \/>\n<p><span style=\"color: #000099\"><b>MathCAD Output &#8211; 1<br \/>\n<\/b><\/span><span style=\"color: #cc0000\">ADJUSTABLE STROKE PUMP<\/span><\/p>\n<p>Determination of the output for a given crank angle.<\/p>\n<p>conv := \u03c0\/180 (to convert degrees into radians)<\/p>\n<p>A<sub>0<\/sub>A = a<sub>2<\/sub>\u00a0:= 55 <sub><b>\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<\/b><\/sub>\u00a0AB = a<sub>3<\/sub> := 240 <sub><b>\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<\/b><\/sub>\u00a0BB<sub>0<\/sub>\u00a0= a<sub>4<\/sub>\u00a0:= 165 <sub><b>\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<\/b><\/sub>\u00a0b<sub>1<\/sub>\u00a0:= 185 <sub><b>\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<\/b><\/sub>\u00a0s<sub>1<\/sub>\u00a0:= 90\u00a0<sub><b>\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<\/b><\/sub>\u00a0c<sub>4<\/sub>\u00a0:= 70<\/p>\n<p>\u03b8\u00a0:= 320 <sub><b>\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<\/b><\/sub> \u03b8\u00a0:= \u03b8<sup>.<\/sup>conv\u00a0<sub><b>\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<\/b><\/sub><\/p>\n<p>x<sub>a<\/sub>\u00a0:= b<sub>1<\/sub>\u00a0+ a<sub>2<\/sub><sup>.<\/sup>cos(\u03b8)\u00a0<sub><b>\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<\/b><\/sub>\u00a0y<sub>a<\/sub> := \u2212s<sub>1<\/sub>\u00a0+ a<sub>2<\/sub><sup>.<\/sup>sin(\u03b8)\u00a0<sub><b> \u00a0\u00a0\u00a0\u00a0\u00a0<\/b><\/sub><\/p>\n<p>s := <span class=\"katex-eq\" data-katex-display=\"false\"> \\displaystyle\\sqrt{{{{{{{{\\text{x}}_{\\text{a}}}}}}^{2}}+{{{{{{\\text{y}}_{\\text{a}}}}}}^{2}}}} <\/span><\/p>\n<p>\u03d5 := angle(x<sub>a<\/sub>, y<sub>a<\/sub>)<\/p>\n<p>\u03c8 := acos<span class=\"katex-eq\" data-katex-display=\"false\"> \\displaystyle \\left[ {\\frac{{{{\\text{a}}_{4}}^{2}+{{\\text{s}}^{2}}-{{\\text{a}}_{3}}^{2}}}{{2{{\\text{a}}_{4}}{{\\text{s}}}}}} \\right] <\/span><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone wp-image-1009\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/eqn218-13.gif\" alt=\"\" width=\"575\" height=\"246\" \/><\/p>\n<p>If we want to analyse the motion for the whole cycle, say for every 5\u00ba, then we must use indices for the variables so that the computer can allocate several cells to that variable. Secondly, besides obtaining numerical output we can also ask MathCAD to plot the result. In such a case, the solution will be as follows:<\/p>\n<hr size=\"2\" \/>\n<p><span style=\"color: #cc0000\"><span style=\"color: #000099\"><b>MathCAD Output &#8211; 2<\/b><\/span><\/span><\/p>\n<p><span style=\"color: #cc0000\">ADJUSTABLE STROKE PUMP<br \/>\n<\/span><\/p>\n<p>Complete cycle analysis.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone wp-image-1010\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/eqn218-14.gif\" alt=\"\" width=\"528\" height=\"234\" \/> \u00a0\u00a0<sub><b>\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<\/b><\/sub><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone wp-image-1011\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/eqn218-15.gif\" alt=\"\" width=\"675\" height=\"444\" \/><\/p>\n<p align=\"center\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-993\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-13.gif\" alt=\"\" width=\"846\" height=\"548\" \/><\/p>\n<p style=\"text-align: left\" align=\"center\">In this case i is not a subscript but an index that tells the computer to allocate an array of cells for that variable with index i. \u03b8<sub>i<\/sub>\u00a0is written as <span style=\"font-family: Symbol\"><strong>q<\/strong><\/span><strong>[i<\/strong>\u00a0not<strong>\u00a0<span style=\"font-family: Symbol\">q<\/span>.i<\/strong>\u00a0and refers to the ith element of the array\u00a0<span style=\"font-family: Symbol\">q<\/span><sub>i<\/sub> as in the case of a subscript. If we want to determine the effect of the distance s<sub>1<\/sub> which is adjusted by the screw, to the displacement of the pump, then we must perform the same analysis for different values of s. In this case we must use another index, k. The result is as follows:<\/p>\n<hr size=\"2\" \/>\n<p><span style=\"color: #cc0000\"><span style=\"color: #000099\"><b>MathCAD Output &#8211; 3<\/b><\/span><\/span><\/p>\n<p><span style=\"color: #cc0000\">ADJUSTABLE STROKE PUMP<br \/>\n<\/span><\/p>\n<p>Analysis for complete cycle and for different strokes:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone wp-image-1012\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/eqn218-16.gif\" alt=\"\" width=\"629\" height=\"582\" \/><\/p>\n<p align=\"center\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-994\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/image218-14.gif\" alt=\"\" width=\"1060\" height=\"618\" \/><\/p>\n<p align=\"left\">As it can be seen from the figure, as the distance s<sub>1<\/sub> is increased from 50 mm (k = 0) to 150 mm (k = 5), the stroke changes 63.8 mm (k = 0) to 155.4 mm.<\/p>\n<p align=\"left\"><a href=\"https:\/\/ocw.metu.edu.tr\/pluginfile.php\/3961\/mod_resource\/content\/12\/ch3\/AdjustableStrokePump.mcd\">You can download the Mathcad file by clicking here<\/a><\/p>\n<p><span style=\"color: #ff0000\"><strong>Example 4:<\/strong><\/span><\/p>\n<p>The example solved in MathCad can also be solved in Excel. Since the number of steps in the calculation has increased, if we use the Excel sheet only, the number of columns and rows to be used will be large. Instead of using the standard build in functions, we can write a function that will evaluate s<sub>16<\/sub>\u00a0for a given value of the input angle\u00a0<span style=\"font-family: Symbol\">q<\/span><sub>12<\/sub>\u00a0and adjustment position s<sub>1<\/sub>Hence the output parameter value can be evaluated in a single step. From &#8220;Tools&#8221; menu select &#8220;Macro&#8221; and then &#8220;Visual Basic Editor&#8221; . Else you can use Alt+F11 keys. When you are in VBA editor, select &#8220;module&#8221; from the insert menu. On the module you can write ant VBA program. Let us write a program named AdjustablePump() which will utize the formulas that we have derived. This program may look like:<\/p>\n<p style=\"padding-left: 40px\">Global Const PI = 3.1415926<br \/>\n<b><strong>Function AdjustablePump(Crank;Coupler;Rocker;FixedLink_X;Eccentricity;B$7;$A8)<\/strong><\/b><br \/>\nDim Xa, Ya As Double<br \/>\nDim S, fi, si, Theta2, Theta4 As Double<br \/>\nDim S_4 As Double<br \/>\nTheta2 = ThetaDegrees * PI \/ 180<br \/>\nXa = FixedLink_X + Crank * Cos(Theta2)<br \/>\nYa = -AdjustmentLength + Crank * Sin(Theta2)<br \/>\nS = Mag(Xa, Ya) \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<span style=\"color: #008080\">&#8216;BoA distance<\/span><br \/>\nfi = Ang(Xa, Ya) \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<span style=\"color: #008080\">\u00a0&#8216;Angle AB0 makes with the horizontal<\/span><br \/>\nSi = AngCos(Rocker, S, Coupler) \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<span style=\"color: #008080\">\u00a0&#8216;Angle BB0A<\/span><br \/>\nTheta4 = fi &#8211; si &#8211; PI \/ 2 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0<span style=\"color: #008080\">&#8216;Rocker angle<\/span><br \/>\nS4 = -(AdjustmentLength + Eccentricity * Sin(Theta4)) \/ Cos(Theta4)<br \/>\nAdjustablePump = FixedLink_X &#8211; Eccentricity * Cos(Theta4) + S4 * Sin(Theta4)<br \/>\n<b>End Function<\/b><\/p>\n<p>&nbsp;<\/p>\n<p>You can easily follow this function routine. However this function uses three functions which are not available in Excel Library. These are\u00a0<strong>Mag(X,Y)<\/strong>\u00a0which determines the hypotenuse of a right angled triangle whose two other sides are X and Y.\u00a0<strong>Ang(X,Y)<\/strong>\u00a0is the same as\u00a0<strong>ATAN2(X,Y)<\/strong>\u00a0function. The reason we use it is that ATAN2() is not available in visual Basic.<strong>\u00a0AngCos(U1,U2, Opposite)<\/strong> applies the cosine theorem to determine the angle opposite to the side called &#8220;Opposite&#8221; in a triangle whose other two sides are U1 and U2. Detailed explanation for these functions and more are given in <strong><a href=\"https:\/\/blog.metu.edu.tr\/eresmech\/mechanisms\/ek2\/\">Appendix-2<\/a><\/strong><b>.<\/b>\u00a0You must insert the Basic.Bas file into your module to use these functions .<\/p>\n<p>Let us place the values of the link lengths into cells B1 to B5 and name these cells as\u00a0<b>Crank<\/b>,\u00a0<b>Coupler<\/b>,\u00a0<b>Rocker<\/b>,\u00a0<b>FixedLink_X<\/b>, and\u00a0<b>Eccentricity<\/b>. Into cell B7let us place an adjustment length s<sub>1<\/sub>\u00a0(s<sub>1<\/sub>\u00a0= 20mm).<\/p>\n<p>Now click on to the cell B8 and into this cell write.\u00a0<b>= AdjustablePump(Crank; Coupler; Rocker; FixedLink_X; Eccentricity; B$7; $A8 )(you can also use &#8220;Insert-Function-User Defined&#8221; commands to reach this function and fill in the required parameters)<\/b>\u00a0Once you enter this function the position of the plunger (s<sub>15<\/sub>) will be determined for the given link lengths, the crank angle in A8 and the adjustmant length in cell B7. If you copy This cell (B8) up to B26, you will be able to determine the plunger position for every crank angle and for the given adjustmant length B7, since when during vertical copying the cell B\/ does not change since we have fixed the row as B$7.<\/p>\n<p>The distance s\u00a0<sub>1<\/sub>\u00a0is changed by an adjustment screw shown. If we want to determine the effect of this adjustment, s\u00a0<sub>1<\/sub>, on the displacement of the plunger, s<sub>15<\/sub>\u00a0, into cells B7, C7, D7, E7, F7 and G7 we place different values of s\u00a0<sub>1<\/sub>. Since the value of s\u00a0<sub>1<\/sub>\u00a0in cell B7 is written as B$7 row 7 does not change (whereas in column C B$7 will be C$7). Now copy cell B8 horizontally from B8 to G8 and then verticall up to G26. Each column will refer to a particilar adjustmant s\u00a0<sub>1<\/sub>\u00a0and each row will refer to a particular crank angle\u00a0<span style=\"font-family: Symbol\">q<\/span><sub>12<\/sub>\u00a0The result is given as a Excell sheet and sketch below.<\/p>\n<p><b><a href=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/09\/ayarlipompa.xls\">You can download the Excel file<\/a><\/b>.<\/p>\n<p>If your company is producing such a pump with different size, everybody in the company can make use of such a file very easily!!!<\/p>\n<table border=\"1\" width=\"100%\">\n<tbody>\n<tr bgcolor=\"#CCCCCC\">\n<td><\/td>\n<td align=\"center\">A<\/td>\n<td align=\"center\">B<\/td>\n<td align=\"center\">C<\/td>\n<td align=\"center\">D<\/td>\n<td align=\"center\">E<\/td>\n<td align=\"center\">F<\/td>\n<td align=\"center\">G<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">1<\/td>\n<td>A2<\/td>\n<td>55<\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">2<\/td>\n<td>A3<\/td>\n<td>240<\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">3<\/td>\n<td>A4<\/td>\n<td>165<\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">4<\/td>\n<td>B1<\/td>\n<td>185<\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">5<\/td>\n<td>C4<\/td>\n<td>70<\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">6<\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<td><\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">7<\/td>\n<td>\u03b8<sub>12<\/sub> \/ s<sub>1<\/sub><\/td>\n<td>20<\/td>\n<td>50<\/td>\n<td>100<\/td>\n<td>150<\/td>\n<td>200<\/td>\n<td>250<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">8<\/td>\n<td>0<\/td>\n<td>252.0937<\/td>\n<td>247.4078<\/td>\n<td>249.6989<\/td>\n<td>254.5738<\/td>\n<td>251.8945<\/td>\n<td>231.2036<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">9<\/td>\n<td>20<\/td>\n<td>252.1405<\/td>\n<td>245.7816<\/td>\n<td>248.1027<\/td>\n<td>255.8408<\/td>\n<td>258.3663<\/td>\n<td>245.9161<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">10<\/td>\n<td>40<\/td>\n<td>252.2511<\/td>\n<td>245.3691<\/td>\n<td>249.6373<\/td>\n<td>261.9564<\/td>\n<td>270.8585<\/td>\n<td>266.6617<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">11<\/td>\n<td>60<\/td>\n<td>252.2215<\/td>\n<td>246.2394<\/td>\n<td>254.6551<\/td>\n<td>273.3395<\/td>\n<td>289.4999<\/td>\n<td>292.9677<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">12<\/td>\n<td>80<\/td>\n<td>252.0927<\/td>\n<td>248.7137<\/td>\n<td>263.8473<\/td>\n<td>290.7546<\/td>\n<td>314.6072<\/td>\n<td>324.3855<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">13<\/td>\n<td>100<\/td>\n<td>252.2614<\/td>\n<td>253.48<\/td>\n<td>278.3194<\/td>\n<td>315.0693<\/td>\n<td>345.947<\/td>\n<td>359.3292<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">14<\/td>\n<td>120<\/td>\n<td>253.5854<\/td>\n<td>261.6704<\/td>\n<td>299.246<\/td>\n<td>346.0074<\/td>\n<td>380.8984<\/td>\n<td>393.3704<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">15<\/td>\n<td>140<\/td>\n<td>257.3465<\/td>\n<td>274.4691<\/td>\n<td>325.9924<\/td>\n<td>379.1378<\/td>\n<td>412.1455<\/td>\n<td>418.5263<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">16<\/td>\n<td>160<\/td>\n<td>264.4316<\/td>\n<td>291.0171<\/td>\n<td>352.1259<\/td>\n<td>403.739<\/td>\n<td>429.014<\/td>\n<td>426.5381<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">17<\/td>\n<td>180<\/td>\n<td>272.9428<\/td>\n<td>305.118<\/td>\n<td>365.8015<\/td>\n<td>409.1448<\/td>\n<td>425.2503<\/td>\n<td>414.8999<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">18<\/td>\n<td>200<\/td>\n<td>277.6904<\/td>\n<td>308.9373<\/td>\n<td>361.3043<\/td>\n<td>395.2032<\/td>\n<td>404.2678<\/td>\n<td>388.3058<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">19<\/td>\n<td>220<\/td>\n<td>276.148<\/td>\n<td>302.3755<\/td>\n<td>344.4052<\/td>\n<td>370.4372<\/td>\n<td>374.5207<\/td>\n<td>354.1352<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">20<\/td>\n<td>240<\/td>\n<td>270.7803<\/td>\n<td>291.0105<\/td>\n<td>323.5058<\/td>\n<td>343.1006<\/td>\n<td>343.1258<\/td>\n<td>318.3708<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">21<\/td>\n<td>260<\/td>\n<td>264.7198<\/td>\n<td>279.2832<\/td>\n<td>303.5609<\/td>\n<td>317.7134<\/td>\n<td>314.1433<\/td>\n<td>284.8884<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">22<\/td>\n<td>280<\/td>\n<td>259.5959<\/td>\n<td>269.1044<\/td>\n<td>286.4366<\/td>\n<td>296.0851<\/td>\n<td>289.5916<\/td>\n<td>256.4869<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">23<\/td>\n<td>300<\/td>\n<td>255.902<\/td>\n<td>260.9577<\/td>\n<td>272.5591<\/td>\n<td>278.8124<\/td>\n<td>270.5826<\/td>\n<td>235.7459<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">24<\/td>\n<td>320<\/td>\n<td>253.6054<\/td>\n<td>254.7785<\/td>\n<td>261.8858<\/td>\n<td>266.0913<\/td>\n<td>257.8185<\/td>\n<td>224.6681<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">25<\/td>\n<td>340<\/td>\n<td>252.4534<\/td>\n<td>250.3422<\/td>\n<td>254.2934<\/td>\n<td>258.0019<\/td>\n<td>251.617<\/td>\n<td>223.5825<\/td>\n<\/tr>\n<tr>\n<td align=\"center\" bgcolor=\"#CCCCCC\">26<\/td>\n<td>360<\/td>\n<td>252.0937<\/td>\n<td>247.4078<\/td>\n<td>249.6989<\/td>\n<td>254.5738<\/td>\n<td>251.8945<\/td>\n<td>231.2036<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p align=\"center\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-2145 size-full aligncenter\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2022\/03\/image218-16e.gif\" alt=\"\" width=\"621\" height=\"412\" \/><\/p>\n<p><strong>A note on the accuracy of computation<\/strong>:\u00a0In the computers the numbers are stored using finite number of digits. For example in MathCAD, all the computations are performed in 15 significant digits. In certain other programs or compilers the number of digits may vary from 6 to 32. In mechanism analysis the number of significant digits used by the hardware or software is more than sufficient and very precise results can be obtained. The exceptions are when the mechanism is at or near a dead centre position, in which case the precision may decrease or an error command may be displayed (such as division by zero or negative value in the square root). This will also indicate a misbehaving mechanism in practice.<\/p>\n<p>In engineering computations<strong>\u00a0\u201cNo result can be more accurate than the accuracy of the input parameters\u201d<\/strong>. Therefore, your result must not contain more significant digits than the number of significant digit for the\u00a0 input parameters. For example, if you are to manufacture the links using classical tools in a machine shop and if the accuracy of the link lengths are within \u00b10.1 mm, an x displacement result cannot be x = 26.231456789 mm. This must be rounded to x = 26.2 mm.<\/p>\n<\/div>\n<\/div><\/div><\/div><\/div><\/div>\n\n\n<p>  <a href=\"https:\/\/blog.metu.edu.tr\/eresmech\/mechanisms\/ch3\/3-7\/\" data-type=\"page\"><img loading=\"lazy\" decoding=\"async\" width=\"38\" height=\"38\" class=\"wp-image-16\" style=\"width: 38px\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/04\/back_button.gif\" alt=\"\"><\/a><a href=\"https:\/\/blog.metu.edu.tr\/eresmech\/mechanisms\/ch3\/\" data-type=\"page\" data-id=\"52\"><img loading=\"lazy\" decoding=\"async\" width=\"38\" height=\"38\" class=\"wp-image-17\" style=\"width: 38px\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/04\/contents_button.gif\" alt=\"\"><\/a><a href=\"https:\/\/blog.metu.edu.tr\/eresmech\/mechanisms\/\" data-type=\"page\" data-id=\"47\"><img loading=\"lazy\" decoding=\"async\" width=\"38\" height=\"38\" class=\"wp-image-18\" style=\"width: 38px\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/04\/home_button.gif\" alt=\"\"><\/a><a href=\"https:\/\/blog.metu.edu.tr\/eresmech\/mechanisms\/ch3\/3-9\/\" data-type=\"page\" data-id=\"92\"><img loading=\"lazy\" decoding=\"async\" width=\"38\" height=\"38\" class=\"wp-image-20\" style=\"width: 38px\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/04\/next_button.gif\" alt=\"\"><\/a><img loading=\"lazy\" decoding=\"async\" width=\"119\" height=\"40\" class=\"wp-image-15\" style=\"width: 119px\" src=\"https:\/\/blog.metu.edu.tr\/eresmech\/files\/2021\/04\/ceres.gif\" alt=\"\">       <\/p>\n","protected":false},"excerpt":{"rendered":"<p>3.8 Numerical Solution of the Loop Closure Equations The equations derived in the previous sections for position analysis can easily be used to obtain numerical results when the values of the link lengths are known. One can use a scientific &hellip;<\/p>\n<p class=\"read-more\"> <a class=\"more-link\" href=\"https:\/\/blog.metu.edu.tr\/eresmech\/mechanisms\/ch3\/3-8\/\"> <span class=\"screen-reader-text\">3-8<\/span> Devam\u0131n\u0131 Oku &raquo;<\/a><\/p>\n","protected":false},"author":7747,"featured_media":0,"parent":1950,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"full-width-page.php","meta":{"footnotes":"","_links_to":"","_links_to_target":""},"class_list":["post-2139","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/blog.metu.edu.tr\/eresmech\/wp-json\/wp\/v2\/pages\/2139","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blog.metu.edu.tr\/eresmech\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/blog.metu.edu.tr\/eresmech\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/blog.metu.edu.tr\/eresmech\/wp-json\/wp\/v2\/users\/7747"}],"replies":[{"embeddable":true,"href":"https:\/\/blog.metu.edu.tr\/eresmech\/wp-json\/wp\/v2\/comments?post=2139"}],"version-history":[{"count":0,"href":"https:\/\/blog.metu.edu.tr\/eresmech\/wp-json\/wp\/v2\/pages\/2139\/revisions"}],"up":[{"embeddable":true,"href":"https:\/\/blog.metu.edu.tr\/eresmech\/wp-json\/wp\/v2\/pages\/1950"}],"wp:attachment":[{"href":"https:\/\/blog.metu.edu.tr\/eresmech\/wp-json\/wp\/v2\/media?parent=2139"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}