水平岩层软弱基座危岩体演化过程研究

以重庆市甄子岩29号危岩体为研究对象，基于现场调查和分析解剖，建立UDEC离散元数值计算模型，对其演化变形过程进行模拟分析。研究表明，在巨大的自重作用下，由于基座岩体岩质软，岩体较破碎，且基座外侧为直立的陡崖，具有完整的临空面，从而使基座岩体易发生压缩流变及剪切流变，进而使危岩体存在滑移垮塌的危险。     

Questioning numerical integration methods for microsphere (and microplane) constitutive equations

In the last few years, more and more complex microsphere models have been proposed to predict the mechanical response of various polymers. Similarly than for microplane models, they consist in deriving a one-dimensional force vs. stretch equation and to integrate it over the unit sphere to obtain a three-dimensional constitutive equation. In this context, the focus of authors is laid on the physics of the one-dimensional relationship, but in most of the case the influence of the integration method on the prediction is not investigated.Here we compare three numerical integration schemes: a classical Gaussian scheme, a method based on a regular geometric meshing of the sphere, and an approach based on spherical harmonics. Depending on the method, the number of integration points may vary from 4 to 983,040! Considering simple quantities, i.e. principal (large) strain invariants, it is shown that the integration method must be carefully chosen. Depending on the quantities retained to described the one-dimensional equation and the required error, the performances of the three methods are discussed. Consequences on stress–strain prediction are illustrated with a directional version of the classical Mooney–Rivlin hyperelastic model. Finally, the paper closes with some advices for the development of new microsphere constitutive equations.     

亲子乐园空调送风方式的数值模拟研究

大型综合室内亲子乐园属于高大空间,设有游乐设施和游戏的特殊性使得对空间的舒适性要求一致,但是送风气流遇阻严重,室内存在较多气流死角,影响室内空气质量和儿童健康。因此其空调设计不仅需要考虑温度、风速的空间均匀度,还要考虑各点的空气龄和PMV-PPD指标。以天津某亲子乐园为研究对象,利用scSTREAM软件对适用于该房间的辐射供冷加新风、置换通风、混合通风三种空调方式的送风效果进行数值模拟分析,从流场的均匀性、人员的热舒适性等方面对模拟结果进行探讨,研究结果表明辐射供冷加新风方式的空间均匀性和PMV指标最佳,混合通风方式的空气龄最小。     

Static and free vibration analysis of functionally graded plates based on a new quasi-3D and 2D shear deformation theories

In this study, two dimensional (2D) and quasi three-dimensional (quasi-3D) shear deformation theories are presented for static and free vibration analysis of single-layer functionally graded (FG) plates using a new hyperbolic shape function. The material of the plate is inhomogeneous and the material properties assumed to vary continuously in the thickness direction by three different distributions; power-law, exponential and Mori–Tanaka model, in terms of the volume fractions of the constituents. The fundamental governing equations which take into account the effects of both transverse shear and normal stresses are derived through the Hamilton's principle. The closed form solutions are obtained by using Navier technique and then fundamental frequencies are found by solving the results of eigenvalue problems. In-plane stress components have been obtained by the constitutive equations of composite plates. The transverse stress components have been obtained by integrating the three-dimensional stress equilibrium equations in the thickness direction of the plate. The accuracy of the present method is demonstrated by comparisons with the different 2D, 3D and quasi-3D solutions available in the literature.     

Modeling the effect of deep traps on the capacitance–voltage characteristics of p-type Si-doped GaAs Schottky diodes grown on high index GaAs substrates

Numerical simulation, using SILVACO-TCAD, is carried out to explain experimentally observed effects of different types of deep levels on the capacitance–voltage characteristics of p-type Si-doped GaAs Schottky diodes grown on high index GaAs substrates. Two diodes were grown on (311)A and (211)A oriented GaAs substrates using Molecular Beam Epitaxy (MBE). Although, deep levels were observed in both structures, the measured capacitance–voltage characteristics show a negative differential capacitance (NDC) for the (311)A diodes, while the (211)A devices display a usual behaviour. The NDC is related to the nature and spatial distribution of the deep levels, which are characterized by the Deep Level Transient Spectroscopy (DLTS) technique. In the (311)A structure only majority deep levels (hole traps) were observed while both majority and minority deep levels were present in the (211)A diodes. The simulation, which calculates the capacitance–voltage characteristics in the absence and presence of different types of deep levels, agrees well with the experimentally observed behaviour.     

A new metaheuristic for optimization: Optics inspired optimization (OIO)

Due to the law of reflection, a concave reflecting surface/mirror causes the incident light rays to converge and a convex surface/mirror causes the light rays to reflect away so that they all appear to be diverging. These converging and diverging behaviors cause that the curved mirrors show different image types depending on the distance between the object and the mirror. We model such optical phenomena metaphorically into the searching process of numerical optimization by a new algorithm called optics inspired optimization (OIO). OIO treats the surface of the numerical function to be optimized as a reflecting surface in which each peak is assumed to reflect as a convex mirror and each valley to reflect as a concave one. Each individual is assumed to be an artificial object (or light point) that its artificially glittered ray is reflected back by the function surface, given that the surface is convex or concave, and the artificial image is formed (a candidate solution is generated within the search domain) based on the mirror equations adopted from physics of optics. Besides OIO, we introduce different variants of it, called ROIO (Rotation based OIO), and COIO (Convex combination based OIO) algorithms and conduct an extensive computational effort to find out the merit of the new algorithms. Our comparisons on benchmark test functions and a real world engineering design application (i.e., optimization of a centrifuge pump) demonstrate that the new algorithms are efficient and compete better than or similar to most of state of the art optimization algorithms with the advantage of accepting few input parameters.     

Development and validation of a micro-fin tubes evaporator model using R134a and R1234yf as working fluids

This paper presents a model of shell and tube evaporator with micro-fin tubes using R1234yf and R134a. The model developed for this evaporator uses the ε-NTU method to predict the evaporating pressure, the refrigerant outlet enthalpy and the outlet temperature of the secondary fluid. The model accuracy is evaluated using different two-phase flow boiling correlations for micro-fin tubes and comparing predicted and experimental data. The experimental tests were carried out for a wide range of operating conditions using R134a and R1234yf as working fluids. The predicted parameter with maximum deviations, between the predicted and experimental data, is the evaporating pressure. The correlation of Akhavan– Behabadi et al. was used to predict flow boiling heat transfer, with an error on cooling capacity prediction below 5%. Simulations, carried out with this validated model, show that the overall heat transfer coefficient of R1234yf has a maximum decrease of 10% compared with R134a.     

Neutral axis depth and moment redistribution in FRP and steel reinforced concrete continuous beams

The neutral axis depth is considered the best parameter for quantifying the moment redistribution in continuous concrete beams, as exemplified in various design codes worldwide. It is therefore important to well understand the variation of neutral axis depth against moment redistribution. This paper describes a theoretical investigation into the neutral axis depth and moment redistribution in concrete beams reinforced with fibre reinforced polymer (FRP) and steel bars. A finite element model has been developed. The model predictions are in favourable agreement with experimental results. Three types of reinforcement are considered, namely, glass fibre, carbon fibre and steel. Various levels of reinforcement ratio are used for a parametric evaluation. The results indicate that FRP reinforced concrete continuous beams exhibit significantly different response characteristics regarding the moment redistribution and variation of neutral axis depth from those of steel reinforced ones. In addition, it is found that the code recommendations are generally unsafe for calculating the permissible moment redistribution in FRP reinforced concrete beams, but the neglect of redistribution in such beams may be overconservative.     

Collapse load of composite laminates: Lower bound evaluation by stress field analytical approximation

The paper proposes a limit analysis approach to define the ultimate load capacity of orthotropic composite laminates under biaxial loading and plane stress conditions. A lower bound to the collapse load multiplier is computed by solving a maximization nonlinear problem, according to the static theorem of limit analysis. To set up the optimization problem a stress field distribution is hypothesized at lamina level, moreover inter-lamina stresses are also considered. The effectiveness and validity of the proposed approach is shown by comparing the obtained numerical predictions both with available experimental data and with other numerical results carried out by means of a different numerical lower bound approach.