The use of object oriented modelling and simulation tools in constraint satisfaction for embodiment design

The software package Dymola, which implements the new, vendor-independent standard modelling language Modelica, exemplifies the emerging generation of object-oriented modelling and simulation tools. This paper shows how, in addition to its simulation capabilities, it may be used as an embodiment design tool, to size automatically a design assembled from a library of generic parametric components. The example used is a miniature model aircraft diesel engine. To this end, the component classes contain extra algebraic equations calculating the overload factor (or its reciprocal, the safety factor) for all the different modes of failure, such as buckling or tensile yield. Thus the simulation results contain the maximum overload or minimum safety factor for each failure mode along with the critical instant and the device state at which it occurs. The Dymola "Initial Conditions Calculation" function, controlled by a simple software script, may then be used to perform automatic component sizing. Each component is minimised in mass, subject to a chosen safety factor against failure, over a given operating cycle. Whilst the example is in the realm of mechanical design, it must be emphasised that the approach is equally applicable to the electrical or mechatronic domains, indeed to any design problem requiring numerical constraint satisfaction.

Mechanical constraint and loading on ferroelectric memory capacitors

The influence of mechanical constraint imposed by device geometry upon the switching response of a ferroelectric thin film memory capacitor is investigated. The memory capacitor was represented by two-dimensional ferroelectric islands of different aspect ratio, mechanically constrained by surrounding materials. Its ferroelectric non-linear behaviour was modeled by a crystal plasticity constitutive law and calculated using the finite element method. The switching response of the device, in terms of remnant charge storage, was determined as a function of geometry and constraint. The switching response under applied in-plane tensile stress and hydrostatic pressure was also studied experimentally. Our results showed that (1) the capacitor's aspect ratio could significantly affect the clamping behaviour and thus the remnant polarization, (2) it was possible to maximise the switching charge through the optimisation of the device geometry, and (3) it is possible to find a critical switching stress at zero electric field and a critical coercive field at zero residual stress. © 2009 Materials Research Society.

Effects of constraint curvature on structural instability: tensile buckling and multiple bifurcations

Bifurcation of an elastic structure crucially depends on the curvature of the constraints against which the ends of the structure are prescribed to move, an effect which deserves more attention than it has received so far. In fact, we show theoretically and we provide definitive experimental verification that an appropriate curvature of the constraint over which the end of a structure has to slide strongly affects buckling loads and can induce: (i.) tensile buckling; (ii.) decreasing- (softening), increasing- (hardening), or constant-load (null stiffness) postcritical behaviour; (iii.) multiple bifurcations, determining for instance two bifurcation loads (one tensile and one compressive) in a single-degree-of-freedom elastic system. We show how to design a constraint profile to obtain a desired postcritical behaviour and we provide the solution for the elastica constrained to slide along a circle on one end, representing the first example of an inflexional elastica developed from a buckling in tension. These results have important practical implications in the design of compliant mechanisms and may find applications in devices operating in quasi-static or dynamic conditions.

Optimal constraint tightening policies for robust variable horizon model predictive control

This paper develops a technique for improving the region of attraction of a robust variable horizon model predictive controller. It considers a constrained discrete-time linear system acted upon by a bounded, but unknown time-varying state disturbance. Using constraint tightening for robustness, it is shown how the tightening policy, parameterised as direct feedback on the disturbance, can be optimised to increase the volume of an inner approximation to the controller's true region of attraction. Numerical examples demonstrate the benefits of the policy in increasing region of attraction volume and decreasing the maximum prediction horizon length. © 2012 IEEE.

Size effects in the constrained deformation of metallic foams

The constrained deformation of an aluminium alloy foam sandwiched between steel substrates has been investigated. The sandwich plates are subjected to through-thickness shear and normal loading, and it is found that the face sheets constrain the foam against plastic deformation and result in a size effect: the yield strength increases with diminishing thickness of foam layer. The strain distribution across the foam core has been measured by a visual strain mapping technique, and a boundary layer of reduced straining was observed adjacent to the face sheets. The deformation response of the aluminium foam layer was modelled by the elastic-plastic finite element analysis of regular and irregular two dimensional honeycombs, bonded to rigid face sheets; in the simulations, the rotation of the boundary nodes of the cell-wall beam elements was set to zero to simulate full constraint from the rigid face sheets. It is found that the regular honeycomb under-estimates the size effect whereas the irregular honeycomb provides a faithful representation of both the observed size effect and the observed strain profile through the foam layer. Additionally, a compressible version of the Fleck-Hutchinson strain gradient theory was used to predict the size effect; by identifying the cell edge length as the relevant microstructural length scale the strain gradient model is able to reproduce the observed strain profiles across the layer and the thickness dependence of strength. © 2002 Elsevier Science Ltd. All rights reserved.

The plastic collapse and energy absorption capacity of egg-box panels

The plastic collapse response of aluminium egg-box panels subjected to out-of-plane compression has been measured and modelled. It is observed that the collapse strength and energy absorption are sensitive to the level of in-plane constraint, with collapse dictated either by plastic buckling or by a travelling plastic knuckle mechanism. Drop weight tests have been performed at speeds of up to 6 m s-1, and an elevation in strength with impact velocity is noted. A 3D finite element shell model is needed in order to reproduce the observed behaviours. Additional calculations using an axisymmetric finite element model give the correct collapse modes but are less accurate than the more sophisticated 3D model. The finite element simulations suggest that the observed velocity dependence of strength is primarily due to strain-rate sensitivity of the aluminium sheet, with material inertia playing a negligible role. Finally, it is shown that the energy absorption capacity of the egg-box material is comparable to that of metallic foams. © 2003 Elsevier Ltd. All rights reserved.

Robust multiplexed model predictive control

This paper extends the recently developed multiplexed model predictive control (MMPC) concept to ensure satisfaction of hard constraints despite the action of persistent, unknown but bounded disturbances. MMPC uses asynchronous control moves on each input channel instead of synchronised moves on all channels. It offers reduced computation, by dividing the online optimisation into a smaller problem for each channel, and potential performance improvements, as the response to a disturbance is quicker, albeit via only one channel. Robustness to disturbances is introduced using the constraint tightening approach, tailored to suit the asynchronous updates of MMPC and the resulting time-varying optimisations. Numerical results are presented, involving a simple mechanical example and an aircraft control example, showing the potential computational and performance benefits of the new robust MMPC.

Using strong sustainability to optimize electricity generation fuel mixes

This work represents a contribution to the field of sustainable electricity system design by using an optimization tool to specify the final mix composition, subject to the constraints of: emissions that are within the biocapacity of the region; a diverse and robust electricity supply system; and supply that at least meets current demand. The 25-country European Union (EU-25) is used as a case study. All the goals, save diversity, can be met by re-structuring the current fuel mix, thus maintaining current consumption levels. The diversity target is only met when consumption is reduced by 10-15% and the constraint on maximum material throughput is relaxed. Re-structuring the mix and reducing consumption is insufficient to achieve a sustainable EU carbon footprint. However, the solution proposed singlehandedly allows the EU to meet its Kyoto emissions target as well as its 2007 policy of a reduction of 20% in greenhouse gas emissions by 2020. © 2007 Elsevier Ltd. All rights reserved.

Designing MPC controllers by reverse-engineering existing LTI controllers

This technical report presents a method for designing a constrained output-feedback model predictive controller (MPC) that behaves in the same way as an existing baseline stabilising linear time invariant output-feedback controller when constraints are inactive. The baseline controller is cast into an observer-compensator form and an inverse-optimal cost function is used as the basis of the MPC controller. The available degrees of design freedom are explored, and some guidelines provided for the selection of an appropriate observer-compensator realisation that will best allow exploitation of the constraint-handling and redundancy management capabilities of MPC. Consideration is given to output setpoint tracking, and the method is demonstrated with three different multivariable plants of varying complexity.

Engineering design and mechatronics: The schemebuilder project

The paper describes the development of a software design aid for use at the conceptual stage of engineering design. It is intended for use in the design of mechatronic products but has wider potential uses. Early approaches were based on function structures and tables of options and the system that evolved allows the assembly of schemes linked by matching their input and output ports. A database of components is provided which can be accessed via different indexes and the designer can easily create and compare alternative schemes at the concept stage. A bond graph approach is used to define the interconnections between components. This allows correct port matching but also provides for future development such as constraint propagation through the design and links to simulation tools. © 1993 Springer-Verlag New York Inc.

The use of syntax and multiple alternatives in the VODIS voice operated database inquiry system

The paper describes the architecture of VODIS, a voice operated database inquiry system, and presents some experiments which investigate the effects on performance of varying the level of a priori syntactic constraints. The VODIS system includes a novel mechanism for incorporating context-free grammatical constraints directly into the word recognition algorithm. This allows the degree of a priori constraint to be smoothly varied and provides for the controlled generation of multiple alternatives. The results show that when the spoken input deviates from the predefined task grammar, a combination of weak a priori syntax rules in conjunction with full a posteriori parsing on a lattice of alternative word matches provides the most robust recognition performance. © 1991.

A comparison of different failure assessment methodologies applied to the NESC-1 test

Predictions for a 75x205mm surface semi-elliptic defect in the NESC-1 spinning cylinder test have been made using BS PD 6493:1991, the R6 procedure, non-linear cracked body finite element analysis techniques and the local approach to fracture. All the techniques agree in predicting ductile tearing near the inner surface of the cylinder followed by cleavage initiation. However they differ in the amount of ductile tearing, and the exact location and time of any cleavage event. The amount of ductile tearing decreases with increasing sophistication in the analysis, due to the drop in peak crack driving force and more explicit consideration of constraint effects. The local approach predicts a high probability of cleavage in both HAZ and base material after 190s, while the other predictions suggest that cleavage is unlikely in the HAZ due to constraint loss, but likely in the underlying base material. The timing of this event varies from ∼150s for R6 predictions to ∼250-300s using non-linear cracked body analysis.