Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

A neural network approach to ship track-keeping control

This paper presents an on-line trained neural net work controller for ship track-keeping problems. Following a brief review of the ship track-keeping control development since the 1980's, an analysis of various existing backpropagation-based neural controllers is provided. We then propose a single-input multioutput (SIMO) neural control strategy for situations where the exact mathematical dynamics of the ship are not available. The aim of this study is to build an autonomous neural controller which uses rudder to regulate both the tracking error and heading error. During the whole control process, the proposed SIMO neural controller adapts itself on-line from a direct evaluation of the control accuracy, and hence the need for a “teacher” or an off-line training process can be removed. With a relatively modest amount of quantitative knowledge of the ship behavior, the design philosophy enables real time control of a nonlinear ship model under random wind disturbances and measurement noise. Three different track-keeping tasks have been simulated to demonstrate the effectiveness of the training method and the robust performance of the proposed neural control strategy     

图象的模型及其统计特性

本文对图象的统计特性及其适用模型作了分析。讨论了以前常用的平稳图象模型的缺陷,对实际图象的统计特性作了分析,认为它是非平稳的,不满足各态历经性,同时在空间上是高度相关的;讨论了相应的随机参数统计模型和描述性统计模型。指出作用于整幅图象的“全局”图象复原算法比每次只独立计算单个象元的“点”图象复原算法更为优越。     

OQ208, a non-ideal VLBI calibrator of flux density

OQ208 is used ore and more frequently as a VLBI calibrator of flux density. A VLBI image of OQ208 at 5 GHz is firstly provided in this paper. Its structure consists of a compact core and a component. The angular distance between the core and the component is about 7 milli-arcsecond; the position angle is −125 degrees. In principle, OQ208 is not an ideal VLBI calibrator of flux density. This paper provides the structural information for correcting the structure effect when OQ208 is to be used as a calibrator.     

Exposure dating and reinterpretation of coarse debris accumulations (‘rock glaciers’) in the Cairngorm Mountains,Scotland

Relict rock glaciers have considerable potential for contributing to palaeoclimatic reconstruction, but this potential is often undermined by lack of dating control and problems of interpretation. Here we reinvestigate and date four proposed ‘rock glaciers’ in the Cairngorm Mountains and show that the morphology of only one of these appears consistent with that of a true rock glacier produced by creep of underlying ice or ice‐rich sediment. All four features comprise rockslide or rock avalanche runout debris, and the possibility that all four represent unmodified runout accumulations cannot be discounted. Surface exposure dating of the four debris accumulations using cosmogenic ¹⁰Be produced uncertainty‐weighted mean ages of 15.4 ± 0.8 ka, 16.2 ± 1.0 ka, 12.1 ± 0.6 ka and 12.7 ± 0.8 ka. All four ages imply emplacement under cold stadial conditions, two prior to the Windermere Interstade of ca. 14.5–12.9 cal. ka BP and two during the Loch Lomond Stade of ca. 12.9–11.5 cal. ka BP. The above ages indicate that paraglacial rock‐slope failure on granite rockwalls occurred within a few millennia after deglaciation. The mean exposure ages obtained for runout debris at two sites – Strath Nethy (16.2 ± 1.0 ka) and Lairig Ghru (15.4 ± 0.8 ka) – are consistent with basal radiocarbon ages from Loch Etteridge, 22 km to the southwest (mean = 15.6 ± 0.3 cal. ka BP) and imply widespread deglaciation of the Cairngorms and adjacent valleys before 15 ka and possibly 16 ka. Copyright © 2008 John Wiley & Sons, Ltd.