一种孔隙介质力学模型及在水泥基材料的应用

基于细观力学和断裂力学的基本理论提出一个新的分析模型, 对孔隙介质的力学性能进行了分析. 依据孔隙介质内部孔隙的几何描述和状态参数,如孔隙率、形状、尺度及分布等,通过等效夹杂理论获得孔隙介质的等效本构方程,其最终变量为应力、应变和孔隙的形态参数. 根据断裂理论中材料承受载荷作用下破坏增长过程中的能量守恒,对孔隙介质变形过程中机械能、弹性应变能和载荷提供的势能进行分析, 根据能量守恒定律建立能量守恒方程,其最终变量也为应力、应变和孔隙的形态参数. 根据等效本构方程和能量守恒方程,获得孔隙介质承受载荷作用下的应力应变关系. 最后将该力学模型应用于水泥基材料,计算水泥基材料的力学性能并与文献中的结果进行对比分析,结果显示模型的计算结果准确有效.      

随机孔隙缺陷对3D编织复合材料力学性能的影响研究

3D编织复合材料经RTM工艺成型将难以避免产生孔隙缺陷,孔隙缺陷主要包括纤维束内干斑和富树脂区孔隙.经分析RTM工艺孔隙缺陷的产生机理,提出了随机孔隙缺陷分布模型,建立了基于细观胞元的含孔隙缺陷有限元分析模型,并应用概率统计方法,分析了缺陷的随机分布和孔隙率对材料等效弹性性能的影响.可为3D编织复合材料力学性能的设计提供依据.     

周期性吸声多孔材料微结构优化设计

多孔材料的吸声性能与材料孔隙率以及材料微结构几何构型存在密切相关.本文采用有限元方法研究了材料微观结构与宏观声学性能参数之间的关系,分析了通孔材料微结构的开孔形状、孔隙率以及孔隙尺寸对材料吸声性能的影响,并建立了在特定频率下具有高声吸收性能的通孔材料微结构几何构型的设计理论和方法,得到了具有较高声能吸收率的多孔材料微结构构型.     

PREDICTIONS OF EFFECTIVE OUT-PLANE SHEAR MODULUS AND SIZE EFFECT OF HEXAGONAL HONEYCOMB

给出了预测六边形蜂窝材料等效剪切模量及其尺寸效应的圆筒扭转力学模型和扭转能量法,建立了等效面外剪切模量G₁₃相对于材料体分比ν、周向单胞数n、圆筒半径r和单胞层数参数m变化的解析表达式;同时将扭转能量法、有限元数值模拟计算和G-A经典细观力学方法进行了比较,从理论上揭示并验证了尺寸效应的存在性. 结果表明,当蜂窝体胞尺寸相对结构尺寸无穷小时,预测结果趋近于细观力学方法的结果. 此外,利用周期性蜂窝材料的结构对称特性,使用体胞子结构有限元计算模型进行等效面外剪切模量及其尺寸效应的预测,在不影响计算结果的前提下极大地提高了计算效率.     

不同孔隙率的碾压混凝土力学性能试验研究

混凝土的破坏受其内部缺陷影响,为了研究孔隙含量对碾压混凝土力学性能的影响,探讨碾压混凝土的力学性能和尺寸效应随孔隙率的变化规律,采用预制孔隙法模拟碾压混凝土中的孔隙开展力学性能试验,分析3种尺寸C20碾压混凝土试件在4种不同孔隙率情况下的抗压强度和劈拉强度的变化情况。结果表明：同一尺寸下,C20碾压混凝土的抗压强度和劈拉强度随着孔隙率增加逐渐降低,碾压混凝土孔隙含量增加对其强度产生了负面影响;碾压混凝土拥有明显的尺寸效应,孔隙率越大,碾压混凝土强度的尺寸效应越显著。     

周期性多孔材料等效剪切模量与尺寸效应研究

提出了以圆筒扭转力学模型为基础, 预测周期性多孔材料等效剪切模量及其研究尺寸效 应的一种简单有效计算方法. 以方形孔和圆形孔两种典型多孔材料为例进行了数值计算求解; 同时, 建立了几何参数随体胞尺寸缩放因子$n$的解析关系, 证明了两种构型的周期性多孔材 料的等效剪切模量均随尺寸缩放因子$n$的增大而减小. 当$n \to \infty $时, 即体胞尺寸相对整体结构无限小时, 多孔材料的等效剪切模量趋近收敛于一个恒定值; 当体胞的材料体分比增大时, 多孔材料等效剪切模量也随之增大. 此外, 依据周 期性多孔材料的结构对称特性, 使用体胞子结构有限元计算模型进行等效剪切模量及其尺寸 效应的预测, 极大地提高了计算效率.     

试样尺寸对核电厂主蒸汽管道断裂韧性的影响研究

以某核电厂SA335 材料主蒸汽管道为研究对象,首先结合SA335 的断裂阻力曲线（J-R 曲线）试验测量结果,提出了一种方法确定韧性金属材料裂纹扩展模型（G-T 模型）中的主要微观参数：初始孔隙率0 f 和初始孔隙间距D;随后,在有限元计算中引入G-T 模型模拟了SA335 紧凑拉伸试样的断裂过程,讨论了试样尺寸对于J-R 曲线的影响．结果表明：试样厚度一致时,初始裂纹长度大的试样对应偏低的断裂韧性;当试样尺寸整体缩放时,较大尺寸的试样对应偏低的断裂韧性．最后,结合实例说明了试样整体尺寸对于主蒸汽管道临界裂纹长度的影响．     

双相钢板料的单向拉伸断裂失效研究(Ⅱ)——弧长法非线性有限元分析

为了获得docol800DP双相钢板料的等效塑性应变和应力三轴度在单向拉伸断裂失效过程中的变化历程,对其拉伸过程进行了弧长法非线性有限元分析。根据试验得到的材料参数建立了双相钢单向拉伸试件的有限元模型,考虑大变形引起的几何非线性和塑性强化引起的材料非线性,模拟了试件拉伸变形的全过程。计算结果表明:弧长法可较好地模拟试件变形的全程,计算得到的典型变形阶段、集中性失稳带的分布及方向与试验结果吻合很好;发生失稳变形时的应变与理论分析结果一致。最后,给出了起裂点处失效参数的变化历程,为定量化研究双相钢材料的断裂失效模型提供准确的数据支持。     

一类多孔固体的等效偶应力动力学梁模型

一维多孔固体结构可采用等效连续介质梁模型来研究其动力学行为. 当类梁结构的高度尺寸和多孔固体单胞结构尺寸相近时,等效模型的力学行为会产生尺寸效应现象. 等效经典模型由于不包含尺度参数而无法描述尺寸相关特点,而广义连续介质力学模型则可以准确地考虑尺寸效应的影响. 基于偶应力理论,对一类单胞含有圆形孔洞的周期性多孔固体类梁结构,给出了分析其横向自由振动的等效连续介质铁木辛柯梁模型. 通过对单胞分析,在应变能等价和几何平均的意义下,定义了等效偶应力介质的材料常数. 利用已有的材料常数,推导了等效铁木辛柯梁的动力学微分方程. 将实际多孔固体结构进行完全的动力学有限元离散计算,所获得的解作为精确解以检验等效梁模型所获得的频率和振型的精度. 振型的比较借助于模态置信准则矩阵方法. 大量算例表明,等效偶应力铁木辛柯梁模型在频率和振型两方面均具有较高的计算精度. 重点研究了单胞孔径的相对大小、类梁结构高度与单胞尺寸比以及类梁结构长高比对等效梁模型精度的影响. 在此基础上,偏保守地建议了多孔固体类梁结构自振分析方法.      

基于自适应网格的材料渗透系数设计

采用基于单元(结点)密度为设计变量进行结构和材料的拓扑优化设计时,有限元网格的密度对优化设计有很大影响. 在以渗透系数为目标进行材料微结构设计时,为了较好地描述单胞中的流固边界,需要将单胞划分为很小的网格,进一步增加了有限元计算和优化分析的规模. 为了降低计算规模, 研究了基于自适应网格的逆均匀化方法,以最大化各向同性等效渗透系数为目标,进行材料微结构设计. 优化迭代过程中,对单胞中流固界面处的网格进行自适应加密,降低优化问题的计算规模. 采用这一算法,对不同初始密度分布得到的单胞优化结果虽然不同,但具有相同的材料微结构,一定程度上说明了该方法的有效性.      

THEORETICAL MODEL OF EFFECTIVE STRESS COEFFICIENT FOR ROCK/SOIL-LIKE POROUS MATERIALS

Physical mechanisms and influencing factors on the effective stress coefficient for rock/soil-like porous materials are investigated, based on which equivalent connectivity index is proposed. The equivalent connectivity index, relying on the meso-scale structure of porous material and the property of liquid, denotes the connectivity of pores in Representative Element Area （REA）. If the conductivity of the porous material is anisotropic, the equivalent connectivity index is a second order tensor. Based on the basic theories of continuous mechanics and tensor analysis, relationship between area porosity and volumetric porosity of porous materials is deduced. Then a generalized expression, describing the relation between effective stress coefficient tensor and equivalent connectivity tensor of pores, is proposed, and the expression can be applied to isotropic media and also to anisotropic materials. Furthermore, evolution of porosity and equivalent connectivity index of the pore are studied in the strain space, and the method to determine the corresponding functions in expressions above is proposed using genetic algorithm and genetic programming. Two applications show that the results obtained by the method in this paper perfectly agree with the test data. This paper provides an important theoretical support to the coupled hydro-mechanical research.     

Tracer transport in a fractured chalk: X-ray CT characterization and digital-image-based (DIB) simulation

A digital-image-based simulation methodology is applied to evaluate the influence of heterogeneous porosity on the evolution of tracer concentrations in imaged tracer tests. Maps of computed tomography (CT)-number are calibrated relative to average porosity, and then thresholded to define porosity maps. These data are then used to automate the distribution of parameters within a finite element representation of the geometry. The technique is applied to characterize the variability of the porosity, the hydraulic conductivity, and the diffusivity for an artificially fractured chalk core (30 × 5 cm). X-ray CT was used both to characterize the initial condition of the core, and then to concurrently monitor the transport of an NaI tracer within the fracture and into the surrounding matrix. The X-ray CT imaging is used to characterize the heterogeneous rock porosity, based on which the hydraulic conductivity, and diffusivity of the chalk were defined and were directly imported into our newly developed three-dimensional FEMLAB-based multiple physics simulator. Numerical simulations have confirmed the observed tracer transport behaviors: (1) The different tracer-penetration distances imaged in the matrix above and below the horizontal fracture are indicative of a greater tracer mass penetrating into the lower matrix; and (2) Transport in the matrix below the fracture was enhanced. The computer simulated tracer concentration distributions compare favorably with those monitored by X-ray CT.     

Finite element modeling of porous titanium

The porous microstructures of metallic foams cause microscopic stress and strain localization under deformation which reduces the damage tolerance and therefore limits application of the materials. In this paper, the deformation of a relatively low porosity porous titanium is examined using two-dimensional (2D) plane strain and three-dimensional (3D) finite element models to identify the accuracy and limitations of such simulations. To generate the finite element models, a simulated microstructure was created based on micrographs of an experimental material. Compared to the 2D models, the 3D models require smaller model size to obtain convergent results. The macroscopic responses predicted by the 3D models are in reasonable agreement with experimental results while the 2D models underestimated the response. In addition, 3D models predicted more uniform microscopic field variable distributions. 2D models predicted higher probability of Von Mises stress and equivalent plastic strain exceeding a certain value and therefore overestimate the failure probability of the material.     

Pore-Scale Simulation of Shear Thinning Fluid Flow Using Lattice Boltzmann Method

The present work attempts to identify the roles of flow and geometric variables on the scaling factor which is a necessary parameter for modeling the apparent viscosity of non-Newtonian fluid in porous media. While idealizing the porous media microstructure as arrays of circular and square cylinders, the present study uses multi-relaxation time lattice Boltzmann method to conduct pore-scale simulation of shear thinning non-Newtonian fluid flow. Variation in the size and inclusion ratio of the solid cylinders generates wide range of porous media with varying porosity and permeability. The present study also used stochastic reconstruction technique to generate realistic, random porous microstructures. For each case, pore-scale fluid flow simulation enables the calculation of equivalent viscosity based on the computed shear rate within the pores. It is observed that the scaling factor has strong dependence on porosity, permeability, tortuosity and the percolation threshold, while approaching the maximum value at the percolation threshold porosity. The present investigation quantifies and proposes meaningful correlations between the scaling factor and the macroscopic properties of the porous media.     

An analytical homogenization method for heterogeneous porous materials

The objective of this work is to develop an analytical homogenization method to estimate the effective mechanical properties of fluid-filled porous media with periodic microstructure. The method is based on the equivalent inclusion concept of homogenization applied earlier for solid–solid mixture. It is assumed that porous media are described by the poroelastic constitutive law developed by Biot where porosity is a material parameter. By solving the governing equations of poroelasticity in Fourier transformed domain, the relation between periodic strain and eigenstrain in porous media is established. This relation is subsequently used in an average consistency condition involving both solid and fluid phase stresses and strains. The geometry of the porous microstructure is captured in the g-integral. This homogenization method can also be applied to estimate the equivalent properties of solid–fluid mixture where a pure solid and fluid can be modeled by assuming very low and high porosity, respectively. Several examples are considered to establish this new method by comparing with other existing analytical and numerical methods of homogenization. As an application, poroelastic properties of cortical bone fibril are estimated and compared with previously computed values.     

Steady water flow through unsaturated aggregated porous materials

The three-dimensional steady water flow through unsaturated aggregated porous materials composed of simple cubic open packed and tetrahedral close packed assemblies of uniform porous spheres is investigated with electric analogues and numerical computations. Water around the spheres is considered to be discontinuous with flow restricted through isolated annular water lenses held by surface tension forces around the contact points between spheres. It is found that the conductance of individual spheres depends only on the size of the water lenses and is independent of the radius of the sphere. For a simple cubic packing the conductance for small lens radii is given by Weber’s formula for flow from an electrified disc into an infinite medium. It follows that the bulk hydraulic conductivity of these assemblies of porous spheres is also independent of aggregate size over a range of water contents. This independence is also shown in measurements of hydraulic conductivity of aggregates of diatomaceous earth that show a convergence to a single relationship between conductivity and water content when there is no longer continuity of water in the macropore space. The effect of the three-dimensional flow through aggregates on solute leaching is demonstrated by considering the numerical results of the stream-tube pattern in a sphere.     

The Effect of Pore Structure on the Electrical Conductivity of Ti

Open-pore Ti foam samples with porosity in the range of 10–70% and average pore size of 150–400 μm was fabricated by powder metallurgy method using polymethyl methacrylate (PMMA) as space holder initially. The resulting foam is anisotropic: the pores are spheroidal, being shorter along the pressing direction than in the pressing plane. The pore anisotropy decreases as the size of the polymethyl methacrylate (PMMA) particles used increases and hence with pore size, which leads to a higher conductivity in the plane of the pressing. As the porosity increases, the conductivity of porous Ti decreases dramatically. The porosity e{\varepsilon} dependence of the electrical conductivity σ could be well described by Maxwell approximation, while the differential effective medium approximation is only suitable to porous Ti with finite size of 400 μm in the porosity range of 40–70%, i.e., high porosity metal with randomly oriented spheroids.     

含界面相Voronoi单元的等效弹性常数的计算

采用含界面相Voronoi单元有限元法,根据广义胡克定律,计算了在给定边界条件下,颗粒增强复合材料的等效弹性常数。建立了含多个随机分布的椭圆形夹杂及界面相的VCFEM模型,分析了夹杂体分比,界面相厚度和界面相弹性模量等因素对颗粒增强复合材料等效弹性常数的影响,并利用普通有限元方法对比验证。结果表明,当界面相弹性模量小于基体与夹杂时,材料的等效弹性模量会随着界面相厚度的增大而减小,随着夹杂体分比的增大而减小,并且界面过薄时,材料的等效弹性模量会随着夹杂体分比的增大而增大;当界面相弹性模量大于基体或夹杂时,材料的等效弹性模量会随着夹杂体分比和界面相厚度的增大而增大。而界面相的厚度和弹性模量对材料的等效泊松比的影响较小,材料的等效泊松比主要受夹杂体分比的影响,与其呈反比关系。     

Free Convection Heat Transfer From A Sphere In A Porous Medium Using A Thermal Non-equilibrium Model

This work studies the free convection heat transfer from a sphere with constant wall temperature embedded in a fluid-saturated porous medium using a thermal non-equilibrium model. The governing equations are transformed into boundary-layer partial differential equations by the coordinate transform, and the obtained governing equations are then solved by the cubic spline collocation method. The temperature distributions for fluid and solid phases are shown for different values of the porosity scaled thermal conductivity ratio, the interphase heat transfer parameter, and the streamwise coordinate. The effects of the porosity scaled thermal conductivity ratio and the interphase heat transfer parameter between solid and fluid phases on the local Nusselt numbers for fluid and solid phases are examined. Results show the local Nusset number for the porous medium can be increased by increasing the porosity scaled thermal conductivity ratio. Moreover, the thermal non-equilibrium effect is more significant for low values of the porosity scaled thermal conductivity ratio or the interphase heat transfer parameter.     

Shock Propagation in a One-Dimensional Flow Through Deformable Porous Media

We consider a one-dimensional incompressible flow through a porous medium undergoing deformations such that the porosity and the hydraulic conductivity can be considered to be functions of the flux intensity. The medium is initially dry and we neglect capillarity, so that a sharp wetting front proceeds into the medium. We consider the open problem of the continuation of the solution possibly in the case of onset of singularities, which can be interpreted as a local collapse of the medium. In particular we analyze the case in which the boundary pressure has a piecewise constant derivative.