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.

Gain saturation in 60-fs mode-locked semiconductor laser

A passively mode-locked optically-pumped InGaAs/GaAs quantum well laser with an intracavity semiconductor saturable absorber mirror emits sub-100-fs pulses. Pulse energy declines steeply as pulse duration is reduced below 100 fs due to gain saturation. © 2010 Optical Society of America.

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.

Two-dimensional Cartesian air-gap element (CAGE) for dynamic finite-element modeling of electrical machines with a flat air gap

This paper extends the air-gap element (AGE) to enable the modeling of flat air gaps. AGE is a macroelement originally proposed by Abdel-Razek et al.for modeling annular air gaps in electrical machines. The paper presents the theory of the new macroelement and explains its implementation within a time-stepped finite-element (FE) code. It validates the solution produced by the new macroelement by comparing it with that obtained by using an FE mesh with a discretized air gap. It then applies the model to determine the open-circuit electromotive force of an axial-flux permanent-magnet machine and compares the results with measurements.

Low-temperature (≤600 °C) semi-insulating oxygen-doped silicon films by the PECVD technique for large-area power applications

This work describes the deposition, annealing and characterisation of semi-insulating oxygen-doped silicon films at temperatures compatible with polysilicon circuitry on glass. The semi-insulating layers are deposited by the plasma enhanced chemical vapour deposition technique from silane (SiH4), nitrous oxide (N2O) and helium (He) gas mixtures at a temperature of 350 °C. The as-deposited films are then furnace annealed at 600 °C which is the maximum process temperature. Raman analysis shows the as-deposited and annealed films to be completely amorphous. The most important deposition variable is the N2O SiH4 gas ratio. By varying the N2O SiH4 ratio the conductivity of the annealed films can be accurately controlled, for the first time, down to a minimum of ≈10-7Ω-1cm-1 where they exhibit a T -1 4 temperature dependence indicative of a hopping conduction mechanism. Helium dilution of the reactant gases is shown to improve both film uniformity and reproducibility. A model for the microstructure of these semi-insulating amorphous oxygen-doped silicon films is proposed to explain the observed physical and electrical properties. © 1995.

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.