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.

Effect of dentin tubules on the mechanical properties of dentin - part III: Numerical, analysis

Imitating a real tooth and the periodontal supporting tissues, we have established a 2D finite element model and carried out a numerical analysis based on the inhomogeneous and anisotropic (IA) stress-strain relation and strength model of dentin proposed in the preceding Parts I and II, and the conventional homogeneous and isotropic (III) model, respectively. Quite a few cases of loadings for a non-defected and a defected tooth are considered. The numerical results show that the stress level predicted by the IA model is remarkably higher than that by the III model, revealing that the effect of the dentin tubules should be taken into a serious consideration from the viewpoint of biomechanics.

Analysis of critical excavation depth for a jointed rock slope using a face-to-face discrete element method

The critical excavation depth of a jointed rock slope is an important problem in rock engineering. This paper studies the critical excavation depth for two idealized jointed rock slopes by employing a face-to-face discrete element method (DEM). The DEM is based on the discontinuity analysis which can consider anisotropic and discontinuous deformations due to joints and their orientations. It uses four lump-points at each surface of rock blocks to describe their interactions. The relationship between the critical excavation depth D-s and the natural slope angle alpha, the joint inclination angle theta as well as the strength parameters of the joints c(r) ,phi(r) is analyzed, and the critical excavation depth obtained with this DEM and the limit equilibrium method (LEM) is compared. Furthermore, effects of joints on the failure modes are compared between DEM simulations and experimental observations. It is found that the DEM predicts a lower critical excavation depth than the LEM if the joint structures in the rock mass are not ignored.

A Finite Element Analysis of a Crack Penetrating or Deflecting into an Interface in a Thin Laminate

The crack tip driving force of a crack growing from a pre-crack that is perpendicular to and terminating at an interface between two materials is investigated using a linear fracture mechanics theory. The analysis is performed both for a crack penetrating the interface, growing straight ahead, and for a crack deflecting into the interface. The results from finite element calculations are compared with asymptotic solutions for infinitesimally small crack extensions. The solution is found to be accurate even for fairly large amounts of crack growth. Further, by comparing the crack tip driving force of the deflected crack with that of the penetrating crack, it is shown how to control the path of the crack by choosing the adhesion of the interface relative to the material toughness.

Influence of Inter-Fiber Spacing and Interfacial Adhesion on Failure of Multi-Fiber Model Composites: Experiment and Numerical Analysis

The uniaxial tension experiments on glass-fiber-reinforced epoxy matrix composites reveal that the fragmentations of fibers display vertically aligned fracture, clustered fracture, coordinated fracture, and random fracture with the increase of inter-fiber spacing. The finite element analysis indicates that the fragmentations of fibers displaying different phenomena are due to the stress concentration as well as the inherent randomness of fiber defects, which is the dominant factor. The experimental results show that matrices adjacent to the fiber breakpoints all exhibit birefringent-whitening patterns for the composites with different interfacial adhesion strengths. The larger the extent of the interfacial debonding, the less the domain of the birefringent-whitening patterns. The numerical analysis indicates that the orientation of the matrix adjacent to a fiber breakpoint is caused by the interfacial shear stress, resulting in the birefringent-whitening patterns. The area of shear stress concentrations decides on the domain of the birefringent-whitening patterns.

Numerical Analysis of Theoretical Model of the Rf MEMs Switches

An improved electromechanical model of the RF MEMS (radio frequency microelectromechanical systems) switches is introduced, in which the effects of intrinsic residual stress from fabrication processes, axial stress due to stretching of beam, and fringing field are taken into account. Four dimensionless numbers are derived from the governing equation of the developed model. A semi-analytical method is developed to calculate the behavior of the RF MEMS switches. Subsequently the influence of the material and geometry parameters on the behavior of the structure is analyzed and compared, and the corresponding analysis with the dimensionless numbers is conducted too. The quantitative relationship between the presented parameters and the critical pull-in voltage is obtained, and the relative importance of those parameters is given.

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

Triple-region structure for turbulent-flow in a square duct - a finite-element approach

The Reynolds-averaged Navier-Stokes equations for describing the turbulent flow in a straight square duct are formulated with two different turbulence models. The governing equations are then expanded as a multi-deck structure in a plane perpendicular to the streamwise direction, with each deck characterized by its dominant physical forces as commonly carried out in analytical work using triple-deck expansion. The resulting equations are numerically integrated using higher polynomial (H-P) finite element technique for each cross-sectional plane to be followed by finite difference representation in the streamwise direction until a fully developed state is reached. The computed results using the two different turbulence models show fair agreement with each other, and concur with the vast body of available experimental data. There is also general agreement between our results and the recent numerical works anisotropic k-epsilon turbulence model.

A surface energy model and application to mechanical behavior analysis of single crystals at sub-micron scale

Size effects of mechanical behaviors of materials are referred to the variation of the mechanical behavior due to the sample sizes changing from macroscale to micro-/nanoscales. At the micro-/nanoscale, since sample has a relatively high specific surface area (SSA) (ratio of surface area to volume), the surface although it is often neglected at the macroscale, becomes prominent in governing the energy effect, although it is often neglected at the macroscale, becomes prominent in governing the mechanical behavior. In the present research, a continuum model considering the surface energy effect is developed through introducing the surface energy to total potential energy. Simultaneously, a corresponding finite element method is developed. The model is used to analyze the axial equilibrium strain problem for a Cu nanowire at the external loading-free state. As another application of the model, from dimensional analysis, the size effects of uniform compression tests on the microscale cylinder specimens for Ni and Au single crystals are analyzed and compared with experiments in literatures. (C) 2009 Elsevier B.V. All rights reserved.

Numerical analysis and comparison of three metal-oxide-type super-resolution near field structures

Based on the Fresnel-Kirchkoff diffraction theory, we build up a Gaussian diffraction model of metal-oxide-type super-resolution near field structure (super-RENS), which can describe far field optical properties. The spectral contrast induced by refractive index and the structural changes in AgOx, PtOx and PdOx thin films, which are the key functional layers in super-RENS, are studied by using this model. Comparison results indicate that the spectral contrast intensively on laser-induced distribution and change of the refractive index in the metal-oxide films. The readout mechanism of the metal-oxide-type super-RENS optical disc is further clarified. This Gaussian diffraction model can be used as a simple and effective method for choosing proper active materials in super-RENS.

Numerical analysis for four-wave mixing based wavelength conversion of differential phase-shift keying signals

We numerically investigate the main constrains for high efficiency wavelength conversion of differential phase-shift keying (DPSK) signals based on four-wave mixing (FWM) in highly nonlinear fiber (HNLF). Using multi-tone pump phase modulation techniques, high efficiency wavelength conversion of DPSK signals is achieved with the stimulated Brillouin scattering (SBS) effects effectively suppressed. Our analysis shows that there is a compromise between conversion efficiency and converted idler degradation. By optimizing the pump phase modulation configuration, the converted DPSK idler's degradation can be dramatically decreased through balancing SBS suppression and pump phase modulation degradation. Our simulation results also show that these multi-tone pump phase modulation techniques are more appropriate for the future high bit rate systems.

Fabrication and analysis of 2x2 thermo-optic SOI waveguide switch with low power consumption and fast response by anisotropy chemical etching

A low power consumption 2 x 2 thermo-optic switch with fast response was fabricated on silicon-on-insulator by anisotropy chemical etching. Blocking trenches were etched on both sides of the phase-shifting arms to shorten device length and reduce power consumption. Thin top cladding layer was grown to reduce power consumption and switching time. The device showed good characteristics, including a low switching power of 145 mW and a fast switching speed of 8 +/- 1 mus, respectively. Two-dimensional finite element method was applied to simulate temperature field in the phase-shifting arm instead of conventional one-dimensional method. According to the simulated result, a new two-dimensional index distribution of phase-shifting arm was determined. Consequently finite-difference beam propagation method was employed to simulate the light propagation in the switch, and calculate the power consumption as well as the switching speed. The experimental results were in good agreement with the theoretical estimations. (C) 2004 Elsevier B.V. All rights reserved.

Thermal stress analysis for GaInAsP multiple quantum well wafer chemically bonded to Si (100)

An n-InP-based InGaAsP multiple-quantum-well wafer was bonded with p-Si by chemical surface activated bonding at 70 degrees C, and then annealed at 450 degrees C. Different thermal expansion coefficients between InP and Si will induce thermal stresses in the bonded wafer. Planar and cross-sectional distributions of thermal stress in the bonded InP-Si pairs were analyzed by a two-dimensional finite element method. In addition, the normal, peeling, and shear stresses were calculated by an analytic method. Furthermore, x-ray double crystalline diffraction was applied to measure the thermal strain and the strain caused by the mismatching of the crystalline orientation between InP (100) and Si (100). The wavelength redshift of the photoluminescence (PL) spectrum due to thermal strain was investigated via the calculation of the band structure, which is in agreement with the measured PL spectra.

Effects and numerical analysis of argon gas flow on the oxygen concentration in Czochralski silicon single crystal growth

Argon gas, as a protective environment and carrier of latent heat, has an important effect on the temperature distribution in crystals and melts. Numeric simulation is a potent tool for solving engineering problems. In this paper, the relationship between argon gas flow and oxygen concentration in silicon crystals was studied systematically. A flowing stream of argon gas is described by numeric simulation for the first time. Therefore, the results of experiments can be explained, and the optimum argon flow with the lowest oxygen concentration can be achieved. (C) 2002 Elsevier Science B.V. All rights reserved.