Improvement of FEM's dynamic property

The discretization size is limited by the sampling theorem, and the limit is one half of the wavelength of the highest frequency of the problem. However, one half of the wavelength is an ideal value. In general, the discretization size that can ensure the accuracy of the simulation is much smaller than this value in the traditional finite element method. The possible reason of this phenomenon is analyzed in this paper, and an efficient method is given to improve the simulation accuracy.

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.

Improving Nucleation in the Fabrication of High-Quality 3D Macro-Porous Copper Film through The Surface-Modification of A Polystyrene Colloid-Assembled Template

A new approach is developed to the fabrication of high-quality three-dimensional macro-porous copper films. A highly-ordered macroporous copper film is successfully produced on a polystyrene sphere (PS) template that has been modified by sodium dodecyl sulfate (SDS). It is shown that this procedure can change a hydrophobic surface of PS template into a hydrophilic surface. The present study is devoted to the influence of the electrolyte solution transport on the nucleation process. It is demonstrated that the permeability of the electrolyte solution in the nanochannels of the PS template plays an important role in the chemical electrodeposition of high-quality macroporous copper film. The permeability is drastically enhanced in our experiment through the surface modi. cation of the PS templates. The method could be used to homogeneously produce a large number of nucleations on a substrate, which is a key factor for the fabrication of the high-quality macroporous copper film.

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.

Micro- and macro-tribological study on a self-assembled dual-layer film

A novel self-assembled dual-layer film as apotential excellent lubricant for micromachines was successfully prepared on single-crystal silicon substrate by chemical adsorption of stearic acid (STA) molecules on self-assembled monolayer of 3-aminopropyltri