神农架地区米心水青冈林和锐齿槲栎林的格局与过程

本文对神农架地区广泛分布的米心水青冈林和锐齿槲栎林的种群和群落学特征、干扰历史、更新策略、生物量、生产量及元素循环特征进行了研究。得到如下结论： 1 米心水青冈林是神农架地区山地垂直分布的地带性植被类型，主要群落学特征为：（1）建群种明显，该区域主要有2种类型，即米心水青冈林和米心水青冈、锐齿槲栎林;（2）群落结构简单，但物种组成丰富，在6600m~2样地中出现高等植物（不含苔藓植物）77科150属271种，组成种类以蔷薇科、百合科、忍冬科、虎耳草科、樟科、杜鹃花科和壳斗科为主;（3）群落乔木层（占重要值的12％）和灌木层（占盖度的15％）中含有一定比例的常绿树种;（4）群落生活型以高位芽植物（70.89%）占绝对优势，其次为地面芽植物（15.50%）和地下芽植物(12.92%)。 2 米心水青冈是多主干的树种，萌枝现象普遍，但萌枝数量不同地点差异较大。通过萌枝产生的枝群体平均密度为257 ± 99.3n•hm~(-2)。枝群体的年龄结构表现为“幼龄个体数目较多型”和“中国年龄阶段数目较多型”，并且有较多的枝群体表现出一致的年龄结构。从整个群落米心水青冈的年龄结构来看，表现出发展型种群的特点。枝群体的分布格局为随机分布。9丛米心水青冈完整的年轮分析结果表明，它们萌枝的时间不是边疆的，而与森林的受干扰有关。根据83个圆盘和生长锥芯资料，米心水青冈在萌枝后成长为乔木层或林冠层的过程中，径向生长表现为5种模式。这是丛株内竞争的结果。萌枝在米心水青冈林的维持和发展过程中，具有重要的生态学作用。 3 锐齿槲栎林是神农架地区山地垂地分布的地带性植被类型，主要群落学特征为：（1）建群种明显，该区域主要有2种类型，即锐齿槲栎林和锐齿槲栎、米心水青冈林;（2）群落乔木层和灌木层中含有一定比例的常绿树种，和暖温带的落叶栎林有较大差异;（3）群落物种组成丰富，不仅具有典型的温带科属，还有典型亚热带分布的科属，组成种类主要以蔷薇科、百合科、忍冬科、虎耳草科、山茱萸科、杜鹃花科、壳斗科和樟科;（4）生活型以高位芽植物（66.32%）占绝对优势，其次为地面芽（23.51%）和地下芽（9.47%）植物。 4 通过样地调查、树干解析及直径分析法，对米心水青冈林和锐齿槲栎林受压和释压历史及更新策略进行了研究。米心水青冈直径生长表现为5种模式。而锐齿槲栎只表现为2种模式。85.9 ± 6.9%的米心水青冈有过受压过程，平均受压2.1 ± 0.8次，平均受压时间为47 ± 24.1a，最长受压时间73a，平均释压次数为1.6 ± 0.7次，平均释压时间为23 ± 21.5a，而60.83%的锐齿槲栎都均有1次受压。平均受压时间为19 ± 14a，受压后没有表现出释压过程。结合高生长和径向生长，认为米心水青冈是耐阴树种，它的更新策略是在林下形成苗性萌枝，在有林窗形成时释压生长进入乔木层;而锐齿槲栎是不耐阴树种，它的更新策略是通过产生大量种子，当有大的林窗时，幼苗在林窗内生长逐步进入乔木层。 5 神农架地区102-130a成熟米心水青冈林的生物量在251.31-358.63T•hm~(-2)之间，平均为288.70 ± 48.30T•~(-2)，20-60a锐齿槲栎林群落生物量在134.85-301.20T•hm~(-2)，平均为231.60 ± 78.10T•hm~(-2)。虽然米心水青冈林和锐齿槲栎林灌木层草本层及藤本植物组成很丰富，但二种类型森林生物量的95％以上集中在乔木层。乔木层生物量主要集中在少数优势种中。在米心水青冈林生物量从大到小的序列中，前5种植物分别占乔木层总生物量58.67%-96.37%不等，同样锐齿槲栎林前5种植物占群落生物量的68.13%-95.26%。常绿植物占乔木层生物量的比例变化较大，米心水青冈林中占2.85-18.70%，锐齿槲栎林中一般常绿植物占0.8-9.98%，只有1个锐齿槲栎林样地常绿植物（主要是粉白杜鹃）占乔木层生物量的44.04%。米心水青冈林生物量根冠比为0.27 ± 0.05，锐齿槲栎林为0.21 ± 0.06。神农架地区米心水青冈林的生物量，在成熟的欧洲水青冈林及日本的水青冈林生物量范围之内，而锐齿槲栎林生物量远远大于我国温带落叶栎林的生物量。 6 神农架地区102-130a米心水青冈林生产量范围在1857-2786g•m~(-2)•a~(-1)，平均为2330 ± 397 g•m~(-2)•a~(-1)。20-60a锐齿槲栎林的生产量范围在1319-2521 g•m~(-2)•a~(-1)，平均为1930 ± 498 g•m~(-2)•a~(-1)。米心水青冈林和锐齿槲栎林乔木层生产量占群落总生产量的95％以上，乔木层各器官生产量大小顺序为叶> 树干> 枝和根，其中叶生产量占乔木层的一半以上，达53.87 ± 2.72%（米心水青冈林）和57.31 ± 6.23%（锐齿槲栎林）。在乔木层生产量从大到小的序列中，前5种植物平均占乔木层总生产量的81.03 ± 13.94%（米心水青冈林，范围在62.75%-92.66%）和84.23 ± 9.68%（锐齿槲栎林，范围在68.54-95.11%）。米心水青冈林和锐齿槲栎林群落地下部分生产量占总生产量的比例分别为11.29 ± 1.02%和9.22 ± 2.72。和我国其它地区地带性植被类型相比。米心水青冈林和锐齿槲栎生产量是较高的，和亚热带绿阔叶林生产量接近，但在器官分配上两者差异较大。 7 米心水青冈林和锐齿槲栎林土壤均呈酸性。其中锐齿槲栎林地土壤酸性更强。土壤元素特征表现为Al>C>K>Mg>Ca>N>S、P的特点，富铝化作用明显。8种元素在群落优势植物不同部位含量差异较大，N、P、K、Ca、Mg基本上是以叶片含量最高，树干或根中最低。仅从叶片来看，元素特征表现为C>Ca、N>K>Mg>S>P、Al。优势植物的C/N和C/P显著高于暖温带落叶阔叶林优势植物。8种元素在米心水青冈林和锐齿槲栎林中积累量分别为147.09 ± 25.60和116.00 ± 37.63 Mg hm~(-2)a~(-1)，其中97％以上积累在群落乔木层。两种森林类型各元素的积累量都表现为C>Ca, N> K> Mg> P> S> Al的特点。米心水青风林和锐齿槲栎林8种元素的年存留量分别为6263 ± 90.8和5946 ±246 kg hm~(-2)a~(-1),其中N、P、K、Ca、Mg 5种主要营养元素的存留量分别为179.7 ± 18.2和169.4 ± 23.5kg hm~(-2) a~(-1)。两种森林类型各元素的存留量都表现为C> N> Ca> K> Mg> S> P> Al。

Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

三峡永久船闸高边坡开挖三维离散元数值模拟:3-D Simulation of the Excavation of High Steep Slope of Three-Gorges Permanent Lock by Distinct Element Method

采用面一面接触的三维离散元刚性块体模型,从实测节理面中取出其中的三组,按照其倾向、倾角和节理间距将三峡永久船闸未开挖的区域划分为10~5个离散单元,通过施加力边界条件,给出了与实测初始地应力场接近的数值模拟结果；然后,分4步模拟了永久船闸的开挖过程。计算结果表明：开挖过程会引起节理面出现张开趋势,个别岩体还会沿着节理面滑移。岩体位移的不对称现象较为自然地说明了由节理引起的岩体各向异性特征。

Determination of Proper Processing Conditions of Material Surface Heat Treatment Through Finite Element Method

A two-dimensional model has been developed based on the experimental results of stainless steel remelting with the laminar plasma technology to investigate the transient thermo-physical characteristics of the melt pool liquids. The influence of the temperature field, temperature gradient, solidification rate and cooling rate on the processing conditions has been investigated numerically. Not only have the appropriate processing conditions been determined according to the calculations, but also they have been predicted with a criterion established based on the concept of equivalent temperature area density (ETAD) that is actually a function of the processing parameters and material properties. The comparison between the resulting conditions shows that the ETAD method can better predict the optimum condition.

A New Finite Element Method for Strain Gradient Theories and Applications to Fracture Analyses

A new compatible finite element method for strain gradient theories is presented. In the new finite element method, pure displacement derivatives are taken as the fundamental variables. The new numerical method is successfully used to analyze the simple strain gradient problems – the fundamental fracture problems. Through comparing the numerical solutions with the existed exact solutions, the effectiveness of the new finite element method is tested and confirmed. Additionally, an application of the Zienkiewicz–Taylor C1 finite element method to the strain gradient problem is discussed. By using the new finite element method, plane-strain mode I and mode II crack tip fields are calculated based on a constitutive law which is a simple generalization of the conventional J2 deformation plasticity theory to include strain gradient effects. Three new constitutive parameters enter to characterize the scale over which strain gradient effects become important. During the analysis the general compressible version of Fleck–Hutchinson strain gradient plasticity is adopted. Crack tip solutions, the traction distributions along the plane ahead of the crack tip are calculated. The solutions display the considerable elevation of traction within the zone near the crack tip.

Distinct Element Analysis on Propagation Characteristics of P-Wave in Rock Pillar with Finite length

以节理岩体等效刚度的概念为基础,讨论了离散元刚性块体模型中节理刚度的选取问题。采用面-面接触模型模拟了纵波在一维岩体中的传播,给出了纵波波形；研究了阻尼比、软弱夹层以及节理间是否可拉对波传播规律的影响。

A finite element method for cracked components of structures

In this paper, a method is developed for determining the effective stiffness of the cracked component. The stiffness matrix of the cracked component is integrated into the global stiffness matrix of the finite element model of the global platform for the FE calculation of the structure in any environmental conditions. The stiffness matrix equation of the cracked component is derived by use of the finite variation principle and fracture mechanics. The equivalent parameters defining the element that simulates the cracked component are mathematically presented, and can be easily used for the FE calculation of large scale cracked structures together with any finite element program. The theories developed are validated by both lab tests and numerical calculations, and applied to the evaluation of crack effect on the strength of a fixed platform and a self-elevating drilling rig.

Iterative method using consistent mass matrix in axisymmetrical finite element analysis of hypervelocity impact

We present in this paper an iterative method using consistent mass matrix in axisymmetrical finite element analysis of hypervelocity impact. To retain the advantage of integration on an element-by-element basis which is at the heart of modern hydrocodes, we suggest that the first step should be to solve for accelerations at an advanced time step by using the lumped mass approach, then iterate using a consistent mass matrix to improve the estimate. Examples are given to show the improved resolution with the new method.

Finite Element Analysis of Deformed Legs of Offshore Platform Structures

The element stiffness matrix of the equivalent beam or pipe element of the deformed leg of the platform is derived by the finite element method. The stresses and displacements of some damaged components are calculated, and the numeri-cal solutions agree well with those obtained by the fine mesh finite element method. Finally, as an application of this method, the stresses of some platform structures are calculated and analyzed.

Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

The application of B-P constitutive equations in finite element analysis of high velocity impact

We present in this paper the application of B-P constitutive equations in finite element analysis of high velocity impact. The impact process carries out in so quick time that the heat-conducting can be neglected and meanwhile, the functions of temperature in equations need to be replaced by functions of plastic work. The material constants in the revised equations can be determined by comparison of the one-dimensional calculations with the experiments of Hopkinson bar. It can be seen from the comparison of the calculation with the experiment of a tungsten alloy projectile impacting a three-layer plate that the B-P constitutive equations in that the functions of temperature were replaced by the functions of plastic work can be used to analysis of high velocity impact.

Stochastic Structural Model of Rock and Soil Aggregates by Continuum-Based Discrete Element Method

This paper first presents a stochastic structural model to describe the random geometrical features of rock and soil aggregates. The stochastic structural model uses mixture ratio, rock size and rock shape to construct the microstructures of aggregates,and introduces two types of structural elements (block element and jointed element) and three types of material elements (rock element, soil element, and weaker jointed element)for this microstructure. Then, continuum-based discrete element method is used to study the deformation and failure mechanism of rock and soil aggregate through a series of loading tests. It is found that the stress-strain curve of rock and soil aggregates is nonlinear, and the failure is usually initialized from weaker jointed elements. Finally, some factors such as mixture ratio, rock size and rock shape are studied in detail. The numerical results are in good agreement with in situ test. Therefore, current model is effective for simulating the mechanical behaviors of rock and soil aggregates.

Application of Three-Dimensional Discrete Element Face-to-Face Contact Model with Fissure Water Pressure to Stability Analysis of Landslide in Panluo Iron Mine

Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical phi numerically calculated is less than the one calculated by use of the limit equilibrium method for the same C. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.