DDA与流形元耦合计算方法及其应用.pdfVIP

第一届中国水利水电岩土力学与工程学术讨论会论文集 DDA与流形元耦合计算方法及其应用 朱爱军1，邓安福2，颜昌武’ (1．贵州大学土木建筑1程学院。贵州责用550003：2蕈庆大学土木工程学院．四川重庆400045) 摘要；在数值流形方洁的体系下，建立DDA与流形元耦合作用模型，推导了DDA块体与流彤元的接触方程，实现填石等散体系统和连续 介质共同作用的数值模拟计算。编制了DDA与流形元耦合分析的大型计算程序，通过算例表明，采片jDDA与流形元耦音计算方法模拟散体 系统和连续介质之间的若同作用足可行的，其对该类问题的模拟更接近分析对象的实蕲情况，有助于从根本上揭示分析对象变形的细观机理 和规律．井能考察更多因素对工程问题的嚣响。 关键词：流形元；DDA：散体系统；填石；连续介质{菇同作用 作者简介：朱爱军(1975一)．男，博士．副教授．主要从事岩土工程数值方法及建筑地基基础研究。 O引 言 Eh (J，Y)的0，l，2阶完全多项式·如果流形单元fq 堆石坝或填石路堤等问题属于大型的散体系统问题， 个物理覆盖仉l，Ue2、…、以。的交集形成，以f的覆盖函 对其进行数值计算，由于传统的数值方法如有限元法等是 数为 建立于连续介质模型之上的，所以一般将坝体或路堤材科 等效为均匀连续介质进行分析…。随着离散元及DDA等不 阱秆∥厶“o啦。‘}， m 连续变形分析方法的发展．从细观上对散体系统材料进行 式中 下标e表示流形单元-工dx·y)为完全多项式的基 分析已成为现实【2I“。从细观上对散体系统进行数值模拟， 本级数。如孕l及以“4为物理覆盖的覆盖自由度，对应了 符合材料本身的物理模型．有助于从根本上揭示散体材料 完全多项式的各系数。物理覆盖的覆盖函数可取0阶到高 的变形机理及力学特性。 阶的完全多项式，对应0、I、2阶情况，式(1)中卅分 对于修建于土质地基上的堆石坝(如紫坪铺高面板堆 别为1、3，6。 石坝【列)或填石路堤．利坝体或路堤的变彤必然受地基变 则流形单元的位移函数为 形的影响，所以这类工程问题的分析必须考虑填石和地基 阱喜W．,阱ua 的共同作用。在力学模型上。对路堤填石采用完全非连续 散体系统模型，地基采用连续介质的力学模型，要考虑填 卦“紫㈤％“蚺0“船)川， 石与地基的共同作用，必须实现非连续体力学模型与连续 介质力学模型的耦合计算，而目前这方面的研究很少。 式中 w,i(x，y)为物理覆盖玑．的权函数 数值流形方法是包含流形元、有限元及DDA 令物理覆盖仉。的变形矩阵为 (DiscontinuousDeformation Analysis)在内的数值方法体 【兄】_【L．毛：…7k】。 系，在这个体系中，可蚪找到流形元法和DDA共同的基 上式中 础．如它们都可以用有限覆盖的思想去阐述．适用相同的 ‰]；殴之]寸“訾㈤％“最“』。c。， 接触处理方法，可以有相同的位移模式，及相同的体系平 衡方程建立方法，因此在数值流形方法的体系下，两种模 令 型的耦合变得自然和方便。所以可以采用DDA来对堆石 {嘎)；k．．吒．。…噍。．以。]r， (4) 体建立计算模型，而地基则可以采用流形元进行模拟，建 由此。式(2)可表达为 立DDA与流形元的接触方程．就可以在数值流形方法体 ㈦=舡№)，jd一Ⅵ。 (5) 系中实现散体系统与连续介质的共同作用分析。 1．2流形元的控制方程 1 流形元的基本理论”1