A study on time schemes for DRBEM analysis of elastic impact wave

 The precise integration and differential quadrature methods are two new unconditionally stable numerical schemes to approximate time derivative with more than the second order accuracy. Recent studies showed that compared with the Houbolt and Newmark methods, they produced more accurate solutions with large time step for the problems where response is primarily dominated by low and intermediate frequency modes. This paper aims to investigate these time schemes in the context of the dual reciprocity BEM (DRBEM) formulation of various shock-excited scalar elastic wave problems, where high modes have important affect on traction response. The Houbolt method was widely recommended in many literatures for such DRBEM dynamic formulations. However, this study found that the damped Newmark algorithm was the most efficient and accurate for impact traction analysis in conjunction with the DRBEM. The precise integration and differential quadrature methods are shown inapplicable for such shock-excited problems due to the absence of numerical damping. On the other hand, we also found that to achieve the same order of accuracy, the differential quadrature method required much less computing effort than the precise integration method due to the use of the Bartels–Stewart algorithm solving the resulting Lyapunov matrix analogization equation. Received 6 November 2000     

转移法色交换

将图Ｇ的着色由一种变为中一种，通常用Ｋｅｍｐｅ法以交换。但是，对于某些情况，用此法无效。针对这个问题，本文提出了一种转移法色交换，它适用于平面图着色，方法直观，清晰且有效。     

A review of error estimation and adaptivity in the boundary element method

When using the boundary element method, the accuracy of the numerical solution depends critically on the discretization of the boundary into elements (panels). The distribution of the panels is one of the most important decisions taken when analyzing a problem, but still the vast majority of users employ empirical guidelines to distribute the panels. This paper reviews the various adaptive schemes that have been proposed for boundary elements. Numerical results are obtained for infinite fluid flow problems and free surface problems and are used to assess the reliability and effectiveness of each method.     

Non‐reflecting artificial boundaries for transient scalar wave propagation in a two‐dimensional infinite homogeneous layer

This paper presents an exact non‐reflecting boundary condition for dealing with transient scalar wave propagation problems in a two‐dimensional infinite homogeneous layer. In order to model the complicated geometry and material properties in the near field, two vertical artificial boundaries are considered in the infinite layer so as to truncate the infinite domain into a finite domain. This treatment requires the appropriate boundary conditions, which are often referred to as the artificial boundary conditions, to be applied on the truncated boundaries. Since the infinite extension direction is different for these two truncated vertical boundaries, namely one extends toward x →∞ and another extends toward x→‐ ∞, the non‐reflecting boundary condition needs to be derived on these two boundaries. Applying the variable separation method to the wave equation results in a reduction in spatial variables by one. The reduced wave equation, which is a time‐dependent partial differential equation with only one spatial variable, can be further changed into a linear first‐order ordinary differential equation by using both the operator splitting method and the modal radiation function concept simultaneously. As a result, the non‐reflecting artificial boundary condition can be obtained by solving the ordinary differential equation whose stability is ensured. Some numerical examples have demonstrated that the non‐reflecting boundary condition is of high accuracy in dealing with scalar wave propagation problems in infinite and semi‐infinite media. Copyright © 2003 John Wiley & Sons, Ltd.     

Neutron Stress Measurement of W‐Fiber Reinforced Cu Composite

Stress measurement methods using neutron and X‐ray diffraction were examined by comparing the surface stresses with internal stresses in the continuous tungsten‐fiber reinforced copper‐matrix composite. Surface stresses were measured by X‐ray stress measurement with the sin²ψ method. Furthermore, the sin²ψ method and the most common triaxal measurement method using Hooke's equation were employed for internal stress measurement by neutron diffraction. On the other hand, microstress distributions developed by the difference in the thermal expansion coefficients between these two phases were calculated by FEM. The weighted average strains and stresses were compared with the experimental results. The FEM results agreed with the experimental results qualitatively and confirmed the importance of the triaxial stress analysis in the neutron stress measurement.     

水平管气力输送数值模拟

本文采用粒子追踪理论对水平管中粒子流动状态进行了数值模拟.文中提出了一种新的粒子碰撞模型,并忽略粒子的旋转和扬力,仅考虑气流阻力和重力对粒子的影响,在输送量为9.13×10~3kg/s,混合比约0.3的条件下,通过数值计算得知:由于粒子间的碰撞,使得异径粒子的速度趋向平均化;粒子间的碰撞是粒子悬浮的原因之一;粒子的浓度沿管轴向是随时间而发生变化的.     

On the design of multiple key hashing files for concurrent orthogonal range retrieval between two disks

This paper concerns the following problem: given a set of multi-attribute records, a fixed number of buckets and a two-disk system, arrange the records into the buckets and then store the buckets between the disks in such a way that, over all possible orthogonal range queries (ORQs), the disk access concurrency is maximized. We shall adopt the multiple key hashing (MKH) method for arranging records into buckets and use the disk modulo (DM) allocation method for storing buckets onto disks. Since the DM allocation method has been shown to be superior to any other allocation methods for allocating an MKH file onto a two-disk system for answering ORQs, the real issue is knowing how to determine an optimal way for organizing the records into buckets based upon the MKH concept.

A performance formula that can be used to evaluate the average response time, over all possible ORQs, of an MKH file in a two-disk system using the DM allocation method is first presented. Based upon this formula, it is shown that our design problem is related to a notoriously difficult problem, namely the Prime Number Problem. Then a performance lower bound and an efficient algorithm for designing optimal MKH files in certain cases are presented. It is pointed out that in some cases the optimal MKH file for ORQs in a two-disk system using the DM allocation method is identical to the optimal MKH file for ORQs in a single-disk system and the optimal average response time in a two-disk system is slightly greater than one half of that in a single-disk system.     

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.     

Homogenization of hardness distribution and distortion in high pressure gas quenching

Quenching with gases rather than oil or other liquid media has the advantages of reducing the risks concerning health and environment, while simultaneously homogenizing the quenching results and minimizing distortion due to a wide range of possible process parameter variations and the pure convective heat transfer. In this contribution, a coupled solution for increasing homogenization of quenching results within high pressure gas quenching will be presented. In the first stage, an experimental test facility was set up for flow investigations and in the second stage a numerical simulation model was generated. The numerical and experimental results of the flow through the chamber were compared for several boundary conditions. Finally, after complete verification of the simulation, the model may be used to assist in parameter variation for optimization of homogeneous high pressure gas quenching.     

Application of boundary element method to 3-D problems of coupled thermoelasticity

This paper deals with a new boundary element method for analysis of the quasistatic problems in coupled thermoelasticity. Through some mathematical manipulation of the Navier equation in elasticity, the heat conduction equation is transformed into a simpler form, similar to the uncoupled-type equation with the modified thermal conductivity which shows the coupling effects. This procedure enables us to treat the coupled thermoelastic problems as an uncoupled one, A few examples are computed by the proposed BEM, and the results obtained are compared with the analytical ones available in the literature, whereby the accuracy and versatility of the proposed method are demonstrated.