当前位置：文档库 › 聚能射流形成及破甲过程的数值模拟分析

聚能射流形成及破甲过程的数值模拟分析

数值模拟报告(DOC)

第一部分：数值模拟技术研究文献综述浅析数值模拟技术 1．引言近年来，随着我国大规模地进行“西部大开发”和“南水北调”等巨型工程，越来越多的岩土工程难题摆在我们面前，单纯依靠经验、解析法显然已不能有效指导工程问题的解决，迫切需要更强有力的分析手段来进行这些问题的研究和分析。自R.W. Clough 上世纪60年代末首次将有限元引入某土石坝的稳定性分析以来，数值模拟技术在岩土工程领域取得了巨大的进步，并成功解决了许多重大工程问题。特别是个人电脑的普及及计算性能的不断提高，使得分析人员在室内进行岩土工程数值模拟成为可能。在这样的背景下，数值模拟特别是三维数值模拟技术逐渐成为当前中国岩土工程研究和设计的主流方法之一，也使得岩土工程数值模拟技术成为当今高校和科研院所岩土工程专业学生学习的一个热点。采用大型通用软件对岩土工程进行数值模拟计算，在目前已成为项目科研、工程设计、风险评估等岩土类项目的必须，学习和掌握Ansys、FLAC3D、UDEC 等数值计算软件已成为学校、科研院所对工程从业人员的基本要求。数值模拟方法主要有限元法、边界元法、加权余量法、半解析元法、刚体元法、非连续变形分析法、离散元法、无界元法和流形元法等，各种方法都有其对应的软件。 2.数值模拟的发展趋势可以说, 继理论分析和科学试验之后, 数值模拟已成为科学技术发展的主要手段之一。随着软件技术和计算机技术的发展, 目前国际上数值模拟软件发展呈现出以下一些趋势：（1）. 由二维扩展为三维。早期计算机的能力十分有限，受计算费用和计算机储存能力的限制，数值模拟程序大多是一维或二维的，只能计算垂直碰撞或球形爆炸等特定问题。随着第三代、第四代计算机的出现, 才开始研制和发展更多的三维计算程序。现在，计算程序一般都由二维扩展到了三维，如LS-DYNA2D 和LS - DYNA3D、AUTODYN2D 和AUTO-DYN3D。（2）.从单纯的结构力学计算发展到求解许多物理场问题。数值模拟分析方法最早是从结构化矩阵分析发展而来，逐步推广到板、壳和实体等连续体固体力学分析，实践证明这是一种非常有效的数值模拟方法。近年来数值模拟方法已发展到流体力学、温度场、电传导、磁场、渗流等求解计算，最近又发展到求解几个交叉学科的问题。例如内爆炸时，空气冲击波使墙、板、柱产生变形，而墙、板、柱的变形又反过来影响到空气冲击波的传播，这就需要用固体力学和流体动力学的数值模拟结果交叉迭代求解。（3）.由求解线性问题进展到分析非线性问题。随着科学技术的发展，线性理论已经远远不能满足设计的要求。诸如岩石、土壤、混凝土等，仅靠线性计算理论就不足以解决遇到的问题，只有采用非线性数值算法才能解决。众所周知，非线性的数值模拟是很复杂的，它涉及到很多专门的数学问题和运算技巧，很难为一般工程技术人员所掌握。为此，近年来国外一些公司花费了大量的人力和资金，开发了诸如LS- DYNA3D、ABAQUS和AU-TODYN等专长求解非线性问题的有限元分析软件，并广泛应用于工程实践。这些软件的共同特点是具有高效

数值模拟

有限元软件在焊接热应力分析中的应用和发展随着焊接温度场、应力场和变形的深入研究，有限元技术的发展与应用，以及近年来由于计算机技术的突飞猛进，目前在进行有限元分析时所用的软件方面已经有了不少优秀的计算分析软件，如:ANSYS,ABAQUS,ADINA,NASTRAN,MARC,SYSWBLD等可供焊接工作者选用[31。我国目前尚不具各开发大型通用有限元软件的条件，没有自主版权的商品化有限元软件，所以我国的有限元发展途径主要是使用、扩充和改进从国外引进的某些有限元软件。这些现有的有限元软件具有自动划分有限元网格和自动整理计算结果，并使之形成可视化图形的前后处理功能。因而，焊接工作者己经无需自己从头编制分析软件，可以利用上述商品化软件，必要时加上二次开发，即可以得到需要的结果，这就明显地加速了焊接模拟技术发展的进程。在国内还很少利用通用有限元软件分析焊接结构应力场的例子。中科院的颜抬霞[231等利用ANSYS对球壳焊接瞬态温度场、应力场进行模拟取得较好的结果。清华大学的鹿安理等利用MARC软件，开发专用用户子程序，使网格自适应技术更趋完善，并用于厚板焊接过程的三维数值模拟，取得了很好的效果，并在模型上利用相似原理及简化热源模型等技术问题进行探讨，提出未来焊接数值模拟应重点研究的几个问题。清华大学的蔡志鹏等人，利用MARC软件，简化热源模型，用串热源模型代替高斯热源进行焊接应力和变形的分析，但其实际分析例子只是进行切割变形分析，其方法的有效性还须进一步验证[241[251[26][271 SYSWELD2．1SYSWELD的特点 SYSWELD的开发最初源于核工业领域的焊接工艺模拟，当时核丁业需要揭示焊接工艺中的复杂物理现象，以便提前预测裂纹等重大危险。随着应用的发展，SYSWELD逐渐扩大了其应用范围，并迅速被汽车工业、航空航天、国防和重型工业所采用。 SYSWEI。D完伞实现了机械、热传导和金属冶金的耦合计算，允许考虑晶相转变及同一时间晶相转变潜热和晶相组织对温度的影响。在具体计算中分两步进行。首先实现温度和晶相组织的计算，然后进行机械力的计算。在机械力计算中，已经充分考虑了第一步计算的结果，如残余应力和应变的影响。 SYSWELD的电磁模型允许模拟点焊和感应加热，并可实现能量损失和热源加载的计算模拟。SYSWELD扩散与析出模型可实现渗碳、渗氮、碳氮共渗模拟，先计算化学元素的扩散和沉积，然后再考虑对热和机械性能的影响。SYS肌IJ)的氢扩散模型能计算模拟氢的浓度，预测冷裂纹的严重危害旧。 2．2 SYSWELD的应用焊接残余应力是焊接过程中影响构件强度和寿命的主要囚素之一，通过计算机仿真分析可准确分析焊接时温度场、应力场的变化规律，焊接时构件的变形情况[Io-11J。利用SYSWELD 对某零件进行激光焊接仿真分析可得出一些结果云图，由这些云图能够判断构件在焊接过程中瞬态温度场的变化情况、构件中的应力以及焊接完成后残余应力的分布情况。最重要的是利用SYSWELD软件能够方便准确地分析焊接过程中材料金相组织的转化情况，为激光动态焊接过程数值仿真提供理论基础。在焊接中，热膨胀与收缩伴随着金相转变而发生，从而导致焊接过程和焊后焊件的结构变形。在SYSWELD中这些囚素都可以进行模拟；焊接过程的热效应引发热应力(结构膨胀与收缩)，同样可以利用SYSWELD进行评估；通过SYSWELD，可以对焊接的内应力和金相结构进行预测，并将这些结果直接应用于产品寿命计算与分析。在许多工业上，经常将等厚或不等厚的材料焊接在一起后进行冲压，这些材料可以相同也可以不同。采用SYSWELD把焊接模拟的内应力和冶金特性作为初始条件加载到冲压模拟软件中，如PAM—STAMP就是专业的冲压模拟软件。

工程材料与材料成型技术教案

教案（理论课） 2010～2011学年第2学期课程名称工程材料与成形技术基础教学系机械工程系授课班级焊接091 主讲教师晏丽琴职称讲师

培黎工程技术学院二○一一年二月课程基本情况

系主任：年月日目录第一章绪论第一节材料加工概述一、材料加工概述二、材料加工的基本要素和流程第二节材料成形的一些基本问题和发展概况一、凝固成形的基本问题和发展概况二、塑性成形的基本问题和发展概况三、焊接成形的基本问题和发展概况四、表面成形的基本问题和发展概况第三节本课程的性质和任务绪论学习思考问题 ·材料加工的基本要素和流程是什么？ ·材料成形存在的基本问题是什么? ·本课程的性质和基本任务是什么? 一、材料加工概述任何机器或设备，都是由许许多多的零件装配而成的。这些零件所用材料有金属材料，也有非金属材料。零件或材料的加工方法多种多样，归纳起来有以下4类： (1)成形加工：用来改变材料的形状尺寸，或兼有改变材料的性能。主要有凝固成形、塑性成形、焊接成形、粉末压制和塑料成形等。 (2)切除加工：用于改变材料的形状尺寸，主要有车、铣、刨、钻、磨等传统的切削加工，以及直接利用电能、化学能、声能、光能进行的特殊加工，如电火花加：[、电解加工、超声加工和激光加工等。 (3)表面成形加工：用来改变零件的表面状态和(或)性能，如表面形变及淬火强化、化学热处理、表面涂(镀)层和气相沉积镀膜等。

(4)热处理加工：用来改变材料或零件的性能，如退火、正火、淬火和回火等。根据零件的形状尺寸特征、工作条件及使用要求、生产批量和制造成本等多种因素，选择零件的加工方法，以达到技术上可行、质量可靠和经济上合理。零件制成后再经过检验、装配、调试，最终得到整机产品。二、材料加工的基本要素和流程材料加工方法的种类虽然繁多，但通过对每种材料加工方法的过程分析表明，它们都可以用建立在少数几个基本参数基础上的统一模式来描述。该模式便于对各种加工方法进行综合分析和横向比较。任何一种材料的加工过程，都是为了达到材料的形状尺寸或性能的变化。而为了产生这种变化，必须具备三个基本要素：材料、能量和信息（图1.2）。因而材料的加工过程，可以用相关材料流程、能量流程和信息流程来描述。三大流程： 1.材料流程表征加工过程特点的类型；要改变形状尺寸和性能的材料状态；能够用来实现这种形状尺寸和性能变化的基本过程； 2.能量流程包括机械过程的能量流程，热过程能量：电能、化学能、机械能 3.信息流程形状信息、性能信息

数值模拟分析实例

华中科技大学体育馆数值模拟分析 6.1分析模型的建立采用有限元软件ANSYS建立该网壳结构有限元分析模型。整体屋盖结构共计1481个节点，4430个单元，16种截面类型。建模时，网壳结构主体结构部分（包括主桁架、次桁架、水平支撑和檩条）采用ANSYS的LINK8杆单元建模，两侧翼的主梁、次梁和支承钢管柱均采用BEAM4梁单元，网壳结构屋面下部混凝土支承结构亦采用BEAM4梁单元。分析时，屋面板、设备管线等荷载等效为节点荷载，施加在结构节点上。在网壳结构有限元分析中，对于杆件采用的LINK8 3-D Spar单元为三维单元，假设材料为均质等直杆，且在轴向上施加载荷，可以承受单向的拉伸或者压缩，每个节点上具有三个自由度，即沿X、丫和Z坐标轴方向。该单元具有塑性、蠕变、应力硬化和大变形等功能，能较好的模拟三维空间桁架单元。对于两侧翼结构和下部支撑体系的柱、梁等结构采用的BEAM4单元是一个轴向拉压、扭转和弯曲单元，每个节点有三个平动自由度和三个转动自由度，具有应力刚化和大变形功能。施工过程模拟分析时考虑时，同时考虑温度效应影响，计算时材料假定为理想弹塑性材料。

6.2分析工况选取按照实际施工顺序，将网壳结构屋盖施工过程划分为 5个工况进行施工数值模拟，计算温度取为该阶段施工完成时的环境温度。工况1: 7榀拱形主桁架安装完毕，但临时支撑未撤除，计算温度为温度15C ； (a )短轴立面 1' 裁灌 1::“ 麻 MM ■■号 -? ■?■ 水 * S-ES-B- Kir *￡靈曲 r-…；： 1 ;u * 图6-1有限元分析模型

图6-4工况2中屋盖结构平面图图6-5工况2中屋盖结构立面图工况3：次桁架、水平支撑及楼梯安装完毕，临时支撑拆除，计算温度为29C; 图6-6工况2中屋盖结构平面图 (b )长轴立面图6-7工况2中屋盖结构立面图工况4：檩条及设备管线安装完毕，计算温度为41C; 完成后拆除其临时支撑，计算温度为8C； (a)短轴立面 (b)长轴立面 (a)短轴立面

材料成型数值模拟设计实验

学生学号实验课成绩学生实验报告书实验课程名称材料成型数值模拟设计实验开课学院材料学院指导教师姓名学生姓名学生专业班级成型1001班 2012-- 2013学年第二学期

实验教学管理基本规范实验是培养学生动手能力、分析解决问题能力的重要环节；实验报告是反映实验教学水平与质量的重要依据。为加强实验过程管理，改革实验成绩考核方法，改善实验教学效果，提高学生质量，特制定实验教学管理基本规范。 1、本规范适用于理工科类专业实验课程，文、经、管、计算机类实验课程可根据具体情况参照执行或暂不执行。 2、每门实验课程一般会包括许多实验项目，除非常简单的验证演示性实验项目可以不写实验报告外，其他实验项目均应按本格式完成实验报告。 3、实验报告应由实验预习、实验过程、结果分析三大部分组成。每部分均在实验成绩中占一定比例。各部分成绩的观测点、考核目标、所占比例可参考附表执行。各专业也可以根据具体情况，调整考核内容和评分标准。 4、学生必须在完成实验预习内容的前提下进行实验。教师要在实验过程中抽查学生预习情况，在学生离开实验室前，检查学生实验操作和记录情况，并在实验报告第二部分教师签字栏签名，以确保实验记录的真实性。 5、教师应及时评阅学生的实验报告并给出各实验项目成绩，完整保存实验报告。在完成所有实验项目后，教师应按学生姓名将批改好的各实验项目实验报告装订成册，构成该实验课程总报告，按班级交课程承担单位（实验中心或实验室）保管存档。 6、实验课程成绩按其类型采取百分制或优、良、中、及格和不及格五级评定。附表：实验考核参考内容及标准观测点考核目标成绩组成实验预习1．预习报告 2．提问 3．对于设计型实验，着重考查设计方案的科学性、可行性和创新性对实验目的和基本原理的认识程度，对实验方案的设计能力 20% 实验过程1．是否按时参加实验 2．对实验过程的熟悉程度 3．对基本操作的规范程度 4．对突发事件的应急处理能力 5．实验原始记录的完整程度 6．同学之间的团结协作精神着重考查学生的实验态度、基本操作技能；严谨的治学态度、团结协作精神 30% 结果分析1．所分析结果是否用原始记录数据 2．计算结果是否正确 3．实验结果分析是否合理 4．对于综合实验，各项内容之间是否有分析、比较与判断等考查学生对实验数据处理和现象分析的能力；对专业知识的综合应用能力；事实求实的精神 50%

梳齿受力数值模拟与分析

第五章梳齿受力数值模拟与分析梳齿是梳齿式采棉机的核心部件，是梳齿式采棉机主要的研究内容，梳齿的好坏直接影响棉花的采摘效果，由于梳齿在采摘过程中受力复杂，配合影响因素较多，因此本文主要选取梳齿在单个棉杆节点处的受力为梳齿的受力状况，主要分析梳齿在节点处的受力变形问题。通过ANSYS分析，为梳齿式采棉机采摘部分的优化设计提供理论依据和实验指导。 5.1 基于ANSYS梳齿模型的建立 5.1.1 几何模型构建及网格划分 ANSYS软件可以用于简单模型的建立，它和其他的二维、三维软件有很好的接口，因此对于复杂的模型一般都是先用PRO/E、SolidWorks、UG等三维造型软件进行建模，对于简单的模型，可以运用CAD、CAXA等二维绘图软件进行建模，然后导入ANSYS进行分析，本文中涉及的梳齿受力的结构比较简单，可以直接用ANSYS软件本身提供的建模工具进行建模。对于梳齿受力的模拟可以是二维模拟也可以是三维模拟，由于本研究对象结构简单，形状均匀，受力简单，因此可以做些简化。在建模时采用二维有限元模型并将出棉口按密闭处理，网状输送管道做静态分析，简化后的影响将在第六章分析验证。本研究主要关心的是梳齿装置受力变形的分布情况，试验证明这样的简化对梳齿影响不大。在模型生成后，要对梳齿受力区域内进行离散，即有限元网格的生成，网格划分可分为结构化网格和非结构化网格两种。采用结构化网格易于生成物面附近的边界层网格，且有较多算法和成熟的受力分析模型，但需要较长的物面离散时间，对复杂的外形必须构造不同的网络拓扑结构，无法实现网格的自动生成，费时费力。非结构化网格消除了结构化网格中节点的结构性限制，节点和单元的分布可控性好，能够很好的处理边界问题，适宜于模拟复杂外形[23]，近年来发展迅速，应用较为广泛。在进行网格划分时要考虑到分离室内的流动特性，进而对网格做适当的调整。本文对梳齿式采棉机采摘部件的分析面采用非结构化网格划分，并适当对一些边界面进行加密，尽可能保持计算区域和实际流动区域的一致性，通过加密可以捕捉到在某些截面产生的巨大的受力梯度对梳齿造成的显著影响，以保证模拟计算的精度。建立梳齿的几何模型，整个梳齿的计算区域被划分单元格，生成的模型网格。 5.1.2 主要边界条件本研究中，主要考虑梳齿在受不同种力对梳齿变形的影响，因此，设计不同

顶盖驱动流数值模拟分析

《数值传热学》作业：顶盖驱动流数值模拟分析

西安科技大学能源学院安全技术及工程申敬杰201112612

顶盖驱动流数值模拟分析顶盖驱动流作为经典的数值计算模型，常常用来考核源程序和计算思想的正确性。这种流动边界条件简单，而且不涉及模型的影响，便于直接评价差分格式的性能。 1.引言数值传热学，又称计算传热学，是指对描写流动与传热问题的控制方程采用数值方法，通过计算机求解的一门传热学与数值方法相结合的交叉学科。数值传热学的基本思想是把原来在空间与时间坐标中连续的物理量的场（如速度场，温度场，浓度场等），用一系列有限个离散点上的值的集合来代替，通过一定的原则建立起这些离散点变量值之间关系的代数方程（称为离散方程）。求解所建立起来的代数方程已获得求解变量的近似值。由于实验方法或分析方法在处理复杂的流动与换热问题时，受到较大的限制，例如问题的复杂性，即无法做分析解，也因为费用的昂贵而无力进行实验测定，而数值计算的方法正具有成本较低和能模拟复杂或较理想的过程等优点，数值传热学得到了飞速的发展。特别是近年来，计算机硬件工业的发展更为数值传热学提供了坚实的物质基础，使数值模拟对流动与传热过程的研究发挥了重要的作用。目前，比较著名的数值模拟分析应用软件有FLUENT、CFX、STAR-CD、和PHOENICS等，而FLUENT是国内外比较流行的商用CFD软件包，该软件以其市场占有率高、计算准确、界面友好、使用简单、应用领域广、物理模型多而获得较高的市场占有率和用户的肯定。 2.物理模型在一个正方形的二维空腔中充满等密度的空气，方腔每边长为0.12m，取雷诺数为Re=12000，由Re=vd/υ，方腔的当量直径d ，计算知d=0.12m，又υ=15.7 ×10 ﹣6m2/s，则顶盖驱动流的速度v=1.57m/s，即其顶板以1.57m/s的速度向右移动，同时带动方腔内流体的流动，流场内的流体为紊流。计算区域示意图如图1所示。 v=1.57m/s L=0.12m 图1 计算区域示意图

焊接接头延性损伤与断裂的数值模拟与分析

焊接接头延性损伤与断裂的数值模拟与分析目录第一章绪论 (1) 1.1 研究背景 (1) 1.2 理论基础及研究现状 (2) 1.2.1 延性损伤力学的发展 (2) 1.2.2 细观损伤力学的现状 (3) 1.2.3 延性裂纹扩展模拟研究现状 (5) 1.3 本文研究内容 (8) 第二章30Cr2Ni4MoV转子钢焊接接头GTN模型材料参数的确定 (9) 2.1 引言 (9) 2.2 实验方法 (9) 2.3 有限元计算模型 (11) 2.4 GTN模型参数的确定 (12) 2.5 结果分析 (15) 2.6 小结 (17) 第三章裂纹位置对焊接接头延性裂纹起裂和扩展的影响 (18) 3.1 概述 (18) 3.2 有限元计算模型 (18) 3.3 有限元计算结果及讨论 (19) 3.3.1 模拟的载荷位移曲线和J阻力曲线 (19) 3.3.2 延性裂纹起裂和扩展的力学分析 (21) 3.3.3 延性裂纹扩展路径 (26) 3.4 结论 (28) 第四章裂尖局部残余应力对焊接接头延性裂纹起裂和扩展的影响 (29) 4.1 概述 (29) 4.2 有限元模型 (30) 4.3 残余应力分布 (31) 4.3.1 母材、焊缝、热影响区中裂尖残余压应力分布 (31) 4.3.2 界面处残余压应力分布 (33) 4.3.3 裂尖残余拉应力分布 (33) 4.4 裂尖局部残余应力对载荷位移曲线和J阻力曲线的影响 (35)

4.4.1 局部残余压应力对载荷位移曲线和J阻力曲线的影响 (35) 4.4.2 残余拉应力对载荷位移曲线的影响 (38) 4.5 残余应力对应力、应变及损伤场的影响 (39) 4.5.1 母材裂尖局部残余压应力对裂尖前应力、应变、损伤场的影响 (39) 4.5.2 焊缝裂尖残余压应力对应力、应变、损伤场的影响 (43) 4.5.3裂尖残余压应力对焊缝和热影响区界面材料裂尖应力、应变、损伤场的影响 (46) 4.5.4母材裂尖局部残余拉应力对裂尖前应力、应变、损伤场的影响 (49) 4.6 小结 (52) 第五章残余应力对三维延性裂纹起裂和扩展的影响 (53) 5.1 概述 (53) 5.2 有限元模型 (54) 5.3 残余应力分布 (55) 5.4 残余应力对载荷位移曲线和J阻力曲线的影响 (58) 5.5 残余应力对裂纹扩展路径的影响 (60) 5.6 残余应力对应力、应变、损伤场的影响 (61) 5.7 小结 (65) 第六章总结 (67)

FLUENT算例 (3)三维圆管紊流流动状况的数值模拟分析

三维圆管紊流流动状况的数值模拟分析在工程和生活中，圆管内的流动是最常见也是最简单的一种流动，圆管流动有层流和紊流两种流动状况。层流，即液体质点作有序的线状运动，彼此互不混掺的流动；紊流，即液体质点流动的轨迹极为紊乱，质点相互掺混、碰撞的流动。雷诺数是判别流体流动状态的准则数。本研究用CFD 软件来模拟研究三维圆管的紊流流动状况，主要对流速分布和压强分布作出分析。 1 物理模型三维圆管长2000mm l =，直径100mm d =。流体介质：水，其运动粘度系数6 2 110m /s ν-=?。 Inlet ：流速入口，10.005m /s υ=，20.1m /s υ= Outlet ：压强出口 Wall ：光滑壁面，无滑移 2 在ICEM CFD 中建立模型 2.1 首先建立三维圆管的几何模型Geometry 2.2 做Blocking 因为截面为圆形，故需做“O ”型网格。

2.3 划分网格mesh 注意检查网格质量。在未加密的情况下，网格质量不是很好，如下图因管流存在边界层，故需对边界进行加密，网格质量有所提升，如下图

2.4 生成非结构化网格，输出fluent.msh等相关文件 3 数值模拟原理紊流流动

当以水流以流速20.1m /s υ=，从Inlet 方向流入圆管，可计算出雷诺数10000υd Re ν ==，故圆管内流动为紊流。假设水的粘性为常数（运动粘度系数62 110m /s ν-=?）、不可压流体，圆管光滑，则流动的控制方程如下： ①质量守恒方程： ()()()0u v w t x y z ρρρρ????+++=???? (0-1) ②动量守恒方程： 2()()()()()()()()()()[]u uu uv uw u u u t x y z x x y y z z u u v u w p x y z x ρρρρμμμρρρ??????????+++=++??????????'''''????+---- ???? (0-2) 2 ()()()()()()()()()()[]v vu vv vw v v v t x y z x x y y z z u v v v w p x y z y ρρρρμμμρρρ??????????+++=++??????????'''''????+- ---???? (0-3) 2 ()()()()()()()()()()[]w wu wv ww w w w t x y z x x y y z z u w v w w p x y z z ρρρρμμμρρρ??????????+++=++??????????'''''????+- ---???? (0-4) ③湍动能方程： ()()()()[())][())][())]t t k k t k k k ku kv kw k k t x y z x x y y k G z z μμρρρρμμσσμμρεσ????????+++=+++????????? ?+ ++-?? (0-5) ④湍能耗散率方程： 212()()()()[())][())][())]t t k k t k k u v w t x y z x x y y C G C z z k k εεμμρερερερεεεμμσσμεεεμρσ??????? ?+++=+++??????????+++-?? (0-6) 式中，ρ为密度，u 、ν、w 是流速矢量在x 、y 和z 方向的分量，p 为流体微元体上的压强。方程求解：采用双精度求解器，定常流动，标准ε-k 模型，SIMPLEC 算法。 4 在FLUENT 中求解计算紊流流动 4.1 FLUENT 设置除以下设置为紊流所必须设置的外，其余选项和层流相同，不再详述。

公路隧道数值模拟(分析讲课版)

广东省交通厅科技项目复杂地质条件下隧道施工安全保障技术研究茶林顶公路隧道初始应力状态及施工力学数值模拟目录 1 工程概况 (1) 2 工程地质条件 (1) 2.1地形地貌 (1) 2.2地质构造 (1) 2.2.1褶皱 (1) 2.2.2断层 (1) 2.3地层岩性 (1) 3 MIDAS/GTS简介 (2) 4隧道岩体应力场的数值模拟 (3) 4.1数值分析模型的建立 (3) 4.2数值模拟结果分析 (4) 4.2.1 最大主应力特征 (4) 4.2.2 最小主应力特征 (7) 4.2.3 最大剪应力特征 (9)

4.3主要结论 (12) 5隧道典型横断面施工力学数值模拟 (12) 5.1计算参数的选取 (12) 5.2数值分析模型的建立 (13) 5.3施工过程控制 (14) 5.4数值分析结果及其分析 (14) 5.3.1围岩位移特征 (14) 5.3.2围岩应力特征 (21) 5.3.3围岩屈服接近度特征 (32) 5.3.4断层带位移特征 (35) 5.3.5断层带应力特征 (41) 5.3.6断层带屈服接近度特征 (50) 5.3.7隧道初期支护结构内力及应力特征 (53) 5.5主要结论 (67) 6 结论和建议 (67)

1 工程概况广梧高速公路茶林顶公路隧道左线起点里程LK71+566，终点里程LK74+261，全长2695m；右线起点里程RK71+632，终点里程LK74+246，全长2614m。为双洞四车道，左、右线隧道分离布设，设计行车速度为80km/h。 2 工程地质条件 2.1地形地貌隧道地处茶林顶重丘山岭区，山体走向总体呈近北东或北西向，地势总体呈南高北低，隧道线路经过最大高程约为355m，隧道进出口丘山体呈缓坡状，自然坡度为10°~20°，隧道中部山顶及山凹两侧山坡坡度较大，约30°~35°，山体植被茂密，主要生长松树和杂草，山体地表发育有数条小沟谷，部分沟谷内有长年流水，地表水量较小，隧道中部为一较大沟谷（分水凹），呈北东方向，平时无水流，但大雨时水量较大。 2.2地质构造 2.2.1褶皱根据地质填图岩性组合分析对比，隧道区存在一背斜褶皱构造，其轴部为泥盆系东岗岭组地层，两翼为泥盆系榴江组地层，为一向北西倾覆背斜构造。 2.2.2断层分布于郁南茶林顶F7断层，地貌上表现为沟谷，推荐线路大致在右线K72＋980处遇该断层，影响带宽20～30m。断层走向NE40?～45?，倾向NW，倾角75?~80?。受断层影响，中泥盆统东岗岭组白云质灰岩破碎，形成构造角砾岩和密集节理带，地表沟谷中有泉水溢出。钻孔岩心显示断层角砾呈棱角、次棱角状，为方解石脉胶结，脉中晶洞及自形方解石发育，反映其晚期活动为张性和正断层特征。 2.3地层岩性第四系覆盖层主要为亚粘土，下伏基岩主要为上泥盆统榴江组砂岩和中泥盆统东岗岭组白云质灰岩，其褶皱相对发育，从地表看，LK71+730(RK71+710)～LK72+950(RK72+940)，LK73+620(RK73+605)～LK74+000(RK73+960)为上泥盆统榴江组砂岩、页岩、石英砂岩，其余为中泥盆统东岗岭组白云质灰岩、灰岩。 1、耕植土（Q pd）：灰褐色，湿，可塑状，主要由亚粘土组成，局部可见有植物根系，偶

FLAC-数值模拟分析

1工程概况研究此段为中条山隧道K9+450～K10+560段，此处隧道最大埋深约540m，主要由太古界涑水群表壳岩组合解州片麻岩（Hgn）地层组成，构成中条山隧道分水岭北侧的主体；此段组成隧道的围岩岩性主要为变粒岩、花岗片麻岩等。该套地层岩性复杂，组合无规律。岩层产状整体倾向南东，倾角一般在50ｏ～70ｏ间变化。在AK9+900～AK10+000段为区域性断层影响段，此断层为破碎岩石组成，将为基岩裂隙水下渗提供通道，隧道开挖必将引起涌水，同时此段围岩稍差，施工时易引起坍塌。此段同时也是中条山北侧泉水主要涌出段，水文情况复杂。总体评价，本段工程地质条件差。在此处，具体运用FLAC3D进行模拟的区段均取洞身YK10+100～YK10+180段。该区段为V级围岩区域，埋深为505～512m，为断层，附近太古界涑水群花岗片麻岩、黑云斜长片麻岩、岩石破碎。隧道断面为SVc型，如图2-2所示。图2-2SVc型隧道断面图隧道衬砌按新奥法原理设计，采用SVc型复合式衬砌，该衬砌适用于隧道洞身V级断层影响带及软弱破碎围岩段的初期支护及衬砌，超前支护各环采用42×4mm注浆小导管超前预加固围岩，长4.5m，环向间距35cm，搭接长度1.3m，斜插角10ｏ～15ｏ，每环37根；初衬以喷、锚、网为主要支护手段：钢拱架为I20a型钢，纵向间距75cm，每榀钢拱架之间采用φ22钢筋连接，环形间距1.0m；锚杆采用D25中空注浆锚杆，长3.5m，间距75cm（纵）×100cm（环），与钢拱架交错布置；喷C25早强混凝土26cm。二次衬砌和仰拱均为C30钢筋混凝土结构，厚50cm。 1.2数值计算模型根据中条山隧道工程的实际状况，为提高计算速度，在保证计算精度的前提下，取桩号YK10+100～YK10+160段采用大型有限差分软件FLAC3D进行建模分析。对于全断面法、预留核心土法、台阶法，由于整个隧道模型左右对称，为减少计算量，可取隧道模型的一半计算，隧道的计算模型I如图2-3所示。图2-3计算模型I示意图 x、y、z各方向的长度分别为60m、60m和140m。模型的表面为自由边界，底部为固定边界条件，四周为法向位移约束边界条件。模型共有10912个节点，9360个实体单元。对于单侧壁导坑法，由于整个隧道模型左右不对称，可取隧道模型的一半计算，隧道的计算模型II如图2-4所示。 x、y、z各方向的长度分别为120m、60m和140m。模型的表面为自由边界，底部为固定边界条件，四周为法向位移约束边界条件。模型共有25637个节点，23220个实体单元。图2-4计算模型II示意图 FLAC3D提供了丰富的单元库供用户选择，这样可以对各种材料进行模拟。计算中V级围岩均采用摩尔—库伦材料，初期支护体系中的喷层混凝土视为线弹性体，二次衬砌混凝土及仰拱均采用

数值模拟步骤演示教学

数值模拟步骤

数值模拟 1、CFD方法简介利用CFD方法，采用流体力学分析软件Fluent对三相分离器的流场进行了研究和分析，为实验研究提供理论支持。 CFD是英文Computational Fluid Dynamics(计算流体动力学) 的缩写,是一门用数值计算方法求解流动主控方程以发现各种流动现象规律的学科]。用CFD 技术进行数值求解的基本思想是: 把原来在空间与时间坐标中连续的物理量的场, 用一系列有限个离散点上的值的集合来代替, 通过一定的原则来建立离散点上变量值之间关系的代数方程, 求解代数方程以获得所求解变量的近似值。其主要用途是对流态进行数值仿真模拟计算，因此，CFD技术的用途十分广泛，可用于传质、传热、动量传递及燃烧等方面的研究。流体机械的研究中多用CFD方法对分离器进行仿真模拟，其基本应用步骤如下： 1) 利用Gimbit进行前处理 a. 根据分离的形状、结构及尺寸建立几何模型； b. 对所建立的几何模型进行网格划分； 2) 利用Fluent进行求解 a. 确定计算模型及材料属性； b. 对研究模型设置边界条件； c. 对前期设置进行初始化，选择监视器，进行迭代计算； 3)利用Fluent进行后续处理，实现计算结果可视化及动画处理。

上述迭代求解后的结果是离散后的各网格节点上的数值,这样的结果不直观。因此需要将求解结果的速度场、温度场或浓度场等用计算机表示出来，这也是CFD 技术应用的必要组成部分。利用CFD方法进行仿真模拟可以对分离器的结构设计及参数选择作出指导，保证设计的准确度，也可以为分离器样机的试验提供理论参考。由于CFD仿真模拟的广泛使用及其重要性，国内外很多学者，如Mark D Turrell、M.Narasimha、师奇威等都对其进行了研究，尤其是A.F. Nowakowski及Daniel J.SUASNABAR等人]对CFD 技术在旋流器模拟方面的应用做了详细的介绍，这些工作对CFD技术的发展起到了积极的促进作用。 2、控制方程流体流动要受物理定律的支配，基本的守恒定律包括：质量守恒定律，动量守恒定律、能量守恒定律。如果流动处于湍流状态，系统还应要遵守附加的湍流输运方程。 1、基本假设（1）鉴于我国各主力油田采出液含水已达到80％以上，故以水代替采出液进行分析计算；天然气的主要成分是是甲烷，故采用甲烷的替代天然气的性质；（2）在工作状态下，流动不随时间变化，流动为稳态；（3）水在管内的流动可以简化成二维流动；（4）不考虑温度的影响，服从绝热流动基本方程。 2、基本控制方程（1）连续性方程在直角坐标系下的质量守恒方程又称连续性方程：

数值模拟软件大全

数值模拟软件大全 GEO-SLOPE Offical WebSite: www. geo-slope. com SLOPE/W: 专业的边坡稳定性分析软件, 全球岩土工程界首选的稳定性分析软件 SEEP/W: 专业的地下渗流分析软件, 第一款全面处理非饱和土体渗流问题的商业化软件 SIGMA/W: 专业的岩土工程应力应变分析软件, 完全基于土(岩)体本构关系建立的专业有限元软件 QUAKE/W: 专业的地震应力应变分析软件, 线性、非线性土体的水平向与竖向耦合动态响应分析软件 TEMP/W: 专业的温度场改变分析软件, 首款最具权威、涵盖范围广泛的地热分析软件 CTRAN/W: 专业的污染物扩散过程分析软件, 超值实用、最具性价比的地下水环境土工软件 AIR/W:专业的空气流动分析软件, 首款处理地下水-空气-热相互作用的专业岩土软件 VADOSE/W: 专业的模拟环境变化、蒸发、地表水、渗流及地下水对某个区或对象的影响分析软件, 设计理论相当完善和全面的环境土工设计软件 Seep3D(三维渗流分析软件)是GeoStudio2007专门针对工程结构中的真实三维渗流问题, 而开发的一个专业软件, Seep3D软件将强大的交互式三维设计引入饱和、非饱和地下水的建模中, 使用户可以迅速分析各种各样的地下水渗流问题. 特点:GeoStudio其实就是从鼎鼎大名的GEO-SLOPE发展起来的, 以边坡分析出名, 扩展到整个岩土工程范围, 基于. NET平台开发的新一代岩土工程仿真分析软件, 尤其是VADOSE/W模块是极具前瞻性的, 环境岩土工程分析的利器. 遗憾的是其模块几乎都只提供平面分析功能. Rocscience Offical WebSite: www. rocscience. com Rocscience 软件的二维和三维分析主要应用在岩土工程和采矿领域, 该软件使岩土工程师可以对岩质和土质的地表和地下结构进行快速、准确地分析, 提高了工程的安全性并减少设计成本. Rocscience 软件对于岩土工程分析和设计都很方便, 可以帮助工程师们得到快速、正确的解答. Rocscience 软件对于用户最新的项目都有高效的解算结果, 软件操作界面是基于WINDOWS 系统的交互式界面. Rocscience 软件自带了基于CAD 的绘图操作界面, 可以随意输入多种格式的数据进行建模, 用户可以快速定义模型的材料属性、边界条件等, 进行计算得到自己期望的结果. Rocscience 软件包括以下十三种专业分析模块: Slide 二维边坡稳定分析模块

基于Fluent的三通管数值模拟及分析

第40卷第2期当代化工 Vol.40，No. 2 2011年2月 Contemporary Chemical Industry February，2011 收稿日期： 2010-08-17 作者简介：魏显达（1983-），男，硕士，黑龙江北安人，2007年毕业于大庆石油学院电子信息工程，研究方向：塔顶流出系统的腐蚀与防基于 Fluent 的三通管数值模拟及分析魏显达，王为民，徐建普（辽宁石油化工大学石油天然气工程学院，辽宁抚顺 113001）摘要：Fluent 软件作为流体力学中通用性较强的一种商业CFD 软件应用范围很广。通过利用Fluent 计算流体动力学（CFD）的软件，对石油工业系统中常见的三通管内部流体进行了模拟分析，得到了三通管内在流体流动时的速度、压力和温度场分布图，为石油管道中的流体输送提供了理论依据。关键词：Fluent；三通管；模拟分析；分布图中图分类号： TQ 018 文献标识码： A 文章编号： 1671-0460（2011）02-0165-03 Numerical Simulation and Analysis of Fluid in Three-way Connection Pipe Based on Fluent Software WEI Xian-da ，WANG Wei-min ，XU Jian-pu (Institute of Petroleum and gas engineering , Liaoning Shihua University, Liaoning Fushun 113001，China ） Abstract : As a commercial CFD software with good universality, the Fluent software has been used extensively. In this paper, Simulation analysis on fluid in the three-way connection pipe of the oil industry was carried out by the software of fluid mechanics computation .Then distribution graphs of velocity , pressure and temperature of fluid in the three-way pipe were gained ，which can offer theoretical basis on fluid transportation in the petroleum pipeline. Key words : Fluent three-way ；Connection pipe ；Simulation analysis ；Distribution graphs Fluent 是目前国际上比较流行的商用CFD 软件包，在美国的市场占有率为60%，广泛应用于流体、热传热和各种化学反应等有关工业。软件包括前处理器（利用Gambit 进行物理建模、网格划分和划定边界层条件）、求解器（根据专业条件不同，采用不同的求解器，并规定物性、外部工作环境和进行数值迭代）和后处理器（把一些数据可视化，满足用户的特定要求）。三通管在石油工业中应用广泛，采用传统的设计开发方法，存在经济成本高，研发周期长等缺陷，耗费大量的人力、物力 [1-2] 。应用CFD 软件，能够在相对较短的设计周期内，较低的成本运行下，准确模拟流动具体过程，如速度场、压力场和温度场等的时变特性等。CFD 技术已经成为不可缺少的设计手段。本文利用Fluent 的超强数值计算和分析能力对三通管道内原油流动时的速度、压强和温度场进行了数值模拟和分析，为石油管道中的流体输送提供了可靠的理论依据。 1 数学模型的建立和分析输油管道管中，原油在三通管内的流动属于湍流，简化方程管道内的流体流动满足质量守恒、动量守恒、能量守恒、状态方程等。连续性方程（连续性方程式质量守恒定律在流体力学中的表现形式）在直角坐标系下表示为（（1）方程） [3-5] ： 0)()()(=??+??+??+??z y x t z y x νννρρρρ （1）式中：V x ，V y ，V z 是速度矢量ν在x 、y 和z 轴方向的分量，t 是时间，ρ是密度。最常用的湍流求解模型是标准k -ε湍流模型。它需要求解湍动能k ((2)方程)和耗散率ε((3)方程)，具体如下所示： Y G G x x M b k i t i k t k ?+++??+??=ρεσμρ μ)[(d d （2） K K k t C G C G C x x b K i t i εμρεσμερεεε2 231)(])[(d d ?++??+??= （3）

塑性成形过程的数值模拟汇总

实验报告塑性成型过程的数值模拟班级：机自07 姓名：欧阳罗辉学号：10011170 2012年12月

一、实验目的：通过本实验的教学，使学生基本掌握有限元技术在板料塑性成形领域的应用情况，拓宽学生的知识面，开阔视野，使学生对塑性成形过程的数值模拟技术有深刻的理解，预测板料弯曲成形的性能。二、教学基本要求：学会使用Dynaform数值模拟软件进行板料弯曲成形过程的仿真模拟，对模拟结果具有一定的分析和处理能力。三、实验内容提要：掌握前处理的关键参数设置，如零件定义、网格划分、模型检查、工具定义、坯料定义、工具定位和移动、工具动画、运行分析。了解后处理模块对模拟结果的分析，如读入d3plot 文件、动画显示变形和生成动画文件、成形极限图分析、坯料厚度变化分析等。四、软件操作过程: 1. 导入压边圈、板料、下模板、上模板图形文件点击File —Import，出现Import File 对话框，找到“ L型弯曲零件图” 选中binder.igs,点击Import，如此，依次导入四个模型文件，最后点击“确定”确认

四个模型导入后，结果如图 2. 重命名文件点击PartLEdit ，出现Edit Part 对话框，这里便要依次更改文件名首先选用红色文件名“ cOOIvOOO 1 ”，在上面的Name 对话框中输入binder ，然后点击 Modify ，以此类推输入 banker 、die 、punch 。 Edit Part 3. 对各图形文件划分有限元网格 1. Binder 零件网格划分

n 点击口图标，出现Part Turn。。。对话框，依次单击banker 2, die 3, punch 4,它们都会被取消选中，只留下binder 1被显示，点击0K确定。然后点击右下角的Current Part，弹出Current Part对话框，选择binder 1，点击OK确定。点击Preproces L Element，弹出Element对话框，选择Part Turn On/Off Select by Cursor 第四个图标（自动模式），将Max Size改为10，点击Select Surfaces点击Displayed Surf，点击OK，点击Apply，点击Yes,点击Exit，点击OK，于是第一个零件网格划分完成。 2. Banker零件网格划分 n 点击i□图标，取消Binder 1零件的显示，添加Banker 2 Select Dy Name All On |AllOfi OK Undo F Only SeledOn 零件的显示，点击OK确定。然后点击右下角的Current Part图标，将当前零件选成Banker 2, 点击OK确定。点击Preproces A Element,弹出Element对话框，选择第三个图标，弹出Control Keysto ne对话框，点击POINTS/NODES，弹出In put Coo。。。对话框，选中Poi nt,然后在绘图区沿顺时针或者逆时针方向依次选中Ba nker零件的四个顶点，如下图所示