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.     

Relationship between footwall composition,crustal contamination,and fluid–rock interaction in the Platreef,Bushveld Complex,South Africa

In the northern limb of the 2.06-Ga Bushveld Complex, the Platreef is a platinum group elements (PGE)-, Cu-, and Ni-mineralized zone of pyroxenite that developed at the intrusion margin. From north to south, the footwall rocks of the Platreef change from Archaean granite to dolomite, hornfels, and quartzite. Where the footwall is granite, the Sr-isotope system is more strongly perturbed than where the footwall is Sr-poor dolomite, in which samples show an approximate isochron relationship. The Nd-isotope system for samples of pyroxenite and hanging wall norite shows an approximate isochron relationship with an implied age of 2.17 ± 0.2 Ga and initial Nd-isotope ratio of 0.5095. Assuming an age of 2.06 Ga, the ɛNd values range from −6.2 to −9.6 (ave. −7.8, n = 17) and on average are slightly more negative than the Main Zone of the Bushveld. These data are consistent with local contamination of an already contaminated magma of Main Zone composition. The similarity in isotope composition between the Platreef pyroxenites and the hanging wall norites suggests a common origin. Where the country rock is dolomite, the Platreef has generally higher plagioclase and pyroxene δ ¹⁸O values, and this indicates assimilation of the immediate footwall. Throughout the Platreef, there is considerable petrographic evidence for sub-solidus interaction with fluids, and the Δ _{plagioclase–pyroxene} values range from −2 to +6, which indicates interaction at both high and low temperatures. Whole-rock and mineral δD values suggest that the Platreef interacted with both magmatic and meteoric water, and the lack of disturbance to the Sr-isotope system suggests that fluid–rock interaction took place soon after emplacement. Where the footwall is granite, less negative δD values suggest a greater involvement of meteoric water. Consistently higher values of Δ _{plagioclase–pyroxene} in the Platreef pyroxenites and hanging wall norites in contact with dolomite suggest prolonged interaction with CO₂-rich fluid derived from decarbonation of the footwall rocks. The overprint of post crystallization fluid–rock interaction is the probable cause of the previously documented lack of correlation between PGE and sulfide content on the small scale. The Platreef in contact with dolomite is the focus of the highest PGE grades, and this suggests that dolomite contamination played a role in PGE concentration and deposition, but the exact link remains obscure. It is a possibility that the CO₂ produced by decarbonation of assimilated dolomite enhanced the process of PGE scavenging by sulfide precipitation.     

三峡工程对河口及邻近海域渔业影响的初步探讨

河口只是河流与海洋之间的一小部分水域,但对人类社会的发展和生物的进化起着重要的作用。它不仅是生物种群繁殖、育幼和栖息的场所,而且是溯河和降海鱼类洄游的必经之路,因此,河口对生物种群的繁衍延续和资源补充,以及保持生态平衡具有十分重要的意义。 然而,河口又是生态环境十分脆弱和敏感的水域,同时也是受人类干扰最为严重的区域。因此,河口的生态环境问题愈来愈引起人们的关注,其中包括长江三峡水库的兴建对长江口及其附近海域生态、环境的影响,特别是对位于河口两侧的我国著名两大渔场——舟山渔场和吕四渔场的影响问题。影响渔业资源的因子是多方面而且错综复杂的。本文主要通过长江径流的变化与环境因子变量的相关分析,并根据三峡水库不同蓄水方案,对河口近海渔业的影响问题进行初步探讨。     

南海北部陆缘地壳结构探测结果分析

深部地震和重力资料反演揭示了南海北部陆缘地壳结构在总体上由北部的华南沿海（厚约30km）向南部的洋盆（5──8km）逐渐减薄。南海的近SN向拉张不仅造成南北方向地壳结构的巨大变化，也造成东西向的明显变化。在南海北部陆缘的西部，局部拉张产生了一系列裂谷构造。西沙海槽作为一条狭窄的陆内裂谷向西延伸，海槽南北两侧地壳厚度超过25km，海槽中部地壳减薄至不足10km。西端的莺歌海盆地地壳厚仅5km，缺少明显的壳内反射－折射。在珠江口盆地中部，地壳厚度在下陆坡明显减薄，地壳下部存在较薄的（3──4km）高速层（地震波速7．2──7．5km·s－1）；在珠江口盆地东部，地壳底部存在约 10km厚、300km宽的高速层。在台湾地区，由于弧陆碰撞，曾经减薄的陆壳在碰撞带增厚，莫霍面深度超过30km。南海北部陆缘在裂谷拉张和海底扩张期间岩浆活动平静，表明南海北部陆缘为非火山型陆缘。     

浅析省级基础测绘数据在城市建模中的应用

利用基础测绘数据中房屋数据获取和编辑的方法,探讨1:10 000基础测绘数据中的建筑物在城市地理信息服务和建模中的应用以及目前基础测绘数据与城市建模数据间的差距.     

基于免疫多智能体的网络入侵主动防御模型

构建了一种基于免疫多智能体的网络入侵主动防御模型ADNIIMA(the active defense model for net-work intrusion based on immune multi-agent),提出免疫智能体概念,建立免疫智能体的逻辑结构及其运行机制,实现了对网络入侵的多层次、分布式主动防御机制,为网络安全保障提供一种新的思路.     

Neoproterozoic alkaline magmatism at Sivamalai, southern India

The Sivamalai alkaline complex lies at the southern margin of the Cauvery Shear System that separates the Archaean and Proterozoic domains of the Southern Granulite Terrain in India. U–Pb TIMS dating of zircon from a pegmatitic syenite sample in the complex yields a concordant age of 590.2 ± 1.3 (2σ) Ma which is interpreted to date the intrusion of the alkaline rocks. A lower concordia intercept at 168 ± 210 Ma defined by two grains with high common lead may indicate post-magmatic disturbances due to recrystallisation which is also evident in the CL images of the zircons. EPMA dating of monazite from a post-kinematic pegmatite which intrudes the crystalline basement hosting the alkaline rocks yields an age of 478 ± 29 (2σ) Ma and provides a lower bracket for the main phase of tectonism in this part of the Southern Granulite Terrain. The Pan-African high-grade metamorphism and ductile deformation has thus most likely affected the alkaline rocks. This is supported by the presence of a metamorphic foliation and extensive recrystallisation textures seen in the rocks. The major and trace element concentrations measured on selected samples reveals the presence of both enriched and depleted rock types. The enriched group includes ferrosyenite and nepheline syenite while the depleted group has only nepheline syenites. The trace element depletion of some nepheline syenites is interpreted to be a result of fractional crystallization involving the removal of accessory phases like zircon, titanite, apatite and allanite.     

淄博市洪山、寨里煤矿地下水串层污染形成原因及防治

淄博市洪山、寨里煤矿区隐伏的奥陶系灰岩裂隙岩溶水污染严重。水质污染突发于矿坑全面闭坑以后，污染范围与矿坑水区一致，地下水污染组份与矿坑水的高含量组份相同。水质恶化主要由矿坑水串层污染所致。在此针对串层污染形成原因及影响因素，提出了防治对策。     

二维TIM中的传播矩阵与面波相速度

本文利用横向各向同性介质（ＴＩＭ）中的传播矩阵和面波形成条件，导出了面波波慢度方程。对瑞利波给出了各向同性及横向各向同性弹性介质中其相速度的解析解；对于勒夫波，给出了横向各向同性介质中的波慢度方程。     

Stochastic Modelling of the Impact of Flood Protection Measures Along the River Waal in the Netherlands

River flooding is a problem of international interest. In the past few years many countries suffered from severe floods. A large part of the Netherlands is below sea level and river levels. The Dutch flood defences along the river Rhine are designed for water levels with a probability of exceedance of 1/1250 per year. These water levels are computed with a hydrodynamic model using a deterministic bed level and a deterministic design discharge. Traditionally, the safety against flooding in the Netherlands is obtained by building and reinforcing dikes. Recently, a new policy was proposed to cope with increasing design discharges in the Rhine and Meuse rivers. This policy is known as the Room for the River (RfR) policy, in which a reduction of flood levels is achieved by measures creating space for the river, such as dike replacement, side channels and floodplain lowering. As compared with dike reinforcement, these measures may have a stronger impact on flow and sediment transport fields, probably leading to stronger morphological effects. As a result of the latter the flood conveyance capacity may decrease over time. An a priori judgement of safety against flooding on the basis of an increased conveyance capacity of the river can be quite misleading. Therefore, the determination of design water levels using a fixed-bed hydrodynamic model may not be justified and the use of a mobile-bed approach may be more appropriate. This problem is addressed in this paper, using a case study of the river Waal (one of the Rhine branches in the Netherlands). The morphological response of the river Waal to a flood protection measure (floodplain lowering in combination with summer levee removal) is analysed. The effect of this measure is subject to various sources of uncertainty. Monte Carlo simulations are applied to calculate the impact of uncertainties in the river discharge on the bed levels. The impact of the “uncertain” morphological response on design flood level predictions is analysed for three phenomena, viz. the impact of the spatial morphological variation over years, the impact of the seasonal morphological variation and the impact of the morphological variability around bifurcation points. The impact of seasonal morphological variations turns out to be negligible, but the other two phenomena appear to have each an appreciable impact (order of magnitude 0.05–0.1 m) on the computed design water levels. We have to note however, that other sources of uncertainty (e.g. uncertainty in hydraulic roughness predictor), which may be of influence, are not taken into consideration. In fact, the present investigation is limited to the sensitivity of the design water levels to uncertainties in the predicted bed level.