An improved return‐map stress update algorithm for finite deformation analysis of general isotropic elastoplastic geomaterials

This paper develops a novel return mapping algorithm for the numerical integration of general isotropic finite strain elastoplastic constitutive models for geomaterials. The constitutive formulation is founded on multiplicative decomposition of the deformation gradient. The logarithmic strain measure as well as the exponential approximation of the plastic flow rule is utilized to restore the standard infinitesimal format return mapping algorithm. Central to the algorithm is the exploitation of a set of three mutually orthogonal unit base tensors for the representation of constitutive relations and the corresponding integration of the rate form of the constitutive equations. The base tensors constitute a local cylindrical coordinate system in the principal space, which allows to formulate the return mapping algorithm in the three‐dimensional space and reduce the dimension of the problem to be analyzed from six down to three. With the proposed approach, direct determination of the principal axes and the transformation procedure between the general space and the principal space, as required in traditional spectral decomposition, are avoided. Furthermore, the matrices that are involved in the inversion evaluation take simple forms, leading to extremely easy inverse computation. As a result, the consistent tangent operator can be streamlined into a form simpler and more compact than those by conventional integration methods. Following the formulation of the integration procedure, a numerical experiment is performed to assess the accuracy and efficiency of the proposed algorithm. Copyright © 2014 John Wiley & Sons, Ltd.     

Capturing the complete stress–strain behaviour of jointed rock using a numerical approach

This paper presents the results of a series of numerical experiments using the synthetic rock mass (SRM) approach to quantify the behaviour of jointed rock masses. Field data from a massive sulphide rock mass, at the Brunswick mine, were used to develop a discrete fracture network (DFN). The constructed DFN model was subsequently subjected to random sampling whereby 40 cubic samples, of height to width ratio of two, and of varying widths (0.05 to 10 m) were isolated. The discrete fracture samples were linked to 3D bonded particle models to generate representative SRM models for each sample size. This approach simulated the jointed rock mass as an assembly of fractures embedded into the rock matrix. The SRM samples were submitted to uniaxial loading, and the complete stress–strain behaviour of each specimen was recorded. This approach provided a way to determine the complex constitutive behaviour of large‐scale rock mass samples. This is often difficult or not possible to achieve in the laboratory. The numerical experiments suggested that higher post‐peak modulus values were obtained for smaller samples and lower values for larger sample sizes. Furthermore, the observed deviation of the recorded post‐peak modulus values decreased with sample size. The ratio of residual strength of rock mass samples per uniaxial compressive strength intact increases moderately with sample size. Consequently, for the investigated massive sulphide rock mass, the pre‐peak and post‐peak representative elemental volume size was found to be the same (7 × 7 × 14 m). Copyright © 2015 John Wiley & Sons, Ltd.     

Using MODFLOW code to approach transient hydraulic head with a sharp‐interface solution

Saltwater intrusion problems have been usually tackled through analytical models because of its simplicity, easy implementation and low computational cost. Most of these models are based on the sharp‐interface approximation and the Ghyben–Herzberg relation, which neglects mixing of fresh water and seawater and implicitly assumes that salt water remains static. This paper provides insight into the validity of a sharp‐interface approximation defined from a steady state solution when applied to transient seawater intrusion problems. The validation tests have been performed on a 3D unconfined synthetic aquifer, which include spatial and temporal distribution of recharge and pumping wells. Using a change of variable, the governing equation of the steady state sharp‐interface problem can be written with the same structure of the steady confined groundwater flow equation as a function of a single potential variable (?). We propose to approach also the transient problem solving a single potential equation (using also the ? variable) with the same structure of the confined groundwater flow equation. It will allow solving the problem by using the classical MODFLOW code. We have used the parameter estimation model PEST to calibrate the parameters of the transient sharp‐interface equation. We show how after the calibration process, the sharp‐interface approach may provide accurate enough results when applied to transient problems and improve the steady state results, thus avoiding the need of implementing a density‐dependent model and reducing the computational cost. This has been proved by comparing results with those obtained using the finite difference numerical code SEAWAT for solving the coupled partial differential equations of flow and density‐dependent transport. The comparison was performed in terms of piezometric heads, seawater penetration, transition zone width and critical pumping rates. Copyright © 2014 John Wiley & Sons, Ltd.     

基于水下光纤应变监测的导管架结构总冰力测量方法

基于测量结构水下光纤应变获取冰力的方法首次应用于渤海JZ20-2NW加锥导管架平台,该方法方便进行零点标定,可获取结构总冰力的绝对信息(包括均值和波动值)。首先介绍了JZ20-2NW平台的现场监测系统,主要包括结构水下应变响应记录冰力信息,甲板上部视频记录同步冰情(冰厚、冰速和来冰方向)信息和拾振器记录结构冰振响应;重点分析了由测点应变向结构总冰力的转化方法,并对总冰力进行了初步分析。将无量纲化实测冰力与5种典型冰力计算模型进行了比较,分析结果表明,锥体宽度与海冰厚度比值(简称"宽厚比")对冰力大小及变化起决定作用。     

�����Ͷ�Ms5.1����ǰ�˲�̨����Ӧ��۲����ϱ仯ԭ�����

???ó????????????С???任????????Ms5.1????????????????????????仯??????з???????????????????2????????????????????????????????????????????????????к??????????????????????????仯???????????????????????????????С??????????????????仯???????????????????????????仯?????????????????????????????仯??     

��������ؿ�Ӧ�����ܶȱ仯�ʼ��乹���˶�����

?????????????????仯????????????е????????????????й????????????????????????????GPS????????????й??????????????????????????????????仯??????????????2001~2007???????????????????????к????????????????????????????????????????????????????????????????????????????2007~2009?????????????????????????????????????????????????????ν??????????????????????????????????2009~2011????????????????????????????????????????????????2007~2009?????     

菲律宾蛤仔不同品系斑马蛤的选择反应及现实遗传力

对菲律宾蛤仔中的斑马蛤、黑斑马、白斑马三品系进行混合选择,测量了各实验组的壳长,计算了不同品系斑马蛤的选择反应和现实遗传力。结果表明:3个品系的上选组壳长均大于对照组,说明单向上选是行之有效的遗传改良方法。在不同的生长发育阶段,斑马蛤品系的选择反应(R)和现实遗传力(2Rh)出现阶段性变化,从平均水平上看,R在幼虫期为0.688±0.184,稚贝期为0.605±0.214,养成期为0.670±0.120;2Rh在幼虫期为0.392±0.081,稚贝期为0.344±088,养成期为0.380±0.058。从一周年的生长情况来看,三个品系总的R为0.654±0.179,斑马蛤、黑斑马蛤和白斑马蛤品系的R分别为0.655±0.155、0.525±0.116、0.782±0.170;三个品系总的2Rh为0.373±0.070,斑马蛤、黑斑马蛤和白斑马蛤品系的2Rh分别为0.372±0.042、0.300±0.028、0.445±0.031。不同品系间R和2Rh次序为白斑马斑马蛤黑斑马,说明各斑马蛤品系的选择效果是不同的,品系间存在着一定程度的遗传差异。     

Time preference of small-scale fishers in open access and traditionally managed reef fisheries

Individuals with high discount rates are likely not partial to conservation because they are unwilling to sacrifice short term benefits for potentially higher gains in the future. Many reef fisheries worldwide are open access, and fishers under open access systems are theorized to discount the future at an infinite rate. In contrast, fishers in a customary managed fishery can be expected to be more long term oriented, and thus possibly have lower discount rates. The present study tests this hypothesis by eliciting the discount rates of fishers in an open access small-scale reef fishery, and compares these rates to those of fishers in a customary managed reef fishery. Results indicate that fishers in both open access and traditionally managed reef fisheries have high annual discount rates that are on average over 200%. Contrary to expectations, fishers under an open access system are not associated with higher discount rates compared to customary management. It also appears that a larger proportion of open access fishers are more long-term oriented than those in the customary managed fishery, which is encouraging for the future conservation and sustainability of open access fisheries resources.     

Post‐1850 changes in Glacier Benito,North Patagonian Icefield,Chile

Glacier Benito is a temperate outlet glacier on the west side of the North Patagonian Icefield. Rates of thinning and ablation were obtained using data collected by the British Joint Services Expedition in 1972/73 and subsequent data collected in 2007 and 2011. Ice‐front recession rates were based on dendrochronological dating for the terminal moraines and aerial and satellite imagery of the ice front in 1944, 1998 and 2002. Between the first Benito survey in 1973 and 2007, the lower glacier thinned by nearly 150 m at an average rate of 4.3 m yr^?1 with the rate increasing to 6.1 m yr^?1 between 2007 and 2011, a 28.7% increase during the latter period. Increases in ice movement and ablation were negligible: ice movement for 1973 and 2007 averaged 0.45 m day^?1 and ablation averaged 0.05 m day^?1. Ice front recession along the glacier's centre line from 1886 to 2002 was approximately 1860 m. Retreat rates between 1886 and 1944 averaged 8.9 m yr^?1. Thereafter glacier asymmetry makes measurement along the glacier centre line unrepresentative of areal change until 1998 when symmetry was restored; retreat between 1944 and 1998 was 15.4 m yr^?1. From 1998 to 2002 the rate increased dramatically to 127.2 m yr^?1. Recession from the southern end of Benito's terminal moraine in the 1850s supports an early date for initial retreat of the Icefield's glaciers.     

玉树MS7.1级地震前的形变应变场变化特征

基于1999～2007年和2007～2009年的GPS站点观测数据计算的速度场,分析了青藏亚板块的运动特征.利用块体(旋转与线性应变)弹性运动公式计算可可西里-巴颜喀拉块体和羌唐块体的应变率、可可西里-玉树断裂带的运动速度和应变率,研究了玉树7.1级地震前的形变应变特征.结果表明:可可西里-巴颜喀拉块体处于青藏亚板块运...