刚塑性有限元模拟中体积变化问题处理方法

分析了在塑性有限元分析金属体积成形问题中引起体积变化的原因 ,指出在采用罚函数方法引入体积不可压缩条件的刚塑性有限元模拟中惩罚因子、增量步长、节点调整、网格密度和重划次数等是引起体积变化的主要原因。针对所分析原因 ,在现有的方法基础上提出基于三维有限元的节点调整和体积补偿新方法 ,并在实际模拟中应用 ,得到较好的效果 ,提高了模拟的精度和可信度。     

VISCOMETRIC FLOWS OF NEUTRALLY BUOYANT SUSPENSIONS IN SECOND-ORDER FLUIDS

The local volume averages of the equations of motion as well as the appropriate boundary conditions for neutrally buoyant suspensions of spheres in second-order fluids (Li and Slattery, 1989a) are employed to analyze both the cone-plate viscometer and the parallel-plate viscometer. The predicted results are compared with measurements of apparent viscosity and of the primary normal stress difference. Predictions of measured apparent viscosities are in excellent agreement with the experimental data up to 20% solids by volume; the average error is no more than 5%. Predictions of apparent viscosities for the limiting case of a very dilute suspension are restricted to less than 10% solids by volume. The primary normal stress difference is predicted to be independent of the volume fraction of solids, in approximate agreement with previously reported data.

Both theory and experiment have demonstrated that a neutrally buoyant suspension of uniform spheres in an incompressible second-order fluid behaves like a single-phase second-order fluid in these two viscometers.     

基于改进自由体积和无规因子贡献的聚合物热力学模型

给出了一个基于改进的自由体积和无规因子的聚合物溶液热力学模型，其无规因子项(NRF)用来描述聚合物溶液中聚合物与溶剂的局部交互作用对剩余Gibbs自由能的贡献，改进的自由体积项(Entropic-FV)用来描述聚合物溶液的混合效应对剩余Gibbs自由能的贡献。文中对45组二元聚合物与溶剂体系进行了关联计算，并与Zafarani-Moattar模型，Flory-Huggins模型和UNIFAC-FV模型的计算结果进行了比较，其平均绝对误差分别为1．910％，2．465％，3．738％和10．149％。     

Development of the meshless finite volume particle method with exact and efficient calculation of interparticle area

The Finite Volume Particle Method (FVPM) is a meshless method based on a definition of interparticle area which is closely analogous to cell face area in the classical finite volume method. In previous work, the interparticle area has been computed by numerical integration, which is a source of error and is extremely expensive. We show that if the particle weight or kernel function is defined as a discontinuous top-hat function, the particle interaction vectors may be evaluated exactly and efficiently. The new formulation reduces overall computational time by a factor between 6.4 and 8.2. In numerical experiments on a viscous flow with an analytical solution, the method converges under all conditions. Significantly, in contrast with standard FVPM and SPH, error depends on particle size but not on particle overlap (as long as the computational domain is completely covered by particles). The new method is shown to be superior to standard FVPM for shock tube flow and inviscid steady transonic flow. In benchmarking on a viscous multiphase flow application, FVPM with exact interparticle area is shown to be competitive with a mesh-based volume-of-fluid solver in terms of computational time required to resolve the structure of an interface.     

Dynamics of fluid circulation in coupled free and heterogeneous porous domains

Fluid flow in coupled free and porous domain, particularly when the porous medium is heterogeneous, is encountered in many hydro-environmental conditions, e.g., leakage from underground pipe, combined groundwater lake-subsurface interactions. One of the most difficult problems in the study of coupled flow behaviour has been the development of a universally applicable modelling scheme for combining the flow regimes. This is because the free/porous interfacial properties (e.g., shear-stress; velocity slip) that govern the coupled flow behaviour are difficult to determine experimentally under hydro-environmental conditions. On the other hand, the implications of various forms of heterogeneity in the porous media properties can be very different on the fluid-flow behaviour. Difficulties may also arise in direct coupling of the model equations that govern the fluid flow in the individual regions (e.g., Navier-Stokes for free-flow region and the Darcy's equation for the porous flow region). Consequently, models of coupled free and porous flow for hydro-environmental conditions are not very well developed at the moment. While there are some indications that fluids in coupled free and porous domains may circulate (i.e., development of flow cells), there is a lack of appropriate 3D analysis on how heterogeneities in porous media may affect such flow patterns. In this paper, we aim to analyse how porous media heterogeneity affects the dynamics of flow circulation in the porous side of a coupled free and porous domain. For this purpose, we analyse flow patterns in several model domains made up of two porous layers with differing permeabilities. The governing model equations are discretised and solved using the standard finite volume method on a staggered cell-centred mesh. The temporal discretisation is done using the explicit method. An in-house graphical user interface (GUI) has been created specifically to aid in the visualisation of otherwise complex flow patterns. The GUI contains many post-processing options and provides a comprehensive tool for the analysis of hydrodynamics and contaminant motion (not discussed in this paper) in coupled free and porous flow domains. This GUI is described in this paper briefly. The effects of altering the aspect ratio (i.e., multi-scale) of the domain on the coupled flow pattern have also been discussed.     

基于GF_1卫星的丹江口水库水面面积-蓄水量-水位相关性研究

基于GF_1卫星16m分辨率多光谱遥感数据及水文观测数据,通过计算NDWI水体指数来提取丹江口水库水面面积,并研究水面面积、水位与蓄水量三者之间的相关性,实现动态分析水库蓄水量变化情况,以期为"南水北调"提供决策支撑。研究发现:(1)南水北调中线首个调度年前后丹江口水库水面面积-水位,水位-蓄水量以及水面面积-蓄水量之间均为显著相关关系,且三对变量依次呈现明显的1次、2次、3次曲线关系,拟合曲线方程的决定系数分别达到0.976、0.999和0.980,因此可以根据其中一个已知变量预测其他两个未知量;(2)水面面积自2014年10月5日至2016年5月5日最大变化幅度相差180km~2,是影响周边环境变化的重要因素;(3)同时蓄水量在首个调度年内,于2014年11月初达到最大值204.10亿m~3,可调水量在16.53~78.93亿m~3范围内浮动,且年内整体变化平缓稳定,保证了2015调度年内南水北调中线工程的供水需求。     

无伸缩缝混凝土防渗衬砌渠道跳仓浇筑施工技术

现浇混凝土防渗渠道普遍采用设置伸缩缝的方法防止混凝土产生裂缝,而伸缩缝的止水易因老化和长草等遭受破坏,导致渠道渗漏。针对上述弊端提出了无伸缩缝混凝土防渗衬砌渠道的隔板跳仓浇筑施工技术。经过在多个大中型灌区节水项目中推广应用,总结出该施工技术的操作步骤和技术要点,包括基面处理措施,浇筑方法,操作要点,浇筑,振捣、抹面等的质量要求及适用范围。多年的运行实践表明,该技术可明显节约投资、缩短工期、延长使用寿命、有效改善防渗效果,且施工简便、工程质量易于保证等。     

导电复合材料的研究

以铝粉，乙炔炭黑，石墨等导电填料，不饱和聚酯树脂为基休一项玻璃纤维布为增强材料制成复合型导电复合。研究了导电填料不同种类，添加量对复合材料性能的影响。     

定容条件下火药实际燃烧规律的数值模拟

为研究火药的实际燃烧规律,建立了定容条件下火药燃烧的修正数学模型,分析了火药实际弧厚的分布及变化、点火不同步和燃速系数变化等因素对火药实际燃烧规律的影响。在此基础上,用考虑综合因素的修正数学模型对定容条件下火药的实际燃烧过程进行了模拟,计算值与实验结果具有较好的一致性,说明修正模型所建立的假设及处理方案是合理的。     

考虑土体强度参数深度依赖性的非平稳随机场边坡滑动深度和体积分析

基于随机有限元方法（RFEM）和Monte-Carlo模拟,建立了能够模拟具有深度依赖性的土体强度参数边坡或具有软弱夹层边坡的非平稳随机场模型。通过辐射扫描方法获取最危险滑动面,得到非平稳随机场边坡滑动深度和滑动体积。将该方法与均质边坡和平稳随机场边坡获得的滑动深度和滑动体积进行对比,研究土体强度参数深度依赖性的非平稳随机场边坡的不同坡度和各向异性对滑动深度和滑动体积的影响。结果表明：基于RFEM的非平稳随机场滑动深度和滑动体积结果为地质灾害风险及后果评估提供了一个参考分析的新角度;参数的深度依赖性程度越高(变化率b越大),则边坡平均滑动深度和滑动体积越小;在大坡度滑坡风险后果评估中,可以基于均质边坡分析结果进行初步预估;在非平稳随机场中,随着各向异性的增大,所有坡度边坡的滑动体积变异系数增大;不同坡度下,水平方向的相关程度越大,则土坡内软弱夹层的不确定性越强,滑动体积的变异系数也不断增大。