ME 3010 Project 1 Consider the planar linkage shown below

 ME 3010 Project 1  Consider the planar linkage shown below. The crank (link with length R) will turn counter-clockwise with constant angular velocity. The slider will experience a force of 100 N directed to the left. Your task is to write a program that calculates the required input torque and the magnitude of the pin reaction forces for two revolutions of the input link. You will then use your program to determine how changing the dimension H changes the average magnitude of the pin reaction forces at pins A, B and C (refer to the free-body diagram for the location of these pins). The first task is to determine the coordinates ,s given  . The constraint equations that must be satisfied are: R L s cos( ) cos( )     (1) R L H sin( ) sin( )     (2) You are required to solve these equations using Newton’s method. The next task is to solve for the time derivatives of ,s given  (constant). Equations (1) and (2) may be differentiated to obtain the following linear systems of equation in the unknowns. You may use the MATLAB “\” operator to solve these.         sin 1 sin cos 0 cos L R L s R                                  (3) Slider R L  H g s R = 0.5m L = 1.25m H = 0.25m = 2π rad/sec             2 2 2 2 sin 1 cos cos cos 0 sin sin L R L L s R L                                       (4) The next step is to draw the free-body diagrams and write the corresponding dynamic equations. Following is the free-body diagram. Note that we are neglecting friction on the slider. Ax Ay u T Bx By m1g Bx By Cx Cy m2g Cy Cx m3g N M PoA The resulting dynamics equations may be cast in matrix form as 1 1 1 1 2 2 2 2 2 3 3 1 0 1 0 0 0 0 0 1 0 1 0 0 0 ( ) 0 0 1 2 2 2 2 0 0 0 1 0 1 0 0 0 0 0 1 0 1 0 ( ) 0 0 0 2 2 2 2 0 0 0 0 1 0 0 x y x y x Ax m a Ay m a g R R R R S C S C Bx By m a Cx m a g L L L L Cy I S C S C T PoA m a                                                                         (5) where m1 = 10 kg, m2 = 25 kg, m3 = 10 kg and I2 = m2 L 2 /12         , cos( ),sin( ) , cos( ),sin( ) C S C S           (6) and 1 1 1 1 2 1 2 2 2 2 2 3 2 ˆ a R a a R R a s i                (7) with     1 2 ˆ ˆ cos( ) sin( ) 2 ˆ ˆ cos( ) sin( ) 2 R R i j L R i j         (8) 1 2 2 ˆ ˆ ˆ k k k          (9) You may use the MATLAB cross function to obtain the above cross products. Note that the accelerations are known functions of the derivatives of the three coordinates. Turn in a brief report containing the following: (i) title page (ii) introduction and objectives (iii) program design (choices for functions) (iv) results including plots (labeled) of simulation results for the base case  (H = 0.25m). 1.  and s vs.  2.  and s vs.  3.  and s vs.  4. A , B , C and T vs.  (v) include a plot displaying the average value of the magnitude of the pin  reaction forces over the complete cycle vs. H for the range    0 7 H 0 7 . . (vi) discussion of results including conclusion for part (v)