Lower bound limit analysis of cohesive‐frictional materials using second‐order cone programming

The formulation of limit analysis by means of the finite element method leads to an optimization problem with a large number of variables and constraints. Here we present a method for obtaining strict lower bound solutions using second‐order cone programming (SOCP), for which efficient primal‐dual interior‐point algorithms have recently been developed. Following a review of previous work, we provide a brief introduction to SOCP and describe how lower bound limit analysis can be formulated in this way. Some methods for exploiting the data structure of the problem are also described, including an efficient strategy for detecting and removing linearly dependent constraints at the assembly stage. The benefits of employing SOCP are then illustrated with numerical examples. Through the use of an effective algorithm/software, very large optimization problems with up to 700 000 variables are solved in minutes on a desktop machine. The numerical examples concern plane strain conditions and the Mohr–Coulomb criterion, however we show that SOCP can also be applied to any other problem of lower bound limit analysis involving a yield function with a conic quadratic form (notable examples being the Drucker–Prager criterion in 2D or 3D, and Nielsen's criterion for plates). Copyright © 2005 John Wiley & Sons, Ltd.     

Coupled mechanics and finite element for non‐linear laminated piezoelectric shallow shells undergoing large displacements and rotations

A theoretical framework is presented for analysing the coupled non‐linear response of shallow doubly curved adaptive laminated piezoelectric shells undergoing large displacements and rotations. The formulated mechanics incorporate coupling between in‐plane and flexural stiffness terms due to geometric curvature, coupling between mechanical and electric fields, and encompass geometric non‐linearity effects due to large displacements and rotations. The governing equations are formulated explicitly in orthogonal curvilinear co‐ordinates and are combined with the kinematic assumptions of a mixed‐field shear‐layerwise shell laminate theory. Based on the above formulation, a finite element methodology together with an incremental‐iterative technique, based on Newton–Raphson method is formulated. An eight‐node coupled non‐linear shell element is also developed. Various evaluation cases on laminated curved beams and cylindrical panels illustrate the capability of the shell finite element to predict the complex non‐linear behaviour of active shell structures including buckling, which is not captured by linear shell models. The numerical results also show the inherent capability of piezoelectric shell structures to actively induce large displacements through piezoelectric actuators, by jumping between multiple equilibrium states. Copyright © 2005 John Wiley & Sons, Ltd.     

Fourier methods for quasi‐periodic oscillations

Quasi‐periodic oscillations and invariant tori play an important role in the study of forced or coupled oscillators. This paper presents two new numerical methods for the investigation of quasi‐periodic oscillations. Both algorithms can be regarded as generalizations of the averaging and the harmonic (spectral) balance methods. The algorithms are easy to implement and require only minimal a priori knowledge of the system. Most importantly, the methods do not depend on an a priori co‐ordinate transformation. The methods are applied to a number of illustrative examples from non‐linear electrical engineering and the results show that the methods are efficient and reliable. In addition, these examples show that the presented algorithms can also continue through regions of sub‐harmonic (phase‐locked) resonance even though they are designed only for the quasi‐periodic case. Copyright © 2006 John Wiley & Sons, Ltd.     

High‐order finite elements applied to the discrete Boltzmann equation

A discontinuous Galerkin approach for solving the discrete Boltzmann equation is presented, allowing to compute approximate solutions for fluid flow problems. Based on a two‐dimensional high‐order finite element and an explicit Euler time stepping scheme, the D2Q9 model is discretized and the results are compared to the exact solution of the Navier–Stokes equation. Four numerical examples are considered, including stationary and instationary problems with curved boundaries. It is demonstrated that the proposed method allows to obtain the desired, highly efficient exponential convergence. Copyright © 2006 John Wiley & Sons, Ltd.     

Two‐phase thermo‐mechanical and macrosegregation modelling of binary alloys solidification with emphasis on the secondary cooling stage of steel slab continuous casting processes

As an approach towards a better modelling of solidification problems, we introduce a thermo‐mechanical and macrosegregation model that considers a solidifying alloy as a binary mixture made of a liquid and a solid phase. Macroscopic conservation laws for mass, momentum and solute are obtained by spatial averaging of the respective microscopic conservation equations. Assuming local thermal equilibrium, a single equation for the conservation of the mixture energy is then written. A single equation can be obtained for the solute as well by invoking a proper microsegregation rule. The numerical implementation in a two‐dimensional finite element code is then detailed. Lastly, some examples of simulations of academic tests as well as industrial applications for continuous casting of steel slabs are discussed. They particularly enlighten the ability of the formulation to describe the formation of central macrosegregation during the secondary cooling of slab continuous casting processes. Copyright © 2006 John Wiley & Sons, Ltd.     

Preservation of strawberries by modified atmosphere packages with other treatments

The impact of modified atmosphere packaging (MAP), in combination with ozone treatment and an edible film coating, on the preservation of strawberries was studied. A combination of 2.5% O₂ with 15% CO₂ was the optimum gas composition for the MAP of strawberries and prolonged shelf‐life by 4–6 days. When the strawberries were pre‐treated with a 4.3 mg/m³ ozone water dip and an edible coating prior to packaging, the shelf‐life was extended to 8–10 days. The results showed that two treatments [No.1: polyvinyl alcohol 134 (2%); monostearatacylglycerol (0.5%); phytic acid (0.05%); sorbitol (0.05%); sodium alginate (0.1%); absolute alcohol (8%): No.2: polyvinyl alcohol 134(1%); soluble starch (1%); glucose (1%); sucrose (1%); sodium alginate (0.1%); sorbitol (0.05%)] gave the best effects on strawberry quality and shelf‐life. All treatments (Nos 1–4) were obviously superior to check experiments for soluble sugar, ascorbic acid, acidity and anthocyanin. Furthermore, the treatments gave superior results for browning index, commodity ratio and permeability of the cell membrane. Copyright © 2006 John Wiley & Sons, Ltd.     

Classes of Split‐plot Response Surface Designs for Equivalent Estimation

When planning an experimental investigation, we are frequently faced with factors that are difficult or time consuming to manipulate, thereby making complete randomization impractical. A split‐plot structure differentiates between the experimental units associated with these hard‐to‐change factors and those that are relatively easy‐to‐change. Furthermore, it provides an efficient strategy that integrates the restrictions imposed by the experimental apparatus into the design structure. In this paper, several industrial and scientific examples are presented to highlight design considerations when a restriction on randomization is encountered. We propose classes of split‐plot response designs that provide an intuitive and natural extension from the completely randomized context. For these designs, the ordinary least‐squares estimates of the model are equivalent to the generalized least‐squares estimates. This property provides best linear unbiased estimators and simplifies model estimation. The design conditions that provide equivalent estimation are presented and lead to design construction strategies to transform completely randomized Box–Behnken, equiradial and small composite designs into a split‐plot structure. Published in 2006 by John Wiley & Sons, Ltd.