《FLAC3D_前处理程序：》.docVIP

FLAC3D 前处理程序： 强大的处理功能 FLAC3D是用于工程力学计算的三维解析有限差分软件。该程序是FLAC二维计算程序在三维空间的扩展，用于模拟三维土体、岩体或其它材料体力学行为，尤其是达到屈服极限时的塑性流变行为。 物体由多面体单元所表示，可以通过调整三维网格的方法，以适应研究体真实的形状， 每个单元力学行为是对应力-应变法则和边界力、约束条件的响应。材料能产生屈服和流动，而且网格也能变形（in large-strain mode），并随材料移动。 FLAC3D采用了解析Lagrangian计算方法和混合离散元技术，以确保塑性破坏和流动的模拟效果非常准确。 FLAC3D 提供十种材料模型，即： (1) null; (2) elastic, isotropic; (3) elastic, orthotropic; (4) elastic, transversely isotropic; (5) Drucker-Prager plasticity; (6) Mohr-Coulomb plasticity; (7) strain-hardening / softening Mohr-Coulomb plasticity; (8) ubiquitous-joint plasticity; (9) bilinear strain-hardening / softening ubiquitous-joint plasticity; (10)modified Cam-clay plasticity。 每一种模型代表特定地质材料的基本力学行为。 总之，FLAC3D为地质工程技术提供了理想的分析工具。 建立计算模型中的困难 FLAC3D为用户提供了初始单元网（PRIMITIVE MESH），共十二种。即：Brick、Degenerate brick、 Wedge、 Pyramid、 Tetrahedron 、Cylinder、 Radial Brick、 Radial Tunnel、 Radial Cylinder、Cylindrical Shell、Cylinder Intersection、Tunnel Intersection。这些初始单元网对于建立三维规整的工程地质体模型具有快速、方便的功效。但是，对于较复杂的地质体，如整个坝区的地质计算模型，FLAC3D表现得比3D-δ等软件更加费力、耗时。 前处理程序简介 为了解决工程中的实际问题，更准确地反映研究区的地质地形条件，实现更快速、便捷地建立较复杂的三维FLAC计算模型，编写了此程序。 程序设计的基本思路 程序采用FORTRAN语言编写。见结构框图。 数据的准备 研究区的地形网格化数据和地层的高程、产状。 网格化数据可以直接从地形图上读取，经Surfer 软件网格化处理后，生成ASIC代码的数据文件（*.DAT）。 地层的高程是指模型左下脚点（X=0，Y=0）对应的地层高程。产状分别用地层在X、Y轴方向上的坡降来表示。 输出结果 程序直接生成与FLAC3D相匹配的单元体（ZONE）文件格式: gen zone brick p0 0 0 0 p1 10 0 0 p2 0 20 0 p3 0 0 10 p4 10 20 0 p5 0 20 15 p6 10 0 10 p7 10 0 0 size 2 5 5 根据一般地质体的特征及单元体的空间组合关系，在此程序中仅选用了四种初始单元体，它们是Brick、 Wedge、 Pyramid、 Tetrahedron。 解决的问题 能迅速建立无断层、无风化带划分的较复杂的研究区的计算模型。见下图。 模型1 模型2 FLAC3D微模型及其说明 共建立了三大类五十多种微模型，这三大类概括无地质界限模型、单地质界限模型、双地质界限模型。说明如下： A．无界限模型： 仅一个FLAC3D微模型。 B．单界限模型： 共计6或8点，即1、2、3、4点和5、6点或5、6、7、8点。 其中根据5、6点或5、6、7、8点在四条边的不同分布，划分为7种分布。每一种分建立2个FLAC3D微模型。分布类型以代码表示。 C．双界限模型： 共计8点 即1、2、3、4点和5、6、7、8点。 其中根据后四点在四条边的不同分布，划分为三类19种分布。每一种分布可建立2-4个FLAC3D微模型。分布类