堆石体的本构模型对其应力变形分析尤为重要.在热力学条件下,分析堆石坝筑坝过程中应力比的变化过程,推导出符合能量守恒定律的堆石料本构模型,给出吉布斯自由能函数指数的范围.在闭合应力路径下模型符合热力学定律,能够反映堆石料的非线性、压硬性等主要工程性质,参数少且独立、物理意义明确.以水布垭面板堆石坝为例,采用热力学本构模型进行了有限元计算,与实测资料相比,其计算结果比其他模型更加符合实际,进一步验证了模型的可靠性.

The researches on constitutive relationship for rockfills are significant because of its wide application. Based on the theory of thermodynamics, a four-parameter K-G model is derived from the Gibbs free energy function. The model is theoretically rigorous and can reflect the nonlinearity, pressure rigidity, partly anisotropy and other engineering properties. Most importantly, only four parameters which are independent of each other are needed in this model. The iso-stress ratio tests and numerical simulations of two-dimensional CFRD validate the reasonability of the model. Furthermore, three-dimensional static calculation of Shuibuya CFRD is performed. The results show that the vertical displacement is 2.02 m, which is closer to the measured value than that of other nonlinear models. Therefore, the four-parameter K-G model can be used to calculate the stress and deformation of CFRDs.

In order to study the dynamic and deformation characteristics and the stability of CFRD through earthquake,the seismic response of Pankou CFRD is studied by using finite element method,focusing on the acceleration response characteristics,the post-earthquake residual deformation and panel deformation,and the seismic safety of dam elements and downstream slope.The results show that the overall seismic performance of the dam is good enough to meet the seismic stability requirements under the given earthquake,and meanwhile,the crest and the downstream slope near the crest has a larger acceleration response and the relative large permanent deformation occurs there.

In the numerical calculation of geotechnical project, it’s difficult to confirm the parameters because of the complexity and the uncertainty of them as the time is changing. However, the back-analysis provides us an effective way. Based on the result of the triaxial test on rock-fill of Shui Bu Ya CFRD, the thesis adopts the direct back-analysis method which combines the BP Neural Network and Genetic Algorithm to calculate the Tsinghua non-linear K-G model parameters of the rock-fill. The back-analysis parameters are used to simulate the filling process of Shui Bu Ya CFRD and predict the displacement of the dam. The thesis provides a technical reference for displacement back-analysis of soil parameters for CFRD.

Face slab is as the water retaining deck of CFRDs,which is very crucial in the overall design of CFRDs.General speaking,the design of face slab includes plinth,face slab,perimeter joints,vertical joints and horizontal joints.The development of design of the high CFRD is based on experiences,however,the prevailing specification for CFRDs in China has no detailed guidelines for the dams over 200m.Based on the conclusion of the construction of CFRDs,including Shuibuya Dam(233 m) and Bakun Dam(205 m),the article illustrate the design of face slab and impervious measures for future design reference.

A vertical two dimensional model for thermo stratified flows is proposed,which solves the velocity and pressure by a fractional step method,and Large Eddy Simulation model is used to calculate the Reynolds stress and turbulent diffusion effect.The recirculation cavity flows caused by temperature difference are simulated,and the relations between the flow patterns and the Rayleigh number are analyzed.The density flows in a reservoir model are also calculated,and the calculated results are compared with others.The compared results showed that the calculation results of this model is good.

In this paper the influence of the growing elastic modulus, time dependent creep and temperature change on the expansion of magnesium oxide(MgO) concrete is studied using simulation analysis. The influence of curing temperature and temperature of concrete on the effective expansion and the compensation effects in different dam areas are also probed. Results show that the expansion of MgO at later age is good for compensating the shrinkage induced by temperature drop and the expansion of MgO can improve the stress status in the restricted zones but may worsen the stress status in other zones.

彭水水电站船闸开挖段地形地质条件较复杂,断层分布较多,为此应用弹塑性有限元分析方法,采用接触力学的直接约束法模拟断层的滑移和脱开等不连续行为,对高边坡开挖过程进行了数值仿真.通过研究对比开挖过程中岩体的变形、应力以及塑性区的大小和分布,来评价开挖过程中边坡的稳定性和支护措施的有效性,为高边坡开挖方案的优化和选定提供指导.

This paper presents a dynamic model of temperature rise caused by cementitious materials hydration based on the basic nature of chemical change. There are many models for computation of temperature rise caused by cementitious materials hydration, but the most takes this temperature rise as a function of concrete age only, which does not consider the impact of initial temperature and temperature rise during hydration and thus cannot reflect the real course of temperature rise. This model not only accords with the general law of chemical change, but also coincides well the experimental data and stimulates well the actual course of temperature rise of concrete under any initial temperature.

A direct constraint procedure was used to model the discontinuous mechanical behavior(such as sliding and disjointing) of faults and weak interlays in steep slopes to evaluate the safety of such slopes.The direct constraint procedure is the most robust procedure for contact mechanics.The elasto-plastic FEM was used to model the contact characteristics between the faults and the rock mass to analyze the excavation process for the Pengshui hydropower station.The results show that the magnitude and distribution of the deformation,the stresses and the plastic yield zones are all greatly increased by the sliding and disjoining,which increase the probability of landslides affecting the slope safety.

Numerical simulation analysis was performed on the exravation process of underground powerhouse of Shuihuya hydropower station based on three-dimension nonlinear finite element method. It adopted the MARC software, including the numerial simulation of excavation, nonlinear solution, finite element calculation, displacement analysis on the calculation results, stress analysis, and the analysis of yield area. The study was focused on the stress and deformation of cave and adjacent rock, the distribution of yield area,and the integral stability of the underground cavern group. The research results provides the basis for the construction of the project.

At present there are many models that calculates expansion of MgO, but none of them really reflects the effect of temperature on expansion. This paper analyzes their defects respectively, and deduces expansion model of dynamics for the first time. This model reveals the truth behind chemical changes and simulates experimental data quite well, it solves the problem of calculating MgO expansion enduring any temperature process exactly.

基于一种新型的接触算法--直接约束法,研究水口水电站混凝土坝和坝后式厂房在上下游水压力下的相互作用情况,给出了厂坝间相互传递的推力与联接缝宽度的关系,为厂房的安全性评价提供了依据.

The turbulent concentration fields in the wake of a circular cylinder were investigated by both numerical and experimental methods. Planar Laser-Induced Fluorescence (PLIF) technique was used to non-intrusively measure the simultaneous two-dimensional concentration fields, which is necessary for research on unsteady flows. The fractional step finite element formulation was extended to simulate the turbulent transportation processes in the wake region. This numerical formulation has third order accuracy and an enlarged stability domain which is suitable for convection dominated flows. The characteristics of the mean and instantaneous concentration distributions were analyzed in detail. The properties of the vortex movement and the concentration distributions were also analyzed. The numerical results were verified by the PLIF experiment.

The influence of gaps in hydrostructures on the safety of the hydrostructures can not be easily evaluated. A direct constraint procedure was developed to analyze the effects of gaps based on the kinematics of the system rather than interface elements. The procedure tracks the motion of the bodies and when contact occurs, places direct constraints on the motion using both kinematic constraints on transformed degrees of freedom and nodal forces. The method was used to analyze the effect of gaps between the dam and a workshop behind the Shuikou Hydropower Station to illustrate the advantages of the proposed procedure. The results show that the direct constraint procedure can be used in analyzing hydro contact force conveniently. 

在各项异性渗流模型的基础上,建立了地下水污染物质输移运动的二维数学预测模型,采用有限元法对渗流和污染物分布进行模拟,并且用解析解检验了数学模型的可靠性.对于污染物随地面降雨入渗的二维非恒定问题,计算结果与已有文献进行比较,证明了数值模拟的可信度.进一步实际模拟在非均匀介质情况下的污染物随地下水输移过程,分析了介质的透水特性对污染物浓度的影响.

The mechanical characteristic of the expanded metal and its simplified analysis model, as well as its using in hydrolic engineering are proposed. A plate model is introduced to simulate the expanded metal, which is used in the experiment of covered way's liner, with equivalent in stiffness. With the finite element analysis method, calculation results using the plate model are close to those of experiments finished before. Accordingly, the simulation of plate model is doable and it also provides a simple and credible analysis method for expanded metal.

In order to make the free vibration of design gate far from the main frequency scope of high-energy pulse of water flow to guarantee the safety of gate, on the basis of considering the impacts of the additional mass caused by liquid-solid coupling on the free vibration of gate, we adopt the three-dimensional boundary finite analysis method to gain the free vibration frequency and vibration pattern characteristics of gate. The calculation results on the fixed-wheel plane gate of Deqing gate showed that the impact of water head on gate is to lower the free vibration of gate and change its vibration patterns. In order to increase the free vibration frequency of gate, it is necessary to optimize the design of the structure of gate.

名誉理事长：夏武平　理事长：王祖望 副理事长：马逸清　张知彬　张亚平　胡锦矗　赵新全　徐宏发（按姓氏笔画为） 秘书长：魏辅文　副秘书长：王德华　李　明　刘定震 理事名单（以姓氏笔画为序）： 　　　　　丁　平　于孝臣　马　勇　马逸清　马建章　王　丁　王小明　王岐山 　　　　　王应祥　王金星　王祖望　王德华　王　勇　方盛国　冯祚建　宁振东 　　　　　卢立仁　卢浩泉　全国强　刘定震　刘季科　孙儒泳　李进华　李保国 　　　　　宋延龄　汪　松　汪诚信　张亚平　张　洁　张明海　张树义　张知彬 　　　　　张恩迪　苏彦捷　吴　毅　郑昌琳　杨　光　杨奇森　周开亚　赵新全 　　　　　施大钊　胡锦矗　钟文勤　高行宜　贾竞波　徐宏发　黄乘明　盛和林 　　　　　蒋志刚　蒋学龙　彭基泰　温得启　鲍毅新　樊乃昌　魏万红　魏辅文 　　*台湾（2名），香港（1名）和澳门（1名）理事暂时空缺，以后增补

This paper presents a numerical model of 3D seepage flow for non-uniform medium. The governing equations of saturated-unsaturated flow are introduced and the finite element method is applied to discretize them in space. The techniques of determining the seepage faces and computing the discharge rate through a specified section are proposed. Therefore, this model describes processes in relatively smaller space and time units. It has been tested and implemented in the research of complex 3D seepage flow. The results reveal that both the seepage flow rate and the level of groundwater are affected by the scale of impervious wall.