MATLAB Book 110.5772/46450Dynamic and Quasi-Static Simulation of a Novel Compliant MEMS Force Amplifier by Matlab/SimulinkKosaErgin1TrabzonLevent1SonmezUmit2KizilHuseyin3Mechanical Engineering, Istanbul Technical University, Istanbul,TurkeyMechanical Engineering, Sharjah University, Dubai,United Arab EmiratesMetallurgical and Materials Engineering, Istanbul Technical University, Istanbul,Turkey1. IntroductionMEMS are micro electromechanical systems having component sizes varying from 1 micrometer to 1 millimeter and provide specific engineering operations. MEMS are used as a micro sensor, micro actuator, micro accelerometer etc. and also have tendency to function rapidly due to having low inertia moment and affected less by disturbances coming from environment due to their small size (Hsu, 2002). Compliant mechanisms having an ability to transmit motion and energy via their flexible hinges and/or flexible components instead of joints and rigid components, perform large deflections (Sreekumar et al, 2008). The large deflections of compliant mechanisms instead of rigid-body mechanisms depend on applied force that causes a much more complexity to nonlinear analysis (Ashok, 2000). Moreover, the geometry of several flexure hinges are modeled as torsion springs in its pseudo-rigid-body mechanisms (Howell, 2001). Flexible segments of compliant mechanism store and transfer energy when it is functioning (Howell, 2001; Tantanawat & S. Kota, 2007). Flexible links having small cross sections instead of traditional joints provide acting of mechanism due to its very low moment of inertia (Howell, 2001; Lobontiu et al, 2001).Compliant four-link mechanism is designed as seen in Fig. 1 achieving force or displacement application according to the output spring constant and also, studied on size optimization to achieve maximum mechanical or geometric benefit at specific spring constants (Parkinson et al, 2001). Large displacement amplifier integrated with comb drive achieves 100 times amplifying of comb drive displacement by means of its design is modeled (Li et al, 2005).Figure 1.Schematic view of force or displacement amplifier’ mechanismCompliant MEMS have been used as a force amplifier in micro actuators and micro-mechanisms (Parkinson et al, 2001). They are preferred since there is no need for assembly, no energy loss due to absence of friction, thus requiring no need for lubrication all of which providing high precision (Kosa et al, 2010). Besides, compliant micro mechanisms could be activated by mechanically (Han et al, 2007; Krishnan & Ananthasuresh, 2008), electro statically (Français et al, 2005; Millet et al, 2004), thermally (Lai et al, 2004; Terre & Shkel, 2004) or electrical (Gomm et al, 2002; Huang & Lan, 2006) induced forces.Moreover, compliant MEMS having two or three clear stable states as named bi-stable or tri-stable behavior respectively were used in micro valve, micro switch, micro clasps applications (Chen et al, 2009; Jensen et al, 2001; Jensen & Howell, 2003; Nathan & Howell, 2003; Wilcox & Howell, 2005). For instance, Jensen designed several mechanisms such as double slider crank, slider-rocker mechanisms and explained the theory of bi-stable behavior (Jensen et al, 2004). Recent studies on compliant mechanisms are focused on novel designs (Kosa et al, 2010), new developed methodologies and optimization in topology (Chour & Jyhjei, 2006; Krishnan & Ananthasuresh, 2008; Pedersen & Seshia, 2004), size and shape (Krishnan & Ananthasuresh, 2008) or the use of finite element methods (Jensen et al, 2001). Compliant micro mechanisms enable mechanical or geometric benefit meaning that the ratio of output force to input force and the ratio of output displacement to input displacement, respectively, and both mechanical and geometric advantage (MA and GA, respectively) are formulized as follows;MA=Fout/FinGA=dout/dinThe energy is conserved during the motion transfer of compliant micro mechanism indicating that the increase in the output force causes decrease in the output displacement and vice versa. So, both mechanical and geometric benefits are significant to provide input to the micro actuators in MEMS applications (Kosa et al, 2010). Optimization of compliant mechanisms such as topology and size optimization is a challenging issue. In topology optimization, it is critical to design a suitable functional configuration of the mechanism to provide desired output motion under applied forces while in size optimization, it is important to achieve desired force or displacement amplification so as to operate under maximum loads (Kota et al, 2001).In this study, novel compliant MEMS force amplifier is designed and simulated by modeling its rigid body mechanism by Matlab/Simulink to determine the dynamic and quasi-static behavior. Kinematic approach is investigated and kinematic equations are derived and velocity and acceleration analysis of the micro mechanism are modeled. Dynamic response of MEMS amplifier is validated at a constant angular velocity and it is concluded that force amplification reaches to infinity at zero-crank angle. It is achieved that force amplification ratio reaches 5093, as the first stage crank angle, Θ2 passes from 0° in quasi-static simulation.2. Mechanism designCompliant MEMS force amplifier’s configuration is schematically shown in Fig. 2. Micro amplifier is composed of two slider-crank mechanisms. The two stage slider-crank amplifier provides force amplifying by means of its novel design. Its aim is to perform high output force at point B under low input forces. Two stages provide much more amplification compare to one stage. For both stages, rigid beams are linked by single thin flexible beams having a width of 3 µm. These flexible beams make the micro mechanism motion possible under operating forces. The micro mechanism stores energy and transfers force by elastic deformation of flexible beams linking rigid beams as both stage-slider cranks get close zero degree crank angle. Afterwards, input force is removed and micro amplifier springs back to its original position by means of flexible links having large deflections.The beams in first stage have a length of 100 µm and width of 25 µm as the beams in second stage have a length of 800 µm and width of 25 µm, as all beams have rectangular cross sectional area. The depths of all beams are chosen as 25µm limited by SOI-MUMPs (Silicon on Insulator Multi User MEMS Process) manufacturing technology (Cohen et al, 2009). 2.1. Grashof theoremIn rigid body model of the MEMS amplifier, four-bar configuration is attained after vector loop equations are derived. Grashof theorem becomes significant to demonstrate the act of micro mechanism. Grashof theorem takes three cases into consideration and states that when both of beams are rocked it is called a double-rocker when both of beams are able to revolve, then it is called double-crank, when the short beam is able to rotate as the long one is rocked, then it is called a crank-rocker mechanism. To determine the moving limit of the micro mechanism, the relation between the lengths of beams turns out to be an important issue. Therefore, selecting the length of a beam plays a crucial role for the micro mechanism. Due the fact that, x1, x2 are assumed as length of the shortest beam and length of the longest beam, respectively, as x3, x4 are the mean lengths of the beams. If x1+x2<= x3+x4, at least one of the beams can rotate and If x1+x2= x3+x4, the mechanism is activated and crank has limited rotation this feature enables beams to pass horizontal positions closely to each other achieving a high force amplifying.

Figure 2.Novel compliant MEMS Force Amplifier3. Analysis of quasi-static behavior3.1. Force and moment equation derivationRigid body model of the compliant micro mechanism is considered. Free body diagram of each beam is sketched and a typical beam model is schematically shown in Fig. 3. Forces acting on each beam is broken down into x and y components as follows;

Figure 3.Free body diagram of beam 2The static force and moment equations of beam 2 is typically shown and derived as; Equation derivation of forces acting on beam 2 along x axis;

$$$$Equation derivation of forces acting on beam 2 along y axis;

$$$$Equation derivation of moments acting on beam 2 along z axis;

$$$$Free body diagram of beam 3 is shown in Fig. 4 and equation derivation of forces acting on beam 3 along x axis;

$$$$Equation derivation of forces acting on beam 3 along y axis;

$$$$Equation derivation of moments acting on beam 3 along z axis;

$$$$Figure 4.Free body diagram of beam 3Free body diagram of beam 5 is shown in Fig. 5 and equation derivation of forces acting on beam 5 along x axis;

$$$$Equation derivation of forces acting on beam 5 along y axis;

$$$$Equation derivation of moments acting on beam 5 along z axis;

$$$$Free body diagram of beam 6 is shown in Fig. 6 and equation derivation of forces acting on beam 6 along x axis;

$$$$Equation derivation of forces acting on beam 6 along y axis;

$$$$Equation derivation of moments acting on beam 6 along z axis;

$$$$Figure 5.Free body diagram of beam 5

Figure 6.Free body diagram of beam 6Free body diagram of slider is shown in Fig. 7 and equation derivation of forces acting on slider along x- and y- axes;

$$$$$$$$Thus, 14 force and moment equations are derived. Equations of relation between internal forces of beams;

$$$$$$$$$$$$$$$$8 equations are derived from the relations between internal forces of beams.

Figure 7.Free body diagram of slider3.2. Equations and initial conditionsThe vector loop equations are derived and broken down into x and y components as force and moment equations. It is seen that linear matrix method could not be used to solve the position problem. To analyze the position behavior of the micro mechanism, nonlinear and transcendental equations should be solved by Matlab and in quasi-static run, initial conditions of Θ2 and (Θ6-90°) are chosen as 10° and 20°, respectively. 3.3. Position analysisThe micro mechanism is a single degree of freedom mechanism and position analysis provides to inform the positions of other links and points as one of the links moves or rotates.To find out position problem of the micro mechanism, nonlinear and transcendental vector loop equations that are derived and solved. The vector loops are schematically shown in Fig. 8. There are two vector loop equations such as;First vector loop equation:

$$Figure 8.Vector loops for the force amplifierDeriving equations according to coordinates of x and y:

$$$$Second vector loop equation:

$$Vector loop equations along x--axis

$$Vector loop equations along y--axis

$$By quasi-static analysis, it is claimed that (360°-Θ3) and Θ2 decreases linearly and are equal to each other during both quasi-static and dynamic simulations run by Matlab/Simulink. As seen in Fig. 9, it is calculated that as Θ5 goes from 70° to 74.0248°, Θ6 reduces from 110° to 105.9756°. Thus, as Θ2 rotates 20°, both Θ5 and Θ2 rotates approximately 4.02° and slightly different from each other. The relation both between Θ5 and Θ2, Θ6 and Θ2 are linear.Figure 9.Plot of Θ5 and Θ6 according to first stage crank angle, Θ2Displacement ratio is defined as Uoutput/Uinput. As the micro mechanism operates under an input force along – x direction, the first stage crank angle starts decreasing and pass from 0° and again starts increasing in an opposite direction and the ratio of output displacement to input displacement decreases as shown in Fig. 10. Beams 5 and 6 moves along –x and –y directions and the length of beams 5 and 6 are 8 times of beams 2 and 3. So, the input displacement increases rapidly than output displacement at close to zero degree crank angles. At negative crank angle values defining opposite directions, the slider gets close to initial position on contrast, beams 5 and 6 continue to get close to their vertical positions meaning that input displacement goes on to increase whereas output displacement begin to decrease. Therefore, after zero-crank angle, the displacement ratio continues to decrease according to Θ2.As the micro mechanism displays, both the second stage crank angle, (Θ6-90°) and the first stage crank angle, Θ2 get close to zero degree, the force amplification defined as Foutput/Finput starts increasing and when Θ2 is 0° and (Θ6-90°) is at about small values, the micro mechanism provides high output force and force amplifying sharply increases as seen in Fig. 11 under 1.7*10-7 in [N]. Also, there are two peaks in force amplification by quasi-static run. As, the first crank angle is close to zero but at still positive value, the force amplifying reaches 5093 and after that step first crank angle gets negative value but it is still close to zero, the force amplification ratio is 4830 at negative direction due to the fact that the slider motion begin to move in opposite direction and also, output force is in opposite direction. It is claimed that the toggle position of the micro mechanism is a very crucial issue meaning that if the initial conditions such as crank angles are adjusted properly to enable both crank angle pass 0° at the same time, the ratio of the output to the input force applied to the mechanism goes to infinity at zero degree crank angles. Figure 10.Plot of displacement ratio according to first stage crank angle, Θ2Figure 11.Plot of force amplifying according to first stage crank angle, Θ24. Dynamic behavior of a novel MEMS amplifier4.1. Inertial and geometric parameters:It is assumed that micro mechanism is made up of silicon having a density of 2.33 g/cm3. For short length of beams, lengths are 100 micron, widths and heights are 25 micron. The mass of short beams is;Ms=2.33*100*25*25*10-15=145625*10-18 [kg] For long length of beams, lengths are 800 micron, widths and heights are 25 micron. The mass of short beams is;Ml=2.33*800*25*25*10-15=1165000*10-18 [kg] The mass of the slider is accepted as 145625*10-18 in kilograms.The mass moments of inertia of the beams are calculated as follows; For short beams;Is = Ms *(L2+a2)/12=145625*10-18*(1002+252)/12=128938802.1*10-18 [kg*µm2] For long beams;Il = Ms *(L2+a2)/12=1165000*10-18 *(8002+252)/12=6.219401042*10-8[kg*µm2] 4.2. Kinematic behavior4.2.1. Velocity analysisKinematic simulation is used to calculate and to plot the velocities and acceleration of the beam of the MEMS amplifier.To understand kinematic behavior of the mechanism, first of all, derivatives of vector loop equations derived in position analysis are taken with respect to time and the velocity equations are arranged as follows;

$$$$$$$$The beam 6 are rotated at a constant speed, 0.01 rad/s, in clockwise direction and the initial conditions of w2, w3, w5,

$$are -0.059378175917485 [rad/s], 0.059378175917485 [rad/s], 0.011371580426033 [rad/s ], 2.062182408251533 [µm/s], respectively.Figure 12.Plot of angular velocity of beam 2 and beam 3 versus time for simulation of force amplifierThe angular velocities of beams 2 and 3 in 1 stage are equal to each other in magnitude. As w3 rotate counter clockwise direction, w2 rotate clockwise direction and the absolute values of the changes in w3 and w2 equal to each other according to time as shown in Fig. 12. Slider slows down until the first stage crank angle, Θ2 pass from 0°. When first stage beams are fully open, as having horizontal position, slider velocity is equal to zero. Then the slider moves to along -x direction and angular velocity of beam 5 decreases according to time as in Fig. 13.Figure 13.Plot of angular velocity of beam 5 and slider velocity versus time for simulation of force amplifier4.2.2. Acceleration analysisTo analyze the acceleration of the beams, second derivatives of the terms must be handled. The second derivatives of the vector loop equations for the micro mechanism are as follows;

$$$$$$$$In acceleration simulation by Simulink, the velocities such as w2, w3, w5,

$$ w6 are considered as known. The beam 6 rotates at a constant speed meaning that acceleration of beam 6 is zero. Acceleration of beam 2 and beam 3 are shown in fig. 14. Both acceleration of beams decrease as the micro mechanism operates under constant w6, angular velocity. The magnitude of acceleration of beam 2 and beam 3 are equal to each other during simulation. Also, as seen in Fig. 15, acceleration of beam 5 and slider decrease as function of time. Figure 14.Acceleration of beam 2 and beam 3 under constant angular acceleration, 𝛼2Figure 15.Acceleration of beam 5 and slider under constant angular acceleration, 𝛼24.3. Acceleration vector equations according to center of massThe linear acceleration of the center of mass equations are not present in vector loop equations that are previously derived. So, there must be equations relating to the acceleration of the center of mass of beams. Equation derivation is as follows and schematic representation of the center of mass acceleration in first and second loops is shown in Fig. 16 and Fig. 17, respectively.

Figure 16.The center of mass acceleration in first loopThe center of mass acceleration of beam 2 along x and y direction;

$$$$$$Figure 17.The center of mass acceleration in second loopThe center of mass acceleration of beam 3 along x and y direction;

$$$$$$The center of mass acceleration of beam 6 along x and y direction;

$$$$$$The center of mass acceleration of beam 5 along x and y direction;

$$$$$$4.4. Force and dynamic analysis of the micro mechanismThe micro mechanism operates under constant angular velocity, 0.01 [rad/s], the slider crank starts increasing and reaches to its maximum value, 200 micron, meaning that the first stage slider crank is fully opened at 3.20 sec. then the first crank angle pass from 0° and slider begins to get close to its initial position and R1 decreases as shown Fig. 18. According to both crank angles, Θ2 and (90-Θ6), the output force increases or decreases. In the first section of the Foutput vs. time curve, first, both crank angles decrease, and two slider cranks start to open and at small crank angles, Foutput sharply increase and at 3.20 sec Θ2 is equal to 0.0013° and at 3.25sec. Θ2 is equal to -0.0012° and at Θ2 these values, Foutput goes to its peak values such as -3.07*104 µN at 3.20 sec. and 3.32*104 µN at 3.25 sec. Θ2 decreases until 3.20 sec. and then it increases, whereas 90°-Θ6 decreases and gets close to small values during the simulation. The magnitude of first peak of Foutput at 3.25 sec. is higher than the magnitude of second peak of Foutput at 3.20 sec. due to the fact that (90°-Θ6) at 3.25 sec. is smaller than the value of (90°-Θ6) at 3.20 sec., meaning that small crank angle value of (90°-Θ6) contributes to get much more output force.Figure 18.Displacement of slider and output force versus time5. ConclusionThe MEMS force amplifier designed in this study is shown to provide high output to input ratio.By quasi-static analysis, 5090 force amplifying is achieved as the first crank angle, Θ2, rotates 10° and passes from its horizontal position and (90°-Θ6) rotates 1.85° and continues to decrease.The maximum amplifying ratio changes based on the initial position of the micro mechanism. So, the toggle of the micro mechanism has a crucial role to get high force output and high force amplification. If the mechanism’s initial position is adjusted properly as both crank angles pass 0° at the same time, the force output and consequently force amplification go to infinity.If pseudo rigid body of the compliant MEMS force amplifier having elastic hinges is modeled as a further study, it would provide us to get much more close response to the micro mechanism’s real behavior. This novel MEMS amplifier design achieves high force amplifying due to its geometric design.By dynamic analysis, high output force is achieved as the micro mechanism operates under 0.01 [rad/s] constant angular velocity of beam 6 at 3.20 sec and at about fully open position of first crank angle.By Simulink, the simulation displays dynamic behavior of the micro compliant mechanism and it is claimed that second stage crank angle rotates 4.01°, whereas first stage crank angle, Θ2, rotates 19.92°.5.1. NomenclautureAshokM.HowellL. L.NortonT. W.2000Limit Positions of Compliant Mechanisms Using the Pseudo-Rigid-Body Model Concept. Mechanism and Machine Theory, 3599115ChenG.WilcoxD. L.HowellL. L.2009Fully Compliant Double Tensural Tristable Micromechanisms (DTTM), Journal of Micromechanics and Microengineering, 1918CohenA.HamesG.MonkD.WilcenskiS.HardyB.2009Soimumps Design Handbook, MEMSCAP Revision 5.FrançaisO.RousseauL.BourouinaT.HaussyJ.TissotA.2005MEMS Memory Based on Bi-stable Mechanical Structures, Proceedings of DTIP of MEMS and MOEMS, 28481303571Switzerland, June, 2005GommT.HowellL. L.SelfridgeR. H.2002In-plane Linear Displacement Bistable Microrelay, Journal of Micromechanics and Microengineering, 12257264HanJ. S.MüllerC.WallrabeU.KorvinkJ. G.2007Design, Simulation, and Fabrication of a Quadstable Monolithic Mechanism with X- and Y-Directional Bistable Curved Beams. Journal of Mechanical Design, 12911982003HsuTai.Ran2002MEMS & Microsystems Design and Manufacture, McGraw Hill (1st edition), 007113051HuangS. C.LanG. J.2006Design and Fabrication of a Micro Compliant Amplifier with a Topology Optimal Compliant Mechanism Integrated with a Piezoelectric Microactuator. Journal of Micromechanics and Microengineering, 165315380960-1317JensenB. D.ParkinsonMatthew. B.KurabayashiKatsuo.HowellLarry. L.BakerMichael. S.2001Design Optimization of a Fully-Compliant Bistable Micromechanism, Proceedings of 2001 ASME International Mechanical Engineering Congress and Exposition, New York, November, 2001JensenB. D.HowellL. L.2003Identification of Compliant Pseudo-Rigid-Body Four-Link Mechanism Configurations Resulting in Bistable Behavior. Journal of Mechanical Design, 125701708JensenB. D.HowellL. L.2004Bistable Configurations of Compliant Mechanisms Modeled Using Four Links and Translational Joints. Journal of Mechanical Design. 126657666KosaE.SonmezU.KizilH.TrabzonL.2010The Design and Analysis of a Novel MEMS Force Amplifier. Turkish Journal of Engineering & Environmental Sciences, 34253259KotaS.JooJ.LiZ.RodgersS. M.SniegowskiJ.2001Design of Compliant Mechanisms: Applications to MEMS. Analog Integrated Circuits and Signal Processing, 29715KrishnanG.AnanthasureshG. K.2008Evaluation and Design of Displacement-Amplifying Compliant Mechanisms for Sensor Applications. Journal of Mechanical Design, 13019LaiY.Mc DonaldJ.KujathM.HubbardT.2004Force, Deflection and Power Measurements of Toggled Microthermal Actuators. Journal of Micromechanics and Microengineering, 144956LarryH. L.2001Compliant Mechanisms, John Wiley & Sons, 047138478USA.LiJ.LiuZ. S.LuC.ZhangQ. X.LiuA. Q.2005A Self Limited Large Displacement Ratio Micromechanical Amplifier. Proceedings of The 13th International Conference on Solid-State Sensors, Actuators and Microsystems, 0-78038-952-2Korea, June, 2005LobontiuN.PaineJ. S. N.GarciaE.GoldfarbM.2001Corner-Filleted Flexure Hinges. Journal of Mechanical Design, 123346352MilletO.BernardoniP.RégnierS.BidaudP.TsitsirisE.CollardDominique.BuchaillotLionel.2004Electrostatic Actuated Micro Gripper using an Amplification Mechanism, Sensors and Actuators A, 114371378NathanD.HowellLarry. L.2003A Self-Retracting Fully Compliant Bistable Micromechanism. Journal of Microelectromechanical Systems, 12273280ParkinsonM. B.JensenB. D.KurabayashiK.2001Design of Compliant Force and Displacement Amplification Micro-Mechanisms. Proceedings of DETC’01 ASME, Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Pittsburgh, Pennsylvania, September, 2001PedersenC. B. W.SeshiaA. A.2004On the Optimization of Compliant Force Amplifier Mechanisms for Surface Micromachined Resonant Accelerometers. Journal of Micromechanics and Microengineering, 1412811293SreekumarM.NagarajanT.SingaperumalM.2008Experimental Investigations of the Large Deflection Capabilities of a Compliant Parallel Mechanism Actuated by Shape Memory Alloy Wires. Smart materials and structures, 17112TantanawatT.KotaS.2007Design of Compliant Mechanisms for Minimizing Input Power in Dynamic Applications. Journal of Mechanical Design, 12910641075TerreJ. C.ShkelA.2004Dynamic Analysis of a Snap-Action Micromechanism, IEEE, 12451248WilcoxD. L.HowellL. L.2005Fully Compliant Tensural Bistable Micromechanisms (FTBM). Journal of Microelectromechanical Systems, 1418