deflt-3d
- 格式:doc
- 大小:1.65 MB
- 文档页数:32
DEFORM-3D基本操作入门QianRF前言有限元法是根据变分原理求解数学物理问题的一种数值计算方法。
由于采用类型广泛的边界条件,对工件的几何形状几乎没有什么限制和求解精度高而得到广泛的应用。
有限元法在40年代提出,通过不断完善,从起源于结构理论、发展到连续体力学场问题,从静力分析到动力问题、稳定问题和波动问题。
随着计算机技术的发展与应用,为解决工程技术问题,提供了极大的方便。
现有的计算方法(解析法、滑移线法、上限法、变形功法等)由于材料的本构关系,工具及工件的形状和摩擦条件等复杂性,难以获得精确的解析解。
所以一般采用假设、简化、近似、平面化等处理,结果与实际情况差距较大,因此应用不普及。
有限元数值模拟的目的与意义是为计算变形力、验算工模具强度和制订合理的工艺方案提供依据。
通过数值模拟可以获得金属变形的规律,速度场、应力和应变场的分布规律,以及载荷-行程曲线。
通过对模拟结果的可视化分析,可以在现有的模具设计上预测金属的流动规律,包括缺陷的产生(如角部充不满、折叠、回流和断裂等)。
利用得到的力边界条件对模具进行结构分析,从而改进模具设计,提高模具设计的合理性和模具的使用寿命,减少模具重新试制的次数。
通过模具虚拟设计,充分检验模具设计的合理性,减少新产品模具的开发研制时间,对用户需求做出快速响应,提高市场竞争能力。
一、刚(粘)塑性有限元法基本原理刚(粘)塑性有限元法忽略了金属变形中的弹性效应,依据材料发生塑性变形时应满足的塑性力学基本方程,以速度场为基本量,形成有限元列式。
这种方法虽然无法考虑弹性变形问题和残余应力问题,但可使计算程序大大简化。
在弹性变形较小甚至可以忽略时,采用这种方法可达到较高的计算效率。
刚塑性有限元法的理论基础是Markov变分原理。
根据对体积不变条件处理方法上的不同(如拉格朗日乘子法、罚函数法和体积可压缩法),又可得出不同的有限元列式其中罚函数法应用比较广泛。
根据Markov变分原理,采用罚函数法处理,并用八节点六面体单元离散化,则在满足边界条件、协调方程和体积不变条件的许可速度场中对应于真实速度场的总泛函为:∏≈∑π(m)=∏(1,2,…,m)(1)对上式中的泛函求变分,得:∑=0(2)采用摄动法将式(2)进行线性化:=+Δun(3)将式(3)代入式(2),并考虑外力、摩擦力在局部坐标系中对总体刚度矩阵和载荷列阵,通过迭代的方法,可以求解变形材料的速度场。
Deform-3D(version6.1)使用步骤Deform—3D是对金属体积成形进行模拟分析的优秀软件,最近几年的工业实践证明了其在数值模拟方面的准确性,为实际生产提供了有效的指导。
Deform—3D的高度模块化、友好的操作界面、强大的处理引擎使得它在同类模拟软件中处于领先地位。
以下将分为模拟准备、前处理、求解器、后处理四部分简要介绍Deform—3D的使用步骤。
一、模拟准备模拟准备阶段主要是为模拟时所用的上模、下模、坯料进行实体造型,装配,并生成数据文件。
实体造型可通过UG、Pro-e、Catia、Solidworks等三维作图软件进行设计,并按照成形要求进行装配,最后将装配体保存为STL格式的文件。
该阶段需要注意的是STL格式的文件名不能含有中文字符;另外对于对称坯料,为了节省求解过程的计算时间并在一定程度上提高模拟精度(增加了网格数量),可把装配体剖分为1/4,1/8或更多后再进行保存。
二、前处理前处理是整个数值模拟的重要阶段,整个模拟过程的工艺参数都需要在该阶段设置,各参数设置必须经过合理设置后才能保证模拟过程的高效性和模拟结果的准确性。
首先打开软件,新建(new problem)→选择前处理(Deform-3D preprocessor)→在存放位置(Problem location)选项卡下选择其他(other location)并浏览到想要存放deform 模拟文件的文件夹→下步的problem name可任意填写。
注意:所有路径不能含有中文字符。
simulation controls)→改变单位(units)为SI,接受弹出窗口默认值;选中模式(mode)选项卡下热传导(heat transfer)。
导入坯料、模具并设置参数:导入毛坯:1、general:通常采用刚塑性模型即毛坯定义为塑性(plastic),之后导入的模具定义为刚性(rigid);温度(temperature):根据成形要求设定坯料预热温度(温热成形时一定注意);材料(material):点击load选择毛坯材料,若材料库中没有对应的材料可选择牌号相近的。
Deform-3d热处理模拟操作热处理工艺在机械制造中占有十分重要的地位。
随着机械制造现代化和热处理质量管理现代化的发展,对热处理工艺提出了更高的要求。
热处理工艺过程由于受到加热方式、冷却方式、加热温度、冷却温度、加热时间、冷却时间等影响,金属内部的组织也会发生不同的变化,因此是个十分复杂的过程,同时工艺参数的差异,也会造成热处理加工对象硬度过高过低、硬度不均匀等现象。
Deform-3d 软件提供一种热处理模拟模块,可以帮助热处理工艺员,通过有限元数值模拟来获得正确的热处理参数,从而来指导热处理生产实际。
减少批量报废的质量事故发生。
热处理模拟,涉及到热应力变形、热扩散和相变等方面,因此计算很复杂,软件采用牛顿迭代法,即牛顿-拉夫逊法进行求解。
它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。
多数方程不存在求根公式,因此求精确根非常困难,甚至不可能,从而寻找方程的近似根就显得特别重要。
方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。
牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程f(x) = 0的单根附近具有平方收敛,而且该法还可以用来求方程的重根、复根等。
但由于目前Deform-3d软件的材料库只带有45钢、15NiCr13和GCr15等三种材料模型,而且受到相变模型的局限,因此只能做淬火和渗碳淬火分析,更多分析需要进行二次开发。
本例以45钢热处理淬火工艺的模拟过程为例,通过应用Deform-3d 热处理模块,让读者基本了解热处理工艺过程有限元模拟的基本方法与步骤。
1 、问题设置点击“文档”(File)或“新问题”(New problem),创建新问题。
在弹出的图框中,选择“热处理导向”(heat treatment wizard),见图1。
图1 设置新问题2、初始化设置完成问题设置后,进入前处理设置界面。
首先修改公英制,将默认的英制(English)修改成公制(SI),同时选中“形变”(Deformation)、“扩散”(Diffusion)和“相变”(Phase transformation),见图2。
DEFORM-3D基本操作入门QianRF前言有限元法是根据变分原理求解数学物理问题的一种数值计算方法。
由于采用类型广泛的边界条件,对工件的几何形状几乎没有什么限制和求解精度高而得到广泛的应用。
有限元法在40年代提出,通过不断完善,从起源于结构理论、发展到连续体力学场问题,从静力分析到动力问题、稳定问题和波动问题。
随着计算机技术的发展与应用,为解决工程技术问题,提供了极大的方便。
现有的计算方法(解析法、滑移线法、上限法、变形功法等)由于材料的本构关系,工具及工件的形状和摩擦条件等复杂性,难以获得精确的解析解。
所以一般采用假设、简化、近似、平面化等处理,结果与实际情况差距较大,因此应用不普及。
有限元数值模拟的目的与意义是为计算变形力、验算工模具强度和制订合理的工艺方案提供依据。
通过数值模拟可以获得金属变形的规律,速度场、应力和应变场的分布规律,以及载荷-行程曲线。
通过对模拟结果的可视化分析,可以在现有的模具设计上预测金属的流动规律,包括缺陷的产生(如角部充不满、折叠、回流和断裂等)。
利用得到的力边界条件对模具进行结构分析,从而改进模具设计,提高模具设计的合理性和模具的使用寿命,减少模具重新试制的次数。
通过模具虚拟设计,充分检验模具设计的合理性,减少新产品模具的开发研制时间,对用户需求做出快速响应,提高市场竞争能力。
一、刚(粘)塑性有限元法基本原理刚(粘)塑性有限元法忽略了金属变形中的弹性效应,依据材料发生塑性变形时应满足的塑性力学基本方程,以速度场为基本量,形成有限元列式。
这种方法虽然无法考虑弹性变形问题和残余应力问题,但可使计算程序大大简化。
在弹性变形较小甚至可以忽略时,采用这种方法可达到较高的计算效率。
刚塑性有限元法的理论基础是Markov变分原理。
根据对体积不变条件处理方法上的不同(如拉格朗日乘子法、罚函数法和体积可压缩法),又可得出不同的有限元列式其中罚函数法应用比较广泛。
根据Markov变分原理,采用罚函数法处理,并用八节点六面体单元离散化,则在满足边界条件、协调方程和体积不变条件的许可速度场中对应于真实速度场的总泛函为:∏≈∑π(m)=∏(1,2,…,m)(1)对上式中的泛函求变分,得:∑=0(2)采用摄动法将式(2)进行线性化:=+Δun(3)将式(3)代入式(2),并考虑外力、摩擦力在局部坐标系中对总体刚度矩阵和载荷列阵,通过迭代的方法,可以求解变形材料的速度场。
Deform-3d热处理模拟操作热处理工艺在机械制造中占有十分重要的地位。
随着机械制造现代化和热处理质量管理现代化的发展,对热处理工艺提出了更高的要求。
热处理工艺过程由于受到加热方式、冷却方式、加热温度、冷却温度、加热时间、冷却时间等影响,金属内部的组织也会发生不同的变化,因此是个十分复杂的过程,同时工艺参数的差异,也会造成热处理加工对象硬度过高过低、硬度不均匀等现象.Deform—3d软件提供一种热处理模拟模块,可以帮助热处理工艺员,通过有限元数值模拟来获得正确的热处理参数,从而来指导热处理生产实际.减少批量报废的质量事故发生.热处理模拟,涉及到热应力变形、热扩散和相变等方面,因此计算很复杂,软件采用牛顿迭代法,即牛顿-拉夫逊法进行求解.它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法.多数方程不存在求根公式,因此求精确根非常困难,甚至不可能,从而寻找方程的近似根就显得特别重要。
方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根.牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程f(x) = 0的单根附近具有平方收敛,而且该法还可以用来求方程的重根、复根等.但由于目前Deform-3d软件的材料库只带有45钢、15NiCr13和GCr15等三种材料模型,而且受到相变模型的局限,因此只能做淬火和渗碳淬火分析,更多分析需要进行二次开发.本例以45钢热处理淬火工艺的模拟过程为例,通过应用Deform-3d 热处理模块,让读者基本了解热处理工艺过程有限元模拟的基本方法与步骤.1 、问题设置点击“文档”(File)或“新问题”(New problem),创建新问题。
在弹出的图框中,选择“热处理导向”(heat treatment wizard),见图1。
图1 设置新问题2、初始化设置完成问题设置后,进入前处理设置界面。
首先修改公英制,将默认的英制(English)修改成公制(SI),同时选中“形变”(Deformation)、“扩散”(Diffusion)和“相变”(Phase transformation),见图2。
Deform轧制模拟流程1. 建立工作目录2. 前处理(导入模型,调整模型位置,设置物理性质,划分网格,设置边界条件(换热面、对称面)3. 设置辊和推板的速度,设置接触条件,check并生成DB文件)4. 运行计算(多线程计算)5. 后处理(云图、点追踪)建立工作目录在硬盘建立专属文件夹(必须是英文文件夹),在deform软件中打开进入前处理页面添加objects ,分别导入模型(注:下图红色方框所示为主程序,前处理的五个主要功能选项, 即:主菜单、材料库、装配、接触、检查与生成DB 文件。
)i?2D/J0 Ver 10.2:* 云 r«T' -0 aafn^XuReoivi □Ipbwd - S&%D0ORM 计窮hDKD ・L0工3d 亠皐: □ ri- r.in7^-T :0<3 TT : n 论亠古调整模型位置(Position object选择top die,在method中选择rotational (旋转),选中x 轴,在角度中输入90,点击apply,将平方的轧辊旋转90度。
如果有必要)在method中选择ofset (平移),选择work piece (轧件)根据道次压下量调整轧件与轧辊之间的相对位置, 推板bottom die同理。
划分网格 选中work piece ,一般默认的是相对网格保证计算速度,而绝对网格的计算精度更 o3o1.2广亠亠I 订工第二『0:企嵌I ;设置边界条件:1.对称面,1/4模型有两个对称面;2•换热面剩余面为换热面高但速度较慢,无论哪种网格划分方式都具有在计算过程中自动重新划分网格的功能 显示单个对象N JE 力打 cf Swi-?:- | : T■ K :«3ScLrrQGw» ■ M r□ mMo** ■■*=-'; if DCOF'H I-JJ t-r'x l :2 :J EW lj2d :F*in<-it*Ml iAKiWi? » * :」佩-口.7二二_£*1I 「■■ ■■■l ・ | 二■■>1 f -dJ!旧tar :coal ,:yiKt肌 toraatiwiFiwiJiiiilg MPNH .加 Pi •也订driE_U2ifl— □ 壬Pl fsl* SWT HlEMrr u 好M 审帆 J«ih 初MAS!邯油-l*J XJ>JB A Cd||l ■ 00fi ■[ h * 却 UGO 舍仓舍 «■&:<4.■上2陆巴J.…r..-“」Q f^v* W 気- 1JC Vtih j7勺:i 咎si 4母<KAQ6®St卫■ -h / * 4规fsl* Sw'B Vi Wflerr 曲"・《7 Mdn. H M >〕・ A L ” CqiTiMt *Hi■1 --------------------------------------------- G^JzETEfla! 5TKILA-ZDS「皿•常:养帝亡加:3:s>-J@ 0臥範J 超■仙皿就FE- 24M RET •讦Uz ・ xa1□:Ph.j Tup 31*T»B 1 cb^scc: ■'■ 3佥毬■•轟从 0»JM ;IdJlASMfa圭]33」rJ*12Oio2Sirruldtan Cartr^hOtJMK ||>V iM-libhAStJ 町S卜 Feint I 吨詢価 h Ml ; 0£t 心0 titc wr*« 15 ■ 7v*-o«i£tjiHPPi F] If EVIT 冊 EMh D*俱* Mdfil 7x>h Virr ^TlSn?肝:Q!>ja A nd 呂屮少 &也廿门 吕 戶|比却Q 追。
Deflt-3D FLOW模块简明操作手册河海大学第一章引言1.1 特征wL 一Delft HydrauliCS 开发一套功能强大的软件包:Delft3D , 主要应用于自由地表水环境。
该软件具有灵活的框架,能模拟二维(水平面或竖直面)和三维的水流、波浪、水质、生态、泥沙输移及床底地貌,及各个过程之间的相互作用。
Delft3D 的主要特征为:所有程序模块都具有高度的整合性和互操作性;能直接应用最新过程知识,为久誉盛名的水力学研究所的研发成果;采用了市场上最为友好的图形用户界面。
1.2 模块简介Delft3D 具体模拟六种现象的时间、空间变化以及其相互联系。
该软件包原则上广泛适用于多种情况。
运用最为普遍的是海岸、河流以及河口地区。
Delft3D 由一系列经过全面测试和验证的模块组成,是相互联系的有机整体。
有关模块包括:一、水动力模块(Delft3D 一FLOW )该模块主要用于浅水非恒定流模拟。
综合考虑了潮汐、风、气压、密度差(由盐度和温度引起)、波浪、紊流(从简单常量到k-模型)以及潮滩的干湿交替。
本模块集成了热量及物质传输方程求解,并在wL-Delft Hydraulics 有关分层水动力学等前沿理论研究基础上开发而成。
Delft3D 的其它模块均可采用了该模块的输出结果。
二、波浪模块(Delft3D 一WAVE )波浪模块主要计算短波在非平整床底上的非稳定传播,考虑风力、底部摩阻力造成的能量消散、波浪破碎、波浪折射(由于床底地形、水位及流场)、浅水变形及方向分布。
三、水质模块(Delft3D 一WAQ)该模块通过考虑一系列泥沙输移和水质过程来模拟远一中水域的水质及泥沙。
该模块包含了若干对流扩散方程求解工具和一个庞大的标准化过程方程库,其方程组对应用户所选择的物质类型。
四、颗粒跟踪模块(Delft3D 一以RT )颗粒跟踪模块为短期的、邻近水域水质模块,通过即时跟踪个体颗粒轨迹来估算其动态、空间(子网格尺度下)密度分布。
Deform-3d热处理模拟操作热处理工艺在机械制造中占有十分重要的地位。
随着机械制造现代化和热处理质量管理现代化的发展,对热处理工艺提出了更高的要求。
热处理工艺过程由于受到加热方式、冷却方式、加热温度、冷却温度、加热时间、冷却时间等影响,金属内部的组织也会发生不同的变化,因此是个十分复杂的过程,同时工艺参数的差异,也会造成热处理加工对象硬度过高过低、硬度不均匀等现象。
Deform-3d软件提供一种热处理模拟模块,可以帮助热处理工艺员,通过有限元数值模拟来获得正确的热处理参数,从而来指导热处理生产实际。
减少批量报废的质量事故发生。
热处理模拟,涉及到热应力变形、热扩散和相变等方面,因此计算很复杂,软件采用牛顿迭代法,即牛顿-拉夫逊法进行求解。
它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法。
多数方程不存在求根公式,因此求精确根非常困难,甚至不可能,从而寻找方程的近似根就显得特别重要。
方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。
牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程f(x) = 0的单根附近具有平方收敛,而且该法还可以用来求方程的重根、复根等。
但由于目前Deform-3d软件的材料库只带有45钢、15NiCr13和GCr15等三种材料模型,而且受到相变模型的局限,因此只能做淬火和渗碳淬火分析,更多分析需要进行二次开发。
本例以45钢热处理淬火工艺的模拟过程为例,通过应用Deform-3d 热处理模块,让读者基本了解热处理工艺过程有限元模拟的基本方法与步骤。
1 、问题设置点击“文档”(File)或“新问题”(New problem),创建新问题。
在弹出的图框中,选择“热处理导向”(heat treatment wizard),见图1。
sgniht图1 设置新问题2、初始化设置完成问题设置后,进入前处理设置界面。
首先修改公英制,将默认的英制t h i ng si nt he i rb ei n ga re go od fo rs (English )修改成公制(SI ),同时选中“形变”(Deformation )、“扩散”(Diffusion)和“相变”(Phase transformation),见图2。
Deflt-3D FLOW模块简明操作手册河海大学张乃东第一章引言1.1 特征wL 一Delft HydrauliCS 开发一套功能强大的软件包:Delft3D , 主要应用于自由地表水环境。
该软件具有灵活的框架,能模拟二维(水平面或竖直面)和三维的水流、波浪、水质、生态、泥沙输移及床底地貌,及各个过程之间的相互作用。
Delft3D 的主要特征为:所有程序模块都具有高度的整合性和互操作性;能直接应用最新过程知识,为久誉盛名的水力学研究所的研发成果;采用了市场上最为友好的图形用户界面。
1.2 模块简介Delft3D 具体模拟六种现象的时间、空间变化以及其相互联系。
该软件包原则上广泛适用于多种情况。
运用最为普遍的是海岸、河流以及河口地区。
Delft3D 由一系列经过全面测试和验证的模块组成,是相互联系的有机整体。
有关模块包括:一、水动力模块(Delft3D 一FLOW )该模块主要用于浅水非恒定流模拟。
综合考虑了潮汐、风、气压、密度差(由盐度和温度引起)、波浪、紊流(从简单常量到k 一模型)以及潮滩的干湿交替。
本模块集成了热量及物质传输方程求解,并在wL 一Delft Hydraulics 有关分层水动力学等前沿理论研究基础上开发而成。
Delft3D 的其它模块均可采用了该模块的输出结果。
二、波浪模块(Delft3D 一WAVE )波浪模块主要计算短波在非平整床底上的非稳定传播,考虑风力、底部摩阻力造成的能量消散、波浪破碎、波浪折射(由于床底地形、水位及流场)、浅水变形及方向分布。
三、水质模块(Delft3D 一WAQ )该模块通过考虑一系列泥沙输移和水质过程来模拟远一中水域的水质及泥沙。
该模块包含了若干对流扩散方程求解工具和一个庞大的标准化过程方程库,其方程组对应用户所选择的物质类型。
四、颗粒跟踪模块(Delft3D 一以RT )颗粒跟踪模块为短期的、邻近水域水质模块,通过即时跟踪个体颗粒轨迹来估算其动态、空间(子网格尺度下)密度分布。
污染物可以是难降解的,也可以遵循简单的一阶降解过程。
该模块也可用于滨岸水域疏浚/浊漏等灾害事件模拟。
五、生态模块(Delft3D 一ECO )Delft 3D 系统采用了不同的藻类生长和营养动力学模块。
例如,研究富营养化现象时,过程库里嵌入了基本控制过程模块,描述生物及非生物生态系统及其相互作用。
除Delft3D 一wAQ 模块里所有和藻类相关的水质变化过程之外,生态模块还包括一些更为细化的水质过程。
六、泥沙输移模块(Delft3D 一SED )该Delft3D 一WAQ 的子模块用来模拟粘性或非粘性、有机或无机、悬移质或推移质泥沙的输移、侵蚀和沉降过程。
该模块包括若干标准运动方程,单独考虑不同的泥沙粒径。
由于忽略床底地貌变化的影响,该模块仅适用于评估短期的泥沙输移过程。
七、动力地貌模块(Delft3D 一MoR )该模块用于计算床底地形的变化,其结果取决于泥沙输移梯度以及用户定义的、和时间有关的边界条件。
模块中包含风和波浪驱动力,以及一系列的运输方程。
该模块的突出特点,是与Delft3D FLOW 和WAVE 模块的动态回馈。
由此,水流和波浪能够根据当地水下地形自行调整,可以给出任意时间范围的预报成果。
本手册以DEFLT-3D(3.23)为版本,重点介绍水动力模块(Delft3D 一FLOW )操作平台。
由于是简明操作手册,因此对于程序的实际原理不做详细介绍。
第二章水动力模块(FLOW)简介FLOW 模块是一个多维(二维或三维)水动力学(和物质输运)模拟程序,该模块综合考虑潮流、气象作用,采用边界拟合较好的曲线网格离散格式,可计算非稳定流和物质输运现象。
在三维模拟过程中,垂向网格采用。
坐标离散。
这样,整个计算场的垂面层数保持不变,从而大大地提高了计算效率。
2.1应用领域1、河川径流模拟2、河流淡水注入海湾流量过程3、咸水入侵4、湖泊、海洋、水库中的温度分层现象5、冷却水取水口及污水排放口6、溶解及污染物运移7、潮流和风生流(如:风暴潮)8、分层流和密度流9、波生流水流模块FLOW 计算得出的水动力学结果(流速、水位、密度、盐度等)都可以用作Delft3D 其他模块的输入数据。
2.3水流模块FLOW操作平台简介应用水流模块建立一套平面水流数学模型总体上需要四个步骤:原体数据的采集和整理,数模的前处理,数值计算过程和数模的后处理。
原体数据的采集和整理是:所建数模的边界文件数据和地形散点数据的采集和整理。
边界文件数据指的是确定模型的范围,如在CAD中读取数据,一般用一条闭合的样条曲线或多段线将模型边界勾画出来,然后要适当调整线段,使所勾画的边界要尽量光滑,其线段上的点要适当密一点。
然后从CAD中读取样条曲线或多段线上点的X、Y坐标,最后建立以后缀名为.ldb的边界文件。
边界文件格式见附表1。
地形散点数据的用来确定模型的地形,如在CAD中读取数据,要尽量使所采集的地形散点反映原体地形,地形点一般需要三个坐标:X、Y和Z。
在CAD 中直接提取地形点坐标比较麻烦,可以编一些小程序提取数据。
最后建立以后缀名为.xyz的地形散点文件。
地形散点文件格式见附表2。
我们最好先建立一个目标文件夹,将上面两个文件放入,然后正式进入建模的环节。
数模的前处理有:网格生成模块(RGFGRID)和模型初始场数据生成模块(QUICKIN)。
数模的数值计算:水流模块(FLOW)数模的后处理:可视化模块(QP和GPP)第三章FLOW模块的前处理3.1 网格生成Delft3D 的Delft3D 一RGFGRID 网格生成模块,可以生成用于Delft3D 各类模块的尺寸可变的正交曲线网格。
网格尺寸可变,便于用户在重点模型区域布置较高密度的网格,而在与之远离的模型边界区域网格密度采用较低密度的网格,以此减小计算量。
此外,网格线可以沿陆地边界和渠道弯曲,能达到和边界的光滑嵌合, Delft3D 一RGFGRID 网格生成模块允许分步生成网格。
先将网格进行大致的样条划分,而后把样条转化成粗疏的网格,然后再采取平滑加密。
在整个过程当中,可以随时生成正交网格。
曲线网格可以用于笛卡尔坐标系和球面坐标系。
下面是网格生成的操作步骤。
打开DEFLT-3D,进入3D界面,界面右下方会有chang working directory 指令,点击进入将路径设定为你的目标夹(即建模的文件夹)。
然后退出,点击Grid-RGFGRID进入网格生成模块。
图3.1为RGFGRID界面。
图3.1 RGFGRID界面点击界面中的File,然后点击open land boundary指令,出现文件框,从目标文件夹中打开后缀名为.ldb的边界文件,如图3.2图3.2点击界面上的指令,界面会出现模型的边界轮廓。
点击界面上的Edit-spline创建spline曲线:一种以节点控制弯曲程度的顺滑的自由曲线,通过编辑(移动\增加\删除)节点可以很容易的调节曲线的曲率和走向,对于刻画不规则的轮廓非常方便.。
是删除键和移动键。
见图3.3。
图3.3 spline曲线例子运用spline曲线将模型分割成如上的形状,这是生成网格的前提。
Spline 曲线分割模块时,要考虑边界形状和网格的大小。
指令Option-change parameters 是用来设定网格各项参数。
点击operations-change splings into grid指令,则会按照设定的网格参数生成网格。
如图3.4。
图3.4 生成网格点击指令view-splines-nosplines,不显示splines曲线,此时如有需要,可以更改网格参数,使之变密,见图3.5。
图3.5 隐藏splines曲线点击edit-line to land boundary使得网格与边界尽量贴近。
具体操作是鼠标左键选中边界网格上一点,在右键网格上另一点,使之与边界吻合。
见图3.6。
图3.6 line to land boundary指令前后网格变化点击view-grid properties-orthogonality会检验网格正交性,见图3.7。
一般来说,正交性由网格相交线的余弦值衡量,一般来说余弦值要小于0.02,非常靠近边界的一些网格,余弦值高于这个标准,可以接受。
图3.7 网格正交性如果网格正交性不符合要求,可点击operations-orthogonalise grid来提高网格正交性。
如图3.8。
图3.8 修正后的网格正交性对于边界网格,可以点击edit-block orthogonalise来提高边界网格的正交性。
见图3.9。
图3.9 边界网格的正交化网格正交化完成后,一个简单的网格就生成了。
网格一些其他属性,如平滑性、疏密度和纵横比等具体请看其应用手册。
网格生成后,点击file-save grid保存网格数据。
3.2 模型初始场的生成Delft3D-QUICKIN模块用于生成和操作模型的水下地形,以及Delft3D-FLOW 水动力模块、Delft3D-WAVE 波浪模块和Delft3D-WAQ 水质模块的初始条件。
采集的深度数据(原始数据)的来源经常各不相同,包括日期不同,质量各异,疏密度差别等等。
为了避免低质量样本“污染”高质量样本,Delft3D-QU 工CK 工N 可以实现数据组的依次加载功能。
下面介绍QUICKIN模块的操作步骤。
进入FLOW-QUICKIN模块,出现QUICKIN模块的界面,见图3.10。
图3.10 QUICKIN模块界面点击file,从中点击open grid和open land boundary打开已做好的网格。
如图3.11。
3.11 边界和网格点击open sample打开.xyz文件。
如图3.12。
图3.12 地形散点下面进行散点的插值,将散点的数值转化为网格节点的数值。
QUICKIN模块中有两种插值方法:三角插值和平均插值。
一般来说,如果网格节点数多于散点数的话,应用三角插值,反之应用平均插值。
从3.12可以看出,有部分节点数多于散点数,其余部分,节点数少于散点数,因此,将分别进行插值。
先进行三角插值,点击edit-polygon编辑一多边形,如图3.13。
图3.13 部分三角插值点击operations-triangular interpolation,进行部分三角插值,见上图。
点击view-siample points-no siample points隐藏散点,见图3.14。
图3.14 隐藏散点点击operations-gridcell averaging。
点击view-siample points-coloured dots使散点呈现。
下面对散点数多于节点数的部分进行平均插值。
点击addsubdel-delete polygon去掉编辑的多边形,在另一部分重新编辑多边形,见图3.15。