数值模拟分析实例

数值模拟分析实例
数值模拟分析实例

华中科技大学体育馆数值模拟分析

6.1分析模型的建立

采用有限元软件ANSYS建立该网壳结构有限元分析模型。整体屋盖结构共计1481个节点,4430个单元,16种截面类型。建模时,网壳结构主体结构部分(包括主桁架、次桁架、水平支撑和檩条)采用ANSYS的LINK8杆单元建模,两侧翼的主梁、次梁和支承钢管柱均采用BEAM4梁单元,网壳结构屋面下部混凝土支承结构亦采用BEAM4梁单元。分析时,屋面板、设备管线等荷载等效为节点荷载,施加在结构节点上。

在网壳结构有限元分析中,对于杆件采用的LINK8 3-D Spar单元为三维单元,假设材料为均质等直杆,且在轴向上施加载荷,可以承受单向的拉伸或者压缩,每个节点上具有三个自由度,即沿X、Y和Z坐标轴方向。该单元具有塑性、蠕变、应力硬化和大变形等功能,能较好的模拟三维空间桁架单元。

对于两侧翼结构和下部支撑体系的柱、梁等结构采用的BEAM4单元是一个轴向拉压、扭转和弯曲单元,每个节点有三个平动自由度和三个转动自由度,具有应力刚化和大变形功能。

施工过程模拟分析时考虑时,同时考虑温度效应影响,计算时材料假定为理想弹塑性材料。

图6-1 有限元分析模型

6.2分析工况选取

按照实际施工顺序,将网壳结构屋盖施工过程划分为5个工况进行施工数值模拟,计算温度取为该阶段施工完成时的环境温度。

工况1: 7榀拱形主桁架安装完毕,但临时支撑未撤除,计算温度为温度15℃;

(a)短轴立面

(b)长轴立面

图6-2 工况1中屋盖结构平面图图6-3 工况1中屋盖结构立面图工况2: 两侧翼结构安装完毕,完成后拆除其临时支撑,计算温度为8℃;

(a)短轴立面

(b)长轴立面

图6-4 工况2中屋盖结构平面图图6-5 工况2中屋盖结构立面图工况3: 次桁架、水平支撑及楼梯安装完毕,临时支撑拆除,计算温度为29℃;

(a)短轴立面

(b)长轴立面

图6-6 工况2中屋盖结构平面图图6-7 工况2中屋盖结构立面图工况4: 檩条及设备管线安装完毕,计算温度为41℃;

(a)短轴立面

(b)长轴立面

图6-8 工况2中屋盖结构平面图图6-9 工况2中屋盖结构立面图工况5: 屋面板及保温层等安装完毕,计算温度为16℃。

6.3分析结果

图6-10所示为华中科技大学网壳结构屋盖施工过程模拟分析杆件应变云图。从图中可以看出,在网壳结构施工过程中,受力较大杆件都集中在支座附近。

(a)施工工况1网壳结构杆件应变图

(b)施工工况2网壳结构杆件应变图

(c)施工工况3网壳结构杆件应变图

(d)施工工况4网壳结构杆件应变图

图6-10 施工过程中网壳结构杆件应变图

7 武汉体育中心体育馆有限元分析

7.1空心球节点有限元分析

7.1.1分析模型

采用有限元软件ANSYS对空间节点进行受力性能进行分析。节点建模采用ANSYS提供的4节点弹塑性壳单元,每个节点有6个自由度,包括3个线位移自由度和3个转动自由度,该单元具有大挠度小应变的特点,能比较好的模拟球节点的受力特性。划分网格时采用映射划分和自由划分、整体划分和局部加密相结合的方法。先对模型分区,杆与球的主体采用映射网格划分法,相贯线处由于边界复杂采用自由划分。网格划分时通过控制单元边长,并在相贯线附近缩小控制尺寸,从而保证计算精度。这样,整个模型共分为56681个单元,28243个结点。加载时,考虑节点的实际受力情况加载,节点支座底部约束所有自由度,各杆上荷载按实际受力比例轴向加载。计算模型见图43。

图43 空心球节点计算模型

分析时考虑了材料非线性和几何非线性。材料为Q345钢,计算时材料假定为弹塑性材料,服从V on-Mises屈服准则,材料塑性按双线性等向强化考虑,弹性模量E=2.0×105N/mm2,屈服后模量取2%E,fy=370N/mm2。分析采用增量迭代的方法,将牛顿拉斐逊方法(Newton-Raphson)、线性搜索技术(Line Search)、应用预测(Predictor)、自适应下降(Adaptive Descent)等加速收敛技术有机结合建立非线性平衡求解方法。收敛准则为位移及不平衡节点力收敛准则。本文不计残余应力和节点区焊缝对节点极限承载力的影响。

7.1.2分析结果

在整个加载过程中,除7、8、9三根杆外,杆上应力呈弹性状态,从受荷点至球杆交界处,应力呈比例增长,由于是轴向加载,无偏心影响,杆周应力均匀分布,基本是轴向应力,环向应力很小,其中应力最大点位于最大受力杆与球的交界处。从计算可知,在设计荷载和检验荷载作用下,除7、8、9三根杆外,其余各杆均未屈服。而在1.4倍设计荷载下,7、8、9三根杆相贯处和杆底端部分区域屈服。

球面受力状态比较复杂,属于空间受力范畴。从计算得知,球杆交界处(焊趾附近)应力集中现象比较严重。在1.5倍设计荷载作用下,8号、9号杆件与空心球体交界处球体出现了局部屈服。

因此,根据上述有限元分析结果,该节点极限承载力可认为是1.5倍设计荷载。

在设计荷载和检验荷载作用下,对应于各杆轴向、支座肋板及索耳板的单向应变计测点应力分析结果见表8所示,对应于杆与球面交汇区及杆相贯处的应变花测点应力结果如表9所示。节点的等效应力云图如图44~ 45所示。

图44 设计荷载作用下节点的单元等效应力图

7.2铸钢节点有限元分析

7.2.1分析模型

采用有限元软件ANSYS对铸钢空间节点进行受力性能进行分析。节点建模采用ANSYS提供的SOLID45单元,SOLID45单元为3-D固体结构单元,由八个节点组成。在单元每一个节点上有三个自由度,即分别沿着三个坐标轴方向。此单元可以进行塑性、蠕变、应力硬化、大变形以及大应变分析。能比较好的模拟实体铸钢节点的受力特性。划分网格时采用映射划分和自由划分、整体划分和局部加密相结合的方法。这样,整个模型共分为61530个单元,13097个结点。分析时考虑了材料非线性和几何非线性,按双线性等向强化考虑,采用V on-Mises 屈服准则。材料参数取弹性模量E=2.10×105MPa,屈服强度fy=375MPa,屈服后弹性模量取2%E,收敛准则为位移及不平衡节点力收敛准则。加载时,考虑节点的实际受力情况,环向索三个方向的线位移均被约束,节点可以做竖向平面内的微小转动。A1、A2、A3方向上荷载取试验中的荷载值轴向加载,计算模型见图65。

图65 节点计算模型

7.2.2分析结果

在整个加载过程中,索夹板及索耳板的应力均呈弹性状态,由应力分析图可见应力较大区域出现在索耳板上及叉耳孔周围。从计算可知,在设计荷载和检验荷载作用下,节点未屈服。

. 在设计荷载和检验荷载作用下单元的等效应力如图65和图66所示。

图65 设计荷载作用下节点的单元等效应力图

图66 检验设计荷载作用下节点的单元等效应力图云阳县G25公路高架桥

附录一有限元分析结果

工况1下桥梁挠度分析结果

工况1下桥梁应变分析结果

工况2下桥梁挠度分析结果

工况2下桥梁应变分析结果

数值分析之幂法及反幂法C语言程序实例

数值分析之幂法及反幂法C 语言程序实例 1、算法设计方案: ①求1λ、501λ和s λ的值: s λ:s λ表示矩阵的按模最小特征值,为求得s λ直接对待求矩阵A 应用反幂法即可。 1λ、501λ:已知矩阵A 的特征值满足关系 1n λλ<< ,要求1λ、及501λ时,可 按如下方法求解: a . 对矩阵A 用幂法,求得按模最大的特征值1m λ。 b . 按平移量1m λ对矩阵A 进行原点平移得矩阵1m B A I λ=+,对矩阵B 用反幂法 求得B 的按模最小特征值2m λ。 c . 321m m m λλλ=- 则:113min(,)m m λλλ=,13max(,)n m m λλλ=即为所求。 ②求和A 的与数5011 140 k k λλμλ-=+最接近的特征值 ik λ(k=0,1,…39): 求矩阵A 的特征值中与k μ最接近的特征值的大小,采用原点平移的方法: 先求矩阵 B=A-k μI 对应的按模最小特征值k β,则k β+k μ即为矩阵A 与k μ最接近的特征值。 重复以上过程39次即可求得ik λ(k=0,1,…39)的值。 ③求A 的(谱范数)条件数2cond()A 和行列式det A : 在(1)中用反幂法求矩阵A 的按模最小特征值时,要用到Doolittle 分解方法,在Doolittle 分解完成后得到的两个矩阵分别为L 和U ,则A 的行列式可由U 阵求出,即:det(A)=det(U)。 求得det(A)不为0,因此A 为非奇异的实对称矩阵,则: max 2()s cond A λλ= ,max λ和s λ分别为模最大特征值与模最小特征值。

岩土工程数值计算方法作业

中国矿业大学 2 级硕士研究生课程考试试卷 考试科目岩土工程数值计算法 考试时间 学生姓名 学号 所在院系 任课教师 中国矿业大学研究生院培养管理处印制

《岩土工程数值计算法》课程报告课程报告分析的论文是安徽理工大学岩土工程专业乔成的硕士学位论文《深部巷道锚网喷支护结构的数值模拟与优化设计研究》。目前,数值分析方法有很多种,如有限差分法、有限单元法、边界单元法、离散单元法等。有着理论推演和试验分析无法比拟的优越性,更加贴近实际工程运用。但其求解问题的方法也是不同求解方法的近似解,要么是对基本方程和相应定解条件的直接近似求解;要么是求解原问题的等效积分方程的近似解;或者将连续的无限自由度问题变成离散的有限自由度问题再求近似解等等。在实际运用的的时候存在很多局限和不合理性。本报告基于硕士学位论文《深部巷道锚网喷支护结构的数值模拟与优化设计研究》的主要内容及该论文中的数值分析方法。对论文里数值计算与行文中存在的问题进行了分析,概括了文中的创新点,对数值分析的运用做出了总体评价,并提出了自己的一些建议。 0论文主要内容简述 文中探讨了深部巷道开挖过程中及开挖之后围岩的变形与力学特征,(岩体变形具有较强的时间效应,表现为流变或蠕变明显;扩容现象突出;大偏应力下岩体内部节理、裂隙、裂纹张开,出现新裂纹;变形非连续性明显,突然剧烈增加,且具有软岩的力学特性。)讨论了影响巷道变形的主要因素,认为地应力水平和围岩性质是影响巷道稳定的主要因素,并通过对工程实测数据与数值模拟分析对比,讨论了巷道开挖后两种关键因素作用下围岩应力场和位移场的分布情况与变化规律。在此基础上,通过围岩分类法,建立了基于定量指标JV的Hoek-Brown强度参数a和s的线性修正本构关系,并将该强度准则应用于数值模拟之中。在数值模拟分析中,利用FLAC3D对钱营孜煤矿风井巷道的进行了锚喷支护模拟分析,并结合实测数据,提出了风井巷道的锚喷支护参数提出了优化方案。 1文中所用有限差分法软件FLAC3D简介 FLAC3D是美国ITASCA公司在FLAC(Fast Lagrangian Analysis of Continua)基础上开发的三维数值分析软件,并在岩土工程数值计算中得到了广泛应用。其可实现对岩石、土和支护结构等建立高级三维模型,进行复杂的岩土数值分析与设计。 程序采用的是快速拉格朗日方法,基于显式差分来获得模型的全部运动方程(包括内变量)的时间步长解。程序将计算模型划分为若干个不同形状的三维单

数值模拟报告(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就是专业的冲压模拟软件。

演讲稿数值分析应用实例.doc

非线性方程求根 问题:在相距100m的两座建筑物(高度相等的点)之间悬挂一根电缆,仅允许电缆在中间最多下垂1m,试计算所需电缆的长度。 设空中电缆的曲线(悬链线)方程为 ] , [ , ) ( 50 50 2 - ∈ + = - x e e a y a x a x (1) 由题设知曲线的最低点)) ( , (0 0y与最高点)) ( , (50 50y之间的高度差为1m,所以有 1 2 50 50 + = +- a e e a a a) ( (2) 由上述方程解出a后,电缆长度可用下式计算: ) ( ) (a a a x a x L e e a dx e e dx x y ds L 50 50 50 50 50 2 1- - - - = ? ? ? ? ? ? + = ' + = =? ? ?(3) 相关Matlab命令: 1、描绘函数] , [ , ) ( ) (1500 500 1 2 50 50 ∈ - - + = - a a e e a a y a a 的图形;

2、用fzero 命令求方程在1250=a 附近的根的近似值x ,并计算)(x y 的函数值; 3、编写二分法程序,用二分法求0=)(a y 在],[13001200内的根,误差不超过310-,并给出对分次数; 4、编写Newton 迭代法程序,并求0=)(a y 在],[13001200内的根,误差不超过310-,并给出迭代次数。 5、编写Newton 割线法程序,并求0=)(a y 在],[13001200内的根,误差不超过310-,并给出迭代次数。

线性方程组求解应用实例 问题:投入产出分析 国民经济各个部门之间存在相互依存的关系,每个部门在运转中将其他部门的产品或半成品(称为投入)经过加工变为自己的产品(称为产出),如何根据各部门间的投入产出关系,确定各部门的产出水平,以满足社会需求,是投入产出分析中研究的课题。考虑下面的例子: 设国民经济由农业、制造业和服务业三个部门构成,已知某年它们之间的投入产出关系、外部需求、初始投入等如表1所示(数字表示产值)。 表1 国民经济三个部门间的关系单位:亿元 假定总投入等于总产出,并且每个部门的产出与它的投入成正比,由上表可以确定三个部门的投入产出表:如表2所示。 表2 三个部门的投入产出表

数值模拟分析实例

华中科技大学体育馆数值模拟分析 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)短轴立面

数值分析在生活中的应用举例及Matlab实现

Matlab 实验报告 学院:数学与信息科学学院班级:信息班 学号:20135034027 姓名:马永杉

最小二乘法,用MATLAB实现 1.数值实例 下面给定的是郑州最近1个月早晨7:00左右的天气预报所得到的温度,按照数据找出任意次曲线拟合方程和它的图像。下面用MATLAB编程对上述数据进行最小二乘拟合。 2、程序代码 x=[1:1:30]; y=[9,10,11,12,13,14,13,12,11,9,10,11,12,13,14,12,11,10,9,8,7,8,9,11,9 ,7,6,5,3,1]; a1=polyfit(x,y,3) %三次多项式拟合% a2= polyfit(x,y,9) %九次多项式拟合% a3= polyfit(x,y,15) %十五次多项式拟合% b1=polyval(a1,x) b2=polyval(a2,x) b3=polyval(a3,x) r1= sum((y-b1).^2) %三次多项式误差平方和% r2= sum((y-b2).^2) %九次次多项式误差平方和% r3= sum((y-b3).^2) %十五次多项式误差平方和% plot(x,y,'*') %用*画出x,y图像% hold on plot(x,b1, 'r') %用红色线画出x,b1图像% hold on plot(x,b2, 'g') %用绿色线画出x,b2图像% hold on plot(x,b3, 'b:o') %用蓝色o线画出x,b3图像% 2.流程图

4.数值结果分析 不同次数多项式拟合误差平方和为: r1=67.6659 r2=20.1060 r3=3.7952 r1、r2、r3分别表示三次、九次、十五次多项式误差平方和。 5、拟合曲线如下图

数值分析课程设计学生题目

《数值分析》课程设计

本课程设计的内容为:每个小组的同学均应完成以下五个案例; 目标:能将数值分析课程中所学的算法知识熟练应用于实际问题中。 案例1 土木工程和环境工程师在设计一条排水渠道时必须考虑渠道的各种参数(如宽度,深度,渠道内壁光滑度)及水流速度、流量、水深等物理量之间的关系。 假设修一条横断面为矩形的水渠,其宽度为B ,假定水流是定常的,也就是说水流速度不随时间而变化。 根据质量守恒定律可以得到 Q=UBH (1.1) 其中Q 是水的流量(s m /3 ),U 是流速(s m /),H 是水的深度(m )。 在水工学中应用的有关流速的公式是 3 /23 /22/1)2()(1H B BH S n U += (1.2) 这里n 是Manning 粗糙系数,它是一个与水渠内壁材料的光滑性有关的无量纲量;S 是水渠 的斜度系数,也是一个无量纲量,它代表水渠底每米内的落差。 把(1.2)代入(1.1)就得到 3 /23 /52/1)2()(1H B BH S n U += (1.3) 为了不同的工业目的(比如说要把污染物稀释到一定的浓度以下,或者为某工厂输入一定量 的水),需要指定流量Q 和B ,求出水的深度。这样,就需要求解 0) 2()(1)(3 /23 /52/1=-+=Q H B BH S n H f (1.4) 一个具体的案例是 s m Q S n m B /5 ,0002.0 ,03.0 ,203==== 求出渠道中水的深度H 。 所涉及的知识——非线性方程解法。 案例2 在化学工程中常常研究在一个封闭系统中同时进行的两种可逆反应 C D A C B A ?+?+2 其中A ,B ,C 和D 代表不同的物质。反应达到平衡是有如下的平衡关系: d a c b a c C C C k C C C k == 22 1 , 其中2 24 1107.3 ,104--?=?=k k 称为平衡常数,),,,(d c b a n C n =代表平衡状态时该物质的浓度。假定反应开始时各种物质的浓度为:

梳齿受力数值模拟与分析

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

岩土工程计算原理和方法

岩土工程数值计算原理与方法 随着计算机的计算速度和存储能力的飞速发展以及计算方法的日益完善,数值模拟方法已经成为研究未知领域的强有力的工具。在岩土工程计算与分析中数值计算原理与方法也发展很快。特别是有限元的发展,促进了岩土工程研究、工程预测、优化设计和计算机辅助设计等的发展。但在工程实际中使用数值计算原理与方法却存在一些问题:例如有些人因缺乏对有限元和工程性质的深入了解,而有限元的迅速发展给他们造成一种假象,认为它是万能的,可以处理几乎所有的岩土工程问题;同时他们又被有限元计算结果的精度所迷惑,不了解这些精确结果后面所隐藏的不确定性,也不了解这些数值方法所采用本构模型的局限性以及相应参数的不确定性;因这些不确定性导致数值计算原理与方法的预测结果与实际情况和实际经验相差很大,又由于部分人计算偏于保守,使得岩土工程师难以接受现代数值计算原理与方法。 1. 岩土工程数值计算原理与方法也具有两面性。 有些人偏向于用其进行岩土工程的分析计算的原因在于: (1)数值计算原理与方法能够做任何传统的分析方法所能做到的分析与计算,而且做得更多、更好。 (2)数值计算原理与方法能够给出复杂数学模型的解。因而能够从机理上预测工程性质,而不是统计和经验性的描述,这是一大优点;而简化或经验分析方法有时只能描述其表面或形式上(统计)的关系,缺乏物理机制的描述和探讨。 (3)该方法既能处理简单问题,也能处理复杂问题。 数值计算原理与方法难以被其他人接受的原因在于: (1)使用复杂,难以被很好的掌握。 (2)数值计算原理与方法本身的不确定性(指与精确的解析方法相比所产生的不确定性,特别是在岩土动力非线性问题中这种不确定性会很大)导致预测结果与工程实际不符。 (3)数值计算原理与方法所使用的物理模型或本构模型有局限性,难以反映实际情况,导致预测结果与工程实际不符。 (4)采用复杂模型要求较多的参数,而这些参数难以用简单试验获得。 (5)既然数值计算原理与方法和传统的分析方法都具有很大的不确定性,还不如采用传统的分析方法,因为传统的方法简单、实用。 (6)精确的数值分析结果会误导使用者迷信这些结果的精确性,而没有认识到其后面隐

顶盖驱动流数值模拟分析

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

西安科技大学能源学院安全技术及工程 申敬杰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):灰褐色,湿,可塑状,主要由亚粘土组成,局部可见有植物根系,偶

数值模拟软件大全

数值模拟软件大全 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 二维边坡稳定分析模块

FLAC-数值模拟分析

1工程概况 研究此段为中条山隧道K9+450~K10+560段,此处隧道最大埋深约540m,主要由太古界涑水群表壳岩组合解州片麻岩(Hgn)地层组成,构成中条山隧道分水岭北侧的主体;此段组成隧道的围岩岩性主要为变粒岩、花岗片麻岩等。该套地层岩性复杂,组合无规律。岩层产状整体倾向南东,倾角一般在50o~70o间变化。在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,斜插角10o~15o,每环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)连续性方程 在直角坐标系下的质量守恒方程又称连续性方程:

相关文档
最新文档