结构抗震计算不同方法之对比研究

利用有限元软件ANSYS对3种不同层高的钢框架模型进行同等烈度条件下抗震分析, 并进行底部剪力法、反应谱法和时程分析法的对比研究, 得到相应位移应力曲线、3种抗震计算方法各自的特点及适用条件. 通过观察时程曲线和震动动画过程, 可使学生对抗震计算有直观认识和深入理解.     

Flow and heat transfer over a generalized stretching/shrinking wall problem—Exact solutions of the Navier-Stokes equations

In this paper, we investigate the steady momentum and heat transfer of a viscous fluid flow over a stretching/shrinking sheet. Exact solutions are presented for the Navier-Stokes equations. The new solutions provide a more general formulation including the linearly stretching and shrinking wall problems as well as the asymptotic suction velocity profiles over a moving plate. Interesting non-linear phenomena are observed in the current results including both exponentially decaying solution and algebraically decaying solution, multiple solutions with infinite number of solutions for the flow field, and velocity overshoot. The energy equation ignoring viscous dissipation is solved exactly and the effects of the mass transfer parameter, the Prandtl number, and the wall stretching/shrinking strength on the temperature profiles and wall heat flux are also presented and discussed. The exact solution of this general flow configuration is a rare case for the Navier-Stokes equation.     

Fracture dynamics problem on mode I semi-infinite crack

By means of the theory of complex functions, fracture dynamics problems concerning mode I semi-infinite crack were studied. Analytical solutions of stress, displacement and dynamic stress intensity factor under the action of moving increasing loads Pt ³/x ³, Px ³/t ², respectively, are very easily obtained using the ways of self-similar functions. The correlative closed solutions are attained based on the Riemann–Hilbert problems.     

固定接触界面切向静弹性刚度问题研究

根据两球体单峰同时受法向、切向载荷时微滑切向应力的分布以及MB模型,给出了界面的总切向接触静弹性条件刚度、总条件法向载荷的解析解。将切向接触静弹性条件刚度的解析解嵌入到有限元软件中,获得了整机的理论模态。以一款八四七厂华中工学院XHK5140型自动换刀计算机数控立式镗铣床上的结合部为研究对象,通过实验对解析解进行了定量验证。研究结果表明：在理论振型与实验振型一致的条件下,界面模型的相对误差在-19．2％～16．8％之间。     

Energy density and release rate study of an elliptic-arc interface crack in piezoelectric solid under anti-plane shear

The anti-plane problem of an elliptical inhomogeneity with an interfacial crack in piezoelectric materials is investigated. The system is subjected to arbitrary singularity loads (point charge and anti-plane concentrated force) and remote anti-plane mechanical and in-plane electrical loads. Using the complex variable method, the explicit series form solutions for the complex potentials in the matrix and the inclusion regions are derived. The electroelastic field intensity factors, the corresponding energy release rates and the generalized strain energy density at the cracks tips are then provided. The influence of the aspect ratio of the ellipse, the crack geometry and the electromechanical coupling coefficient on the energy release rate and the strain energy density is discussed and shown in graphs. The results indicate that the energy release rate increases with increment of the aspect ratio of the ellipse and the influence of electromechanical coupling coefficient on the energy release rate is significant. The strain energy density decreases with increment of the aspect radio of the ellipse and it is always positive for the cases discussed. The energy release rate, however, can be negative when both mechanical and fields are applied.     

Convective boiling in a parallel microchannel heat sink with a diverging cross section and artificial nucleation sites

To develop a highly stable microchannel heat sink for boiling heat transfer, three types of diverging microchannels (Type 1, Type 2 and Type 3) were designed to experimentally investigate the effect of different distributions of artificial nucleation sites (ANS) on the enhancement of flow boiling heat transfer, in 10 parallel diverging microchannels with a mean hydraulic diameter of 120 μm. Water was used as the working fluid with mass flux, based on the mean cross section area, ranging from 99 to 297 kg/m² s. The Type-1 system did not contain any ANS; the Type-2 system contained ANS distributed uniformly along the downstream half of the channel; and the Type-3 system contained ANS distributed uniformly along the entire channel. The ANS are laser-etched pits on the bottom wall of the channel and have a mouth diameter of approximately 20-22 μm, as indicted by the heterogeneous nucleation theory. The results of the present study reveal that the presence of ANS for flow boiling in parallel diverging microchannels significantly reduces the wall superheat and enhances the boiling heat transfer performance. The Type-3 system shows the best boiling heat transfer performance.     

Development of a multi‐scale ocean model by using particle Laplacian method for anisotropic mass transfer

The direct injection of CO₂ into the deep ocean is one of the feasible ways for the mitigation of the global warming, although there is a concern about its environmental impact near the injection point. To minimize its biological impact, it is necessary to make CO₂ disperse as quickly as possible, and it is said that injection with a pipe towed by a moving ship is effective for this purpose. Because the injection ship moves over a spatial scale of O(10²km), a mesoscale model is necessary to analyse the dispersion of CO₂. At the same time, since it is important to investigate high CO₂ concentration near the injection point, a small‐scale model is also required. Therefore, in this study, a numerical model was developed to analyse CO₂ dispersion in the deep ocean by using a fixed mesoscale and a moving small‐scale grid systems, the latter of which is nested and moves in the former along the trajectory of the moving ship. To overcome the artificial diffusion of mass concentration at the interface of the two different grid systems and to keep its spatial accuracy almost the same as that in the small‐scale, a particle Laplacian method was adopted and newly modified for anisotropic diffusion in the ocean. Copyright © 2010 John Wiley & Sons, Ltd.     

A direct‐forcing pressure‐based lattice Boltzmann method for solving fluid–particle interaction problems

A direct‐forcing immersed boundary‐lattice Boltzmann method (IB–LBM) is developed to simulate fluid–particle interaction problems. This method uses the pressure‐based LBM to solve the incompressible flow field and the immersed boundary method to handle the fluid–particle interactions. The pressure‐based LBM uses the pressure distribution functions instead of the density distribution functions as the independent dynamic variables. The main idea is to explicitly eliminate the compressible effect due to the density fluctuation. In the IB method, a direct‐forcing method is introduced to capture the particle motion. It directly computes an IB force density at each lattice grid from the differences between the pressure distribution functions obtained by the LBM and the equilibrium pressure distribution functions computed from the particle velocity. By applying this direct‐forcing method, the IB–LBM becomes a purely LBM version. Also, by applying the Gauss theorem, the formulas for computing the force and the torque acting on the particle from the flows are derived from the volume integrals over the particle volume instead of from the surface integrals over the particle surface. The order of accuracy of the IB–LBM is demonstrated on the errors of velocity field, wall stress, and gradients of velocity and pressure. As a demonstration of the efficiency and capabilities of the new method, sedimentation of a large number of spherical particles in an enclosure is simulated. Copyright © 2010 John Wiley & Sons, Ltd.     

Enhanced multi‐objective optimization of a dimpled channel through evolutionary algorithms and multiple surrogate methods

Using a multi‐objective evolutionary algorithm (MOEA) and enhanced surrogate approximations, the present study demonstrates the numerical analysis and optimization of staggered‐dimple channels. Two surrogates, the response surface approximation (RSA) model and the Kriging (KRG) model, are applied in light of the surrogate fidelity of the approximate analysis. An enhanced Pareto‐optimal front is obtained by performing local resampling of the Pareto‐optimal front, which provides relatively more accurate Pareto‐optimal solutions in the design space for each surrogate model. Three dimensionless design variables are selected, which are related to geometric parameters, namely, the channel height, dimple print diameter, dimple spacing, and dimple depth. Two objective functions are selected that are related to the heat transfer and pressure loss, respectively. The objective‐function values are numerically evaluated through Reynolds‐averaged Navier–Stokes analysis at the design points that are selected through the Latin hypercube sampling method. Using these numerical simulations two surrogates, viz, the RSA and Kriging models, are constructed for each objective function and a hybrid MOEA is applied to obtain the Pareto‐optimal front. For the particular implementation of surrogate models, it is observed that Pareto‐optimal predictions of the RSA model are better than those of the KRG model, whereas the KRG model predicts equally well at the off‐Pareto‐region (region away from the Pareto‐optimal solutions), which is not the case with the RSA model. The local resampling of the Pareto‐optimal front increases the fidelity of the approximate solutions near the Pareto‐optimal region. The ratios of the channel height to the dimple print diameter and of the dimple print diameter to the dimple pitch are found to be more sensitive along the Pareto‐optimal front than the ratio of the dimple depth to the print diameter. The decrease of the ratio of the channel height to the dimple diameter and the increase of the ratio of the dimple print diameter to the pitch lead to greater heat transfer at the expense of the pressure loss, whereas the ratio of the dimple depth to the print diameter is rather insensitive to Pareto‐optimal solutions. Pareto‐optimal solutions at higher values of the Nusselt number are associated with higher values of the pressure loss due to the increased recirculation, mixing of fluid and vorticity generation. Copyright © 2010 John Wiley & Sons, Ltd.     

An efficient curvilinear non‐hydrostatic model for simulating surface water waves

An efficient curvilinear non‐hydrostatic free surface model is developed to simulate surface water waves in horizontally curved boundaries. The generalized curvilinear governing equations are solved by a fractional step method on a rectangular transformed domain. Of importance is to employ a higher order (either quadratic or cubic spline function) integral method for the top‐layer non‐hydrostatic pressure under a staggered grid framework. Model accuracy and efficiency, in terms of required vertical layers, are critically examined on a linear progressive wave case. The model is then applied to simulate waves propagating in a canal with variable widths, cnoidal wave runup around a circular cylinder, and three‐dimensional wave transformation in a circular channel. Overall the results show that the curvilinear non‐hydrostatic model using a few, e.g. 2–4, vertical layers is capable of simulating wave dispersion, diffraction, and reflection due to curved sidewalls. Copyright © 2010 John Wiley & Sons, Ltd.