下载文档原格式
科技资讯科技资讯S I N &T NOLOGY I NFO RM TI ON 2008NO .28SC I EN CE &TECH NO LOG Y I N FOR M A TI O N 信息技术本文利用大型有限元软件F L AC-3D ,对上海市宜山路站9号线一期工程基坑开挖过程,建立了基坑工程开挖、支护模拟的有限元模型,模拟分析了该工程开挖、支护全过程及在此过程中支护结构后面土体的变形和沉降分析。
在分析中,采用空气单元(N ul l 单元)模拟挖掉土体,通过时间步来定义基坑工程开挖、支护过程,从而可以计算分析每一工况下基坑支护结构的内力、变形和地表沉降。
1显式有限差分软件FLA C -3DF L A C3D 的输入和一般的数值分析程序不同,它可以用交互的方式,从键盘输入各种命令,也可以写成命令(集)文件,类似于批处理,由文件来驱动。
因此,采用F L A C 程序进行计算,必须了解各种命令关键词的功能,然后,按照计算顺序,将命令按先后,依次排列,形成可以完成一定计算任务的命令文件。
F L AC3D 是二维的有限差分程序F L AC2D 的护展,能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型,在外力作用下,当材料发生屈服流动后,网格能够相应发变形和移动(大变形模式)。
F L A C3D 采用的显式拉格朗日算法和混合—离散分区技术能够非常准确的模拟材料的塑性破坏和流动。
由于无需形成刚度矩阵,因此,基于较小内存空间就能够求解大范围的三维问题。
FL AC3D 采用ANSI C++语言编写的。
1.1基坑开挖模拟的方法在城市地区进行深基坑工程开挖,初始应力场为重力场。
当用有限元法进行模拟基坑开挖时,一般采用反转应力释放法和空单元法,但使用有限差分法进行模拟时,一般使用空单元法。
空单元法的开挖效果是通过被挖掉单元的“空单元化”,即将要挖掉单元的刚度矩阵乘以一个很小的比例因子,使其刚度贡献变得很小可忽略不计,同时使其质量、荷载等效果的值也设为零来实现的,故称为空单元法。
随着国家铁路网建设推进，多地高铁站拔地而起，与之对应的还有集合交通客运中心、地下停车场、地下接送客和公交停车场等多功能一体的枢纽配套工程。
实际项目建设时序基本上是高铁站先于配套建成，为实现旅客换乘无缝衔接，枢纽配套往往会临近高铁站，枢纽一般为地下两层，这导致地下空间的基坑开挖施工必然会对已建高铁站产生影响。
国内学者对此进行了一定的研究，王菲[1]通过数值模拟深基坑开挖的影响，从累积沉降、桥墩差异沉降和水平变形等指标判断基坑设计的安全合理性，朱一康[2]采用工程类比和有限元分析相结合的方法，提出设计、施工及监测等方面的建议。
该文以某综合交通枢纽地下空间工程为例，采用MIDAS GTS有限元软件数值模拟分析深大基坑开挖对临近敏感建筑的影响。
1 工程概况该项目整体位于高铁站的东侧，站前广场地下室大部分分为地下2层，基坑开挖深度约12m，局部地下3层紧邻高铁站房，基坑开挖深度约17m，基坑开挖面积约9万m2，地下广场基坑内部还存在已建的地铁站结构，车站底部位于地下广场基坑坑底以下约12m，项目位置关系如图1所示。
该项目地处三角洲平原地貌，地势平坦开阔，场区内主要分布农田和苇塘，土层分布均匀，自上而下分别为①素填土、②粉土、③粉质黏土与粉土、④粉土与粉砂、⑤粉砂、⑥粉土与粉砂、⑦粉土与粉砂、⑧粉土与粉质黏土、⑨粉砂与粉土、⑩粉质黏土、11粉土粉砂与粉质黏土、12粉土与粉砂和13中粗砂。
场区内地下水丰富，周边分布有河沟，常年平均水位在地面下1m左右，粉土、粉砂层渗透系数较大，基坑底部以上及以下4m范围均为孔隙潜水含水层，坑底以下为4m~12m微承压含水层和16m~30m承压含水层，各土层主要计算参数见表1。
高铁站房为大跨钢筋混凝土框架结构形式，线侧站房两层，局部地下一层，线下站房一层，两侧局部设夹层，建筑高度约24m。
结构采用桩承台基础形式，桩基采用钻孔灌注桩，桩径800mm，桩长36~40m，桩端持力层在承压含水层以下的中粗砂层。
7数值模拟7.1数值模拟方法简介数值模拟技术作为一种研究手段，已经被广泛的应用于各行各业领域的研究中。
目前，数值分析方法主要分为二大类：一类是以有限差分法为代表，其特点是直接求解基本方程和相应的定解条件的近似解；另一类数值分析方法是首先建立和原问题基本方程及相应定解条件等效的积分方法，然后据之建立近似解法。
LS-DYNA作为世界上最著名的通用显示动力分析程序，能够模拟真实世界的各种复杂问题，特别适合求解各种二维三维非线性结构的高速碰撞，爆炸和金属成型等非线性动力冲击问题，同时可以求解传热，流体及流固耦合问题，在工程应用如汽车安全设计，武器系统设计，金属成型，跌落仿真等领域被广泛应用。
本次采用ANSYS/LS-DYNA，进行混凝土支撑梁结构爆破拆除数值模拟研究。
在ANSYS/LS-DYNA环境下，数值模拟的实现总体上分为两个过程：在ANSYS中建立结构实体模型，完成有限元网格的划分，输出有限元模型信息即输出关键字文件；编辑关键字文件，在DYNA环境下完成对结构倒塌过程的数值模拟计算。
对结构有限元模型的建立过程，数值模拟中采用的钢筋和混凝土材料模型、接触方式等各种计算控制项进行了阐述。
LS-DYNA程序中主要提供如下几种计算方法：（1）Lagrange算法坐标固定在物质上或者说随物质一起运动和变形，处理自由面和物质界面非常直观，由于网格始终对应物质，因此能够精确的跟踪材料边界和描述物质之间的界面，这是Lagrange算法的主要优点。
但是，由于网格随材料流动而变形，一旦网格变形严重，就会引起数值计算的不稳定，甚至使得计算无法继续进行（如发生负体积或复杂声速等问题）。
因此，Lagrange算法在处理大变形大位移问题时，有其无法克服的弊端。
（2）Euler算法网格被固定在空间，是不变形的。
物质通过网格边界流进流出，物质的大变形不直接影响时间步长的计算。
因此，欧拉算法在处理大变形问题方面具有优势。
欧拉方法通过输运项计算体积、质量、动量和能量的流动。
欧拉计算可以直接通过在离散化格式中包括迁移导数项进行，或通过二步操作完成。
二步法操作的第一步主要是拉格朗日计算，第二步输运阶段是重分计算网格相当于回到它的原来状态。
LS-DYNA程序采用后一种方法。
欧拉算法的缺点是网格中物质边界不清晰，难以捕捉各物质界面。
（3）ALE方法吸取了欧拉法和拉格朗日法两种方法的优点。
ALE算法能够进行自动重分网格操作。
它包括拉格朗日时间步，然后是一个输运步。
输运步可以采用三种方法：1. 发生合理的网格变形时空间网格不再重分（拉格朗日）；2. 发生严重的网格变形时重分成原始形状（欧拉）；3. 发生严重的网格变形时重分为合理的形状，因此允许网格拓扑（拉格朗日和欧拉）。
混凝土是土木工程结构中应用极为广泛的材料，其最本质的特点是材料组成的不均匀性，并且存在初始微裂缝。
从混凝土受单轴压力时的应力应变关系来看，混凝土卸载时有残余变形，不符合弹性关系；如果对其应用弹塑性本构关系，又很难精确定义屈服条件。
此外，混凝土在到达应力顶峰后，其应力-应变关系曲线有一下降段，即存在应变软化现象，所有这些都给建立混凝土的本构关系带来困难。
多年以来，众多学者进行了大量的试验和理论研究，提出了各种各样的混凝土本构模型。
第一：混凝土本构模型可以分为下面几种：（1）线弹性类本构模型。
线弹性类本构模型是以弹性力学为基础的模型，当混凝土无裂缝时，将混凝土看成线弹性匀质材料而采用线弹性本构模型。
虽然混凝土的变形特性是非线性的，但在一些特定的情况下(比如描述混凝土受拉时的工作性能)，采用线弹性类本构模型进行分析还是有足够的精度的，其线弹性本构关系可用广义虎克定律来表示：kl ijkl ij C εσ=（1）式中，ijkl C 为材料弹性常数，为四阶张量，共有81个常数。
按照材料假设的不同，又可分为各向异性本构模型、正交各向异性本构模型、各向同性本构模型等，其中ijkl C 根据材料的不同而变换。
（2）塑性理论类本构模型：塑性理论类本构模型是以塑性流动理论为基础，代表性的模型主要有：Mises 条件的模型、理想弹塑性脆性断裂模型、应变或工作硬化塑性理论模型等。
模型中考虑了混凝土加载路径和混凝土的硬化，在混凝土的应力-应变全曲线中，有上升段和下降段。
自从Drucke公设和Пyushin公设出现之后，经典塑性力学得到飞速发展，混凝土塑性力学模型也是基于这些公设建立的。
以塑性理论为基础的混凝土本构模型，在对其加载面，包括初始屈服面，后续加载面和破坏包络面等特征面的研究中，这些特征面若以应力空间来表示时，当应力达到屈服后，材料发生应力松弛；若以应变空间表示时，当应变达到松弛面后，材料发生应变松弛。
基于应力状态屈服面或破坏包络面的塑性理论类型的本构模型有弹性-全塑性模型、线弹性-硬化塑性-断裂模型等；基于松弛面的塑性理论类型的本构模型有塑性模型、塑性断裂模型、硬化断裂模型等。
所有这些模型所做的假设与混凝土的实际性能还存在很大的差别，而且模型的表达式和计算均较复杂，目前还不便于应用。
（3）其他力学理论类本构模型，许多学者还以新型交叉的力学分支的理论为基础研究混凝土的本构模型。
内时理论模型：内时理论模型最初由Valanis于1971年提出，其基本概念为：塑性和粘塑性材料内任一点的现时应力状态是该点整个变形领域内和温度历史的泛函，而特别重要的是该历史是用一个取决于变形中的材料特性和变形程度的内时(Intrinsic Time)来量度的。
这种模型采用了非弹性应变能逐渐积累的方法而不需考虑塑性理论中的屈服面和流动法则，所以该理论尤其适合没有屈服面的混凝土材料。
由于内时理论能描述混凝土的复杂变形的历史，因而为各国学者所重视。
但由于表达式过多，确定参数又不容易，所以对其推广和应用仍有大量工作要做。
损伤理论模型：针对混凝土材料本身是一种具有固有缺陷－微裂纹的特点，很多学者将损伤力学引用到混凝土本构模型的建立中。
损伤力学研究材料或构件从原生缺陷到形成客观微裂纹直至断裂的过程。
也就是通常指的微裂纹的萌生、扩展或演变、宏观微裂纹的形成、裂纹的稳定扩展和失稳扩展全过程。
损伤力学主要是在连续介质力学和热力学的基础上，用固体力学方法研究材料的宏观力学性能的演变直到破坏的全过程。
20世纪70年代末期，损伤力学局限于研究材料在宏观裂纹出现以前的阶段，当宏观裂纹出现以后，则采用断裂力学的理论和方法进行研究，这是无耦合的分析方法。
实际上，在宏观裂纹出现后，材料的损伤对裂纹尖端附近及其它区域的应力和应变均有影响。
因此，合理的方法应该将损伤耦合到本构方程中进行分析和计算。
这样由于本构方程中将有关的力学参数和损伤进行了耦合，所以分析和计算就变得更为复杂。
非线弹性类本构模型：为了克服线弹性模型的缺点，然后产生了以割线模型形式的非线性弹性类本构模型。
这类模型中，具有代表性的是超弹性模型和亚弹性模型。
除了上述模型外，还有以许多种理论相结合到处的混凝土的本构模型。
第二：LS-DYNA程序中提供的混凝土材料模型：在材料模型方面，LS-DYNA目前拥有近150余种金属和非金属材料模型，涵盖了弹性、弹塑性、超弹、泡沫、玻璃、土壤、混凝土、流体、炸药等材料模型以及多种状态方程，可以考虑材料的失效、损伤、粘性、蠕变、与温度相关、与应变率相关等材料性质。
此外，程序还支持用户自定义材料功能。
LS-DYNA程序中提供的混凝土材料模型主要有以下几种：塑性随动模型*MAT-PLASTIC-KINEMATIC、H_J_C模型、*MAT_DRUCKER_PRAGER 模型和*MAT_SOIL_AND_FOAM模型。
本论文有限元模型中的混凝土本构模型选用LS-DYNA中的塑性随动模型*MAT-PLASTIC-KINEMATIC。
塑性随动模型（*MAT_PLASTIC_KINEMATIC）。
这是一种各向同性、随动硬化或各向同性和随动硬化混合模型，与应变率相关，可考虑实效，参数简单，较容易确定。
7.2 钢筋混凝土支撑梁结构爆破拆除有限元模型的建立7.2.1爆破拆除数值模拟基本假设考虑到建立整体钢筋混凝土支撑梁结构模型过程的复杂性，计算结果的可行性及准确性，数值模拟中对模型进行了简化，模型的简化出于以下基本假设：（1）不考虑混凝土支撑梁内的炸药爆炸对整体结构的作用。
如果模拟中按照实际情况建立模型，单元数目巨大，受到当前计算机能力的限制，无法实现数值计算。
由于爆炸对钢筋梁的变形影响较小，所以在模拟中可以直接删除爆破药包附近的混凝土，并没有考虑钢筋梁结构在炸药爆炸下的响应问题。
（2）钢筋混凝土梁在爆破后，剩余部分的表面是光滑的，由于上面假设是直接删除爆破部分，所以无法模拟混凝土爆破后的形状。
（3）钢筋混凝土的支柱地面没有相对位移，简化为完全固结，即不考虑钢筋梁的坍塌与支撑柱和地面之间的相互作用。
7.2.2钢筋混凝土有限元模型建立过程采用ANSYS建立有限元模型时，采用的数值单位均为国际单位，即长度单位为m，时间单位为s，质量的为kg，压力的单位为Pa，速度单位为m/ s。
为了由实际情况对比实验模拟的准确情况，我们选用工程实例作为有限元模拟的对象。
综合考虑有限元计算的精确性和时效性，来确定单元尺寸和单元数量。
7.2.3钢筋混凝土模型的选取采用数值模拟方法进行钢筋混凝土结构爆破拆除分析前，首先要就研究的具体问题选择好用于模拟钢筋混凝土材料的物理模型。
当前，数值计算中，用于模拟混凝土钢筋材料的有限元模型主要有两类：第一种是把钢筋和混凝土单元各自划分为足够小的单元，分别考虑钢筋和混凝土的贡献，称为分离式模型。
第二种也是把钢筋和混凝土单元包含在一个单元之中，统一考虑钢筋和混凝土的作用，称为整体式模型。
（1）整体式模型在整体式模型中，将钢筋弥散于整个单元中，并把单元视为连续均匀的材料。
钢筋对整个结构的贡献，可以通过调整单元的材料力学性能参数来体现，例如提高材料的屈服强度、材料的弹性模量等。
其优点是建模方便，分析效率高，但是缺点是不适用于钢筋分布较不均匀的区域，且得到钢筋内力状态比较困难。
主要用于有大量钢筋且钢筋分布较均匀的构件中。
(2)分离式模型，位移协调利用空间梁单元beam161建立钢筋模型，和混凝土单元共用节点。
其优点是建模方便，可以任意布置钢筋并可直观获得钢筋的内力。
缺点是建模比整体式模型要复杂，需要考虑共用节点的位置，且容易出现应力集中问题。
（3）分离式模型，界面单元前两种混凝土和钢筋组合方法假设钢筋和混凝土之间位移完全协调，没有考虑钢筋和混凝土之间的位移，而通过加入界面单元的方法，可以进一步提高分析的精度。
同样利用空间梁单元beam161建立钢筋模型。
不同的是混凝土单元和钢筋单元之间利用弹簧模型来建立连接。
不过一般钢筋混凝土结构中钢筋和混凝土之间都有比较良好的锚固，钢筋和混凝土之间滑移带来的问题不是很严重，一般不必考虑。
