A finite deformation analysis of polycrystal by considering the influence of grain boundary

By combining grain boundary (GB) and its influence zone, a micromechanic model for polycrystal is established for considering the influence of GB. By using the crystal plasticity theory and the finite element method for finite deformation, numerical simulation is carried out by the model. Calculated results display the microscopic characteristic of deformation fields of grains and are in qualitative agreement with experimental results.

Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

Numerical Simulation of Heat Transfer and Kinetics in the Bulk Growth of Silicon Carbide

The growth process of 2-inch silicon carbide (SiC) single crystals by the physical vapor transport method (or modified Lely method) has been modeled and simulated. The comprehensive process model incorporates the calculations of radio frequency (RF) induction heating, heat and mass transfer and growth kinetics. The transport equations for electromagnetic field, heat transfer, and species transport are solved using a finite volume-based numerical scheme called MASTRAPP (Multizone Adaptive Scheme for Transport and Phase Change Process). Temperature distribution for a 2-inch growth system is calculated, and the effects of induction heating frequency and current on the temperature distribution and growth rate are investigated. The predicted results have been compared with the experimental data.

Research on Bearing Capacity Calculation Method of Large Diameter Cylinder in Soft Clay

A large diameter cylinder inserted in soils is a new type of engineering structures used in offshore and port engineering. The mechanism of its bearing capacity and the analysis of its stability are important to its design and applications. In this paper, the finite element method is used to analyze the reacting forces of the soft soil foundation on the structure under the wave action. A simplified method is proposed, based on the plastic limit method, for the safety and stability analysis. Our analysis shows that the assumptions made in this paper and the mechanism used are reasonable, and the results obtained are appropriate. The calculation method is very efficient and can be used to evaluate main parameters of the structure in its preliminary designs.

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.

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.

An Instrumented Indentation Method for Young's Modulus Measurement with Accuracy Estimation

The relationships between indentation responses and Young's modulus of an indented material were investigated by employing dimensional analysis and finite element method. Three representative tip bluntness geometries were introduced to describe the shape of a real Berkovich indenter. It was demonstrated that for each of these bluntness geometries, a set of approximate indentation relationships correlating the ratio of nominal hardness/reduced Young's modulus H (n) /E (r) and the ratio of elastic work/total work W (e)/W can be derived. Consequently, a method for Young's modulus measurement combined with its accuracy estimation was established on basis of these relationships. The effectiveness of this approach was verified by performing nanoindentation tests on S45C carbon steel and 6061 aluminum alloy and microindentation tests on aluminum single crystal, GCr15 bearing steel and fused silica.

Finite element analysis of stress and strain distributions in InAs/GaAs quantum dots

In this paper, we perform systematic calculations of the stress and strain distributions in InAs/GaAs truncated pyramidal quantum dots (QDs) with different wetting layer (WL) thickness, using the finite element method (FEM). The stresses and strains are concentrated at the boundaries of the WL and QDs, are reduced gradually from the boundaries to the interior, and tend to a uniform state for the positions away from the boundaries. The maximal strain energy density occurs at the vicinity of the interface between the WL and the substrate. The stresses, strains and released strain energy are reduced gradually with increasing WL thickness. The above results show that a critical WL thickness may exist, and the stress and strain distributions can make the growth of QDs a growth of strained three-dimensional island when the WL thickness is above the critical value, and FEM can be applied to investigate such nanosystems, QDs, and the relevant results are supported by the experiments.

Reactive Mold Filling in Resin Transfer Molding Processes with Edge Effects

Reactive mold filling is one of the important stages in resin transfer molding processes, in which resin curing and edge effects are important characteristics. On the basis of previous work, volume-averaging momentum equations involving viscous and inertia terms were adopted to describe the resin flow in fiber preform, and modified governing equations derived from the Navier-Stokes equations are introduced to describe the resin flow in the edge channel. A dual-Arrhenius viscosity model is newly introduced to describe the chemorheological behavior of a modified bismaleimide resin. The influence of the curing reaction and processing parameters on the resin flow patterns was investigated.

Computer simulation of the damage evolution of fiber-reinforced polypropylene-matrix composites with matrix defects

The damage evolution of fiber-reinforced polypropylene-matrix composites with matrix defects was studied via a Monte Carlo technique combined with a finite element method. A finite element model was constructed to predict the effects of various matrix defect shapes on the stress distributions. The results indicated that a small matrix defect had almost no effect on fiber stress distributions other than interfacial shear stress distributions. Then, a finite element model with a statistical distribution of the fiber strength was constructed to investigate the influences of the spatial distribution and the volume fraction of matrix defects on composite failure. The results showed that it was accurate to use the shear-lag models and Green's function methods to predict the tensile strength of composites even though the axial stresses in the matrix were neglected.

Stress transfer and damage evolution simulations of fiber-reinforced polymer-matrix composites

The stress transfer from broken fibers to unbroken fibers in fiber-reinforced thermosetting polymer-matrix composites and thermoplastic polymer-matrix composites was studied using a detailed finite element model. In order to check the validity of this approach, an epoxy-matrix monolayer composite was used as thermosetting polymer-matrix composite and a polypropylene (PP)-matrix monolayer composite was used as thermoplastic polymer-matrix composite, respectively. It is found that the stress concentrations near the broken fiber element cause damage to the neighboring epoxy matrix prior to the breakage of other fibers, whereas in the case of PP-matrix composites the fibers nearest to the broken fiber break prior to the PP matrix damage, because the PP matrix around the broken fiber element yields. In order to simulate composite damage evolution, a Monte Carlo technique based on a finite element method has been developed in the paper. The finite element code coupled with statistical model of fiber strength specifically written for this problem was used to determine the stress redistribution. Five hundred samples of numerical simulation were carried out to obtain statistical deformation and failure process of composites with fixed fiber volume fraction.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

滑坡型泥石流转化机理研究—以四川省达县青宁乡滑坡为例

Slide-debris flow is debris flow which is transformed from landslide consecutively in a short time, it comprises of two phases: First, Landslide starts to slide; Second, landslide changes to debris flow. Slide-debris flow which brings great property and life loss happens frequently at home and abroad. In order to forecast the happening possibility and scope of slide-debris flow, transfromation mechanism of Slide-debris flow must be studied. Research on transformation mechanism of slide-debris flow is intersectant science of landslide kinetics and debris flow starting theory, It is a fringe problem as well as front problem of geological hazard. This paper takes Qingning slide-debris flow in Da County, Sichuan Province for example and has studied the mechanism of its instability and transfromation into debris flow through indoor test (including usual soil test and ring shear test) and digital modeling method.The research gets the following conclusions. Qingning Landslide took place mainly because of confined water head arising from rainfall infiltration. Before Landslide occurring, it rained continuously for 22 days, accumulated precipitation arrived at 521.6mm.Investigation shows that strata of Qingning Landslide contains quaternary loose accumulation, slip soil and highly weathered bedrock, which is a good condition for formation of confined water in the slope. Further more, groundwater seepage in the slope body and corresponding slope safety factor before landslide occurring have been computed through finite element method. The result shows that because of infiltration of rainfall, confined water head in the slope arose sharply, accordingly, the safety factor of the slope declined quickly. The result also shows that force put on the slide body by the rock mass detached from Dazhaiyan mountain was the direct factor for landslide occurring. Qingning slide-debris transformation mode has been summarized, the process the landslide changed into debris flow is divided into three phases in the prospective of macroscopic geological condition: landslide occurring, transformation and debris flow. Landslide occurring phase is from slope’ local creeping slide to Landslide occurring; transformation phase contains slide body sliding on the slide bed after slide occurring and sliding on the slope after shearing opening; debris flow phase is that slide body breaks up completely and flows downward into the ditches. The transformation mechanism of Qingning slide-debris flow has been studied through indoor ring shear test of slip soil. The result shows that transformation mechanism contains two points: first, during slide body sliding on the slide bed and slope after shearing opening, shearing shrinkage, grain crushing and grain layering brought about declining of its volume and produced excess pore water pressure, and because producing velocity of excess pore water pressure is much greater than its dissipating velocity, shear strength of slide body decreased sharply because of accumulated pore water pressure. Second, grains crushing and grains layering during slide body sliding brought about thick liquefied layer at the bottom of the slidebody, liquefied layer contained high water content and its shear strength was very low, its thickness increased as the sliding displacement increasing. Liquefied layer makes slide body sliding fast and easily break down to debris flow. Excess pore water pressure and liquefied layer made shear strength of slidebody became very low, furthermore, water in the pit of slope joining in the slidebody was also a facter that made slidebody accelerate the transformation. Influence of slide body thickness and fine grains content to transformation of slide-debris flow has been studied through ring shear test. The result reaches two conclusions. First, thickness of slide body affects transformation of slide-debris flow by two ways, porewater pressure and effect of “soft base” increases as thickness of slide body increasing.so the thicker slide body is ,the easier transformation is. Second, actual dissipating velocity of porewater pressure should be considered when studying the influence of fine grains content to tranformation of slide-debris flow. There should be a critical content of fine grains which makes the difference of producing and dissipating velocity of water pore pressre greatest, this value is the best for slide-debris transformation. The whole process of slide-debris flow transformation is reproduced through discrete element method. Transformation mechanism of slide-debris flow is studied through monitoring various parameters including pore water pressure, grain crushing and grain layering in the slide body during the transformation. The result confirms and supplements the transformation mechanism of slide-debris flow got from ring shear test well.

同位网格摄动有限体积格式求解浮力驱动方腔流

利用对流扩散方程的摄动有限体积格式,在Rayleigh数从10~3到10~8的范围内对浮力驱动方腔流动问题作了数值模拟.对流扩散方程的摄动有限体积格式具有一阶迎风格式的简洁形式,使用相同的基点,重构近似精度高,特别是两相邻控制体中心到公共界面的距离相等或不相等,PFV格式公式相同等优点.在数值模拟中,无论均匀网格还是非均匀网格均获得与DSC方法、自适应有限元法、多重网格法等Benchmark解相符较好的数值结果,证明UPFV格式对高Rayleigh数对流传热问题的适用性和有效性.