沖壓機(jī)床液壓控制系統(tǒng)設(shè)計(jì)

沖壓機(jī)床液壓控制系統(tǒng)設(shè)計(jì),沖壓,機(jī)床,液壓,控制系統(tǒng),設(shè)計(jì) 南京理工大學(xué)泰州科技學(xué)院 畢業(yè)設(shè)計(jì)（論文）任務(wù)書 系　　　部： 機(jī)械工程系 專 業(yè)： 機(jī)械工程及自動化 學(xué) 生 姓 名： 陸輝 學(xué) 號： 05010127 設(shè)計(jì)(論文)題目： 沖壓機(jī)床液壓控制系統(tǒng)設(shè)計(jì) 起 迄 日 期: 2009年3 月 09 日 ~ 6 月 14日 設(shè)計(jì)(論文)地點(diǎn): 南京理工大學(xué)泰州科技學(xué)院 指 導(dǎo) 教 師: 張衛(wèi) 專業(yè)負(fù)責(zé)人: 龔光容 發(fā)任務(wù)書日期: 2009 年 2 月 26 日  任務(wù)書填寫要求 1．畢業(yè)設(shè)計(jì)（論文）任務(wù)書由指導(dǎo)教師根據(jù)各課題的具體情況填寫，經(jīng)學(xué)生所在專業(yè)的負(fù)責(zé)人審查、系部領(lǐng)導(dǎo)簽字后生效。此任務(wù)書應(yīng)在第七學(xué)期結(jié)束前填好并發(fā)給學(xué)生； 2．任務(wù)書內(nèi)容必須用黑墨水筆工整書寫或按教務(wù)處統(tǒng)一設(shè)計(jì)的電子文檔標(biāo)準(zhǔn)格式（可從教務(wù)處網(wǎng)頁上下載）打印，不得隨便涂改或潦草書寫，禁止打印在其它紙上后剪貼； 3．任務(wù)書內(nèi)填寫的內(nèi)容，必須和學(xué)生畢業(yè)設(shè)計(jì)（論文）完成的情況相一致，若有變更，應(yīng)當(dāng)經(jīng)過所在專業(yè)及系部主管領(lǐng)導(dǎo)審批后方可重新填寫； 4．任務(wù)書內(nèi)有關(guān)“系部”、“專業(yè)”等名稱的填寫，應(yīng)寫中文全稱，不能寫數(shù)字代碼。學(xué)生的“學(xué)號”要寫全號； 5．任務(wù)書內(nèi)“主要參考文獻(xiàn)”的填寫，應(yīng)按照國標(biāo)GB 7714—2005《文后參考文獻(xiàn)著錄規(guī)則》的要求書寫，不能有隨意性； 6．有關(guān)年月日等日期的填寫，應(yīng)當(dāng)按照國標(biāo)GB/T 7408—2005《數(shù)據(jù)元和交換格式、信息交換、日期和時(shí)間表示法》規(guī)定的要求，一律用阿拉伯?dāng)?shù)字書寫。如“2009年3月15日”或“2009-03-15”。 畢 業(yè) 設(shè) 計(jì)（論 文）任 務(wù) 書 1．本畢業(yè)設(shè)計(jì)（論文）課題應(yīng)達(dá)到的目的： 通過本課題旨在讓學(xué)生了解沖壓機(jī)床液壓控制系統(tǒng)的設(shè)計(jì)方法，并運(yùn)用AMEsim軟件對液壓控制系統(tǒng)進(jìn)行數(shù)字建模與動態(tài)仿真，幫助企業(yè)實(shí)現(xiàn)對沖壓機(jī)床液壓系統(tǒng)的優(yōu)化。 通過該課題可以培養(yǎng)學(xué)生解決實(shí)際問題的能力，提高學(xué)生綜合應(yīng)用能力和獨(dú)立工作能力，為學(xué)生最終走向工作崗位打下基礎(chǔ)。 2．本畢業(yè)設(shè)計(jì)（論文）課題任務(wù)的內(nèi)容和要求（包括原始數(shù)據(jù)、技術(shù)要求、工作要求等）： 本課題針對企業(yè)需求，首先研究3150KN四柱式?jīng)_壓機(jī)床液壓控制系統(tǒng)的組成，進(jìn)而基于AMESim進(jìn)行沖壓機(jī)床液壓控制系統(tǒng)設(shè)計(jì)并對關(guān)鍵部件進(jìn)行數(shù)字建模仿真，并完成動態(tài)仿真實(shí)驗(yàn)，實(shí)現(xiàn)幫助企業(yè)實(shí)現(xiàn)沖壓機(jī)床液壓系統(tǒng)優(yōu)化的目的。 原始數(shù)據(jù)具體要求如下： 3150KN液壓機(jī)基本參數(shù)為：公稱力：3150KN，頂出力：630KN，液體最大工作壓力：25MPa，滑塊行程：800mm，頂出行程：300mm，開口高度：1250mm，工作臺面有效尺寸左右×前后：基型1120×1120mm2 技術(shù)要求具體要求如下： 滑塊行程速度：快速下行V1 150mm/s，慢速下行V2 25mm/s，工作V310mm/s，回程V4 80mm/s，頂缸頂出及退回速度V5 50mm/s 工作要求具體要求如下： 1) 對3150KN四柱式液壓控制系統(tǒng)進(jìn)行需求分析并進(jìn)行總體設(shè)計(jì)； 2) 進(jìn)行機(jī)床液壓控制系統(tǒng)系統(tǒng)的詳細(xì)設(shè)計(jì)，計(jì)算并確定液壓系統(tǒng)各部分的尺寸，其中主要包括油缸尺寸和液壓泵規(guī)格的確定和各個(gè)閥的選型； 3) 運(yùn)用AMESim軟件對所設(shè)計(jì)的機(jī)床液壓控制系統(tǒng)進(jìn)行仿真調(diào)試，解決換向壓力沖擊大和回程速度慢的問題； 4) 按要求提供設(shè)計(jì)文檔； 5) 畢業(yè)設(shè)計(jì)論文。 畢 業(yè) 設(shè) 計(jì)（論 文）任 務(wù) 書 3．對本畢業(yè)設(shè)計(jì)（論文）課題成果的要求〔包括畢業(yè)設(shè)計(jì)論文、圖表、實(shí)物樣品等〕： 1. 機(jī)床液壓控制系統(tǒng)圖及仿真程序 2. 畢業(yè)設(shè)計(jì)論文 4．主要參考文獻(xiàn)： [1] 郭成，儲家佑，等. 現(xiàn)代沖壓技術(shù)手冊[M]. 北京: 中國標(biāo)準(zhǔn)出版社，2005. [2] 程燕軍，柳舟通，等. 沖壓與塑料成型設(shè)備[M]. 北京: 科學(xué)出版社，2005. [3] 李永堂，雷步芳，高雨茁. 液壓系統(tǒng)建模與仿真[M]. 北京: 冶金工業(yè)出版社，2003. [4] 張利平. 液壓控制系統(tǒng)及設(shè)計(jì)[M]. 北京：化學(xué)工業(yè)出版社，2006. [5] 李華聰，李吉. 機(jī)械/液壓系統(tǒng)建模仿真軟件AMEsim[J]. 計(jì)算機(jī)仿真, 2006，(12): 294-298. [6] 王平, 葉曉葦. 沖床加工設(shè)備及自動化[M]. 武漢：華中科技大學(xué)出版社，2006. [7] 俞新陸，楊津光. 液壓機(jī)的結(jié)構(gòu)與控制[M]. 北京：機(jī)械工程出版社，1997. [8] 俞新陸. 液壓機(jī)的設(shè)計(jì)與應(yīng)用[M]. 北京: 機(jī)械工業(yè)出版社, 2006. [9] 付永領(lǐng)，祁曉野. AMEsim系統(tǒng)建模和仿真—從入門到精通[M]. 北京：北京航天航空大學(xué)出版社，2006. [10] 蘇東海，于江華.液壓仿真新技術(shù)AMEsim及應(yīng)用[J].機(jī)械.2006，(11):35-37 [11] 盧長耿, 李金良. 液壓控制系統(tǒng)的分析與設(shè)計(jì)[M]. 北京:煤炭工業(yè)出版社, 1991. [12] 劉海麗，李華聰.液壓機(jī)械系統(tǒng)建模仿真軟件AMEsim及其應(yīng)用[J].機(jī)床與液壓:2006，(6):124-126. [13] 邊海岸, 戴雙獻(xiàn). 四柱式萬能液壓機(jī)液壓系統(tǒng)原理研究[J].中國機(jī)械工程師, 2007, (7): 92-93. 畢 業(yè) 設(shè) 計(jì)（論 文）任 務(wù) 書 5．本畢業(yè)設(shè)計(jì)（論文）課題工作進(jìn)度計(jì)劃： 起 迄 日 期 工 作 內(nèi) 容 2009年 3月9 日 ～3 月29日 3月30日～4 月12 日 4月13日～ 4 月 27日 4月28日～5 月 16日 5月17日～5 月 31日 6月1 日～6 月 5 日 完成外文資料翻譯、文獻(xiàn)綜述和開題報(bào)告 完成3150KN四柱式?jīng)_壓機(jī)床液壓控制系統(tǒng)總體方案，設(shè)計(jì)計(jì)算并確定液壓系統(tǒng)各部分的尺寸 基于AMESim進(jìn)行沖壓機(jī)床液壓控制系統(tǒng)數(shù)字建模 完成動態(tài)仿真實(shí)驗(yàn)，解決換向壓力沖擊大和回程速度慢的問題 完成畢業(yè)論文，準(zhǔn)備畢業(yè)答辯 論文答辯 所在專業(yè)審查意見： 負(fù)責(zé)人： 年 月 日 系部意見： 系部主任： 年 月 日 mented in the AMESim simulation tool. Body and joint components are the basic components 1569-190X/$ - see front matter C211 2005 Elsevier B.V. All rights reserved. * Corresponding author. Tel.: +33 4 72 43 85 58. E-mail address: wilfrid.marquis-favre@insa-lyon.fr (W. Marquis-Favre). Simulation Modelling Practice and Theory 14 (2006) 25–46 of this library. Due to the library philosophy requirements, the mathematical models of the components have required a generic vector calculus based formulation of the constraint equa- tions. This formulation uses a set of dependent generalized coordinates. The dynamics equa- tions are obtained from the application of JourdainC213s principle combined with the Lagrange multiplier method. The body component mathematical models consist of di?erential equations in terms of the dependent generalized coordinates. The joint component mathematical models are based on the Baumgarte stabilization schemes applied to the geometrical, kine- matic and acceleration constraint equations. The Lagrange multipliers are the implicit solution of these Baumgarte stabilization schemes. The first main contribution of this paper is the expression of geometrical constraints in terms of vectors and their exploitation in this form. The second important contribution is the adaptation of existing formulations to the AMESim philosophy. C211 2005 Elsevier B.V. All rights reserved. A planar mechanical library in the AMESim simulation software. Part I: Formulation of dynamics equations Wilfrid Marquis-Favre * , Eric Bideaux, Serge Scavarda Laboratoire d’Automatique Industrielle, Institut National des Sciences Applique′es de Lyon, Ba?t. St Exupe′ry, 25, avenue Jean Capelle, F-69621 Villeurbanne Cedex, France Received 25 March 2003; received in revised form 17 December 2004; accepted 8 February 2005 Available online 17 March 2005 Abstract This paper presents the mathematical developments of a planar mechanical library imple- doi:10.1016/j.simpat.2005.02.006 domains. One can now carry out modeling, analysis and simulation for systems con- sisting of pneumatic, powertrain, hydraulic resistance, thermal, electromagnetic and 26 W. Marquis-Favre et al. / Simulation Modelling Practice and Theory 14 (2006) 25–46 cooling components for instance. The restriction to only one-dimensional motion for the mechanical components has motivated the development of a two-dimensional mechanical library. Keywords: AMESim; Planar mechanics; Dynamics equations; Constraint equations; Lagrange multi- pliers; Baumgarte stabilization 1. Introduction This paper, organized in two parts, presents a new library for the simulation tool AMESim [2]. The first part is dedicated to the theoretical developments of the library. The second part shows the composition of the library as it was primarily implemented in AMESim and illustrates it with an application example of a seven-body mechanism. This library proposes components belonging to the planar mechanical domain. The objective with this library was not to compete with multi- body system software tools that are better adapted to this domain. The objective was more to enlarge the range of industrial applications capable of being treated by AMESim. From a theoretical point of view the challenge of implementing this library was to fit existing mechanical formulations to the inherent requirements of AMESim philosophy. The solution has been found by adapting the dynamic equa- tions expressed from JourdainC213s principle and the Lagrange multiplier method together with BaumgarteC213s stabilization. Also a generic feature of the formulation has been researched over the library components (bodies and joints) and one key contribution of this paper is concerned with this generic feature. Basically the formu- lation consists of expressing the geometric constraints associated with joints in terms of vectors and carrying out the developments of this form. The result is the set up, for kinematic and acceleration constraints, of a unique expression that fits every joint presented in the library. The generic feature of the formulation proposed thus enables the derivation of joint contraints to be systematized. One can then imagine a new joint with its corre- sponding vector constraint and derive straightforwardly the corresponding mathe- matical model by applying the proposed formulation. Also, in the context of predefined component models, the given formulation clearly shows the frontiers of the di?erent mathematical models in terms of inputs and outputs. Therefore it also helps to define in which models output equations must be implemented. Also, the formulation proposed intrinsically enables closed loop structures to be dealt with. AMESim (for Advanced Modeling Environment for performing Simulations) is organized in component libraries. The components, represented by symbolically technologically suggested icons, can be interconnected exactly like the system under study. AMESim was first applied to electrohydraulic engineering systems with simple one-dimensional mechanical systems (like inertia, springs, and dampers in transla- tion or in rotation). It recently opened its libraries to a variety of other component ally reduced and a number of constraints are a priori symbolically eliminated. Like- W. Marquis-Favre et al. / Simulation Modelling Practice and Theory 14 (2006) 25–46 27 wise, tools based on bond graph (e.g. 20Sim [1] or MS1 [18]) can deal with multibody systems in a pluridisciplinary context (e.g. [4,7]). The essential feature of bond graph language is its ability to describe the energy topology of a model at an acausal level. This enables all the model variables to be globally assigned and all the equations to be globally organized. This also eliminates superfluous dependencies of the multi- body models. Section 2 presents an overview of some multibody codes and object-oriented tools, as well as the environmental requirements of AMESim. These requirements have some implications on how the 2D library is built. Section 3 details the theoret- ical developments that enabled the mathematical models of the library components to be set up. Section 4 concludes this first part. 2. Constraints of AMESim library philosophy After a brief overview of multibody code principles and some object-oriented tools, a presentation of AMESim requirements is given. Concerning multibody codes a state of the art is given by [23]. Details are not reproduced here and readers are referred to this book for a more profound presen- tation. Although more than a decade has passed and certain tools are no longer developed and others have changed, this state of the art book gives a good idea of the main principles that can be used as a basis for multibody codes. Also this over- view enables the library proposed to be positioned with respect to these codes. There are di?erent approaches for writing dynamic equations. The approaches most repre- sented in multibody codes are, the Newton–Euler equations applied to each body, the Newton–Euler equations applied to sets of bodies, LagrangeC213s equations and KaneC213s equations [13,14]. The variables, in whose terms the dynamic equations are written, are either absolute coordinates or relative coordinates. Also supplementary methods are used for reducing the index of the Di?erential–Algebraic Equations. The principal ones are the coordinate partitioning method, the projection matrix method, the Baumgarte stabilization and the penalty formulation [9]. The first two methods aim at working with a set of independent generalized coordinates while the Baumgarte stabilization enables the constraints, together with the di?erential equations, to be handled and the penalty formulation increases the di?erential sys- tem order by introducing extra dynamics into the model. In the domain of the object-oriented tools to which AMESim may be attached, certain enable multibody systems to be treated with a di?erent approach to the mod- elling. For instance Dymola [21] is, like AMESim, based on well-identified techno- logical components in a pluridisciplinary context but it sets up the mathematical model in a di?erent way. Basically each component model consists of equations not oriented in terms of variable assignments nor organized a priori. Then, at the component connection stage, all the mathematical models are gathered in an implicit form and the compilation carries out the variable assignments and the organization of the equations in a consistent manner. Thus the order of the whole model is glob- It is now important to show the key features of AMESim to justify how the planar mechanical library has been implemented. Its feature oriented towards engineering systems and its user friendliness make AMESim work with well-identified technolog- ical components, symbolically manipulated by means of icons. These icons are inter- connected, one to the other and identically to the engineering system architecture under study. Fig. 1 shows an example of a door locking system using a permanent magnet modelled in AMESim. The icons displayed here belong to the magnetic, mechanical and signal libraries. This simple example shows the coupling between mechanical and magnetic domains where one circuit, fed by a permanent magnet (right-hand side magnetic circuit), is forced to move with respect to another passive circuit (left-hand side circuit). The main components consist of a permanent magnet (rectangle with a compass needle inside), three magnetic circuit parts characterized 28 W. Marquis-Favre et al. / Simulation Modelling Practice and Theory 14 (2006) 25–46 by a certain reluctance (rectangles with C212squareC213 ports with a diagonal cross inside), two variable air-gaps (vertical twin rectangles), two mechanical nodes (both sides of the air-gap components), a signal generator with a signal-to-displacement converter (in the centre of the right-hand side circuit), and a component for the set of the mag- netic medium characteristics (B–H diagram in a circle). Each component can be associated with one model from a set of component compatible mathematical models. As soon as the model has been chosen the component conserves this mathematical model. Contrary to acausal tools, AMESim works with component models that have equations both a priori oriented in terms of variable assignments and organized. This feature requires implementing new models in a predefined calculus scheme. Also the mathematical formulation of a component model has to be organized in order to fit into other potential component connections. So each component associated with a mathematical model has a predetermined set of input and output variables. It can thus be considered as a causal model. The connection of the components enables the exchange of those variables on the way out a component for those variables that are calculated by its mathematical model (outputs) and, the exchange of those vari- ables on the way in a component for those variables that are calculated by a con- nected component mathematical model (inputs). This causal feature of AMESim philosophy is the main constraint when implementing new components. This di?ers Fig. 1. Example of an AMESim model representation. from other object-oriented tools, based on acausal component models or acausal phenomenon models, like Dymola, or tools with a bond graph input (e.g. 20-Sim or MS1). Fig. 2 gives an example of two components in the mechanical (a mass in transla- tion) and the hydraulic (a two way hydraulic pump) domains respectively. The con- necting ports of the components show the variables exchanged by them and especially the outputs (C212exitingC213 arrows) passed to the connected components and the inputs (C212enteringC213 arrows) received from the connected components. These con- necting ports are intimately associated with power ports since two of the variables W. Marquis-Favre et al. / Simulation Modelling Practice and Theory 14 (2006) 25–46 29 exchanged at these ports are power variables. Fig. 2 examples illustrate two key features of a library oriented simulation tool. The first one is the domain port concept. It shows how AMESim can deal with plu- ridisciplinary systems. The second feature is the connecting port constraints. Since one component mathematical model requires given inputs to then calculate its state and its outputs, not all combinations of the component connections are allowed. For instance the Fig. 2 examples cannot be connected one to the other by any port. How- ever a mass component may be connected to a spring component or a damper component. Another key feature of a library oriented simulation software tool is the modular- ity concept. This often results in symmetrical components with respect to their con- necting port. This symmetry property, though not generalized to all components in AMESim, has been adopted for the planar mechanical library. The reason will ap- pear obvious when components of this library are presented. In the context of planar mechanisms and rigid bodies the library is not restricted to any mechanical domain application. The library also accepts closed loop struc- tures. Although relative coordinates are generally more e?cient for dynamic equa- tion formulation, AMESim philosophy requires the use of absolute coordinates. The absolute coordinates of the mass center have been chosen for each body. Nev- ertheless the planar feature of the library does not require any specific variables for the body orientation. Thus the absolute angular position has been chosen for each body as well. Once again, due to AMESim philosophy, the equations of the compo- nents cannot be globally reorganized when the components are connected. This for- bids the use of the coordinate partitioning method or the projection method to decrease the index of the Di?erential–Algebraic Equation systems. For this reason the Baumgarte stabilization has also been used in the library. Fig. 2. Example of two AMESim components. 3. Theoretical developments of the library components As has been explained in the previous section the library must be organized in 30 W. Marquis-Favre et al. / Simulation Modelling Practice and Theory 14 (2006) 25–46 with m i the body mass and g the gravity acceleration. We consider here ~y 0 as the ascendant vertical axis. The star superscript indicates virtual quantities. The coe?- cients of the virtual velocities in A* are derived from the kinetic coenergy (e.g. [6]) of a body by the equation: A q ? doT dto_q C0 oT oq with q a generalized coordinate e3T 1 A nomenclature is given in Appendix A. well-identified technological components. It has been decided to base the planar mechanical library on a body component and on joint components. The body com- ponent is associated with a supposed rigid material item of a mechanism. Its behav- ior is essentially governed by its kinetic state. The joint components are associated with the abstract items that represent the attachment of bodies in a mechanism. They are supposed to be ideal and their mathematical model is based on the constraints that they impose on the connected bodies. 3.1. Body component mathematical model The mathematical model of the body component is based on JourdainC213s Principle formulation (e.g. [5,23]) 1 : A C3 ? P C3 e1T where A* is the virtual power developed by the acceleration quantities and P* the virtual power developed by the actions on the body. In the library philosophy there is no a priori privileged candidate for the role of the generalized coordinates. For a planar motion, the generalized coordinates, which have been chosen, are the absolute mass center coordinates projected onto the absolute frame of reference ex G i ;y G i T and the absolute angular position h i (Fig. 3). This choice enables the more general case of a body motion to be dealt with. The body motion restriction will be determined by the joint constraints, as shown later. With this choice of generalized coordinates ex G i ;y G i ;h i T Eq. (1) members may now be written A C3 ? A x _x C3 G i tA y _y C3 G i tA h _ h C3 i P C3 ? Q x _x C3 G i teQ y C0m i gT_y C3 G i tQ h _ h C3 i e2T W. Marquis-Favre et al. / Simulation Modelling Practice and Theory 14 (2006) 25–46 31 Applied to Fig. 3 body in planar motion these quantities are written simply: A x ? m i €x G i A y ? m i €y G i with I i the body moment of inertia around eG i ;~z 0 T A h ? I i € h i e4T Q x , Q y , and Q h are the generalized forces including the constraint actions resulting from the fact that x G i , y G i , and h i are not necessarily independent after the connection of a body component to a joint component. From Eq. (1) and by taking a compat- ible virtual transformation with the joints as they exist at time t, we can now write the three identities that constitute the formulation basis for the body components. These three identities are m i €x G i ? Q x m i €y G i ? Q y C0m i g I i € h i ? Q h e5T Fig. 3. Schema of a body in planar motion. This formulation requires that the expression of the three generalized forces Q x , Q y , and Q h be further developed in order to fit any potential connected joint component. First let us inspect the case of a body with only one connecting port at a point M. Let us also consider simply a given action on the body characterized by a wrench about point M (e.g. [17]): fWg : ~ F ? F x ~x 0 tF y ~y 0 force ~ MeMT?C z ~z 0 torque about point M ( e6T The virtual power developed by this action is P C3 ? ~ F C1 ~ V 0 C3 eMTt ~ MeMTC1 ~ X 0 C3 i e7T where ~ V 0 C3 eMT is the virtual absolute velocity of point M and ~ X 0 C3 i is the virtual abso- lute angular velocity of the body. The velocity transport (e.g. [11]) enables Eq. (7) to be written as: 32 W. Marquis-Favre et al. / Simulation Modelling Practice and Theory 14 (2006) 25–46 3.2. Joint component mathematical model First a general formulation is given for the joint component mathematical model. It is then illustrated in the example of a translational joint. Let us consider this time two bodies connected by a joint. By the only fact that both bodies are connected (a joint component between two body components) their generalized coordinates (x G i , y G i and h i for body i and x G j , y G j and h j for body j) are no longer independent. In the library philosophy the constraint equations are ex- pressed in the joint component, which in turn furnishes the constraint actions to the body components. These constraint actions correspond to the variables F x , F y , and C z previously presented and passed to each body component. The general expressions of these variables are now determined. The joints considered in the planar mechanical library generate only geometrical constraints. These constraints may be expressed in a general way in an implicit form by Eq. (10) (e.g. [15]). g k eq 1 ; ...;q n T?0 for k ? 1tom e10T with n the number of generalized coordinates involved in the constraints and m the constraint number. It is supposed here that the constraints are scleronomic [16], which means that time does not explicitly appear in the constraint equations. At the kinematic level these equations become P C3 ? ~ F C1 ~ V 0 C3 eG i Tt ~ F C1 ~ X 0 C3 i C2G i M C131C131! C16C17 t ~ MeMTC1 ~ X 0 C3 i ? ~ F C1e_x C3 G i ~x 0 t _y C3 G i ~y 0 Tt ~ MeMTtG i M C131C131! C2 ~ F C16C17 C1 _ h C3 i ~z 0 e8T From Eq. (8) we can clearly identify the generalized forces used in the dynamic for- mulation of a body component: Q x ? ~ F C1~x 0 ? F x Q y ? ~ F C1~y 0 ? F y Q h ? ~ MeMTtG i M C131C131! C2 ~ F C16C17 C1~z 0 ? C z t ~ F C1e~z 0 C2G i M C131C131! T e9T Since G i M C131C131! is a characteristic vector of the body, the variables F x , F y , and C z , char- acterizing the given force at point M, are the only variables passed to the body at the connecting port. The variables Q x , Q y , and Q h are calculated in the body component model on Eq. (9) basis. It is shown in the next section that th