求解非稳定对流占优水质模型的非标准Galerkin有限元法

在求解非稳定地下水溶质运移模型时,若对流项占优,则模型表现出双曲方程的特性。针对这种特性,采用非标准Galerkin有限元方法进行求解是解决这类问题的有效途径。分别采用Wavelet-Galer-kin有限元方法、迎风有限元方法和特征有限元方法对强对流溶质运移模型进行了求解,并将其结果与标准Galerkin有限元和解析解进行对比。结果表明:标准Galerkin有限元方法会产生强烈的数值振荡;Wavelet-Galerkin有限元方法的时空定位效果好;迎风有限元方法能够有效降低数值振荡现象,但迎风因子对解的影响较大,而且会带来时间延迟;特征有限元方法能够提高解的精度,故可以认为特征有限元方法是求解强对流地下水溶质运移模型的首选方法。     

双排桩-锚杆支护的有限元模拟

基坑支护模型往往在计算假定的基础上建立,而依据计算假定建立的模型往往不能得到和实际相符的结果。但有限元软件可以较为真实的模拟基坑开挖的过程中应力和位移的变化情况,是土工数值分析的重要手段。本文运用大型有限元软件ABAQUS建立双排桩-锚杆结构有限元模型,全面的考虑土体的特性、桩土的共同作用及桩间土对支护结构的影响等因素,分析支护结构在土体开挖荷载作用下的内力和变形,为设计和施工提供了参考。     

软土地下工程抗震数值模拟的若干关键问题

针对我国地下工程建设项目抗震安全性能设计和评估的需求,并基于目前缺少完善的地下结构抗震分析方法和专门的地下结构抗震设计规范的现状,重点阐述了采用有限元方法进行地下工程数值模拟和分析时需要特别关注的三方面问题：计算范围、土一结构动力相互作用和边界条件。着重分析了有限元网格的划分、土一结非线性接触问题和不同边界条件的适用性,以及这些建模条件对计算结果可靠性的影响。对基于数值手段开展地下工程抗震安全性能的计算、分析和评价,制定地下结构抗震设计规范或规程,可提供有益参考。     

浅海导管架平台结构损伤诊断试验研究

基于模型试验研究了浅海导管架平台的结构损伤诊断,采用标量型ARMA模型识别出结构的三阶模态参数,分别采用遗传算法和模态应变能法对试验模型不同模拟损伤工况进行了损伤诊断。基于遗传算法,提出了采用不同自由度一阶模态参数的损伤诊断方法。研究表明,由于海洋平台的前三阶模态分别为不同自由度的一阶模态,采用两阶或三阶模态识别损伤时,损伤杆件影响较小的自由度方向的模态参数对结构的损伤诊断将产生不利的影响,成为识别噪声。因此,可分别采用不同自由度的一阶模态进行结构损伤诊断。     

Qualitative and quantitative analysis of vorticity,strain and area change in general non-isochoric 2D deformation

A scale free representation of a general non-isochoric 2D deformation is presented which is amenable to mathematical analysis. By describing deformation in 2D in terms of polar coordinates the stretching and rotational histories of linear elements separate and are easily analysed both qualitatively and quantitatively. An analysis of finite strain combined with dynamical considerations allows the derivation of equations which may be used to estimate finite strain, area change and kinematic vorticity number. Numerical investigation of method developed here was carried out and it was found to perform well unless large area changes occur in combination with large components of simple shear. A re-analysis of natural data indicates the method is consistent.     

Transtensional folding

Strain modeling shows that folds can form in transtension, particularly in simple shear-dominated transtension. Folds that develop in transtension do not rotate toward the shear zone boundary, as they do in transpression; instead they rotate toward the divergence vector, a useful feature for determining past relative plate motions. Transtension folds can only accumulate a fixed amount of horizontal shortening and tightness that are prescribed by the angle of oblique divergence, regardless of finite strain. Hinge-parallel stretching of transtensional folds always exceeds hinge-perpendicular shortening, causing constrictional fabrics and hinge-parallel boudinage to develop.These theoretical results are applied to structures that developed during oblique continental rifting in the upper crust (seismic/brittle) and the ductile crust. Examples include (1) oblique opening of the Gulf of California, where folds and normal faults developed simultaneously in syn-divergence basins; (2) incipient continental break-up in the Eastern California-Walker Lane shear zone, where earthquake focal mechanisms reflect bulk constrictional strain; and (3) exhumation of the ultrahigh-pressure terrain in SW Norway in which transtensional folds and large magnitude stretching developed in the footwall of detachment shear zones, consistent with constrictional strain. More generally, folds may be misinterpreted as indicating convergence when they can form readily in oblique divergence.     

Strain distribution within a km-scale,mid-crustal shear zone: The Kuckaus Mylonite Zone,Namibia

The subvertical Kuckaus Mylonite Zone (KMZ) is a km-wide, crustal-scale, Proterozoic, dextral strike-slip shear zone in the Aus granulite terrain, SW Namibia. The KMZ was active under retrograde, amphibolite to greenschist facies conditions, and deformed felsic (and minor mafic) gneisses which had previously experienced granulite facies metamorphism during the Namaqua Orogeny. Lenses of pre- to syn-tectonic leucogranite bodies are also deformed in the shear zone. Pre-KMZ deformation (D₁) is preserved as moderately dipping gneissic foliations and tightly folded migmatitic layering. Shear strain within the KMZ is heterogeneous, and the shear zone comprises anastomosing high strain ultramylonite zones wrapping around less deformed to nearly undeformed lozenges. Strain is localized along the edge of leucogranites and between gneissic lozenges preserving D₁ migmatitic foliations. Strain localization appears controlled by pre-existing foliations, grain size, and compositional anisotropy between leucogranite and granulite. The local presence of retrograde minerals indicate that fluid infiltration occurred in places, but most ultramylonite in the KMZ is free of retrograde minerals. In particular, rock composition and D₁ fabric heterogeneity are highlighted as major contributors to the strain distribution in time and space, with deformation localization along planes of rheological contrast and along pre-existing foliations. Therefore, the spatial distribution of strain in crustal-scale ductile shear zones may be highly dependent on lithology and the orientation of pre-existing fabric elements. In addition, foliation development and grain size reduction in high strain zones further localizes strain during progressive shear, maintaining the anastomosing shear zone network established by the pre-existing heterogeneity.     

Strain rates at high temporal resolution from curved inclusion trails in garnet,Passo del Sole,Central Swiss Alps

Quantitative strain rates at outcrop scale are very difficult to obtain, but they may be estimated from crystals with curved inclusion trails by calculating rotation rates from growth rates and corresponding deflections of the internal foliation. Garnet in a quartzose pelite at Passo del Sole in the central Swiss Alps is extraordinarily valuable for calculation of strain rates during Alpine orogenesis, because the unusual zoning patterns clearly define the kinetics of its nucleation and growth. Complex concentric zoning patterns can be correlated from one crystal to another in a hand sample, based on compositional and microstructural similarities; the ubiquity of these features demonstrates that all garnet crystals nucleated at nearly the same time. Compositional bands whose radial widths are proportional to crystal size provide evidence for growth governed by the kinetics of intergranular diffusion of locally sourced nutrients. Together, these constraints increase the reliability of estimates of rates of garnet growth, and the strain‐rate calculations that depend on them. To obtain growth rates, P–T conditions during garnet crystallization were modelled in a series of pseudosections, and compositional evolution was connected to rates of garnet growth by means of an independently determined heating rate. These growth rates, combined with measured amounts of curvature of inclusion trails, indicate that the time‐averaged strain rate at Passo del Sole during Alpine metamorphism was on the order of 10^?14 s^?1. Strain rates calculated using rotational v. non‐rotational models are similar in magnitude. The constraints on crystallization kinetics also allow direct calculation of strain rates during individual stages of garnet growth, revealing short‐term increases to values on the order of 10^?13 s^?1. These higher strain rates are correlated with the growth of concentric high‐Ca or high‐Mn zones in garnet, which implies that strain softening associated with the transient passage of fluids is responsible for acceleration of deformation during these intervals.     

Crust restoration for the western Lachlan Orogen using the strain-reversal,area-balancing technique: implications for crustal components and original thicknesses

Strain reversal of structural/stratigraphic profiles at different scales in the western Lachlan Orogen provides a perspective on original crustal thickness estimates, the former depositional basin width of the proto-western Lachlan Orogen, the original sedimentary-fan thickness, and the possible length extent of lower crust lost by subduction. Retrodeformation using strain-reversal techniques allows basin reconstruction giving an original width of the western Lachlan Orogen basin receptor of between 800 km (minimum) and ～1150 km (maximum), depending on the amount of stratal duplication allowed in the turbidites. Crude area balancing of the regional cross-section, adding in sectional volume lost by erosion and assuming strain compatibility between the upper and lower crust, suggests that the predeformation crustal thickness ranges between 15 km and ～21 km, with a lower crustal thickness (oceanic lithosphere) of ～9 km and a turbidite fan thickness of ～6 km (minimum) and ～12 km (maximum allowable), respectively. Disparity between the calculated fan thickness and that derived from measured stratigraphic sections adjusted for strain (～6 km) indicates that some form of crustal stacking must be important in structural thickening of the turbidite crustal component. By varying shortening due to fault stacking, mass balance dictates the mismatch of the upper crustal (uc) and lower crustal (lc) retrodeformed lengths, and therefore provides an estimate of lower crustal loss by subduction. End members range from: (i) a 12 km-thick fan without fault duplication, a basin width of ～800 km where uc = lc giving no lower crustal loss by subduction; to (ii) a ～6 km fan, requiring duplication by faulting, a basin of ～1150 km where uc > lc, and ～360 km of lower crust length (～30%) lost by subduction. This suggests that the total thickness of underplated igneous material in the western Lachlan Orogen is low, probably < ～2 km.     

钱塘江古海塘塘基防渗加固现场试验研究

目前仍然坚守在防洪御潮第一线的钱塘江明清古海塘,由于年代久远,部分塘段的塘趾部位存在流土现象,因此,选用合理且不破坏古海塘风貌的塘基防渗加固方案来保护古海塘十分必要。以浙江海宁某明清古海塘试验段为依托,通过现场试验,对比研究了4种塘基防渗加固保护方案,通过有限元分析、现场渗透试验、加固后开挖检验、施工期对古海塘的影响监测等手段,分析总结了各加固方案的优缺点,得到了机械适当改制后的“塘趾旋喷注浆”方案为优选方案的结论,研究成果对下一步钱塘江临江古海塘保护工程的塘基防渗加固方案选取具有重要参考价值。