An eight‐node hybrid‐stress solid‐shell element for geometric non‐linear analysis of elastic shells

This paper presents eight‐node solid‐shell elements for geometric non‐linear analysis of elastic shells. To subdue shear, trapezoidal and thickness locking, the assumed natural strain method and an ad hoc modified generalized laminate stiffness matrix are employed. A selectively reduced integrated element is formulated with its membrane and bending shear strain components taken to be constant and equal to the ones evaluated at the element centroid. With the generalized stresses arising from the modified generalized laminate stiffness matrix assumed to be independent from the ones obtained from the displacement, an extended Hellinger–Reissner functional can be derived. By choosing the assumed generalized stresses similar to the assumed stresses of a previous solid element, a hybrid‐stress solid‐shell element is formulated. Commonly employed geometric non‐linear homogeneous and laminated shell problems are attempted and our results are close to those of other state‐of‐the‐art elements. Moreover, the hybrid‐stress element converges more readily than the selectively reduced integrated element in all benchmark problems. Copyright © 2002 John Wiley & Sons, Ltd.     

A Refinement Calculus for Shared-Variable Parallel and Distributed Programming

Parallel computers have not yet had the expected impact on mainstream computing. Parallelism adds a level of complexity to the programming task that makes it very error-prone. Moreover, a large variety of very different parallel architectures exists. Porting an implementation from one machine to another may require substantial changes. This paper addresses some of these problems by developing a formal basis for the design of parallel programs in the form of a refinement calculus. The calculus allows the stepwise formal derivation of an abstract, low-level implementation from a trusted, high-level specification. The calculus thus helps structuring and documenting the development process. Portability is increased, because the introduction of a machine-dependent feature can be located in the refinement tree. Development efforts above this point in the tree are independent of that feature and are thus reusable. Moreover, the discovery of new, possibly more efficient solutions is facilitated. Last but not least, programs are correct by construction, which obviates the need for difficult debugging. Our programming/specification notation supports fair parallelism, shared-variable and message-passing concurrency, local variables and channels. The calculus rests on a compositional trace semantics that treats shared-variable and message-passing concurrency uniformly. The refinement relation combines a context-sensitive notion of trace inclusion and assumption-commitment reasoning to achieve compositionality. The calculus straddles both concurrency paradigms, that is, a shared-variable program can be refined into a distributed, message-passing program and vice versa. Received July 2001 / Accepted in revised form May 2002     

The Logarithmic Image Processing Model: Connections with Human Brightness Perception and Contrast Estimators

The logarithmic image processing (LIP) model is amathematical framework based on abstract linear mathematicswhich provides a set of specific algebraic and functionaloperations that can be applied to the processing of intensityimages valued in a bounded range. The LIP model has been provedto be physically justified in the setting of transmitted lightand to be consistent with several laws and characteristics ofthe human visual system. Successful application examples havealso been reported in several image processing areas, e.g.,image enhancement, image restoration, three-dimensional imagereconstruction, edge detection and image segmentation.The aim of this article is to show that the LIP model is atractable mathematical framework for image processing which isconsistent with several laws and characteristics of humanbrightness perception. This is a survey article in the sensethat it presents (almost) previously published results in arevised, refined and self-contained form. First, an introductionto the LIP model is exposed. Emphasis will be especially placedon the initial motivation and goal, and on the scope of themodel. Then, an introductory summary of mathematicalfundamentals of the LIP model is detailed. Next, the articleaims at surveying the connections of the LIP model with severallaws and characteristics of human brightness perception, namelythe brightness scale inversion, saturation characteristic, Weber'sand Fechner's laws, and the psychophysical contrast notion. Finally,it is shown that the LIP model is a powerful and tractable framework for handling the contrast notion. This is done througha survey of several LIP-model-based contrast estimators associated with special subparts (point, pair of points,boundary, region) of intensity images, that are justified bothfrom a physical and mathematical point of view.     

线性时变系统的区间稳定性与鲁棒稳定性

本文应用向量比较定理研究线性时变系统的区间稳定性和具非线性时变摄动的线性时变系统的鲁棒稳定性，所得的新结果包含文献的一些主要结果作为特例，本文的研究方法说明向量比较方法是分析区间稳定性和鲁棒稳定性的一种自然而有力的工具。     

Best L1 Piecewise Monotonic Data Modelling

Let us consider n data measurements of a univariate process that have been altered by random errors. We assume that an underlying model function has a substantially smaller number of turning points than the observed ones. We propose algorithms that make least the sum of the moduli of the errors by requiring k monotonic sections, alternately increasing and decreasing, in the sequence of the smoothed values. The main difficulty in this calculation is that the optimal positions of the joins of the monotonic sections have to be found automatically among so many combinations that it is impossible to test each one separately. Moreover, the calculation seems to be very intractable to general optimization techniques because O(n^k) local minima can occur. It is shown that dynamic programming can be used for separating the data into optimal disjoint sections of adjacent data, where each section requires a single L₁ monotonic calculation. This procedure is highly efficient, requiring at most O(kn²) computer operations and O(n) best L₁ monotonic calculations to subranges of data for a global minimum.     

A feasibility study of a direct-mapping parallel processing method to solve linear equations in load-flow calculation

With the growing size and complexity of power systems, system analysis—such as transients calculation—takes much time. Hence, fast calculation methods are required. Although parallel processing is a hopeful method, there have been difficulties in the parallel solution of linear equations which appear in power-flow calculations by the Newton-Raphson method. This paper aims at the fast calculation of the power-flow problem by means of parallel processing. In order to improve the suitability to the parallel solution of the differential equation in transients calculation, we assume the use of a direct-mapping parallel processing machine to map directly the network of a power system onto a network of processors. Under this assumption, we propose a new parallel-processing-oriented method in which the linear equation is solved by linear iterations between nodes with Aitken acceleration. We simulate the method on three model power systems and compare this Parallel Iterative Method (PIN) with a Parallel Direct Method (PDM) which uses the banded matrix according to the number of operations required. As a result, we can expect that PIM may solve linear equations faster than PDM with m processors, although the PIM might be inferior to the PDM with m × m processors, where m denotes the half-band width of the banded matrix.     

二阶A型Cooper变换的增强型

对二阶Ａ型Ｃｏｏｐｅｒ变换进行增强，放宽可用性条件，并证明了其正确性，最后举例说明其应用。     

Validation of a Method of Determination of Apparent Diffusivity Versus Composition in Solids

A method is presented for computing the values of apparent diffusivity in solids with respect to the concentration of the diffusing substance (water or sodium chloride). This method does not require any assumption upon the mathematical relationship between diffusivity and concentration. It can be applied to experimental measurements of local concentration versus position within the solid (profiles) with relatively few measurements (circa 10) and a mathematical smoothing of the experimental data by using an artificial neural network model. The method was first validated on simulated data obtained by using a constant diffusivity value and on experimental profiles when the relation between diffusivity and concentration was given. It was then applied to original experimental moisture profiles obtained by putting gelatin gels with different initial moisture contents into contact for up to 14 days. The method was also successfully applied to five sets of experimental moisture and sodium chloride profiles taken from the literature and obtained from different food products. Apparent diffusivities calculated by our method were found in agreement with those obtained by authors using different numerical methods to compute the diffusivity values.     

Solving thermal and phase change problems with the eXtended finite element method

 The application of the eXtended finite element method (X-FEM) to thermal problems with moving heat sources and phase boundaries is presented. Of particular interest is the ability of the method to capture the highly localized, transient solution in the vicinity of a heat source or material interface. This is effected through the use of a time-dependent basis formed from the union of traditional shape functions with a set of evolving enrichment functions. The enrichment is constructed through the partition of unity framework, so that the system of equations remains sparse and the resulting approximation is conforming. In this manner, local solutions and arbitrary discontinuities that cannot be represented by the standard shape functions are captured with the enrichment functions. A standard time-projection algorithm is employed to account for the time-dependence of the enrichment, and an iterative strategy is adopted to satisfy local interface conditions. The separation of the approximation into classical shape functions that remain fixed in time and the evolving enrichment leads to a very efficient solution strategy. The robustness and utility of the method is demonstrated with several benchmark problems involving moving heat sources and phase transformations. Received 20 May 2001 / Accepted 19 December 2001