潮汐河口平面二维非恒定流数值模拟

本文基于三角形网格划分，采用有限元法，建立了平面二维非恒定潮流数学模型。采用了质量集中、压缩存储的处理方法和“预估校正”的时间推进算法，较好地解决了有限元计算存储量和计算速度问题。对长江口南北支河段进行了验证计算，结果表明模型能较好模拟长江口南北支河段的潮流运动情况，且计算稳定性好、速度快、精度高。     

Microstructure evolution during metal forming processes

Recrystallization and grain growth evolutions during metal forming processes are considered. Coupling between the thermo-mechanical and microstructure processes is realized. Die forging of a rear-axle flange is simulated numerically on the base of the finite element method. Material parameters of the models are obtained experimentally. The influence of interpass and holding times on grain size distributions in the end product is shown.     

Detection of multiple cracks using frequency measurements

A method for detection of multiple open cracks in a slender Euler-Bernoulli beams is presented based on frequency measurements. The method is based on the approach given by Hu and Liang [J. Franklin Inst. 330 (5) (1993) 841], transverse vibration modelling through transfer matrix method and representation of a crack by rotational spring. The beam is virtually divided into a number of segments, which can be decided by the analyst, and each of them is considered to be associated with a damage parameter. The procedure gives a linear relationship explicitly between the changes in natural frequencies of the beam and the damage parameters. These parameters are determined from the knowledge of changes in the natural frequencies. After obtaining them, each is treated in turn to exactly pinpoint the crack location in the segment and determine its size. The forward, or natural frequency determination, problems are examined in the passing. The method is approximate, but it can handle segmented beams, any boundary conditions, intermediate spring or rigid supports, etc. It eliminates the need for any symbolic computation which is envisaged by Hu and Liang [J. Franklin Inst. 330 (5) (1993) 841] to obtain mode shapes of the corresponding uncracked beams. The proposed method gives a clear insight into the whole analysis. Case studies (numerical) are presented to demonstrate the method effectiveness for two simultaneous cracks of size 10% and more of section depth. The differences between the actual and predicted crack locations and sizes are less than 10% and 15% respectively. The numbers of segments into which the beam is virtually divided limits the maximum number of cracks that can be handled. The difference in the forward problem is less than 5%.     

The Problem of Estimating the Uncertainty of the Fixed Point Temperature Due to the Presence of Impurities

The main problems in estimating the uncertainty due to the presence of impurities in the material: the absence in the documents of a unique definition of the quantity being measured – the fixed point temperature, the lack of agreed methods of realizing the fixed points and the insufficient justification of the equations used for the estimate, are analyzed. The results of a measurement of the melting point of gallium are presented, which differ in the conditions in which they were realized – with one and two liquid-solid interfaces and different solidification conditions, are presented. The reasons for the considerable difference between the values of the temperatures obtained in these experiments are discussed.     

多层复合膜的成型技术

阐述了用AutoLisp程序模拟塑料熔体在2层复合膜机头内的流动情况，并介绍了复合膜机头的内部结构。     

T-stress in orthotropic functionally graded materials: Lekhnitskii and Stroh formalisms

A new interaction integral formulation is developed for evaluating the elastic T-stress for mixed-mode crack problems with arbitrarily oriented straight or curved cracks in orthotropic nonhomogeneous materials. The development includes both the Lekhnitskii and Stroh formalisms. The former is physical and relatively simple, and the latter is mathematically elegant. The gradation of orthotropic material properties is integrated into the element stiffness matrix using a “generalized isoparametric formulation” and (special) graded elements. The specific types of material gradation considered include exponential and hyperbolic-tangent functions, but micromechanics models can also be considered within the scope of the present formulation. This paper investigates several fracture problems to validate the proposed method and also provides numerical solutions, which can be used as benchmark results (e.g. investigation of fracture specimens). The accuracy of results is verified by comparison with analytical solutions.     

On connections between 3-d microstructures and their 2-d images

The problem of recovery of three-dimensional microstructural parameters from two-dimensional images is addressed in its inverse form. Three-dimensional inhomogeneous solids with known microstructures are generated virtually and images of their cross-sections are analyzed. Three different types of randomly oriented microcracks and microcracks having preferential orientation with some scatter. Somewhat unexpected result is that, in all the cases considered, three-dimensional defect densities are very close to average values of corresponding two-dimensional parameters.     

Accurate higher-order analytical approximate solutions to nonconservative nonlinear oscillators and application to van der Pol damped oscillators

In general, this paper deals with general nonlinear oscillations of a nonconservative and single degree-of-freedom system with odd nonlinearity and, in particular, it presents accurate higher-order analytical approximate solutions to van der Pol damped nonlinear oscillators having odd nonlinearity and the Rayleigh equation. By combining the linearization of the governing equation with harmonic balancing and the method of averaging, we establish accurate analytical approximate solutions for the general weakly damped nonlinear systems. Unlike the classical harmonic balance method, simple linear algebraic equations instead of nonlinear algebraic equations are obtained upon linearization prior to harmonic balancing. The combination of these two methods results in very accurate transient response of the periodic solution. In addition and for the first time, this paper also presents a method for deducing fourth-, fifth- and higher-order linearized governing equations from the lower-order equations without the requirement of formulating the problem from the first principle. Three examples including the van der Pol damped nonlinear oscillator are presented to illustrate the excellent agreement with approximate solution using the exact frequency.     

Origin of the negative wake behind a bubble rising in non-Newtonian fluids

The present work aims at understanding the behavior of individual bubbles in non-Newtonian fluids. By means of a Particle Image Velocimetry (PIV) device, the complete flow field around either a single non-spherical bubble rising in polyacrylamide (PAAm) solutions or a solid sphere settling down in the same fluids shows for the first time the similar coexistence of three distinct zones: a central downward flow behind the bubble or the sphere (negative wake), a conical upward flow surrounding the negative wake zone, and an upward flow zone in front of the bubble or the sphere. This excludes then the possible influence of the interface deformation on the negative wake. A theoretical lattice Boltzmann scheme coupled to a sixth-order Maxwell model was developed for computing the complex flow field around a solid sphere. The good agreement with the experimental measurements provides evidence that the physical mechanism responsible for the negative wake in such fluids could be related to the fluid's viscoelastic properties.