Automatic generation of equations of motion for multibody system in discrete event simulation framework

AbstractIn this paper, the development of a simulation program that can automatically generate equations of motion for mutibody systems in the discrete event simulation framework is presented. The need to analyze the dynamic response of mechanical systems that are under event triggered conditions is increasing. General mechanical systems can be defined as multibody systems that are collections of interconnected rigid bodies, consistent with various types of joints that limit the relative motion of pairs of bodies. For complex multibody systems, a systematic approach is required to efficiently set up the mathematical models. Therefore, a dynamics kernel was developed to automatically generate the equations of motion for multibody systems based on multibody dynamics. The developed dynamics kernel also provides the numerical solver for the dynamic analysis of multibody systems. The general multibody dynamics kernel cannot deal with discontinuous state variables, event triggered conditions, and state triggered conditions, though. To enable it to deal with multibody systems in discontinuous environments, the multibody dynamics kernel was integrated into a discrete event simulation framework, which was developed based on the discrete event system specification (DEVS) formalism. DEVS formalism is a modular and hierarchical formalism for modeling and analyzing systems under event triggered conditions, which are described by discontinuous state variables. To verify the developed program, it was applied to an block-lifting and transport simulation, and dynamic analysis of the system is carried out

An intrinsic Hamiltonian formulation of the dynamics of LC-circuits

First, the dynamics of LC-circuits are formulated as a Hamiltonian system defined with respect to a Poisson bracket which may be degenerate, i.e., nonsymplectic. This Poisson bracket is deduced from the network graph of the circuit and captures the dynamic invariants due to Kirchhoff's laws. Second, the antisymmetric relations defining the Poisson bracket are realized as a physical network using the gyrator element and partially dualizing the network graph constraints. From the network realization of the Poisson bracket, the reduced standard Hamiltonian system as well as the realization of the embedding standard Hamiltonian system are deduce

MIT Space Engineering Research Center

The Space Engineering Research Center (SERC) at MIT, started in Jul. 1988, has completed two years of research. The Center is approaching the operational phase of its first testbed, is midway through the construction of a second testbed, and is in the design phase of a third. We presently have seven participating faculty, four participating staff members, ten graduate students, and numerous undergraduates. This report reviews the testbed programs, individual graduate research, other SERC activities not funded by the Center, interaction with non-MIT organizations, and SERC milestones. Published papers made possible by SERC funding are included at the end of the report

Graph-theoretic Sensitivity Analysis of Dynamic Systems

The main focus of this research is to use graph-theoretic formulations to develop an automated algorithm for the generation of sensitivity equations. The idea is to combine the benefits of direct differentiation with that of graph-theoretic formulation. The primary deliverable of this work is the developed software module which can derive the system equations and the sensitivity equations directly from the linear graph of the system. Sensitivity analysis refers to the study of changes in system behaviour brought about by the changes in model parameters. Due to the rapid increase in the sizes and complexities of the models being analyzed, it is important to extend the capabilities of the current tools of sensitivity analysis, and an automated, efficient, and accurate method for the generation of sensitivity equations is highly desirable. In this work, a graph-theoretic algorithm is developed to generate the sensitivity equations. In the current implementation, the proposed algorithm uses direct differentiation to generate sensitivity equations at the component level and graph-theoretic methods to assemble the equation fragments to form the sensitivity equations. This way certain amount of control can be established over the size and complexity of the generated sensitivity equations. The implementation of the algorithm is based on a commercial software package \verb MapleSim[Multibody] and can generate governing and sensitivity equations for multibody models created in MapleSim. In this thesis, the algorithm is tested on various mechanical, hydraulic, electro-chemical, multibody, and multi-domain systems. The generated sensitivity information are used to perform design optimization and parametric importance studies. The sensitivity results are validated using finite difference formulations. The results demonstrate that graph-theoretic sensitivity analysis is an automated, accurate, algorithmic method of generation for sensitivity equations, which enables the user to have some control over the form and complexity of the generated equations. The results show that the graph-theoretic method is more efficient than the finite difference approach. It is also demonstrated that the efficiency of the generated equations are at par or better than the equation obtained by direct differentiation.1 yea

Micromechanics of Fiber Networks Including Nonlinear Hysteresis and its Application to Multibody Dynamic Modeling of Piano Mechanisms

Many engineering applications make use of fiber assemblies under compression. Unfortunately, this compression behavior is difficult to predict, due to nonlinear compliance, hysteresis, and anelasticity. The main objective of this research is to develop an algorithm which is capable of incorporating the microscale features of the fiber network into macroscopic scale applications, particularly the modeling of contact mechanics in multibody systems. In micromechanical approaches, the response of a fiber assembly to an external force is related to the response of basic fiber units as well as the interactions between these units, i.e. the mechanical properties of the constituent fibers and the architecture of the assembly will both have a significant influence on the overall response of the assembly to compressive load schemes. Probabilistic and statistical principles are used to construct the structure of the uniformly-distributed random network. Different micromechanical approaches in modeling felt, as a nonwoven fiber assembly with unique mechanical properties, are explored to gain insight into the key mechanisms that influence its compressive response. Based on the deformation processes and techniques in estimating the number of fiber contacts, three micromechanical models are introduced: (1) constitutive equations for micromechanics of three-dimensional fiberwebs under small strains, in which elongation of the fibers is the key deformation mechanism, adapted for large deformation ranges; (2) micromechanical model based on the rate theory of granular media, in which bending and torsion of fibers are the predominant elemental deformations used to calculate compliances of a particular contact; and (3) a mechanistic model developed using the general deformation theory of the fiber networks with fiber bending at the micro level and a binomial distribution of fiber contacts. A well-established mechanistic model, based on fiber-to-fiber friction at the micro level, is presented for predicting the hysteresis in compression behavior of wool fiberwebs. A novel algorithm is introduced to incorporate a hysteretic micromechanical model - a combination of the mechanistic model with microstructural fiber bending, which uses a binomial distribution of the number of fiber-to-fiber contacts, and the friction-based hysteresis idea - into the contact mechanics of multibody simulations with felt-lined interacting bodies. Considering the realistic case in which a portion of fibers slides, the fiber network can be treated as two subnetworks: one from the fibers with non-sliding contact points, responsible for the elastic response of the network, and the other consisting of fibers that slide, generating irreversible hysteresis phenomenon in the fiberweb compression. A parameter identification is performed to minimize the error between the micromechanical model and the elastic part of the loading-unloading experimental data for felt, then contribution of friction was added to the obtained mechanistic compression-recovery curves. The theoretical framework for constructing a mechanistic multibody dynamic model of a vertical piano action is developed, and its general validity is established using a prototype model. Dynamic equations of motion are derived symbolically for the piano action using a graph-theoretic formulation. The model fidelity is increased by including hammer-string interaction, backcheck wire and hammer shank flexibility, a sophisticated key pivot model, nonlinear models of bridle strap and butt spring, and a novel mathematical contact model. The developed nonlinear hysteretic micromechanical model is used for the hammer-string interaction to affirm the reliability and applicability of the model in general multibody dynamic simulations. In addition, dynamic modeling of a flexible hub-beam system with an eccentric tip mass including nonlinear hysteretic contact is studied. The model represents the mechanical finger of an actuator for a piano key. Achieving a desired finger-key contact force profile that replicates that of a real pianist's finger requires dynamic and vibration analysis of the actuator device. The governing differential equations for the dynamic behavior of the system are derived using Euler-Bernoulli beam theory along with Lagrange's method. To discretize the distributed parameter flexible beam in the model, the finite element method is utilized. Excessive vibration due to the arm flexibility and also the rigid-body oscillations of the arm, especially during the period of key-felt contact, is eliminated utilizing a simple grounded rotational dashpot and a grounded rotational dashpot with a one-sided relation. The effect on vibration behavior attributed to these additional components is demonstrated using the simulated model

Coarse-grained simulations to predict structure and properties of polymer nanocomposites

textPolymer Nanocomposites (PNC) are a new class of materials characterized by their large interfacial areas between the host polymer and nanofiller. This unique feature, due to the size of the nanofiller, is understood to be the cause of enhanced mechanical, electrical, optical, and barrier properties observed of PNCs, relative to the properties of the unfilled polymer. This interface can determine the miscibility of the nanofiller in the polymer, which, in turn, influences the PNC's properties. In addition, this interface alters the polymer's structure near the surface of the nanofiller resulting in heterogeneity of local properties that can be expressed at the macroscopic level. Considering the polymer-nanoparticle interface significantly influences PNC properties, it is apparent that some atomistic level of detail is required to accurately predict the behavior of PNCs. Though an all-atom simulation of a PNC would be able to accomplish the latter, it is an impractical approach to pursue even with the most advanced computational resources currently available. In this contribution, we develop (1) an equilibrium coarse-graining method to predict nanoparticle dispersion in a polymer melt, (2) a dynamic coarse-graining method to predict rheological properties of polymer-nanoparticle melt mixtures, and (3) a numerical approach that includes interfacial layer effects and polymer rigidity when predicting barrier properties of PNCs. In addition to the above, we study how particle and polymer characteristics affect the interfacial layer thickness as well as how the polymer-nanoparticle interface may influence the entanglement network in a polymer melt. More specifically, we use a mean-field theory approach to discern how the concentration of a semiflexible polymer, its rigidity and the particle's size determine the interfacial layer thickness, and the scaling laws to describe this dependency. We also utilize molecular dynamics and simulation techniques on a model PNC to determine if the polymer-nanoparticle interaction can influence the entanglement network of a polymer melt.Chemical Engineerin