

Concurrent Dynamics delivers high fidelity Multibody Dynamics Simulation programs to assess dynamics and control performances of mechanisms and mobile vehicles to expedite:
Satellite 'Ibuki': courtesy of JAXA MBS simulations and dynamics formulation notes follow. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 1. MBS Simulations We have readytouse simulations for a variety of vehicles and operations. Beside the ones below, more are in the robotics page and in the works. The O(n) equations of motion in these simulations are derived from the NewtonEuler equations of a typical joint connected body. They run on PC's that have Matlab/Simulink or C/C++ Compilers installed. Their main program, XSIM1,
Contact us regarding any of the simulations below by writing to support@concurrentdynamics.com. Click 'pdf' for details on creating the model file using the editor, buildx.exe, and on the mdl file or the control.dll file designed for that application and more. We welcome suggestions on vehicle simulations that you would like to have. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ XSIM1 Simulations on Simulink Platform. The Xsim1seriesdll's are our processors that make it possible to simulate the dynamics and control of a wide variety of multibody systems in Simulink. It's the green block in each of the five examples below. For each application, user would build a model file to define the model parameters for the equations of motion processing and operational input and output parameters. With that, the user's task is to define the mapping from the engine output(ydata) to its input(udata) to close the control loop. Click on the pdf or the Simulink diagram to see the details of each of the simulations presented. SV1sim: Simulates a vehicle with 4 wheels and 8 jets under a nominal rwa ACS and a jet ACS. The control task is to maintain vehicle attitude at LVLH given some initial attitude disturbance. pdf Arraysim: Simulates an expanded version of SV1sim vehicle. It retains the same ACS, but now has two arrays each with a yoke and 3 panels. The control tasks include a coordinated array deployment using the 'gear' and 'lock' constraints and the maintenance of LVLH attitude with the ACS. pdf Sep1sim: Simulates the ejection of a satellite from a spinning mothership. The control tasks include ejecting the satellite, deploying the arrays after separation, followed by nulling the satellite rate and capturing LVLH attiude with jet ACS. The remainder of the simulation is to maintain LVLH attitude using RWA ACS. Examples show the VMASS (variable mass) feature of the simulation engine in managing a vehicle mass property discontinuity at the instant of the separation. pdf Tethersim: Simulates the deployment of a payload with a tether. Tether length and pushoff force are adjustable. Simulation of a case where the vehicletopayload mass ratio of ten is presented. The vehicle flies in an inclined near circular LEO orbit. Simulations with different mass ratio can be designed. The tether stiffness and damping must be adjusted according to the mass ratio to obtain a successful tether deployment, i.e. small steady state libration angle. pdf CMG4sim: Simulates a vehicle with 4 SGCMG's and 6 jets. It has a CMG ACS and it fires jet thrusters to reduce excess vehicle momentum. The vehicle flies in a near circular inclined LEO orbit with an aspherical gravity model. Users can redesign this simulation with different CMG mounting configurations and test the corresponding CMG ACS. LVLH control is presented, sunnadir ACS can be implemented also. pdf ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Simulatons with Cprogrammable XSIM1 Examples in these packages run on PC's that have C/C++ compilers installed. Each package has the program XSIM1.exe to run simulations based on a user specified model file and a user supplied control.dll file for a user conceived application. Our editor, buildx.exe and application library make the programming of that control.dll file easy by creating a near functional template to accomodate user's control system. The user's manual in each package shows how the model file is built and how the associated control.dll is programmed and compiled for each example.
. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 2. Multibody Dynamics Notes Kinematics Kinematics of a mechanism is about computing the motion of bodies in the mechanism given its generalized coordinates and rates. The motion includes the attitude, position and velocity of the bodies in some reference frame over time. As an example, the generalized coordinates of a chained system with 1 dof rotational joints are the joint angles of that system . More complex mechanisms can have 3 dof rotational or translational joints between bodies. In these cases, the generalized coordinates are the relative quaternions and the relative displacement coordinates of the joints in the local body frame accordingly. The kinematics of a generic treestructured system with 1 dof rotational joints is presented in the attached file. pdf Order of MBD Solutions Speed of solving the generalized accelerations in multibody dynamics (MBD) problems is a function of the number of arithmetic operations involved,i.e. additions and multiplications. It can be expressed by a polynomial in N or n. N is the number of bodies in the system and n is the degrees of freedom of the system. This polynomial is the exact computational cost of the MBD solution in relation to N. For example, if the polynomial is kN, for some integer k, then the solution is said to be order(N). If that polynomial is aN+bN^2+cN^3 for some integers (a,b,c),then the solution is said to be order(N^3). The same can be said of the order if the polynomial is in n. Between MBD solutions of the same problem, the lower order solutions execute faster than the higher order ones and progressively more so as N increases. This topic is expanded in the attached file. pdf Generalized Mass Matrix Generalized mass matrix is an important component in the study of the system kinetic energy, momentum and force equations of an MBS. This matrix is a second order partial derivative of the system kinetic energy with respect to the generalized rates with the energy given in terms of total angular rates and velocities of the bodies in the system. Obtaining this partial derivative can be challenging even for a small scale mechanism. A simpler way to compute this matrix for any treestructured MBS using recursive operators is shown in the attached file. pdf NewtonEuler Equations EulerLagrange equations, D'Alembert's work principle, NewtonEuler equations and few other principles/methods are often used in formulating MBD problems. Regardless of the physics principle used, the final eom for the targeted mechanism for the same generalized coordinates are identical. One variation in the eom is that they could be formulated to solve for the system accelerations or for the time derivative of system momenta. Still, solving either form of eom and propagating the motion with the associated kinematics equations should yield the same dynamics response for the same MBD problem. Furthermore, the eom must comply with momentum conservation when no external force or torque is applied. pdf Hamilton's Equations The Hamiltonian is the total energy of the system and it permits a transformation of the Lagrangian equations into a system of 2n first order differential equations in generalized coordinates and generalized momenta, n being the degrees of freedom in the system. These equations produce identical results as those by acceleration based eom. pdf Joint Motion Joints connect bodies in a mechanism and define the relative motion between them. It is the coordinated movements at these joints that creates the trajectory of the robot arm grapple device, the cyclic up and down motion of a piston in a combustion engine, and the precise placement or pointing of a device on a hexapod platform. Most mechanisms are designed with 1 dof rotational joints. Yet, it is not uncommon to see universal joints, ball joints or prismatic joints in mechanisms. Kinematics equations for four types of joints are presented in the attached file. pdf Redundant vs Minimal Coordinates Formulations MBD formulations using redundant coordinates can be done by assigning 6 dof's to each joint in the system. It can also be done by treating each body in the system as a free body thereby giving it 6 dof in the inertial frame. Either way, constraints are defined to keep the values of the redundant coordinates at each joint to zero to form the MBD eom.This approach creates a DAE model that has the right number of dof's fitting the simulated mechanism. In both cases,the systemwide accelerations can be solved by two simple but different O(N) procedures if the joint activation forces and jointconstraint forces are given. However,the required jointconstraint forces for this solution method for treestructured MBS's can only be solved in an O(N^3) manner on sequential processors, i.e. PC's.This effectively makes the said DAE solution an O(N^3) process. The NE and Hamilton's methods of minimal coordinates formulation discussed earlier for treestructured MBS do not have Lagrange multipliers in the eom as required by the redundant coordinates approach. In these cases, the derived eom size equals the degrees of freedom in the system. Furthermore, the generalized accelerations can be solved from these eom's in an O(N) manner given the jointactivation forces. Notably, the joint constraint forces were eliminated during the formulation process. The execution speed and eomsize advantage of the minimal coordinates eom formulations speak to its strength over the formulation with redundant coordinates methods. It must be mentioned that the use of Lagrange multipliers is very important in many MBD problems with closedloop, nonjoint and contact constraints. We will explore this topic later. System Angular Momentum Testing Multibody system angular momentum is a state of motion given the vehicle's generalized coordinates and rates. It is constant in the inertial space when no external torque is exerted on the system. Let P1 test be a comparison of this state computed by the simulation program against the one computed independently based on the body angular rates, velocities and stored angular momenta that are output from the program. The two must be identical. A discrepancy may mean that the implemented mass property is different from the intended one or that the simulator did not compute the system mass matrix properly. Let P2 test be a momentum conservation test. Notably, this test applies to vehicles whose base body has 6 dof wrt the inertial space and it does not apply to ground fixed systems. The reason for the latter is that the constraint force at the interface between the base body and the ground is an external force to the system and it is never zero over time. Passing the P1 and P2 tests means that the dynamics principle underlying the simulation and its implementation are credible and that the associated result is useful. The postsim system angular momentum is as easy to calculate as it is powerful. That calculation by Eq.(9) in the attached file requires no knowlege of the jointtopology of the mechanism or the dynamics formulation method used by the simulation other than that the base body of the system has 6dof wrt the inertial space. pdf References Click here to see selected papers and books published on multibody dynamics formulations. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Contact support@concurrentdynamics.com Interesting Websites:

Concurrent Dynamics International
