udec动力计算需要注意的问题

《udec动力计算需要注意的问题》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、刚度选取模型中块体材料与节理材料的刚度或性质反差较大或者块体尺寸存在较大差异，UDEC将需要运行很长时间才收敛。节理法向和切向刚度的力学计算应小于10倍节理毗邻单元的等效刚度。即和(3.37)式中，和分别是块体材料的体积模量和剪切模量；是毗邻节理单元在法线方向的最小宽度。如果节理刚度大于10倍等效刚度，则模型求解时间将明显大于其刚度比限制在10以内的情况。在另一方面，如果法向刚度，很低，还可能存在问题。应对模型代表性应力所产生的节理法向位移有一个大致的估计（）。该位移与代表性的单元尺寸相比是小的。如果该值大

2、于毗邻单元尺寸的10％左右，则一个单元存在错误或刚度应提高。网格划分UDEC在计算中，网格单元的尺寸对计算结果有很大影响。一般而言，网格单元的尺寸越小，计算精度越高，但会占用更多的内存和消耗更多的计算时间。一般将块体切割成大量的有限差分三角形网格单元，网格单元的尺寸应不大于波长的1/8～/10。边界条件设置UDEC动力计算分析中，计算模型的边界计算模型的区域边界可能引起向外传播波的反射回模型和不允许必要的能量辐射。当模型的尺寸足够大的时候，材料阻尼将吸收大部分从边界处反射回来的能量，从而最小化消除影响。然而

3、，大尺寸的模型会增大计算负担。为了解决这个问题，Lysmer和Kuhlemeyer(1969)开发了粘性边界，通过在边界的法向和切向设置减震器来吸收能量。减震器提供的反向力为：式中，为法向反力(Pa)；为切向反力(Pa)；为介质的密度(kg/m3)；为应力波在介质中传播的纵波波速(m/s)；为应力波在介质中传播的横波波速(m/s)；边界质点的法向震动速度(m/s)；为边界质点的切向振动速度(m/s)。，。K为体积模量(Pa)，G为切向模量(Pa)，为介质的密度(kg/m3)。