6稳定渗流问题的分析.doc

PAGE  PAGE 6 土坝稳定渗流问题的分析 一、问题和条件： 土坝稳定渗流分析，采用生死单元技术，迭代计算浸润线，根据前次计算结果，不断修改单元的浸润逸出点位置，直到满足精度要求。本算例的模型比较简单，对均质矩形土体稳定渗流问题进行计算，上游水头为+6m，下游水头为+1m，坝体宽度为4m，坝体高度为6m。试做渗流分析。 渗透系数k=0.02 二、进入ANSYS界面 单击 开始 → 程序→ansys8.0→configure ANSYS product 然后在File Management中定义Working Directory（工作路径） 如：class； 定义Job Name（工作文件名）。如：calculate seepage of earthdam。 三、定义单元及材料 1 定义程序中需要调用的数组 运行main中parameters→array parameters→define/edit→add，接着在对话框中par选项中输入数组名字:nodefour，在Type中选择Array，在I、J、K中对应输入4、1、1，然后点击OK按钮。 Smart,off Antype,static NROPT,FULL OUTPR,BASIC,ALL 2 新建单元类型 运行主菜单 Preprocessor → Element Type → Add/Edit/Delete (新建/编辑/删除单元类型)命令，接着在对话框中单击“Add”按钮新建单元类型。 3 定义单元类型 先选择单元为Thermal solid，接着选择Quad 4node 55，然后单击“OK”按钮确定，完成单元类型的选择。 4 关闭单元类型对话框 回到单元类型对话框，已经新建了Thermal solid单元，单击对话框中“Close”按钮关闭对话框。 5设置材料属性 运行主菜单 Preprocessor → Material props→ Material Models (材料属性)命令。选择材料属性命令后，系统会显示材料属性设置对话框，接着在对话框Thermal中选择conductivity → Isotropic →在kxx中输入渗透系数0.02。点击OK按钮。 6 关闭材料属性对话框 完成材料属性的设置后，可在对话框右上方单击，关闭材料属性设置。 四、建立框架几何图形 1 生成节点 运行主菜单 Preprocessor → Modeling → Create → Nodes → In Active CS 生成节点于目前坐标系统命令。 接着在Node number 输入1，接着依照顺序输入第一点XYZ的坐标值为0,0,0，完成输入后单击“Apply”按钮继续生成第二点。 如上依次 2 0,6,0 3 4,6,0 4 4,0,0 2 生成线条 运行主菜单 Preprocessor → Modeling →Create →lines →lines →Straight line，依次选择节点1、节点2???点击apply; 选择节点2、节点3，点击apply; 选择节点3、节点4，点击apply; 选择节点4、节点1，点击apply。 3 生成面单元 运行主菜单 Preprocessor → Modeling → Create →Areas →arbitrary → By lines生成单元命令。 依次选择线1、线2、线3、线4，完成输入后单击“OK”按钮生成面单元。 4 按等分每条线为32等份的格式分网格 运行主菜单 Preprocessor → Meshing → Size cntrls →Manualsize →lines → All lines，生成一个对话框，在Ndiv 对应项输入每条线的等分数目，该例题中输入32，点击ok按钮。 5 划分网格 运行主菜单 Preprocessor → Meshing → Meshtool，生成一个对话框依次选择Quad、free，点击Mesh按钮。网格划分完毕。 6 显示单元和节点编号 运行下拉菜单中 PlotCtrls → Numbering 编号显示命令。 接着在编号控制对话框中将选择打开 Node number和 Element number 显示单元编号，完成后单击“OK”按钮。此时系统将会以颜色区分单元编号。 五、施加荷载并进行分析 *get, n_max, node, , num, max !获取最大节点数目，用于计算时循环使用 *get, e_max, elem, , num , max !获取最大单元数目，用于计算时循环使用 *dim, n_temp, array, n_max !将最大节点数目赋值给n_temp变量 *dim, n_pre, array, n_