Collapsed isoparametric element as a singular element for a crack normal to the bi-material interface

It is shown that the variable power singularity of the strain field at the crack tip can be obtained by the simple technique of collapsing quadrilateral isoparametric elements into triangular elements around the crack tip and adequately shifting the side-nodes adjacent to this crack tip. The collapsed isoparametric elements have the desired singularity at crack tip along any ray. The strain expressions for a single element have been derived and in addition to the desired power singularity, additional singularities are revealed. Numerical examples have shown that triangular elements formed by collapsing one side lead to excellent results.

Potential energy model for turbulent entrainment across density interface in a linearly stratified fluid

A potential energy model is developed for turbulent entrainment in the absence of mean shear in a linearly stratified fluid. The relation between the entrainment distance D and the time t and the relation between dimensionless entrainment rate E and the local Richardson number are obtained. An experiment is made for examination. The experimental results are in good agreement with the model, in which the dimensionless entrainment distance D is given by DBAR = A(i)(SBAR)-1/4(fBAR)1/2(tBAR)1/8, where A(i) is the proportional coefficient, S is the dimensionless stroke, fBAR is the dimensionless frequency of the grid oscillation, tBAR the dimensionless time.

Influence of suspended particles and particle cover over the density interface

The effect of the particle cover over the density interface between two layers of fluids and of the suspended solid particles in the upper turbulcnt layer on the turbulent entrainment has been studied experimentally. The entrainment distance D is a function of the time of power: D=kt, where =0.200-0.130p. For suspended particles in the upper layer and pure 2-layer fluid is equal to 0.200, but the value of k for the suspended particles is smaller than that for the pure 2-layer fluid. The non-dimensional entrainment velocity is E=KRiln, where n=1.50+0.93 p. It is shown that the particle cover over the interface changes the power of Ril in the entrainment and hinders the turbulent entrainment. The variation rule of E for the suspended particles is the same as that for the pure 2-layer fluid, but the K value of the former is smaller than that of the latter. The turbulent mixing mechanism has been discussed.

Crack deflection at an interface between dissimilar elastic-materials

A crack intersecting an interface between two dissimilar materials may advance by either penetrating through the interface or deflecting into the interface. The competition between deflection and penetration can be assessed by comparison of two ratios: (i) the ratio of the energy release rates for interface cracking and crack penetration; and (ii) the ratio of interface to material fracture energies. Residual stresses caused by thermal expansion misfit can influence the energy release rates of both the deflected and penetrating crack. This paper analyses the role of residual stresses. The results reveal that expansion misfit can be profoundly important in systems with planar interfaces (such as layered materials, thin film structures, etc.), but generally can be expected to be of little significance in fiber composites. This paper corrects an earlier result for the ratio of the energy release rate for the doubly deflected crack to that for the penetrating crack in the absence of residual stress.

Influence of secondary void damage in the matrix material around voids

The mechanism of ductile damage caused by secondary void damage in the matrix around primary voids is studied by large strain, finite element analysis. A cylinder embedding an initially spherical void, a plane stress cell with a circular void and plane strain cell with a cylindrical or a flat void are analysed under different loading conditions. Secondary voids of smaller scale size nucleate in the strain hardening matrix, according to the requirements of some stress/strain criteria. Their growth and coalescence, handled by the empty element technique, demonstrate distinct mechanisms of damage as circumstances change. The macroscopic stress-strain curves are decomposed and illustrated in the form of the deviatoric and the volumetric parts. Concerning the stress response and the void growth prediction, comparisons are made between the present numerical results and those of previous authors. It is shown that loading condition, void growth history and void shape effect incorporated with the interaction between two generations of voids should be accounted for besides the void volume fraction.

The effect of strain softening in the matrix material during void growth

Furthermore, the compressed flow driven by the piston is discussed. The consistent solution of gasdynamical equations including solar gravity is obtained for the unsteady and two-dimensional configuration, which is applied to the region between the piston and shock wave. This solution may satisfy the jump conditions of shock wave, which separates the region of compressed flow and quiet corona.

颗粒增强复合材料的界面开裂与尺度效应

采用基于Huang等提出的塑性应变梯度传统理论发展的有限元方法，模拟了颗粒增强金属基复合材料的界面开裂与颗粒尺度效应．分别针对考虑颗粒与基体间界面开裂和不开裂两种情况进行分析，并将考虑界面开裂的模拟结果与实验结果进行比较，证明了模型的有效性，同时也获得应变梯度理论中所包含的材料特征尺度参量的取值．

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.

Effects of Thermocapillary Convection on the Growth Process in Industrial 8-inch Silicon Czochralski Growth

Czochralski (CZ) crystal growth process is a widely used technique in manufacturing of silicon crystals and other semiconductor materials. The ultimate goal of the IC industry is to have the highest quality substrates, which are free of point defect, impurities and micro defect clusters. The scale up of silicon wafer size from 200 mm to 300 mm requires large crucible size and more heat power. Transport phenomena in crystal growth processes are quite complex due to melt and gas flows that may be oscillatory and/or turbulent, coupled convection and radiation, impurities and dopant distributions, unsteady kinetics of the growth process, melt crystal interface dynamics, free surface and meniscus, stoichiometry in the case of compound materials. A global model has been developed to simulate the temperature distribution and melt flow in an 8-inch system. The present program features the fluid convection, magnetohydrodynamics, and radiation models. A multi-zone method is used to divide the Cz system into different zones, e.g., the melt, the crystal and the hot zone. For calculation of temperature distribution, the whole system inside the stainless chamber is considered. For the convective flow, only the melt is considered. The widely used zonal method divides the surface of the radiation enclosure into a number of zones, which has a uniform distribution of temperature, radiative properties and composition. The integro-differential equations for the radiative heat transfer are solved using the matrix inversion technique. The zonal method for radiative heat transfer is used in the growth chamber, which is confined by crystal surface, melt surface, heat shield, and pull chamber. Free surface and crystal/melt interface are tracked using adaptive grid generation. The competition between the thermocapillary convection induced by non-uniform temperature distributions on the free surface and the forced convection by the rotation of the crystal determines the interface shape, dopant distribution, and striation pattern. The temperature gradients on the free surface are influenced by the effects of the thermocapillary force on the free surface and the rotation of the crystal and the crucible.

The Effect of the Thickness of Fe2 Al5 Phase Layer at Fe/Al Interface on the Mechanics Behavior of Hot Dip Aluminizing Coating

Hot Dip Aluminized Coatings with different thickness were prepared on Q235 steel in aluminum solutions with different temperature for certain time. Through tensile tests and in-situ SEM observations, the effect of the coating's microstructure on the tensile strength of the samples was studied. It was disclosed at certain aluminum solution temperature,transaction layers mainly composed of Fe2 Al5 phase got thicker with time prolonging, and this changed initial crack's extending direction from parallel with to vertical with stretching direction. The change in crack direction decreased tensile strength of samples, thus made the coating easy to break. It was concluded that the existence of thick Fe2 Al5 phase layer was the basic reason for the lowering of tensile strength of the coating.

Investigations of cohesive zone properties of the interface between metal film and ceramic substrate at nano- and micro-scales

Cohesive zone characterizations of the interface between metal film and ceramic substrate at micro- and nano-scales are performed in the present research. At the nano-scale, a special potential for special material interface (Ag/MgO) is adopted to investigate the interface separation mechanism by using MD simulation, and stress-separation relationship will be obtained. At the micro-scale, peeling experiment is performed for the Al film/Al2O3 substrate system with an adhesive layer at the interface. Adhesive is a mixture of epoxy and polyimide with mass ratio 1:1, by which a brittle cohesive property is obtained. The relationships between energy release rate, the film thickness and the adhesive layer thickness are measured during the steady-state peeling. The experimental result has a similar trend as modeling result for a weak adhesion interface case.