Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Invariance and scaling properties in the distributions of contributing area and energy in drainage basins

The cumulative probability distributions for stream order, stream length, contributing area, and energy dissipation per unit length of channel are derived, for an ordered drainage system, from Horton's laws of network composition. It is shown how these distributions can be related to the fractal nature of single rivers and river networks. Finally, it is shown that the structure proposed here for these probability distributions is able to fit the observed frequency distributions, and their deviations from straight lines in a log-log plot.     

The solubility of iron sulphides in synthetic and natural waters at ambient temperature

A critical evaluation of literature values for the solubility products, K _sp ^NBS = [Fe²⁺][HS^–] Fe²⁺ HS^– (H _NBS ⁺ )^–1, of various iron sulphide phases results in consensus values for the pKs of 2.95 ± 0.1 for amorphous ferrous sulphide, 3.6 ± 0.2 for mackinawite, 4.4 ± 0.1 for greigite, 5.1 ± 0.1 for pyrrhotite, 5.25 ± 0.2 for troilite and 16.4 ± 1.2 for pyrite.Where the analogous ion activity products have been measured in anoxic freshwaters in which there is evidence for the presence of solid phase FeS, the values lie within the range of 2.6–3.22, indicating that amorphous iron sulphide is the controlling phase. The single value for a groundwater of 2.65 (2.98 considering carbonate complexation) agrees. In seawater four values range between 3.85 to 4.2, indicating that mackinawite or greigite may be the controlling phase. The single low value of 2.94 is in a situation where particularly high fluxes of Fe (II) and S (–II) may result in the preferential precipitation of amorphous iron sulphide. Formation of framboidal pyrite in these sulphidic environments may occur in micro-niches and does not appear to influence bulk concentrations. Calculations show that the formation of Fe₂S₂ species probably accounts for very little of the iron or sulphide in most natural waters. Previously reported stability constants for the formation of Fe (HS)₂ and (Fe (HS)₃)^– are shown to be suspect, and these species are also thought to be negligible in natural waters. In completely anoxic pore waters polysulphides also have a negligible effect on speciation, but in tidal sediments they may reach appreciable concentrations and lead to the direct formation of pyrite. Concentrations of iron and sulphide in pore waters can be controlled by the more soluble iron sulphide phase. The change in the IAP with depth within the sediment may reflect ageing of the solid phase or a greater flux of Fe (II) and S (–II) nearer the sediment surface. This possible kinetic influence on the value of IAPs has implications for their use in geochemical studies involving phase formation.     

Quaternary fault segmentation and interaction in the epicentral area of the 1561 earthquake (Mw = 6.4), Vallo di Diano, southern Apennines, Italy

The main structural characteristics of the Caggiano and Polla faults, exposed in the epicentral area of the 1561 earthquake (Mw = 6.4), southern Italy, have been investigated in detail to assess their spatial and temporal properties, and to evaluate their seismogenic potential. These right stepping normal faults show an overlap of about 7 km and an across strike separation of about 4 km. The geometric relationships between the Caggiano and Polla faults, but also the displacement distribution along each fault, demonstrate that they have been strongly interacting throughout the Pleistocene. Nevertheless, geological evidence of Holocene tectonic activity was mainly recognized along the Caggiano Fault (faulted late glacial deposits) and in the southernmost part of the Polla Fault (faulted deposits of probably Late Pleistocene age). This suggests that the Caggiano Fault can be considered as the most tectonically active fault in the Vallo di Diano Fault System. By calculating Coulomb stress changes, we have constrained modes of mechanical interactions between the two faults in a scenario compatible with the 1561 earthquake. This approach allows us to argue that both the Caggiano and the Polla Faults are probably linked at depth, and part of the same seismogenic structure which may be potentially responsible for composite ruptures with magnitude ≥ 6.5.     

Revised age of proximal deposits in the Zagros foreland basin and implications for Cenozoic evolution of the High Zagros

The regionally extensive, coarse-grained Bakhtiyari Formation represents the youngest synorogenic fill in the Zagros foreland basin of Iran. The Bakhtiyari is present throughout the Zagros fold-thrust belt and consists of conglomerate with subordinate sandstone and marl. The formation is up to 3000 m thick and was deposited in foredeep and wedge-top depocenters flanked by fold-thrust structures. Although the Bakhtiyari concordantly overlies Miocene deposits in foreland regions, an angular unconformity above tilted Paleozoic to Miocene rocks is expressed in the hinterland (High Zagros).

The Bakhtiyari Formation has been widely considered to be a regional sheet of Pliocene–Pleistocene conglomerate deposited during and after major late Miocene–Pliocene shortening. It is further believed that rapid fold growth and Bakhtiyari deposition commenced simultaneously across the fold-thrust belt, with limited migration from hinterland (NE) to foreland (SW). Thus, the Bakhtiyari is generally interpreted as an unmistakable time indicator for shortening and surface uplift across the Zagros. However, new structural and stratigraphic data show that the most-proximal Bakhtiyari exposures, in the High Zagros south of Shahr-kord, were deposited during the early Miocene and probably Oligocene. In this locality, a coarse-grained Bakhtiyari succession several hundred meters thick contains gray marl, limestone, and sandstone with diagnostic marine pelecypod, gastropod, coral, and coralline algae fossils. Foraminiferal and palynological species indicate deposition during early Miocene time. However, the lower Miocene marine interval lies in angular unconformity above ~ 150 m of Bakhtiyari conglomerate that, in turn, unconformably caps an Oligocene marine sequence. These relationships attest to syndepositional deformation and suggest that the oldest Bakhtiyari conglomerate could be Oligocene in age.

The new age information constrains the timing of initial foreland-basin development and proximal Bakhtiyari deposition in the Zagros hinterland. These findings reveal that structural evolution of the High Zagros was underway by early Miocene and probably Oligocene time, earlier than commonly envisioned. The age of the Bakhtiyari Formation in the High Zagros contrasts significantly with the Pliocene–Quaternary Bakhtiyari deposits near the modern deformation front, suggesting a long-term (> 20 Myr) advance of deformation toward the foreland.     

伶仃洋沉积动力特点的研究

位于珠江三角洲东侧的伶仃洋，因径流下泄与潮流进退的流向不一，使它各分流口门的出口水道都有主槽和支槽之分，即都有主干水道和分汊水道。陆架高盐海水入侵又使伶仃洋内沉积动力过程在空间分布上发生差异，如沉积物分布有粗－细－稍粗之分；而水体中的密度、速度差异，常常产生锋带，对水下地形的发展有不可忽视的影响。因而全面认识发生在伶仃洋内的沉积动力作用，对深水航道的选线极为重要。     

钦州湾潮流特征分析

根据1994年5月和11—12月两个航次的调查及近年有关海湾的海流观测资料，分析了钦州湾的潮流特征，观测资料来自分布于该海湾东中西三个航道及相关区域的8个测站。钦州湾涨落潮流特征：落潮流速大于涨潮流速：东部最大涨落潮流速小于西部：夏季落潮流总是大于冬季：龙门水道附近的流速最大，而其余区域的流速相对较小。外湾(钦州湾)余流，是气旋式环流：水体东进而西出。从湾内来的泥沙和污染物质，主要从西部进入外海；中、东槽是主要“进水”通道，外海低泥沙含量、少污染物质的清洁水从这里潮流而上。因此，与西部相比，东部盐度高，水质清洁，底质重金属含量少。     

Semi-Probabilistic Evaluation of System Redundancy for Existing Offshore Structures

In this paper, a new definition of structure system redundancy is proposed in view of the various measures for structure redundancy. By introducing the terms of structure system failure at the mechanism level and equivalent reliability index, the safety for existing offshore platforms can be evaluated by the semi-probabilistic method presented in this paper. Some numerical examples are given and satisfactory results have been obtained.     

郯庐断裂带及其周缘中新生代盆地发育特征

郯庐断裂带作为中国东部滨太平洋地区一条巨型走滑构造带，对其周缘中、新生代盆地的发育、演化起着重要的控制作用。随着太平洋板块俯冲方向从NNW向NW到NWW的变化，郯庐断裂带的活动方式逐步从中生代左行走滑-左行斜向滑动过渡到早第三纪以左行斜向-倾向滑动，晚第三纪-第四纪转为倾滑-右行斜向滑动-右行走滑。走滑活动经历了一个循序渐进的周期演变过程。随着郯庐断裂活动方式的演变，其周缘中、新生代盆地的发育逐渐向北迁移，其中南段周缘盆地主要为中生代盆地，中段周缘盆地主要为中、新生代叠加盆地，中北段周缘盆地主要为早第三纪盆地。每个盆地都经历了拉分(伸展)裂陷到挤压反转的演化过程。此外，在同一时期、同一区域剪切应力场作用下，不同区段因其走向变化导致局部应力场变化，在增压弯曲部位发生会聚、挤压、隆升；而释压拉张部位发生离散、伸展、沉降，从而盆地发育。