Automatic Differentiation of Rigid Body Dynamics for Optimal Control and Estimation

Many algorithms for control, optimization and estimation in robotics depend on derivatives of the underlying system dynamics, e.g. to compute linearizations, sensitivities or gradient directions. However, we show that when dealing with Rigid Body Dynamics, these derivatives are difficult to derive analytically and to implement efficiently. To overcome this issue, we extend the modelling tool `RobCoGen' to be compatible with Automatic Differentiation. Additionally, we propose how to automatically obtain the derivatives and generate highly efficient source code. We highlight the flexibility and performance of the approach in two application examples. First, we show a Trajectory Optimization example for the quadrupedal robot HyQ, which employs auto-differentiation on the dynamics including a contact model. Second, we present a hardware experiment in which a 6 DoF robotic arm avoids a randomly moving obstacle in a go-to task by fast, dynamic replanning

Efficient multi-scale modelling of path dependent problems – complas 2017

With growing capabilities of computers use of multi-scale methods for detailed analysis of response with respect to material and geometric nonlinearities is becoming more relevant. In this paper focus is on MIEL (mesh-in-element) multi-scale method and its implementation with AceGen and AceFEM based on analytical sensitivity analysis. Such implementation enables efficient multi-scale modelling, consistency and quadratic convergence also for two-level path following methods for the solution of path dependent problems

Global Stability Analysis of Healthy Situation for a Coupled Model of Healthy and Cancerous Cells Dynamics in Acute Myeloid Leukemia

In this paper we aim to study the global stability of a coupled model of healthy and cancerous cells dynamics in healthy situation of Acute Myeloid Leukemia. We also clarify the effect of interconnection between healthy and cancerous cells dynamics on the global stability The interconnected model is obtained by transforming the PDE-based model into a nonlinear distributed delay system. Using Lyapunov approach, we derive necessary and sufficient conditions for global stability for a selected equilibrium point of particular interest (healthy situation). Simulations are conducted to illustrate the obtained results

Feedback methods for inverse simulation of dynamic models for engineering systems applications

Inverse simulation is a form of inverse modelling in which computer simulation methods are used to find the time histories of input variables that, for a given model, match a set of required output responses. Conventional inverse simulation methods for dynamic models are computationally intensive and can present difficulties for high-speed applications. This paper includes a review of established methods of inverse simulation,giving some emphasis to iterative techniques that were first developed for aeronautical applications. It goes on to discuss the application of a different approach which is based on feedback principles. This feedback method is suitable for a wide range of linear and nonlinear dynamic models and involves two distinct stages. The first stage involves design of a feedback loop around the given simulation model and, in the second stage, that closed-loop system is used for inversion of the model. Issues of robustness within closed-loop systems used in inverse simulation are not significant as there are no plant uncertainties or external disturbances. Thus the process is simpler than that required for the development of a control system of equivalent complexity. Engineering applications of this feedback approach to inverse simulation are described through case studies that put particular emphasis on nonlinear and multi-input multi-output models

The genotype-phenotype relationship in multicellular pattern-generating models - the neglected role of pattern descriptors

Background: A deep understanding of what causes the phenotypic variation arising from biological patterning processes, cannot be claimed before we are able to recreate this variation by mathematical models capable of generating genotype-phenotype maps in a causally cohesive way. However, the concept of pattern in a multicellular context implies that what matters is not the state of every single cell, but certain emergent qualities of the total cell aggregate. Thus, in order to set up a genotype-phenotype map in such a spatiotemporal pattern setting one is actually forced to establish new pattern descriptors and derive their relations to parameters of the original model. A pattern descriptor is a variable that describes and quantifies a certain qualitative feature of the pattern, for example the degree to which certain macroscopic structures are present. There is today no general procedure for how to relate a set of patterns and their characteristic features to the functional relationships, parameter values and initial values of an original pattern-generating model. Here we present a new, generic approach for explorative analysis of complex patterning models which focuses on the essential pattern features and their relations to the model parameters. The approach is illustrated on an existing model for Delta-Notch lateral inhibition over a two-dimensional lattice. Results: By combining computer simulations according to a succession of statistical experimental designs, computer graphics, automatic image analysis, human sensory descriptive analysis and multivariate data modelling, we derive a pattern descriptor model of those macroscopic, emergent aspects of the patterns that we consider of interest. The pattern descriptor model relates the values of the new, dedicated pattern descriptors to the parameter values of the original model, for example by predicting the parameter values leading to particular patterns, and provides insights that would have been hard to obtain by traditional methods. Conclusion: The results suggest that our approach may qualify as a general procedure for how to discover and relate relevant features and characteristics of emergent patterns to the functional relationships, parameter values and initial values of an underlying pattern-generating mathematical model

Airborne photogrammetry and LIDAR for DSM extraction and 3D change detection over an urban area : a comparative study

A digital surface model (DSM) extracted from stereoscopic aerial images, acquired in March 2000, is compared with a DSM derived from airborne light detection and ranging (lidar) data collected in July 2009. Three densely built-up study areas in the city centre of Ghent, Belgium, are selected, each covering approximately 0.4 km(2). The surface models, generated from the two different 3D acquisition methods, are compared qualitatively and quantitatively as to what extent they are suitable in modelling an urban environment, in particular for the 3D reconstruction of buildings. Then the data sets, which are acquired at two different epochs t(1) and t(2), are investigated as to what extent 3D (building) changes can be detected and modelled over the time interval. A difference model, generated by pixel-wise subtracting of both DSMs, indicates changes in elevation. Filters are proposed to differentiate 'real' building changes from false alarms provoked by model noise, outliers, vegetation, etc. A final 3D building change model maps all destructed and newly constructed buildings within the time interval t(2) - t(1). Based on the change model, the surface and volume of the building changes can be quantified