水工隧洞施工及运行的ANSYS模拟

水工隧洞施工及运行的ANSYS模拟4.1数值模拟对象4.2有限元建模4.2.1单元定义和材料定义4.2.2实体建模4.2.3网格划分4.2.4边界条件和初始条件4.3水工隧洞施工过程模拟4.3.1初始状态模拟4.3.2毛洞开挖工况模拟4.3.3毛洞支护工况模拟4.3.4计算结果查看及处理4.4水工隧洞运行期模拟4.4.1运行期内水压力的模拟4.4.2运行期外水压力的模拟水工隧洞施工及运行的ANSYS模拟由于ANSYS在水利工程中应用面广，可以广泛用于水利工程的各个专业领域中，包括水工隧洞、地下厂房、高边坡、重力坝、拱坝、截流堰等水工结构；水轮机组的动力分析；水文预测以及高速水力学等。基于对ANSYS基本操作的进一步熟悉，并建立对水工结构ANSYS分析的概念，本章以一个典型水工隧洞的开挖过程为例，简单介绍ANSYS在水利工程中的应用，并以此作为初学者的入门实例。数值模拟对象对于实际工程而言，对所要数值模拟对象的熟悉程度是进行有效的ANSYS建模和正确进行有限元分析的基础，熟悉的内容主要包括：研究对象地形地质条件（不同的地质分层、断层、节理、裂隙等）、地层及结构的物理力学参数（如果涉及到渗流分析或温度分析，则还需要水力学参数和热力学参数等）、纵横剖面、水文条件以及荷载条件等，以及工程的施工工法，工序安排等，从而为有限元的建模提供前提条件。需要注意的是，作为有限元数值模拟，只是对实际工程的高度近似，换句话说，不可能达到百分之百的相同。因此，对实际工程需要进行一定的简化，否则是无法、也不可能进行数值模拟的。图4-1水工隧道的简单实例问题描述：以一个简单隧洞为例，隧洞内径6m，衬砌厚40cm，地层均质，隧洞进行全断面开挖，开挖后进行一次性衬砌支护。问题抽象：从描述中可以分析，分析为平面应变问题，问题中涉及两种材料（岩石和混凝土衬砌），研究区域根据一定的规则选取为100m×100m（在后面的章节中进行介绍）矩形区域，工程分析过程分为3步，即初始状态>毛洞开挖>支护。有限元建模启动ANSYSProductLaucher，定义好工作目录和文件名称。建议不同的工程建立不同的工作目录，文件名称尽量取与工程名称相关且最好包含日期信息，以便日后对计算过程的回顾和再利用。如目录取为Shuigong，文件名取为Tunnel060824，如图4-2。然后运行Run（如目录不存在，则会弹出对话框提示，直接点击确定，则在相应位置新建，若已存在，则点击Browse去选取，文件名同样如此），进入ANSYS主操作界面。图4-2工作目录和文件名的定义单元定义和材料定义1．单元定义从前面分析中可以看出，本次分析为平面应变问题，单元形式为平面应变单元。则单元的定义过程如下：进入前处理状态Preprocessor（在完成网格划分前都不用退出前处理器）；点击ElementType>Add/edit/delet，弹出ElementType对话框（图4-3）。图4-3ElementType对话框图4-4ElementType对话框点击增加单元按钮Add，弹出单元类型库对话框LibraryofElementType（图4-4）。在左侧A框中选择Solid，对应在左侧B框中选择Quad4node42类型，即为平面四边形单元，单元编号默认为1，点击Ok退出，则在ElementType对话框中显示已定义的单元Plane42（图4-5）。图4-5平面单元的选取图4-6已定义的平面单元Plane42平面应变的定义。一般来说，Plane42默认为平面应力单元，因此需要对其进行修改。点击图4-6中选项按钮Options…，弹出对话框（图4-7），在A框K3的下拉菜单中选择Planestain项，在B框K5项选择Integrationpts，点击Ok确定，退出到图4-6状态。图4-7平面应变的定义点击Close，完成单元定义。2．材料定义从4.1中分析可知，整个分析中涉及到2种材料，即围岩体和衬砌混凝土（均假定为各向同性材料）。因此需在ANSYS中定义2种材料，定义过程如下：点击MaterialProps>MaterialModels，弹出材料模型行为定义框DefineMaterialmodelbehavior（图4-8），左侧列表中默认一种未定义参数的材料，这里定义为岩石材料，编号1。图4-8材料模型行为定义框双击右侧框内的Structural（结构）>Linear（线性）>Elastic（弹性）>Isotropic（各向同性），弹出对于材料1的线性各向同性材料性质的定义框（图4-9）。图4-9材料性质定义在EX框（弹性模量）框中输入4e9（单位为Pa），Prxy（泊松比）框中输入0.25，点击OK确定，退回图4-8状态，则左侧框中出现线性各向同性性质项（图4-10）；图4-10已定义的线性各向同性性质项点击图4-10右侧框中的密度项Density，弹出对话框密度定义对话框，输入2400。图4-11密度的定义（通过以上1-4步则定义了材料1线弹性分析所需的参数，图4-12）图4-12已定义的材料1参数定义材料2（混凝土衬砌）的参数。在图4-12中点击材料Material>Newmodel（新建材料），则左侧框中出现MaterialModelNumber2，重复（1）～（4），定义材料2的参数为Ex=2e10，Prxy=0.167，Dens=2500。图4-13。图4-13材料2（衬砌混凝土）参数的定义在图4-12中点击材料Material>Exit退出材料定义，点击或进行保存。通过以上步骤，则定义了所需的材料，若有更多材料，则重复以上步骤即可，需要注意的是，不同的分析（弹性分析，弹塑性分析，蠕变分析，动力分析）所需的参数是不同的。实体建模所谓实体建模，就是对通过创建点、线、面的方式，在ANSYS实际反映研究对象的体形特征，从而建立网格划分的基础。本实例的最终实体模型包括围岩区域、衬砌区域（开挖后要回填）、开挖区域（不回填，形成隧洞）三个区域，在实体上表现为面Area。建立面有两种方式，一是通过生成点生成线，再由线生成面，二是直接通过ANSYS工具生成面。1．通过点、线生成面建模过程：（1）生成关键点点击Modelling>Create>Keypoints>InActiveCs（通过输入坐标直接生成关键点），弹出对话框（图4-14）。其中第一行为关键点编号，不输则按顺序系统自动进行编号，第二行为三个方向坐标值，平面问题时则Z值不输。按上述方法依次创建点（-50,-50），（50,-50），（50,50）和（-50,50），完成四个关键点的创建，则工作面出现4个关键点（图4-15）。图4-14关键点的创建图4-15已创建的关键点（2）由点生成线点击Modelling>Create>Lines>Lines>StraightLine（通过关键点创建直线）。弹出选点建线对话框CreatStraightLine（图4-16）。先连接关键点1和关键2。鼠标依次点击1点和2点，则生成线1（自动按顺序编号为1），按此方法，可以连接点2和3、3和4，4和1，依次生成线2～4（图4-17）。图4-16创建直线对话框图4-17创建直线（3）由线生成面点击Modeling>Create>Areas>Arbitrary>ByLines（通过线创建面）。弹出CreatAreabyLines（图4-18）。依次点击线1、2、3、4，点击OK，则形成一个矩形面（图4-19）。图4-18创建直线对话框图4-19创建面（4）创建隧洞及衬砌由于隧洞及衬砌为圆形/圆环形，可以通过ANSYS的工具直接创建，再通过布尔运算对圆进行操作，生成开挖区、衬砌及围岩区。具体操作如下：点击Modeling>Create>Areas>Circle（圆形）>SolidCircle（实体圆），弹出实体圆面创建框（图4-20），图中的三个输入框中，WPX和WPY为圆心位置（不输值时，默认圆心位置在整体坐标系的原点；若圆心位置不在整体坐标系原点，则需通过移动工作平面来实现圆心位置的确定）。这里，隧洞内径为6m，外径为6.4m。在图4-20对话框中的Radius中输入6，点击Apply，创建内圆；再输入6.4，点击OK，创建外圆。（点击鼠标右键，点击Replot，可以看到三个有重复区域的面，图4-21）图4-20创建实体圆对话框图4-21创建的三个实体圆通过上面的操作，在工作面上存在3个面Area，但三个面存在重合，因此，必须对三个面进行布尔操作，形成相互独立但边界连接的三个区域，即围岩区、衬砌区和开挖区。点击Modeling>Operate>Booleans>Overlap>Areas（对面进行叠合操作），弹出对话框（图4-22）。选择对话框中的PickAll，对所有的面进行操作，从而形成三个区域。图4-22对面进行Overlap对话框图4-23经过布尔操作的圆面（5）精细实体模型经过上面的操作，初步形成了实体模型，可以进行网格的自由划分，但这样获得的网格往往没有规律性，同时没有突出重点（比如在隧洞周边的网格需要密一些）等等。因此，可以对这些区域进行切割，划分为4个相同且对称的区域。具体做法是：在矩形区域的4个线/边中点创建4个关键点，连接成直线，然后用直线去切割3个面。ANSYS操作：显示所有的线：UtilityMenu>Plot>Lines；在线1的中点创建关键点：Modeling>Create>Keypoints>OnLinew/Ration，对话框（图4-24），选择要创建点的线，点击Apply，弹出对话框（图4-25），在Lineratio(0-1)中输入0.5（表示中点），KeypointNumber默认，点击OK，即在创建了线1上的中点（图4-25），按此方法，创建其他3条边线上的中点关键点。图4-24在线的中点创建关键点连接关键点13和15，14和16，形成2条直线（按图4-16～图4-17的方法创建）；显示所有的元素：UtilityMenu>Plot>Multi-Plots；对3个面进行切割：点击Modelling>Operate>Booleans>Divide>AreabyLine，弹出对话框DivideAreabyLine（图4-25），点击PickAll，弹出对话框，再点击PickAll，完成对面域的划分。为了便于查看面，将面以颜色进行区分，做法为：UtilityMenu>PlotCtrls>Numbering…，弹出对话框，对图4-26中三个框选的部分按图中进行选择，点击OK，即可显示面（图4-27）。图4-25切割面图4-26以颜色显示面设置图4-27以颜色显示的面通过以上操作，即完成了所有实体模型的创建。由于在上述操作过程中，可能存在一些重合的元素，必须对此进行处理。在ANSYS中的做法是：点击Preprocessor>NumberingCtrls>MergeItems，弹出对话框，选择框中的All（图4-28），点击OK；点击Preprocessor>NumberingCtrls>CompressNumbers…，弹出对话框（图4-29），选择All，点击OK结束。图4-28消除重合的点、线、面/体图4-29压缩点、线、面的编号（从1开始且编号连续）同时，为了以后便于选择，可以通过ANSYS的工具创建面的组合，形成围岩组、衬砌组和开挖区组。ANSYS的做法为：图4-30选择面设置图4-31选择面图4-32选择的围岩区域先只选择围岩区域：UtilityMenu>Select>Entities…，弹出对话框（图4-30），按图中框选部分进行选择，点击OK；弹出对话框（图4-31），然后在工作区选择围岩区（图4-32）中的R1、R2、R3和R4，点击OK；显示面，UtilityMenu>Plot>Area，即可看到所选择的区域。创建组件：UtilityMenu>Select>Comp/Assembly>CreateComponent…，弹出对话框（图4-33），按框选部分设定，点击OK，即创建一个以ROCK命名的组件；图4-33创建组件的设置选择所有面：UtilityMenu>Select>Everything；然后显示所有面：UtilityMenu>Plot>Areas；按A中的方法选择衬砌区的面Area（图4-33）；并按B中的方法定义衬砌组Lining；图4-33选择的衬砌区的面并创建Lining组件按C、D中的方法确定开挖组件Kaiwa，图4-34。图4-34选择的开挖区的面并创建Kaiwa组件查看已创建的组件可以通过以下方式进行：UtilityMenu>Select>ComponentManager…，弹出对话框（图4-35），即可看到已定义的3个组件。图4-35已创建的组件2．直接生成面以上是通过面的下级元素（点和线）来逐步生成面，ANSYS中提供了直接创建面的方法。具体做法为：Preprocessor>Modelling>Create>Area>Rectangle（四边形面）>ByDimension，弹出对话框（图4-36），输入相应的尺寸，创建的面与通过上述方法创建的面完全一样。图4-36直接生成面下面的步骤同上面方法一致。点击或进行保存。网格划分建立了经过优化的实体模型后，即可进行平面网格的划分。1．指定面的属性在开始划分网格前，必须对面的属性（单元形式和材料号）进行指定，换句话说，在进行网格划分前，材料和单元类型必须已经定义/存在。在本问题中，隧洞开挖前，所有的面均属于岩石材料（即材料1），因此，指定所有面的属性均为Plane42单元，材料号为1。具体做法为：点击Meshing>MeshAttributes>Allareas…，弹出对话框（图4-37），按框选部分的选择进行设定（材料号1，单元类型Plane42），点击OK，即完成了面的属性指定。图4-37面属性的指定2．网格划分密度的指定尽管可以进行自由划分，但获得的网格很难满足用户的要求，因此对于面域，可以通过指定线的划分密度（长度或数目）来对面域进行映射划分Mapped（对于3边形映射划分，三条边的数目要相同，且最好为偶数；4边形时，对边的划分数目必须相等）。具体做法为：A．显示所有的线：UtilityMenu>Plot>LinesB．指定衬砌及开挖区内线的划分密度点击Meshing>MeshingTool，弹出对话框（图4-38），选择Lines>Set，弹出线选择对话框（图4-39），选择衬砌的内外圈（共8条线）以及开挖区4根线，点击OK，弹出线划分密度设定对话框（图4-40）。线密度的设定可以通过设定长度或数目来进行，这里选择数目，本模型中设定为10，点击OK。图4-38线密度划分对话框图4-39选择要设定的线图4-40线密度的指定图4-41衬砌厚度线（圆圈部分）衬砌厚度的确定也按上述方法，选择4条厚度方向的线，设定数目为2。图4-42过渡线C．过渡线的设定对图4-42中A、B、C、D四条线而言，由于我们关注隧洞周边，因此，网格密度要大一些，即网格以隧洞为圆心，朝外侧网格逐渐过渡，增大。具体做法为：按B中的方法，选择这4条线，按图4-43进行设定，点击OK，完成设定。需要注意的是，由于线的方向的原因，过渡系数可能为小数，也可能为整数，操作时可先试设定，若刚好相反，则取原过渡系数的倒数，如试取5，不合适的话，则取1/5。图4-43过渡线数目的设定需要说明的是，映射划分一般针对3边形或4边形（平面情况），若为5边形，则需要另外的映射方法，可按下面的规则进行设定，即把5边形的其中两边假象为一条边，则“变为”四边形，则按四边形划分规则进行设定（也就是说，假象为“一条边”的两条线的划分数目和要等于这两条线对应的一条边的划分数目，另外两边数目也要相等）。图4-44所有线的划分数目汇总3．网格划分在设定所有线的划分密度后，即可进行网格划分，划分分两种情况（先显示为面）：A．对于3边形或四边形点击Meshing>MeshingTool，弹出对话框（图4-45），按图中框选进行设定，点击Mesh后，弹出面选择框（图4-46），选择开挖区和衬砌区（3边形或四边形）的面后，点击OK，工作面显示为图4-47。图4-45划分工具图4-46面选择框图4-47已划分的网格A．对于5边形选择显示面UtilityMenu>Plot>Area。点击Meshing>MeshingTool，弹出对话框（图4-48），按图中框选进行设定，点击Mesh后，弹出面选择框（图4-49），选择且只能选择其中一个面（5边形），点击OK，再选择5边形的其中4个角点（图4-50，由于假象两条线为“一条线”，因此，选择时不再选择这两条线之间的交点），即可完成5边形面的划分，图4-51。图4-48划分工具图4-49面选择框图4-50选择点的顺序图4-51已划分的网格按照上面方法划分其余3个5边形面的划分，点击UtilityMenu>Plot>Elements，显示划分后的网格如图4-52。图4-52划分后的网格点击或进行保存。点击File>Saveas，在弹出的对话框中，将文件另存为Tunnel-mesh.db（用于后面的计算）。在完成有限元网格的划分后，即完成了有限元的建模，可以进行求解。通过下面进入求解模块。Mainmenu>Solution。边界条件和初始条件1．边界条件对于水工隧洞的计计算区域而言言，所选的区区域一般进行行位移边界约约束，即左右右为水平X法向约束，底底部为铅直YY法向约束。AANSYS的的做法为：DefineLLoads>>Applyy>Struucturaal>Dissplaceement>>OnNoodes（施施加在节点上上），弹出节节点选择对话话框，选择BBox窗选方方式，选择左左右边界的节节点（图4--54），点点击OK，弹出对对话框（图44-55），选选择X方向约束，点点击OK，图4-53节点选选择框图4-54被选择的节节点图4-55指定约约束方向UXX图4-56被约束的左左右边界同样，窗选底部节节点，指定UUY约束，完完成底部边界界条件设定。图4-57指定约约束方向UYY图4-58被约束的底底部边界2．初始条件初始条件包括初始始应力状态，本本问题中主要要是自重应力力情况，即要要设定重力加加速度。DefineLLoads>>Applyy>Struucturaal>Ineertia>>Graviity>Gllobal，弹弹出对话框，在在铅直Y方向输入9..8，点击OK，完成重重力加速度的的设定。图4-59重力加加速度的设定定点击或进行保存。水工隧洞施工过程程模拟在本问题，隧洞的的施工模拟包包括三个过程程，即初始应应力状态（荷荷载步1）→毛洞开开挖（荷载步步2）→毛洞支支护（荷载步步3）。本节中中按这三个步步骤进行计算算模拟。（保保证你处于SSolutiion求解模模式状态）由于在下面的分析析中涉及到单单元的生（开开挖）或死（激激活，如开挖挖回填、施作作衬砌等），因因此必须打开开N－R模式，具体体做法为：AAnalyssisTyype>AnnalysiisOpttions（分分析选项）（有有时在点击AAnalyssisTyype后，不不能发现AnnalysiisOpttions，此此时可点击SSolutiion>UnnabriddgedMMenu，即即可显示出来来），弹出对对话框（图44-60），选选择FulllN-R选选项，点击OOK即完成设设置。图4-60FulllN-RR的设置初始状态模拟初始状态是指研究究对象所在的的地质区域未未受人工扰动动（开挖、支支护等人类活活动）之前就就已客观存在在的状态，包包括自重应力力场和构造应应力场，其中中以自重应力力场为主。初初始状态的求求解过程为：：A．设定初始状态分分析的荷载步步Loadsstep=11，具体做法法为：点击AnalyssisTyype>Sool’Contrrols，弹弹出对话框（图图4-61），按图中框框选设定，点点击OK，完成荷荷载步设置；；图4-61荷载步步的设置B．求解，Solvve>CurrrentLs，点击击OK，开始求求解，当出现现下列界面时时，证明求解解成功，点击击Closee，完成求解解。图4-62经过上面过程，即即完成了初始始状态的计算算。点击或进行保存。毛洞开挖工况模拟拟毛洞成洞过程实际际包括两个部部分，即衬砌砌部分和开挖挖后不回填部部分，在本步步中，衬砌和和开挖区都会会被挖去，在在ANSYSS程序中成为为“杀死Kiill”。具具体过程为（显显示界面为单单元）：A．设定初始状态分分析的荷载步步Loadsstep=22，具体做法法为：点击AnalyssisTyype>Sool’Contrrols，弹弹出对话框（图图4-63），按图中框框选设定，点点击OK，完成荷荷载步设置；；图4-63第2荷荷载步的设置置B．根据前面4.22.2中创建建的三个组件件，选择要开开挖的衬砌区区和开挖区单单元。点击UltiliityMeenu>Seelect>>CompoonentManagger，弹出出下面对话框框，先在下面面的列表中选选择Kaiwwa和Lininng两个组件件，再点击上上面的按钮后后，关闭该对对话框。图4-64显示选择的面UlltilittyMennu>Ploot>Areeas。然后后选择附着于于所选面的单单元。UltilityyMenuu>Seleect>Enntitiees…，在弹出的的对话框中，按按下图设定后后，点击OKK，选择显示示单元UlttilityyMenuu>Plott>Elemment，如下图(图4-64))，即为本步步中要开挖的的单元。图4-64选择要要开挖的单元元C．在命令行中输入入Ekilll,all，杀杀死/开挖所选的的单元。D．选择全部元素：：UtiliitiesMenu>>Selecct>eveerthinng；并显示示右键>Reploot；E．进行毛洞开挖求求解Solvve>CurrrentLs，点击击OK，开始求求解，当出现现下列界面时时，证明求解解成功，点击击Closee，完成求解解。图4-65求解成成功显示为单元UtiilitieesMennu>Ploot>Eleementss，点击或进行保存。毛洞支护工况模拟拟毛洞稳定到一定状态态后，即可进进行衬砌施作作，形成隧洞洞主体结构，有有限元对此过过程的模拟过过程为（显示示为单元UttilitiiesMeenu>Pllot>Ellementts）：A．指定荷载步Loaadstepp＝3；点击AnalyssisTyype>Sool’Contrrols，弹弹出对话框（图图4-66），按图中框框选设定，点点击OK，完成荷荷载步设置；；图4-66第3荷荷载步的设置置B．选择先前设定的衬衬砌区域/面Lininng及相应附附着在面上的的单元；点击UltiliityMeenu>Seelect>>CompoonentManagger，弹出出下面对话框框，先在下面面的列表中选选择Liniing两个组组件，再点击击上面的按钮钮后，关闭该该对话框。图4-64选择衬衬砌区域显示选择的面UlltilittyMennu>Ploot>Areea。然后选选择附着于所所选面的单元元。UltilityyMenuu>Seleect>Enntitiees…，在弹出的的对话框中，按按下图设定后后，点击OKK，选择显示示单元UlttilityyMenuu>Plott>Elemment，如下图(图4-64))，即为本步步中要开挖的的单元。图4-64选择要要开挖的单元元C．激活所选的衬砌单单元，在命令令行输入Eaalive,,all，激激活/施作衬砌区单元；4、由于原来的衬砌砌部分设定为为岩石材料，因因此需要改变变衬砌单元材材料属性为混混凝土属性，具具体为：选定定衬砌单元，Solution>Loadstepopts>Other>Changematprops>Changematnum/弹出菜单，材料编号改为2：图4-67衬砌材材料性质的改改变D．选择全部元素：：UtiliitiesMenu>>Selecct>eveerthinng；并显示示右键>Reploot；E．进行衬砌支护后的的求解Solvve>CurrrentLs，点击击OK，开始求求解，当出现现下列界面时时，证明求解解成功，点击击Closee，完成求解解。图4-68求解成成功显示为单元UtiilitieesMennu>Ploot>Eleementss，点击或进行保存。计算结果查看及处处理通过以上分析，对对水工隧洞施施工期围岩及及衬砌结构的的应力、变形形规律进行了了简单的数值值模拟，下面面就数值模拟拟的结果进行行简单的查看看与处理。首首先通过以下下过程进入通通用后处理模模块：MaiinMennu>GenneralPostpproc。所所有的后处理理操作都必须须在后处理状状态中。1．计算结果的读取取首先查看计算荷载载步的情况：：GenerralPoostProo>ResuultsSSummarry，弹出对对话框。图4-69计算结结果汇总图4-69结果读取方方式计算结果的读取可可以通过多种种方式，以BByPicck为例：点击ReadRResultts>ByPick，弹弹出菜单后，选选择某荷载步步，点击Reead（读取取）→Cloose即可。图4-70通过BByPicck读取结果果2．计算结果的查看看在读取某荷载步的的结果后，就就可以对某荷荷载步下的结结果进行绘图图或文本列表表。计算结果果分两种：节节点解和单元元解。节点解： PlootRessults>>CoutoourPllot（绘制制等值线）>>NodallSoluu（节点解） ListRessults>>NodallSoluution（列列表结果）单元解： PlootRessults>>CoutoourPllot（绘制制等值线）>>ElemeentSoolu（单元元解） ListRessults>>ElemeentSoolutioon（列表结结果）以荷载步2（毛洞洞开挖工况）的的结果为例：：（首先对开挖掉的的单元不选择择显示，先选选择组件，然然后选择附着着在上面的单单元，开挖掉掉的单元就不不显示了）A．绘制节点位移UUY点击PlotRResultts>CouutourPlot>>NodallSoluu，弹出的对对话框中，选选择DofSoolutioon→Y-CommponenntofDispllacemeent，并作作相应设置后后，绘出UYY等值线图（如如图4-711）。图4-71绘制UUY位移分量量等值线图同理可以绘制UXX和合位移。B．绘制节点应力向向量点击PlotRResultts>CouutourPlot>>NodallSoluu，弹出的对对话框中，选选择Streess→Y-CommponenntofStresss，并作相相应设置后，绘绘出SY等值线图图（如图4--72），可可看到两侧洞洞壁出现最大大的竖直方向向压应力。图4-72绘制SSY位移分量量等值线图同理可以绘制其它它应力分量，包包括第1主应力（1sstPriicinpaalStrress）和和第3主应力（3sstPriicinpaalStrress）注意：ANSYSS的应力大小小与方向规定定和材料力学学规定相同，与与我们工程上上的规定刚好好相反，也就就是说ANSSYS中第1主应力实际际为工程上应应用的最小主主应力，第33主应力实际际为工程上应应用的最大主主应力。若要绘制某些特定定部位的计算算结果，则只只需选定关注注区域的单元元进行绘制即即可。若要对比开挖前后后的变形图，则则在对话框中中选择变形后后图形同时显显示变形前边边界。图4-73绘制变变形前后对比比图3．荷载工况组合在前面的分析中，隧隧洞的开挖经经历了3个步骤，在在每一步的分分析中均包含含了自重应力力和自重位移移，对于应力力而言，我们们不仅关注自自重应力，同同时关注开挖挖后的附加应应力和二次应应力；而对于于位移，自重重位移实际是是不存在的（经经过几百万年年的稳定后，位位移基本消失失），因此对对于开挖而言言，我们更关关注于附加位位移，这就涉涉及到荷载工工况之间的逻逻辑运算（加加或减）的问问题，也就是是荷载工况组组合。在本问问题中，评价价开挖引起的的位移或应力力时，就必须须“减掉”自自重位移或位位移，即荷载载步2要减去荷载载步1；而在评价价衬砌施作引引起的附加应应力或位移时时，则必须减减去开挖和自自重引起的位位移或应力，即即荷载步3要减去荷载载步2。荷载工况况的定义和运运算做法为：：A．定义荷载工况：：GenerralPoostprooc>LoaadCasse>CreeatLooadCaase…，弹出对话话框中，选择择从结果文件件ResulltsFiiles创建建荷载工况，图4-74定义荷荷载步1为荷载工况况1同理分别定义荷载载步2为荷载工况况2（图4-755），荷载步步3为荷载工况况3（图4-766）。图4-75定义荷荷载步2为荷载工况况2图4-76定义荷荷载步3为荷载工况况3当要分析开挖引起起的附加位移移和应力或其其它参量时，则则要减去初始始状况的位移移和应力，即即用荷载工况况2“减去”荷荷载工况1，做法为：：A．读取荷载工况22（图4-777）LoadCasse>ReaadLoaadCasse…，在弹出的的对话框中，输输入2，点击OK；图4-77读取荷荷载工况B．减去荷载工况11（图4-788）LoadCasse>Subbtractt…，图4-78减去荷荷载工况重新绘制UY方向向的附加位移移，则图形变变为图4-779，隧洞顶顶部发生向下下位移，而隧隧洞底部发生生向上的回弹弹位移。附加加的最大压应应力变为图44-80，隧隧洞底部和顶顶部发生拉应应力，两侧发发生最大压应应力。图4-79附加竖竖直UY位移图4-80附加最大主主应力4．路径操作对于实际工程而言言，我们可能能关注于应力力或位移在某某些区域或范范围内的变化化趋势，这就就需要用到AANSYS的的路径操作。以以本问题为例例，如果我们们需要知道应应力或围岩沿沿径向的变化化趋势，则可可按下面步骤骤进行操作：：A．定义路径GeneralPostpproc>PPathOOperattions>>DefinnePatth>Bynodess（通过节点点定义）…，弹出节点点选择对话框框，依次选择择路径上的节节点（图4--81），点点击OK，弹出对对话框（图44-81），定定义路径的名名称为XX，图4-81路径的的定义B．映射数据到指定定路径PathOpeeratioons>MaapOnttoPatth…，弹出对话话框，指定需需要映射的参参量，如UXX，UY，SX，SY，S1，S3等，如图图4-82。图4-82映射数数据C．绘制指定参量沿沿路径的变化化曲线PathOpeeratioons>PllotPaathIttem>OnnGrapph…，在弹出的的对话框中，选选择要绘制的的变化曲线（可可以多选，但但最好位移和和应力要分开开），图4--82。图4-82绘制位位移沿路径XX的编号曲曲线图4-83绘制主主应力沿路径径XX的编号曲曲线水工隧洞运行期模模拟前面对水工隧洞施施工期的开挖挖过程进行了了简单的模拟拟，分析了开开挖过程围岩岩及衬砌结构构的应力应变变规律。本节节对水工隧洞洞在运行期，受受内、外水压压力作用下的的衬砌变形及及应力规律。运行期内水压力的的模拟对于隧洞的内水压压力，需要区区别无压隧洞洞和有压隧道道（图4-884）。因此，施施加内水压力力分为有压和和无压两种情情况，有限元元模型采用与与前面开挖支支护后的隧洞洞，并进入求求解状态Soolutioon。A．有压隧洞B．无压隧隧洞图4-84无压和和有压隧洞1．有压隧洞有压隧洞情况下，作作用与隧洞衬衬砌内表面的的水压力可认认为是常量，可可以通过ANNSYS施加加Presssure命令令，直接进行行施加。A．选定衬砌单元且且选择附着于于单元的节点点由于衬砌材料号为为2，因此通过过如下方式进进行选择。UtilitieesMennu>Sellect>EEnitiees…，在弹出的的对话框中，选选择Elemment→ByAtttribuutes→MaterrialNNum，在Min,Maax,Incc空白框中输输入2，点击Appply→Plot→OK。选择附着于衬砌单单元上的节点点。UtilitieesMennu>Sellect>EEnitiees…，在弹出的的对话框中，选选择Nodees→Attacchedtto→Elemeent，点击击Applyy→Plot→OK。图4-85选择衬衬砌单元图4-86选择择附着于衬砌砌单元的节点点B．选择所有衬砌内内表面节点并并施加水压力力Solutionn>DefiineLooads>AApply>>Struccturall>Presssure>>Onnoodes，在在弹出的对话话框中，采用用圆形Cirrcle选择择方式（图44-87），选选择所有的衬衬砌内表面节节点，点击OOK，在弹出出的对话框中中，输入内水水压力1e66，点击OK完成水压压力设定（图图4-88）。图4-87圆形选择择方式图4-88设定内水压压力此时，施加的水压压力表现为一一条线，这里里可以通过设设置，将Prressurre以箭头Arrrow表示。UtilitieesMennu>PlootCtlss>Symbbol…，在弹出的的对话框中，按按图4-899进行设置后后，即可将PPressuure以箭头头表示。模型型上显示的水水压力如图44-90。图4-89水压力表表示方式图4-90模型上显显示的水压力力2．无压隧洞假定水位在半圆高高