Practical Equivalent Continuum Model for Simulation of Jointed Rock Mass Using FLAC3D

A simple practical equivalent continuum numerical model for simulating the behavior of jointed rock mass has been extended to three-dimensional using FLAC3D. This model estimates the properties of jointed rock mass from the properties of intact rock and a joint factor (Jf), which is the integration of the properties of joints to take care of the effects of frequency, orientation, and strength of joint. A new FISH function has been written in FLAC3D specifically for modeling jointed rocks using the Duncan and Chang hyperbolic model. This model has been validated first with simple element tests at different confining pressures for different rocks with different joint configurations. Explicit modeling of the joints has also been done in element tests and results obtained compare well with the results of equivalent continuum model and also with experimental results. Further, this has been applied for a case study of a large underground power house cavern in the Himalayas. The analysis has been done under various stages of excavation, assigning a null model available in FLAC3D for simulating the excavation.     

Two-Dimensional Physical and Numerical Modeling of a Pile-Supported Earth Platform over Soft Soil

This paper focuses on the mechanisms occurring in a granular earth platform over soft ground improved by rigid piles. Two-dimensional physical model experiments were performed using the Schneebeli’s analogical soil to investigate the load transfer mechanisms by arching and the settlement reduction and homogenization. Experimental outputs are compared to results obtained on a numerical model using a plane strain continuum approach. The impact of the constitutive model complexity to simulate the platform material behavior was first assessed, but no large difference was recorded. As far as the proposed model, which takes the main features of the observed behavior satisfactorily into account, the numerical procedure could be validated and the parametric studies extended numerically. Both approaches of this study underlined the main geometrical and geotechnical parameters which should inevitably be taken into account in a simplified design method, namely the capping ratio, the platform height, and the platform material shear strength.     

Applications of the discrete element method in mechanical engineering

Compared to other fields of engineering, in mechanical engineering, the Discrete Element Method (DEM) is not yet a well known method. Nevertheless, there is a variety of simulation problems where the method has obvious advantages due to its meshless nature. For problems where several free bodies can collide and break after having been largely deformed, the DEM is the method of choice. Neighborhood search and collision detection between bodies as well as the separation of large solids into smaller particles are naturally incorporated in the method. The main DEM algorithm consists of a relatively simple loop that basically contains the three substeps contact detection, force computation and integration. However, there exists a large variety of different algorithms to choose the substeps to compose the optimal method for a given problem. In this contribution, we describe the dynamics of particle systems together with appropriate numerical integration schemes and give an overview over different types of particle interactions that can be composed to adapt the method to fit to a given simulation problem. Surface triangulations are used to model complicated, non-convex bodies in contact with particle systems. The capabilities of the method are finally demonstrated by means of application examples. Commemorative Contribution.     

异质谷间转移电子器件的计算机模拟

使用弛豫时间近似求出了描述两能谷间电子交换的方程。在此基础上建立了双能谷电子的连续性方程和油松方程。对具体的ＧａＡｓ／ＡｌｘＧａ１－ｘＡｓ异质谷间转移电子器件结构求解了这些方程组，得到了直流工作时器件内的电场分布、双能谷中载流子的布居以及两个能谷的电流分布。这些结果正好和常规Ｇｕｎｎ器件相反，说明了两种器件的不同工作机理。最后通过对低Ａｌ组份势垒结构的计算说明了Ｘ谷电子注入对器件工作的重要作用。     

熔窑玻璃液流数值模拟的新发展

熔窑玻璃液流数值模拟的新发展赵国昌（清华大学核能技术设计研究院，１０００８４）ＮｅｗＤｅｖｅｌｏｐｍｅｎｔｏｆｔｈｅＮｕｍｅｒｉｃａｌＳｉｍｕｌａｔｉｏｎｆｏｒＦｌｏｗａｎｄＨｅａｔＴｒａｎｓｆｅｒｏｆＧｌａｓｓＭｅｌｔｉｎＴａｎｋＦｕｒｎａｃｅ￥Ｚ...     

CNC中插补和轴伺服的微机仿真

本文级出了在ＰＣ机上对一个自己开发ＣＮＣ系统插祉和轴伺服控制进行仿真的结果。它包含按插补原理对直线、圆和渐开线的插补以及对轴伺服控制中的位置环和前馈控制的仿真。轴伺服马达通过一个二阶连续系统模型来仿真，最后得出对伺服进行前馈控制补偿可显著改善系统的输出。     

Comparison of Displacement and Mixing Diffusers

In order to cumpare the peformance of different supply diffuers of ventilation air, the airflow passern, temperature stratifiation and contaminant dispersion in a furnitured office ventilated by three kinds of air diffuer were numerically investigated. The air diffuers studied in this paper are a quarter-cylinder displacement diffuer on the floor and mixing diffuers (linear and vortex diffuers) on the ceiling. The heat sources in the of-fice are considered to be 50% convective and 50% radiative. The k-? two-equatwn model of turbulence is employed to predict the turbulent diffusion. The results show that the displacement diffuser provides a rather uniform flow field with low velocify in most areas, and the vertical temperature difference from floor to ceiling is as high as 6 K. With the linear diffuser, the air velociry is high, and the temperature is uniform both horizontally and vertically. The air velocity generated by the vortex diffuser is moderate. The distributions of the temperature and the contaminant are rather uniform.     

A numerical model for the dissolution of spherical particles in binary alloys under mixed mode control

A general numerical model is described for the dissolution kinetics of spherical particles in binary systems for any combination of first order reactions at the particle-matrix interface and long distance diffusion within the matrix. The model is applicable to both finite and infinite media and handles both complete and partial particle dissolution. It is shown that interfacial reactions can have a strong effect on the dissolution kinetics, the solute concentration at the particle-matrix interface and the solute concentration profile in the matrix.     

Numerical study on PEMFC’s geometrical parameters under different humidifying conditions

Steady-state and three-dimensional simulations were carried out to study the influences of geometrical parameters on the performance of PEMFC under different hydrating conditions. Flow fields, species transport, transport of water in polymer membrane and movement of liquid water in cathode and anode porous layers were determined, in order to accomplish a complete estimation of ohmic and concentration losses of PEMFC power. The geometrical parameters were thickness of the polymer membrane, cathode catalyst layer as well as gas channel to rib width ratio. Every simulation was made under different relative humidities of inlet flows (50 and 100%) for every change of characteristic length. Results show that the influence of the geometrical parameters on ohmic and concentration losses is of considerable importance. The performance of PEMFC is seriously affected under dehydrating conditions. However, such performance may be considerably improved by using suitable geometrical parameters. Cathode and anode liquid saturation may not only affect the transport of species, but also the polymer electrolyte water content. These results show the importance of simultaneously calculating both the water absorption and desorption through the polymer electrolyte and the liquid saturation in the cathode and anode porous mediums to obtain an actual view of ohmic and concentration losses of the PEMFC performance.     

A polydispersed particle system representation of the porosity for non-saturated cementitious materials

In this paper, the porosity of cementitious materials is described in terms of pore size distribution by means of a 3-dimensional overlapping sphere system with polydispersivity in size. On the basis of results established by Lu and Torquato [B. Lu, S. Torquato, Nearest-surface distribution functions for polydispersed particle systems, Phys. Rev. A 45(8) (1992) 5530-5544] and Torquato [S. Torquato, Random Heterogeneous Media: Microstructure and Macroscopic Properties. Springer-Verlag: New York, 2001] providing relations for nearest-neighbor distribution functions, the volume fraction of pores having a radius larger than a prescribed value is explicitly expressed. By adopting an appropriate size distribution function for the sphere system, it is shown that the pore size distribution of cementitious materials as detected for instance by mercury intrusion porosimetry (MIP), which generally points out several pore classes, can be well approached. On the basis of this porosity representation, the evaluation of the capillary pressure in function of the saturation degree is provided. The model is then applied to the simulation of the saturation degree versus relative humidity adsorption curves. The impact of the pore size distribution, the temperature and the thickness of the adsorbed water layer on these parameters are assessed and analyzed for three model materials having different pore characteristics.