Kinking of an interface crack between 2 dissimilar anisotropic elastic solids

A detailed analysis of kinking of an interface crack between two dissimilar anisotropic elastic solids is presented in this paper. The branched crack is considered as a distributed dislocation. A set of the singular integral equations for the distribution function of the dislocation density is developed. Explicit formulas of the stress intensity factors and the energy release rates for the branched crack are given for orthotropic bimaterials and misoriented orthotropic bicrystals. The role of the stress parallel to the interface, sigma0 is taken into account in these formulas. The interface crack can advance either by continued extension along the interface or by kinking out of the interface into one of the adjoining materials. This competition depends on the ratio of the energy release rates for interface cracking and for kinking out of the interface and the ratio of interface toughness to substrate toughness. Throughout the paper, the influences of the inplane stress sigma0 on the stress intensity factors and the energy release rates for the branched crack, which can significantly alter the conditions for interface cracking, are emphasized.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Numerical study of bluff body flow structures

Our recent progress in numerical studies of bluff body flow structures and a new method for the numerical analysis of near wake flow field for high Reynolds number flow are introduced. The paper consists of three parts. In part one, the evolution of wake vortex structure and variation of forces on a flat plate in harmonic oscillatory flows and in in-line steady-harmonic combined flows are presented by an improved discrete vortex method, as the Keulegan-Carpenter number (KC) varies from 2 to 40 and ratios of U-m to U-0 are of O(10(-1)), O(10) and O(10), respectively. In part 2, a domain decomposition hybrid method, combining the finite-difference and vortex methods for numerical simulation of unsteady viscous separated flow around a bluff body, is introduced. By the new method, some high resolution numerical visualization on near wake evolution behind a circular cylinder at Re = 10(2), 10(3) and 3 x 10(3) are shown. In part 3, the mechanism and the dynamic process for the three-dimensional evolution of the Karman vortex and vortex filaments in braid regions as well as the early features of turbulent structure in the wake behind a circular cylinder are presented numerically by the vortex dynamics method.

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.

Coupling of melting-solidification phase-change convection with floating-zone convection under microgravity

A finite element algorithm is used to analyze the process of floating zone crystal growth under microgravity. The effect of phase change convection coupled with surface tension convection is considered. The results show that the rate of crystal growth is very important. The single-crystal-melt interface is steeper than the feed-melt interface during the process of crystal growth. When the rate exceeds a critical value, the Marangoni vortex near the feed-melt interface will become so large that a secondary vortex will exist.

Preliminary-results of gaas single-crystal growth under high gravity conditions

GaAs single crystals have been grown under high gravity conditions, up to 9g0, by a recrystallization method with decreasing temperature. The impurity striations in GaAs grown under high gravity become weak and indistinct with smaller striation spacings. The dislocation density of surcharge-grown GaAs increases with increase of centrifugal force. The cathodoluminescence results also show worse perfection in the GaAs grown at high gravity than at normal earth gravity.

Crack extension and kinking in laminates and bicrystals

Singular fields at the tip of an interface crack in anisotropic solids are reviewed with emphasis on establishing a framework to quantify fracture resistance under mixed mode conditions. The concepts of mode mixity and surface toughness are unified by using generalized interface traction components. The similarity between the anisotropic theory and existing isotropic theory is shown. Explicit formulae are given for misoriented orthotropic bimaterials with potential applications envisioned including composite laminates and semiconductor crystals. Competition between crack extension along the interface and kinking into the substrate is investigated using a boundary layer formulation. Several case studies reveal the role of anisotropy. An explicit complex variable representation for orthotropic materials and a solution to a dislocation interacting with a crack are presented in two self-contained Appendices.

Mixed-mode fracture mechanisms near the fatigue threshold of aisi-316 stainless-steel

Near threshold, mixed mode (I and II), fatigue crack growth occurs mainly by two mechanisms, coplanar (or shear) mode and branch (or tensile) mode. For a constant ratio of ΔK_I/ΔK_II the shear mode growth shows a self-arrest character and it would only start again when ΔK_I and ΔK_II are increased. Both shear crack growth and the early stages of tensile crack growth, are of a crystallographic nature; the fatigue crack proceeds along slip planes or grain boundaries. The appearance of the fracture surfaces suggest that the mechanism of crack extension is by developing slip band microcracks which join up to form a macrocrack. This process is thought to be assisted by the nature of the plastic deformation within the reversed plastic zone where high back stresses are set up by dislocation pile-ups against grain boundaries. The interaction of the crack tip stress field with that of the dislocation pile-ups leads to the formation of slip band microcracks and subsequent crack extension. The change from shear mode to tensile mode growth probably occurs when the maximum tensile stress and the microcrack density in the maximum tensile plane direction attain critical values.

SiC_p/6151Al复合材料的微结构效应

利用分离式Hopkinson压杆和MTS通用材料试验机研究了SiC_p/6151Al颗粒增强复合材料在不同应变率下的变形行为和增强颗粒的尺寸对复合材料微结构及变形行为的影响。结果表明，对于在不同应变率下的SiC_p/6151Al复合材料，增强颗粒尺寸大小的流动应力高于增强颗粒尺寸的流动应力。根据位错强化理论中的Hall-Petch关系对这个结果进行了解释。首次在实验上观测到增强颗粒对复合材料微损伤-微带形成的影响，并根据微带（microband）形成的双位错墙理论（double dislocation walls）, 分析了增强颗粒对复合材料微带损伤及力学性能影响的微结构效应。

高压扭转铜试样的微观组织与压缩性能

通过高压扭转对铜试样施加不同程度的变形，研究了样品扭转面（ND面）和纵截面（TD面）上微观组织特征．对ND面，在较小的剪应变下，原始晶粒形貌模糊，晶粒内部形成等轴状的位错胞及亚晶结构；随变形量的增大，亚晶间取向差及亚晶内部的位错密度增大，最后形成亚微米尺度的等轴晶粒．对TD面，变形初期原始晶粒被拉长，晶粒内部为位错墙分割成的层状结构，层内为拉长的位错胞；随变形程度的增大，拉长晶粒的宽度减小，与剪切方向的夹角减小，晶内层状组织间距减小，并逐渐演化成拉长的亚晶组织；进一步增大变形，晶粒拉长痕迹消失，变形组织与ND面相似，为等轴状亚微米晶粒．压缩实验表明，经16圈扭转后，整个试样上的压缩性能基本均匀，σ0．2达到385MPa，应变率敏感性指数增大至0．021．

水中悬浮隧道锚索在剪切流中的涡激响应

以Iwan改进的尾流振子模型为基础，给出了剪切流中水中悬浮隧道锚索横向涡激振动的工程分析方法，分析了悬浮隧道重浮比以及剪切流特性对锚索涡激响应的影响．计算结果表明，悬浮隧道重浮比的改变，使锚索各阶模态的频锁区域分布发生变化，从而改变了结构的涡激振动响应；剪切流陡度参数的增大会使锚索涡激响应的幅值减小，若使用均匀流场代替剪切流场计算悬浮隧道锚索的涡激振动，会过高估计其响应幅值．

中等雷诺数圆柱尾流旋涡脱落的控制

介绍通过实验对圆柱尾流旋涡脱落进行抑制的方法及其结果。实验模型的展径比为38，实验的雷诺数范围为3 * 10~2～1.6 * 10~3。抑制方法是在圆柱（直径为D）表面沿展向每隔一定间距伸出一直径0.18D、长度为1.5D的小棒，实验结果表明，当棒间距小于3D，棒与来流夹角在30 °～ 90 °范围内，可有效抑制旋涡脱落。

湍流大小尺度涡量方程组及其封闭

利用流体大小尺度（LSS）方程组推导出湍流大小尺度涡量（LSSV）方程组，给出两个关于湍流大小尺度涡量的命题，从而得到湍流封闭大小尺度涡量（CLSSV）方程组。同时，对近程相互作用命题进行了推广。

Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

Nanocrystalline Formation during Mechanical Attrition of Cobalt

The nanocrystalline (nc) formation was studied in cobalt (a mixture of c (hexagonal close packed) and gamma (face-centered cubic) phases) subjected to surface mechanical attrition treatment. Electron microscopy revealed the operation of {10(1) over bar 0}< 11(2) over bar 0 > prismatic and {0001}< 11(2) over bar 0 > basal slip in the E phase, leading to the successive subdivision of grains to nanoscale. In particular, the dislocation splitting into the stacking faults was observed to occur in ultrafine and nc grains. By contrast, the planar dislocation arrays, twins and martensites were evidenced in the gamma phase. The strain-induced gamma ->epsilon martensitic transformation was found to progress continuously in ultrafine and nc grains as the strain increased. The nc formation in the gamma phase was interpreted in terms of the martensitic transformation and twinning.