基于clump的颗粒形状系数的构建.docxVIP

基于clump的颗粒形状系数的构建 1 颗粒形状对力学性能的影响 j.s.dads等人的研究结果表明，砂化学的力学性质对微观参数（如形状）的影响较大。常规宏观分析方法多以土体作为连续介质研究,且由于测量和量化分析的困难,对颗粒形状影响的独立深入研究仍较少。基于离散元的颗粒流软件的出现,克服了传统连续介质力学模型的宏观连续性假设,可以从细观层面上对砂土的工程特性进行数值模拟。 使用颗粒流方法对颗粒形状影响的研究在国内外已取得一些进展。在国外,B.S.H.Kim和N.Kim研究了颗粒的形状、颗粒正向与切向刚度比等对颗粒组原生各向异性的影响;M.Zeghal用3个颗粒形成一个“cluster”,作固结试验并往复加载,研究了颗粒形状对回弹模量的影响;A.K.Ashmawy等研究了在不排水往复加载试验中,颗粒形状对液化性质的影响;A.A.Lizcano等在离散元中,用任意多边形颗粒取代圆形颗粒,研究了颗粒形状对颗粒组的零界应力水平、配位数、力链形成等的影响;J.C.Santamarina和G.C.Cho研究了颗粒形状对颗粒组固有各向异性以及应力诱导各向异性的影响。在国内,曾远和周健认为,颗粒形状对剪胀效应的影响比颗粒摩擦因数影响更明显;常在研究了颗粒形状对片状黏性土力学性质的影响。 尽管已有研究在一定程度上揭示了颗粒形状对颗粒材料力学特性的影响,但从颗粒变形机制的角度定义颗粒形状系数,以及在常规土工试验的模拟中定量探讨颗粒形状与试验结果关系仍亟待进一步深入研究,如:(1)在颗粒堆积试验中,对颗粒粒径以及粒间摩擦因数的研究较多,对颗粒形状与自然休止角以及天然孔隙率联系的研究较少;(2)在双轴试验中,材料的微观参数,如颗粒摩擦因数、孔隙率、刚度比等对材料宏观力学特性影响的研究较多,但对于颗粒形状,特别是颗粒形状的量化参数与材料的宏观力学特性联系的研究仍较少;(3)在直剪试验中探讨颗粒形状与材料抗剪强度,并从微细观角度分析其原因的研究还少有报道。 基于上述分析,本文在PFC2D中用等截面积法和等密度法生成圆形、类椭圆形、类正方形、类三角形等4种不同形状的颗粒组,并结合颗粒材料变形特征,定义颗粒的“圆度”与“凹凸度”,构建一个能定量描述颗粒形状的形状参数。用PFC2D内置FISH语言开发颗粒堆积试验、双轴试验和直接剪切试验,探讨颗粒形状对试验结果的影响,分析类砂土材料宏观力学特性随颗粒形状变化的规律。 2 生成和描述 2.1 颗粒形状对材料力学特性的影响 图1为在PFC2D中采用clump命令形成4种形状颗粒黏结示意图。 类椭圆形颗粒、类正方形颗粒、类三角形颗粒的生成除需要形心点的坐标外,还需要一个表示其倾斜角度的方位角。在PFC2D模型中方位角是随机指定的,其程序主要是通过urand和pi命令的结合实现。为了考察颗粒形状对材料宏观力学特性的影响,模型颗粒的生成需要保证2个条件: (1)截面积相等,即3种异形颗粒外轮廓线所围成的面积与半径为R1的圆形颗粒面积相同。对类椭圆形有R21=0.707R1,R22=R23=0.69R1,R24=R25=0.577R1,R26=R27=0.424R1,对类正方形有R3=0.443 9R1,对类三角形有R4=0.698R1。 (2)质量密度相等。PFC2D中为满足密度相等条件,需要通过调整黏结颗粒的密度达到,其换算公式如下: 式中:ρd为颗粒黏结体的密度,∑Ad为各黏结颗粒面积之和,Ap为由外轮廓线围成的面积,ρp为颗粒密度(注意:在PFC2D中,颗粒都假设为单位厚度的圆盘)。 本文颗粒生成顺序都是先在模型中生成圆形颗粒,再按照等面积法和等密度法将圆形颗粒替换成其他形状颗粒。 2.2 颗粒外轮廓特征的表征 对颗粒形状影响的研究首先需要对颗粒形状进行恰当的描述。已有的颗粒形状量化分析方法,如涂新斌和王思敬给出的二维颗粒参数以及E.T.Bowman等采用分形理论和傅里叶分析方法对颗粒形状进行描述的方法,或没有将颗粒形状定义与其变形特征相联系,或其数学计算过程复杂,限制了其推广和应用。在工程中,通常用一个数或几个数从不同角度来表述颗粒的形状特征。为此笔者拟在已有研究的基础上,结合颗粒材料变形过程,将形状参数与定义的物理意义相联系。实际颗粒材料的受力变形是由颗粒的移动产生的,在不考虑颗粒的破碎时,颗粒的移动包含颗粒的滑动(平移)及翻滚(自转)2个过程。颗粒自转和平移都受到周围颗粒的阻碍,转动和平移过程都与颗粒的外轮廓特征密切相关。颗粒愈接近于圆形,其自转愈容易。此外,颗粒的表面凹凸程度与其平移过程中受到的阻力联系较紧密,颗粒的凹槽愈深,那么与其他颗粒的咬合作用更加明显,颗粒移动产生的抗阻力就愈大。土体的宏观内摩擦力是颗粒抵抗滑动与翻滚能力的体现(用参数内摩擦角?表达),从细观分析有3部分来源:(1)颗粒抵抗转动能