FLAC3D前后处理及三维模型构建方法资料精要.ppt

Flac3d前后处理和三维模型构建方法 Flac3d前后处理和三维模型构建方法 AutoCAD在DWG文件中以实体名存储各种点、线、面、实体。实体名是AutoCAD中专门设置的一种特殊数据类型，是赋予图形实体的数字标号。通过AutoLisp语言开发的程序可以访问各种点、线、面、实体信息，例如边坡中的地形线，开挖面，地下厂房，引水洞、平硐等 现阶段常用的插值方法有：Kriging法、Delaunay法、Polymesh法、Shepard法、Multiquadric法等等。滑动最小二乘法（Moving Least Squares Method，简称MLSM）是一种精度较高的方法，其克服了经典最小二乘法在拟合中的不足，实质为方差泛函极小化序列与最小二乘法。 Sufer中都是栅格插值方法, 所以如果要按给定任意网格插值的话,需要自己编程. Flac3d前后处理和三维模型构建方法 FLAC3D数值模型分解后的最小单元类型为Brick（六面体）、Wedge（三棱柱体）、Pyramid（金字塔体）和Tetrahedral（四面体）单元。 Flac3d前后处理和三维模型构建方法 赋予初始应力 gen zone brick size 6 6 6 model elas prop bulk 5e9 shear 3e9 ini sxx=-5e6 pause solve 加上应力边界条件 gen zone brick size 6 6 6 model elas prop bulk 5e9 shear 3e9 ini sxx=-5e6 apply sxx=-5e6 range x -0.1 0.1 apply sxx=-5e6 range x 5.9 6.1 pause pl con sxx solve pl con sxx 注意是约束x方向---------变形参数并不影响应力场，它只影响位移场。如果考虑重力的时候，变形参数影响泊松比，从而影响侧压力系数，才会影响水平的应力场分布。 new ;Example 3.24 Initial stress state for an irregular free surface gen zone brick size 15 15 10 p0 0,0,0 edge=100.0 model elastic prop shear 3e8 bulk 5e8 def mountain gp = gp_head loop while gp # null zz = sqrt(gp_xpos(gp)^2 + gp_ypos(gp)^2) dz = 0.06 * sin(0.2 * zz + 100.0) ; Sum Fourier terms for dz = dz + 0.06 * sin(0.22 * zz - 20.3) ; quasi-random surface dz = dz - 0.04 * sin(0.33 * zz + 33.3) ; topology. gp_zpos(gp) = 0.5 * gp_zpos(gp) * (1.0 + dz) gp = gp_next(gp) end_loop end mountain fix x range x -.1 .1 fix x range x 99.9 100.1 fix y range y -.1 .1 fix y range y 99.9 100.1 fix z range z -.1 .1 按照自重求初始地应力，形成的应力场表现为竖直向=自重应力，水平向=自重应力×侧压力系数k0。k0=v/(1-v). K0在这个方法中显然是和泊松比有关,然而土的泊松比变化不大,在0.25-0.35左右, 由此求得的k0在0.3-0.6之间,和phi在45-27度之间. 在上面,泊松比是通过语句prop bulk 5e9 shear 3e9间接来定义的. 这种方法会产生很大的位移,当然FLAC有一个投机取巧的方法来减小位移: 就是在求初始应力的时候,把材料设为弹性,并且用很大的弹性模量确(比如上面的model elas; prop bulk 5e11 shear 3e11).修改变形参数不影响初始应力场的分布，只影响位移场的大小。 在此之后,初始化位移为0,并把材料的参数恢复正常. 但是这种方法存在两个问题： 第一是,由此求得初始侧压力系并不和实际一致,因此初始地应力也和实际不一致. 在这一点上用第二种方法可以得到比较正确的初始应力场. 或者也可以通过正确的k0反算一个不切实际的泊松比v代入来求. 第二也是由于设置材料为弹性的,所以是不会有破坏产生,不管初始产生的剪应力有多大,然而在材料设为弹塑性时,如果剪应力过大