Optimal H2 order-one reduction by solving eigenproblems for polynomial equations

A method is given for solving an optimal H2 approximation problem for SISO linear time-invariant stable systems. The method, based on constructive algebra, guarantees that the global optimum is found; it does not involve any gradient-based search, and hence avoids the usual problems of local minima. We examine mostly the case when the model order is reduced by one, and when the original system has distinct poles. This case exhibits special structure which allows us to provide a complete solution. The problem is converted into linear algebra by exhibiting a finite-dimensional basis for a certain space, and can then be solved by eigenvalue calculations, following the methods developed by Stetter and Moeller. The use of Buchberger's algorithm is avoided by writing the first-order optimality conditions in a special form, from which a Groebner basis is immediately available. Compared with our previous work the method presented here has much smaller time and memory requirements, and can therefore be applied to systems of significantly higher McMillan degree. In addition, some hypotheses which were required in the previous work have been removed. Some examples are included.

Bayesian two-sample tests

In this paper, we present two classes of Bayesian approaches to the two-sample problem. Our first class of methods extends the Bayesian t-test to include all parametric models in the exponential family and their conjugate priors. Our second class of methods uses Dirichlet process mixtures (DPM) of such conjugate-exponential distributions as flexible nonparametric priors over the unknown distributions.

Experimental and theoretical comparison between static and dynamic torsional soil tests

The dynamic properties of dry Leighton Buzzard sand have been investigated using a resonant column test apparatus. These data are compared with very low frequency cyclic tests on identical specimens of sand. The comparison indicates that the properties of dry sand are independent of frequency. A simple one-dimensional model of kinematic hardening plasticity is used to predict the dynamic behaviour of the sand. The input parameters for this model are based on the results of static tests. These may be conducted on standard laboratory equipment with only minor modifications. The predictions are in good agreement with the measured data.

Centrifuge tests investigating inclined grout micro-piles as a method of liquefaction remediation for existing buildings

This paper presents a series of centrifuge tests carried out to investigate the performance of non-structural inclined micro-piles as a potential liquefaction remediation method for existing buildings. Both a single-degree-of-freedom frame structure and a two-storey, two-degree-of-freedom frame structure were used as model buildings in these tests. Centrifuge tests were carried out with and without micro-piles in the foundation soil for each structure. Results primarily from the tests with the SDOF structure are presented in this paper. It is found that the micro-piles have some beneficial effect by increasing shear strains in the soil in their vicinity and hence causing dilation in these zones. However, they also increase structural accelerations by transmitting accelerations from deep in the soil and the beneficial effects from increased dilation are outweighed by the detrimental migration of pore pressures.

Alternative uses of partitioning interwell tracer tests for DNAPL source zone characterisation

The heterogeneous nature of the subsurface and associated DNAPL morphologies often poses the greatest limitation to source zone clean-up strategies. Hence, detailed site characterisation techniques are required. The data presented in this paper has been collected from a series of laboratory 2-D tank experiments and numerical simulations of Partitioning Interwell Tracer Tests (PITT) in a wide range of aquifer conditions and DNAPL morphologies. Alternative uses of tracer breakthrough data have been developed In order to characterise the mass flux generated from the DNAPL source. By combining the laboratory and numerical data, a relationship between normalised mass flux and tracer-based average source zone DNAPL saturation has been established. Knowledge of such a relationship allows remediation targets to be identified, clean-up efficiencies to be evaluated, and increases the accuracy of any risk assessment.

Risk-sensitive optimal feedback control accounts for sensorimotor behavior under uncertainty

Many aspects of human motor behavior can be understood using optimality principles such as optimal feedback control. However, these proposed optimal control models are risk-neutral; that is, they are indifferent to the variability of the movement cost. Here, we propose the use of a risk-sensitive optimal controller that incorporates movement cost variance either as an added cost (risk-averse controller) or as an added value (risk-seeking controller) to model human motor behavior in the face of uncertainty. We use a sensorimotor task to test the hypothesis that subjects are risk-sensitive. Subjects controlled a virtual ball undergoing Brownian motion towards a target. Subjects were required to minimize an explicit cost, in points, that was a combination of the final positional error of the ball and the integrated control cost. By testing subjects on different levels of Brownian motion noise and relative weighting of the position and control cost, we could distinguish between risk-sensitive and risk-neutral control. We show that subjects change their movement strategy pessimistically in the face of increased uncertainty in accord with the predictions of a risk-averse optimal controller. Our results suggest that risk-sensitivity is a fundamental attribute that needs to be incorporated into optimal feedback control models. © 2010 Nagengast et al.

Estimating material elasticity by spherical indentation load-relaxation tests on viscoelastic samples of finite thickness

A two-step viscoelastic spherical indentation method is proposed to compensate for 1) material relaxation and 2) sample thickness. In the first step, the indenter is moved at a constant speed and the reaction force is measured. In the second step, the indenter is held at a constant position and the relaxation response of the material is measured. Then the relaxation response is fit with a multi-exponential function which corresponds to a three-branch general Maxwell model. The relaxation modulus is derived by correcting the finite ramp time introduced in the first step. The proposed model takes into account the sample thickness, which is important for applications in which the sample thickness is less than ten times the indenter radius. The model is validated numerically by finite element simulations. Experiments are carried out on a 10% gelatin phantom and a chicken breast sample with the proposed method. The results for both the gelatin phantom and the chicken breast sample agree with the results obtained from a surface wave method. Both the finite element simulations and experimental results show improved elasticity estimations by incorporating the sample thickness into the model. The measured shear elasticities of the 10% gelatin sample are 6.79 and 6.93 kPa by the proposed finite indentation method at sample thickness of 40 and 20 mm, respectively. The elasticity of the same sample is estimated to be 6.53 kPa by the surface wave method. For the chicken breast sample, the shear elasticity is measured to be 4.51 and 5.17 kPa by the proposed indentation method at sample thickness of 40 and 20 mm, respectively. Its elasticity is measured by the surface wave method to be 4.14 kPa. © 2011 IEEE.