Dynamic function of damage and its implications

Dynamic function of damage is the key to the problem of damage evolution of solids. In order to understand it, one must understand its mesoscopic mechanisms and macroscopic formulation. In terms of evolution equation of microdamage and damage moment, a dynamic function of damage is strictly defined. The mesoscopic mechanism underlying self-closed damage evolution law is investigated by means of double damage moments. Numerical results of damage evolution demonstrate some common features for various microdamage dynamics. Then, the dynamic function of damage is applied to inhomogeneous damage field. In this case, damage evolution rate is no longer equal to the dynamic function of damage. It is found that the criterion for damage localization is closely related to compound damage. Finally, an inversion of damage evolution to the dynamic function of damage is proposed.

Dynamic stress intensity factor of a functionally graded material under antiplane shear loading

The dynamic response of a finite crack in an unbounded Functionally Graded Material (FGM) subjected to an antiplane shear loading is studied in this paper. The variation of the shear modulus of the functionally graded material is modeled by a quadratic increase along the direction perpendicular to the crack surface. The dynamic stress intensity factor is extracted from the asymptotic expansion of the stresses around the crack tip in the Laplace transform plane and obtained in the time domain by a numerical Laplace inversion technique. The influence of graded material property on the dynamic intensity factor is investigated. It is observed that the magnitude of dynamic stress intensity factor for a finite crack in such a functionally graded material is less than in the homogeneous material with a property identical to that of the FGM crack plane.

A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.

非饱和土中镉离子传输模型参数反演

使用自行设计的的真空实验装置,采用一维实验室土柱实验方法在不同实验条件下进行了15组实验,并基于实验所得的穿透曲线,使用梯度正则化方法反演得到了镉离子在8种葡萄牙土样和1种比利时土样中相应的传输模型参数.同时通过数值模拟实验验证了梯度正则化方法的有效性和可靠性,进而证明所得结果是有效的.

Torsional impact response of a penny-shaped interface crack in bonded materials with a graded material interlayer

In this paper, the dynamic response of a penny-shaped interface crack in bonded dissimilar homogeneous half-spaces is studied. It is assumed that the two materials are bonded together with such a inhomogeneous interlayer that makes the elastic modulus in the direction perpendicular to the crack surface is continuous throughout the space. The crack surfaces art assumed to be subjected to torsional impact loading. Laplace and Hankel integral transforms are applied combining with a dislocation density,function to reduce the mixed boundary value problem into a singular integral equation with a generalized Cauchy kernel in Laplace domain. By solving the singular integral equation numerically, and using a numerical Laplace inversion technique, the dynamic stress intensity factors art obtained. The influences of material properties and interlayer thickness on the dynamic stress intensity factor are investigated.

Cauchy singular integral equation method for transient antiplane dynamic problems

In this paper, by use of the boundary integral equation method and the techniques of Green basic solution and singularity analysis, the dynamic problem of antiplane is investigated. The problem is reduced to solving a Cauchy singular integral equation in Laplace transform space. This equation is strictly proved to be equivalent to the dual integral equations obtained by Sih [Mechanics of Fracture, Vol. 4. Noordhoff, Leyden (1977)]. On this basis, the dynamic influence between two parallel cracks is also investigated. By use of the high precision numerical method for the singular integral equation and Laplace numerical inversion, the dynamic stress intensity factors of several typical problems are calculated in this paper. The related numerical results are compared to be consistent with those of Sih. It shows that the method of this paper is successful and can be used to solve more complicated problems. Copyright (C) 1996 Elsevier Science Ltd

3-dimensional transient wave response in a cracked elastic solid

The three-dimensional transient wave response problem is presented for an infinite elastic medium weakened by a plane crack of infinite length and finite width. Tractions are applied suddenly to the crack, which simulates the case of impact loading. The integral transforms are utilized to reduce the problem to a standard Fredholm integral equation in the Laplace transform variable and sequentially invert the Laplace transforms of the stress components by numerical inversion method. The dynamic mode I stress intensity factors at the crack tip are obtained and some numerical results are presented in graphical form.

Non-linear theory of a positive-column in a magnetic-field

A nonlinear theory of an intermediate pressure discharge column in a magnetic field is presented. Motion of the neutral gas is considered. The continuity and momentum transfer equations for charged particles and neutral particles are solved by numerical methods. The main result obtained is that the rotating velocities of ionic gas and neutral gas are approximately equal. Bohm's criterion and potential inversion in the presence of neutral gas motion are also discussed.

YAG激光毛化技术进展

二十世纪八十年代，当比利时冶金研究中心(CRM)开发出CO_2激光毛化冷轧辊技术后，尝试用YAG激光进行轧辊毛化一直吸引着众多的研究者，这是因为YAG(1.06μm)激光波长比CO_2(10.6μm)激光波长短一个量级，材料对YAG激光有更高的吸收率，并用YAG激光可以聚焦到更小的光斑尺寸，同时使用电信号驱动的声光开关技术便于对毛化分布进行可设定控制。但是用传统声光调制的YAG激光虽然可以碇以很高的脉冲频率(>30kHz)，但单脉冲有量仅为10mJ左右，难以达到辊面毛化粗糙度的要求，因此人们认为YAG激光用于毛化的主要困难是脉冲能量太小。

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.

Study On Buoyancy Convection Phenomenon In The Crystal Growth Process

Real-time phase shift Mach-Zehnder interference technique, imaging technique, and computer image processing technique were combined to perform a real-time diagnosis of NaClO3 crystal, which described both the dissolution process and the crystallization process of the NaClO3 crystal in real-time condition. The dissolution fringes and the growth fringes in the process were obtained. Moreover, a distribution of concentration field in this process was obtained by inversion calculation. Finally, the buoyancy convection phenomenon caused by gravity in the crystal growth process was analyzed. The results showed that this convection phenomenon directly influences the growth rate of each crystal face in the crystal.

用于多无人机编队飞行的PIDA+逆飞行控制器

设计了一种结合PIDA控制及逆控制的多无人机编队飞行控制器.分析了多无人机编队飞行中飞行耦合气流的估计范围,针对个体无人机的轨迹控制模型讨论了其飞行控制器需求,结合逆控制的解耦性能和PIDA控制器低模型依赖、简单结构、快速性等优点,设计了一种PIDA+逆的飞行控制器,最后,以控制器输出约束、飞行气动力参数变化以及气动耦合扰动为指标对PID控制和PIDA+逆控制的控制效果进行了对比仿真,仿真结果表明,PIDA十逆控制能够有效的抑制模型参数不确定性带来的影响,并且在控制器输出约束条件下仍能保持较好的动态性、准确性和鲁棒性.

Space-resolved spectral diagnosis of line-shaped laser-produced magnesium plasmas

Space-resolved spectra of line-shaped laser-produced magnesium plasmas in the normal direction of the target have been obtained using a pinhole crystal spectrograph. These spectra are treated by a spectrum analyzing code for obtaining the true spectra and fine structures of overlapped lines. The spatial distributions of electron temperature and density along the normal direction of the target surface have been obtained with different spectral diagnostic techniques. Especially, the electron density plateaus beyond the critical surface in line-shaped magnesium plasmas have been obtained with a fitting technique applied to the Stark-broadened Ly-alpha wings of hydrogenic ions. The difference of plasma parameters between those obtained by different diagnostic techniques is discussed. Other phenomena, such as plasma satellites, population inversion, etc., which are observed in magnesium plasmas, are also presented.

Enhancing Kerr nonlinearity in an asymmetric double quantum well via Fano interference

We investigate the enhancement of Kerr nonlinearity in an asymmetric GaAs double quantum well via Fano interference, which is caused by tunneling from the excited subband to the continuum. In our structure, owing to Fano interference, the Kerr nonlinearity can be enhanced by appropriately choosing the values of the detunings and the intensity of the pump field, while cancel the linear and nonlinear absorptions.

Group velocity of probe light pulse in an open lambda system

The group velocity of the probe light pulse (GVPLP) propagating through an open Lambda-type atomic system with a spontaneously generated coherence is investigated when the weak probe and strong driving light fields have different frequencies. It is found that adjusting the detuning or Rabi frequency of the probe light field can realize switching of the GVPLP from subluminal to superluminal. Changing the relative phase between the probe and driving light. elds or atomic exit and injection rates can lead to GVPLP varying in a wider range, but cannot induce transformation of the property of the GVPLP. The absolute value of the GVPLP always increases with Rabi frequency of the driving light field increasing. For subluminal and superluminal propagation, the system always exhibits the probe absorption, and GVPLP is mainly determined by the slope of the steep dispersion.