# An object oriented approach to modelling of flexible multibody system: focus on joint constrains.

1. INTRODUCTION

Attention of this paper is paid to the object-oriented modelling of the flexible multibody systems with emphasis to modelling of joint constrains. Physical modelling of such complex system is obtained using bond graphs (Karnopp et al. 2000; Damic&Montgomery, 2003). Bond graphs provide modelling of multi-domain systems (electrical, hydraulic, mechanical, etc.) on the same way therefore impose a strong tool in the modelling of mechatronic system. The models of different types of joints, which connect components of the flexible bodies, was developed. This is achieved using object-oriented software environment of BondSim[R] (Damic &Montgomery, 2003). The methods used differ from others in that it is based on a velocity formulation and direct integration of the resulting equations of the motion. The models are computer generated in form of differential-algebraic equations--DAEs and solved using a backward differential formula BDF. Two numerical simulations dealing with flexible multibody systems are performed in order to verify proposed procedure. The results obtained are in good agreement with the reported in the literature and show how model of complex system can be developed efficiently and with good accuracy.

2. BOND GRAPH MODEL OF FLEXIBLE LINKS

The flexible multibody system can be represented as a collection of elementary components. Flexible components are considered as long slender beams whose models are built as aggregation of bond graph 3D beam finite element (FE) component models (Beam3D), Fig.1a. 2D bond graph component has been developed in (Damic,2006; Cohodar et al., 2009.). To develop basic 3D beam component model the co-rotation formulation is applied, (Battini&Pacoste, 2002; Damic&Cohodar,2006). Fig. 1b shows configuration of the basic beam 3D FE. It uses three coordinate frames:

1. Global (inertial) frame, represented by triads (I,J,K),

2. Co-rotation frame (i,j,k) defined on such way that moving of beam finite element with respect to the frame consists only deformation and

3. Section coordinate frame (l,m,n), rigidly attached to the element cross section.

[FIGURE 1 OMITTED]

Two element nodes of beam finite element are represented by two ports. Flows and efforts ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]) at ports are linear and angular velocities and corresponding forces and moment, respectively, defined in the global frame.

3. BOND GRAPH MODEL OF JOINT CONSTRAINS

The articulation between two consecutive links of robot manipulator can be realized by one of two basic joints: a prismatic or a revolute joint. Other type of joints, for instance cylindrical can be modelled as combination of prismatic and revolute joints. Bond graph model of joint can be represented by component Joint, Fig. 2a with two power ports, described by pair vectors--flow and effort. Let joint (i) connect link to link (i), Fig.2b. The flows at left and right ports are represented by vectors composed of linear and angular velocities of link [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and link [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], respectively (expressed in the global frame). If joint i is revolute the angular velocity of link i with respect to link (i-1) can be expressed in the joint coordinate frame ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]) as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

where [[theta].sub.i] is angle of rotation. In the case of the prismatic joint (with joint variable [d.sub.i]) relative linear velocity is given by:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (2)

[FIGURE 2 OMITTED]

Development of revolute joint model is based on the assumption that it is rigidly attached to the previously beam element. Thus, joint orientation can be defined through the orientation of the end right cross-section (denote with index 2 in Fig.1a) using rotation matrix [R.sub.2] = [[l.sub.2] [m.sub.2] [n.sub.2]].

The joint orientation with respect to the right end crosssection does not change and can be represented by a constant matrix [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The joint orientation with respect to inertial frame is defined by:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3) For the revolute joint, relationship between angular velocities at the output and input ports is defined by:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

For prismatic joint linear velocity at the output port can be expressed as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

4. NUMERICAL EXAMPLES

To demonstrate and verify applicability of proposed models two numerical examples are considered in this paper: slider-crank mechanism, and two-component robot arm. In both simulations the time step is 1e-3 s and error tolerance 1e-6.

4.1 Slider-crank mechanism

To validate the models a flexible slider-crank mechanism is adopted as the first numerical example. Bond graph model of the slider-crank mechanism (Fig.3a) consists of a slider block, crankshaft and connected rod that are connected by the revolute joint. Crankshaft performs only a rotation and is connected to the ground by a revolute joint. Motion of the other side of connected rod is restricted to the horizontal siding plane. The slider block is modelled as a massless particle. The lengths of the crankshaft and connected rod are 0.5 m and 1 m, respectively. The dimensions of cross-section area of both components are the same: H=50 mm and W=5 mm, as well as are the material properties: modulus elasticity is E = 2.1e10 Pa, density p = 7.8e3 kg/[m.sup.3]. The mechanism is initially in the horizontally and straighten position and performs the motion under the influence of its own weight. The mid-point deflection of connection rod is shown in Fig.3b and is in a good agreement with results reported in (Gerstmayr&Schoberl, 2006).

[FIGURE 3 OMITTED]

4.2 Two-component robot arm

The second numerical example is the two-component robot arm, taken from (Ibrahimbegovic&Mamouri, 2000) and (Kromer et al., 2004). Fig. 4 presents geometry and physical properties of the robot structure. Links are interconnected by a spherical joint. The left side of the first link is connected to the base by a cylindrical joint and subjected to the prescribed displacement and rotation as shown in Fig.4a. Both links are discretized by ten finite elements components.

[FIGURE 4 OMITTED]

Simulation interval is set to 3 s. Displacement components of free end link are depicted in Fig.5. The results agree well with that obtained in (Ibrahimbegovic&Mamouri, 2000) and (Kromer et al., 2004).

[FIGURE 5 OMITTED]

5. CONCLUSION

In this paper, systematic development of the model of a constrained flexible multibody system has been demonstrated where special attention has been focused on modelling of the joint constrains. Due to component model approach each developed component can be put in the library and taken out of it when developing new systems. Several numerical simulations are carried out to test accuracy of developed model and their results are compared with the results of other authors showing good agreements. The proposed method can be successfully applied to modelling and simulation of other mechatronics systems as well.

6. REFERENCES

Battini, M. & Pacoste, C. (2002). Co-rotational beam elements with warping effects in instability problems, Comput. Methods Appl. Mech. Engrg., 191: 1755-789

Cohodar, M. Borutzky, W. & Damic, V. (2009). Comparison of different formulation of 2D beam elements based on Bond Graph technique, Simulation Modelling Practice and Theory, 17:107-124

Damic V. & Cohodar, M. (2006). Bond Graph Based Modelling and Simulation of Flexible Robotic Manipulators, in: Proc. 20th European Conference on Modelling and Simulation ECMS, Bonn

Damic, V. & Montgomery, V. (2003). Mechatronic by Bond Graphs, Springer--Verlag, ISBN 3-540-42375-3, Berlin Heidelberg

Damic, V. (2006). Modelling flexible body systems: a bond graph component model approach. Mathematical and Computer Modelling of Dynamical Systems, No.12, (175-187)

Gerstmayr, J. & Schoberl, J. (2006). A 3D finite element methods for flexible multibody systems. Multibody Syst. Dyn, 15., (309-324), DOI 10.1007/s11044-006-9009-3

Ibrahimbegovic, A. & Mamouri, S. (2000). On rigid components and joint constraints in nonlinear dynamics of flexible multibody systems employing 3D geometrically exact beam model, Computer methods in applied mechanics and engineering, 188: 805-831

Karnopp, D.C., Margolis, D.L.&Rosenberg, R.C. (2000) System dynamic, Modelling and Simulation of Mechatronic System, John Wiley&Sons, Inc., ISBN:0-471-33301-8

Kromer, V. Dufosse, F. & Gueury, M. (2004) An object-oriented design of a finite element code: application to multibody system analysis. Advanced in Engineering Software, 35, 273-287

Attention of this paper is paid to the object-oriented modelling of the flexible multibody systems with emphasis to modelling of joint constrains. Physical modelling of such complex system is obtained using bond graphs (Karnopp et al. 2000; Damic&Montgomery, 2003). Bond graphs provide modelling of multi-domain systems (electrical, hydraulic, mechanical, etc.) on the same way therefore impose a strong tool in the modelling of mechatronic system. The models of different types of joints, which connect components of the flexible bodies, was developed. This is achieved using object-oriented software environment of BondSim[R] (Damic &Montgomery, 2003). The methods used differ from others in that it is based on a velocity formulation and direct integration of the resulting equations of the motion. The models are computer generated in form of differential-algebraic equations--DAEs and solved using a backward differential formula BDF. Two numerical simulations dealing with flexible multibody systems are performed in order to verify proposed procedure. The results obtained are in good agreement with the reported in the literature and show how model of complex system can be developed efficiently and with good accuracy.

2. BOND GRAPH MODEL OF FLEXIBLE LINKS

The flexible multibody system can be represented as a collection of elementary components. Flexible components are considered as long slender beams whose models are built as aggregation of bond graph 3D beam finite element (FE) component models (Beam3D), Fig.1a. 2D bond graph component has been developed in (Damic,2006; Cohodar et al., 2009.). To develop basic 3D beam component model the co-rotation formulation is applied, (Battini&Pacoste, 2002; Damic&Cohodar,2006). Fig. 1b shows configuration of the basic beam 3D FE. It uses three coordinate frames:

1. Global (inertial) frame, represented by triads (I,J,K),

2. Co-rotation frame (i,j,k) defined on such way that moving of beam finite element with respect to the frame consists only deformation and

3. Section coordinate frame (l,m,n), rigidly attached to the element cross section.

[FIGURE 1 OMITTED]

Two element nodes of beam finite element are represented by two ports. Flows and efforts ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]) at ports are linear and angular velocities and corresponding forces and moment, respectively, defined in the global frame.

3. BOND GRAPH MODEL OF JOINT CONSTRAINS

The articulation between two consecutive links of robot manipulator can be realized by one of two basic joints: a prismatic or a revolute joint. Other type of joints, for instance cylindrical can be modelled as combination of prismatic and revolute joints. Bond graph model of joint can be represented by component Joint, Fig. 2a with two power ports, described by pair vectors--flow and effort. Let joint (i) connect link to link (i), Fig.2b. The flows at left and right ports are represented by vectors composed of linear and angular velocities of link [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and link [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], respectively (expressed in the global frame). If joint i is revolute the angular velocity of link i with respect to link (i-1) can be expressed in the joint coordinate frame ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]) as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

where [[theta].sub.i] is angle of rotation. In the case of the prismatic joint (with joint variable [d.sub.i]) relative linear velocity is given by:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (2)

[FIGURE 2 OMITTED]

Development of revolute joint model is based on the assumption that it is rigidly attached to the previously beam element. Thus, joint orientation can be defined through the orientation of the end right cross-section (denote with index 2 in Fig.1a) using rotation matrix [R.sub.2] = [[l.sub.2] [m.sub.2] [n.sub.2]].

The joint orientation with respect to the right end crosssection does not change and can be represented by a constant matrix [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The joint orientation with respect to inertial frame is defined by:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3) For the revolute joint, relationship between angular velocities at the output and input ports is defined by:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

For prismatic joint linear velocity at the output port can be expressed as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

4. NUMERICAL EXAMPLES

To demonstrate and verify applicability of proposed models two numerical examples are considered in this paper: slider-crank mechanism, and two-component robot arm. In both simulations the time step is 1e-3 s and error tolerance 1e-6.

4.1 Slider-crank mechanism

To validate the models a flexible slider-crank mechanism is adopted as the first numerical example. Bond graph model of the slider-crank mechanism (Fig.3a) consists of a slider block, crankshaft and connected rod that are connected by the revolute joint. Crankshaft performs only a rotation and is connected to the ground by a revolute joint. Motion of the other side of connected rod is restricted to the horizontal siding plane. The slider block is modelled as a massless particle. The lengths of the crankshaft and connected rod are 0.5 m and 1 m, respectively. The dimensions of cross-section area of both components are the same: H=50 mm and W=5 mm, as well as are the material properties: modulus elasticity is E = 2.1e10 Pa, density p = 7.8e3 kg/[m.sup.3]. The mechanism is initially in the horizontally and straighten position and performs the motion under the influence of its own weight. The mid-point deflection of connection rod is shown in Fig.3b and is in a good agreement with results reported in (Gerstmayr&Schoberl, 2006).

[FIGURE 3 OMITTED]

4.2 Two-component robot arm

The second numerical example is the two-component robot arm, taken from (Ibrahimbegovic&Mamouri, 2000) and (Kromer et al., 2004). Fig. 4 presents geometry and physical properties of the robot structure. Links are interconnected by a spherical joint. The left side of the first link is connected to the base by a cylindrical joint and subjected to the prescribed displacement and rotation as shown in Fig.4a. Both links are discretized by ten finite elements components.

[FIGURE 4 OMITTED]

Simulation interval is set to 3 s. Displacement components of free end link are depicted in Fig.5. The results agree well with that obtained in (Ibrahimbegovic&Mamouri, 2000) and (Kromer et al., 2004).

[FIGURE 5 OMITTED]

5. CONCLUSION

In this paper, systematic development of the model of a constrained flexible multibody system has been demonstrated where special attention has been focused on modelling of the joint constrains. Due to component model approach each developed component can be put in the library and taken out of it when developing new systems. Several numerical simulations are carried out to test accuracy of developed model and their results are compared with the results of other authors showing good agreements. The proposed method can be successfully applied to modelling and simulation of other mechatronics systems as well.

6. REFERENCES

Battini, M. & Pacoste, C. (2002). Co-rotational beam elements with warping effects in instability problems, Comput. Methods Appl. Mech. Engrg., 191: 1755-789

Cohodar, M. Borutzky, W. & Damic, V. (2009). Comparison of different formulation of 2D beam elements based on Bond Graph technique, Simulation Modelling Practice and Theory, 17:107-124

Damic V. & Cohodar, M. (2006). Bond Graph Based Modelling and Simulation of Flexible Robotic Manipulators, in: Proc. 20th European Conference on Modelling and Simulation ECMS, Bonn

Damic, V. & Montgomery, V. (2003). Mechatronic by Bond Graphs, Springer--Verlag, ISBN 3-540-42375-3, Berlin Heidelberg

Damic, V. (2006). Modelling flexible body systems: a bond graph component model approach. Mathematical and Computer Modelling of Dynamical Systems, No.12, (175-187)

Gerstmayr, J. & Schoberl, J. (2006). A 3D finite element methods for flexible multibody systems. Multibody Syst. Dyn, 15., (309-324), DOI 10.1007/s11044-006-9009-3

Ibrahimbegovic, A. & Mamouri, S. (2000). On rigid components and joint constraints in nonlinear dynamics of flexible multibody systems employing 3D geometrically exact beam model, Computer methods in applied mechanics and engineering, 188: 805-831

Karnopp, D.C., Margolis, D.L.&Rosenberg, R.C. (2000) System dynamic, Modelling and Simulation of Mechatronic System, John Wiley&Sons, Inc., ISBN:0-471-33301-8

Kromer, V. Dufosse, F. & Gueury, M. (2004) An object-oriented design of a finite element code: application to multibody system analysis. Advanced in Engineering Software, 35, 273-287

Printer friendly Cite/link Email Feedback | |

Author: | Damic, Vjekoslav; Cohodar, Maida; Kulenovic, Malik |
---|---|

Publication: | Annals of DAAAM & Proceedings |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2010 |

Words: | 1370 |

Previous Article: | The oil pressure distribution modifications due to the effect of unaligned shaft. |

Next Article: | Considerations about a 3D matrix based model for a porous scaffold and a cell suspension. |

Topics: |