A multiscale finite element formulation for axisymmetric elastoplasticity with volumetric locking

Similar to plane strain, axisymmetric stress problem is also highly kinematics constrained. Standard displacement‐based finite element exhibits volumetric locking issue in simulating nearly/fully incompressible material or isochoric plasticity under axisymmetric loading conditions, which severely underestimates the deformation and overestimates the bearing capacity for structural/geotechnical engineering problems. The aim of this paper is to apply variational multiscale method to produce a stabilized mixed displacement–pressure formulation, which can effectively alleviate the volumetric locking issue for axisymmetric stress problem. Both nearly incompressible elasticity and isochoric J₂ elastoplasticity are investigated. First‐order 3‐node triangular and 4‐node quadrilateral elements are tested for locking issues. Several representative simulations are provided to demonstrate the performance of the linear elements, which include the convergence study and comparison with closed‐form solutions. A comparative study with pressure Laplacian stabilized formulation is also presented. Copyright © 2009 John Wiley & Sons, Ltd.     

水平层状岩体边坡动力响应中的结构面效应研究

根据应力波传播原理分析了水平层状岩体边坡中应力波传播特征,建立了应力波在该类边坡中传播的模型.利用离散元软件UDEC分析了不同频率垂向压缩应力波作用下边坡动力响应规律中的结构面效应.结果表明:边坡中的水平层面对坡顶的动力响应有明显影响.低频应力波作用下,水平层状岩体边坡坡顶的垂向峰值速度较均质坡体相同部位的峰值速度的增...     

MS8.0汶川地震对MS7.0芦山地震成核失稳的影响

依据速率-状态依赖性摩擦本构关系,并结合中下地壳和上地幔的黏弹性松弛效应,以震后库仑破裂应力变化和同震动态库仑破裂应力变化的计算为基础,模拟了两种应力变化对芦山地震断层的失稳发震时间的影响,研究了2008年M_S8.0汶川地震与2013年M_S7.0芦山地震之间的触发关系.计算得到,汶川地震在芦山地震震源断层面上产生的动态应力变化的峰值为0.127 MPa;此外,经过近5年中下地壳和上地幔的黏弹性松弛效应,芦山地震震源断层面上受到的震后应力变化值为0.025 MPa.结果表明,芦山地震的震源断层在应力积累逐渐接近临界状态的某一特定时期内,受到了汶川地震产生的动态应力变化、静态应力变化以及黏弹性松弛效应造成应力变化的共同触发作用,且动态应力的延迟触发作用可能更为显著.最后对芦山地震之后研究区域的应力变化场进行了初步探讨.     

大理岩破坏的声发射双主频特征及其机制初探

因试样制备、试验操作等限制,获得饱和岩石的直接拉伸强度较为困难。为对饱和硬岩的直接拉伸强度进行合理估算并验证其可靠性,同时分析水对估算误差的影响,设计完成了干燥、饱和状态下两种硬岩的单轴压缩和直接拉伸声发射试验。基于Griffith理论和Hoek-Brown(H-B)强度准则,通过裂纹体积应变法确定岩石起裂应力,估算饱和硬岩的直接拉伸强度和H-B强度参数m_i。结合声发射波形信号频谱特性,分析两者估算结果的准确性。结果表明,基于Griffith准则并采用起裂应力估算饱和硬岩的直接拉伸强度是一种较可靠的方法。对于干燥硬岩,该方法确定的直接拉伸强度估算值比试验值偏小。饱和硬岩直接拉伸强度估算值与试验值基本接近,这是因为硬岩饱和后,延性和蠕变趋向增强,并发生更多的微观拉破裂。饱和硬岩的m_i值增大与增多的微观拉破裂有关,m_i值变化可用于描述饱和硬岩的软化程度。基于Griffith准则并采用起裂应力估算的岩石直接拉伸强度和H-B强度参数m_i可用于岩体工程的早期设计。     

Modeled Sverdrup flow in the North Atlantic from 11 different wind stress climatologies

In studies of large-scale ocean dynamics, often quoted values of Sverdrup transport are computed using the Hellerman–Rosenstein wind stress climatology. The Sverdrup solution varies, however, depending on the wind set used. We examine the differences in the large-scale upper ocean response to different surface momentum forcing fields for the North Atlantic Ocean by comparing the different Sverdrup interior/Munk western boundary layer solutions produced by a 1/16° linear numerical ocean model forced by 11 different wind stress climatologies. Significant differences in the results underscore the importance of careful selection of a wind set for Sverdrup transport calculation and for driving nonlinear models. This high-resolution modeling approach to solving the linear wind-driven ocean circulation problem is a convenient way to discern details of the Sverdrup flow and Munk western boundary layers in areas of complicated geometry such as the Caribbean and Bahamas. In addition, the linear solutions from a large number of wind sets provide a well-understood baseline oceanic response to wind stress forcing and thus, (1) insight into the dynamics of observed circulation features, by themselves and in conjunction with nonlinear models, and (2) insight into nonlinear model sensitivity to the choice of wind-forcing product.The wind stress products are evaluated and insight into the linear dynamics of specific ocean features is obtained by examining wind stress curl patterns in relation to the corresponding high-resolution linear solutions in conjunction with observational knowledge of the ocean circulation. In the Sverdrup/Munk solutions, the Gulf Stream pathway consists of two branches. One separates from the coast at the observed separation point, but penetrates due east in an unrealistic manner. The other, which overshoots the separation point at Cape Hatteras and continues to flow northward along the continental boundary, is required to balance the Sverdrup interior transport. A similar depiction of the Gulf Stream is commonly seen in the mean flow of nonlinear, eddy-resolving basin-scale models of the North Atlantic Ocean. An O(1) change from linear dynamics is required for realistic simulation of the Gulf Stream pathway. Nine of the eleven Sverdrup solutions have a C-shaped subtropical gyre, similar to what is seen in dynamic height contours derived from observations. Three mechanisms are identified that can contribute to this pattern in the Sverdrup transport contours. Along 27°N, several wind sets drive realistic total western boundary current transport (within 10% of observed) when a 14 Sv global thermohaline contribution is added (COADS, ECMWF 10 m re-analysis and operational, Hellerman–Rosenstein and National Centers for Environmental Prediction (NCEP) surface stress re-analysis), a few drive transport that is substantially too high (ECMWF 1000 mb re-analysis and operational and Isemer–Hasse) and Fleet Numerical Meteorology and Oceanography Center (FNMOC) surface stresses give linear transport that is slightly weaker than observed. However, higher order dynamics are required to explain the partitioning of this transport between the Florida Straits and just east of the Bahamas (minimal in the linear solutions vs. 5 Sv observed east of the Bahamas). Part of the Azores Current transport is explained by Sverdrup dynamics. So are the basic path of the North Atlantic Current (NAC) and the circulation features within the Intra-Americas Sea (IAS), when a linear rendition of the northward upper ocean return flow of the global thermohaline circulation is added in the form of a Munk western boundary layer.     

Neotectonics of the Kharaulakh sector of the Laptev Shelf

Seismotectonic deformation and crustal stress pattern have been studied comprehensively in major seismogenic structures of the Kharaulakh sector of the Verkhoyansk fold system and adjacent parts of the Chersky seismotectonic zone. The study focuses on neotectonic structures, deep structure, and systems of active faults, as well as tectonic stress fields inferred by tectonophysical analysis of Late Cenozoic faults and folds. The results, along with geological and geophysical data, reveal main strain directions and structural patterns of crustal stress and strain in the Arctic segment of the Eurasia-North America plate boundary. The area is a junction of mid-ocean and continental structures evolving in a mixed setting of extension, compression, and their various combinations. The rotation pole of the two plates is presumably located near Buor-Khaya Bay. In this case, extension is expected to act currently upon the neotectonic structures north of the bay and compression to control those in the south and southeast. This inference is consistent with the identified zoning of stress and strain in the Kharaulakh sector.     

双壳贝类CK1基因家族的鉴定及表达分析

酪蛋白激酶1 (CK1)是一种重要的蛋白激酶家族,其在DNA损伤应答和修复、细胞增殖和凋亡、胚胎发育和稳态等重要的生物学过程中都具有复杂多样的调控作用。为了理解CK1基因家族在双壳贝类中的特征、进化及生物学功能,实验采用比较基因组学及生物信息学的方法对双壳贝类CK1基因家族进行了鉴定分析,通过虾夷扇贝高温应激实验,研究了CK1基因家族在高温应激时的表达规律。结果显示,在虾夷扇贝、栉孔扇贝、长牡蛎与侏儒蛤中均存在CK1α,CK1αlike,CK1δ,CK1γ3,并发现CK1ε,CK1γ1,CK1γ2在双壳贝类中发生了丢失。时空表达分析发现,双壳贝类CK1基因均在胚胎发育早期集中表达,并以多细胞/囊胚期为界呈现出母源/合子特异性表达模式。CK1为关键基因的基因共表达模块显著富集在错配修复、核苷酸剪切修复等相关通路上,暗示了CK1基因在双壳贝类胚胎发育过程中参与调控DNA损伤修复过程,从而维持早期胚胎发育时的基因组稳定性。双壳贝类CK1基因在成体中呈现组织特异的表达模式,主要在鳃和雄性性腺中具有相对较高的表达量。虾夷扇贝受到高温应激后,其鳃中PyCK1α, PyCK1αlike和PyCK1δ...     

Oxygen radical formation in anoxic transgression and anoxia-reoxygenation: Foe or phantom? Experiments with a hypoxia tolerant bivalve

Intertidal blue mussels, Mytilus edulis, experience hypoxia reoxygenation during tidal emersion and resubmersion cycles, and this is often suggested to represent a major stress for the animals, especially for their respiratory tissues, the gills. We exposed mussels to experimental short and prolonged anoxia and subsequent reoxygenation and analyzed the respiratory response in excised gill tissue and the effects of treatment on reactive oxygen species (mainly ROS: superoxide anion, O₂^·− and hydrogen peroxide, H₂O₂), formation using live imaging techniques and confocal microscopy. Our aim was to understand if this “natural stress” would indeed produce oxidative damage and whether antioxidant defenses are induced under anoxia, to prevent oxidative damage during reoxygenation. Exposure to declining pO₂ in the respiration chamber caused an increase of gill metabolic rate between 21 and 10 kPa, a pO₂ range in which whole animal respiration is reported to be oxyregulating. Exposure of the animals to severe anoxia caused an onset of anaerobiosis (succinate accumulation) and shifted high and low critical p_c values (p_c1: onset of oxyregulation in gills, p_c2: switch from oxyregulation to oxyconformity) to higher pO₂. Concentrations of both ROS decreased strongly during anoxic exposure of the mussels and increased upon reoxygenation. This ROS burst induced lipid peroxidation in the mantle, but neither were protein carbonyl levels increased (oxidative damage in the protein fraction), nor did the tissue glutathione concentration change in the gills. Further, analysis of apoptosis markers indicated no induction of cell death in the gills. To our knowledge, this is the first paper that directly measures ROS formation during anoxia reoxygenation in mussels. We conclude that hypoxia tolerant intertidal mussels do not suffer major oxidative stress in gill and mantle tissues under these experimental conditions.     

Numerical Modeling of Sediment Mass Sliding and Tsunami Generation: The Case of February 7, 1963, in Corinth Gulf,Greece

Numerical modeling of a landslide process at the continental slope by taking into account slope parameters is performed for the tsunami event of February 7, 1963, reported in the Corinth Gulf, central Greece. A layered sediment structure was considered, and an initial external dynamic action at the landslide process was introduced. The results obtained were quite consistent with the observational runup data as well as with the results arrived at by conventional rigid-body and viscous-fluid models. However, this approach permits to describe in detail the formation of several wave groups and particular tsunami characteristics that are strongly dependent on the landslide model.     

Risk Analysis of Breakwater Caisson Under Wave Attack Using Load Surface Approximation

A new load surface based approach to the reliability analysis of caisson-type breakwater is proposed. Uncertainties of the horizontal and vertical wave loads acting on breakwater are considered by using the so-called load surfaces, which can be estimated as functions of wave height, water level, and so on. Then, the first-order reliability method（FORM） can be applied to determine the probability of failure under the wave action. In this way, the reliability analysis of breakwaters with uncertainties both in wave height and in water level is possible. Moreover, the uncertainty in wave breaking can be taken into account by considering a random variable for wave height ratio which relates the significant wave height to the maximum wave height. The proposed approach is applied numerically to the reliability analysis of caisson breakwater under wave attack that may undergo partial or full wave breaking.