保温热力管有限元热固耦合分析

保温热力管有限元热固耦合分析
保温热力管有限元热固耦合分析

ansys流固耦合模态分析

有问题可以发邮件给我一起讨论xw4996@https://www.360docs.net/doc/0217494949.html, FSI流固耦合命令求解流固耦合问题 使用ANSYS计算结构在水中的模态时, FLUID29,FLUID30单元分别用来模拟二维和三维流体部分,相应的结构模型则利用PLANE42单元和SOL ID45等单元来构造,其中,PLANE42和SOL ID45分别是用来构造二维和三维结构模型的单元。FLUID30是流体声单元,主要用于模拟流体介质及流固耦合问题。该单元有8 个节点,每个节点上有4 个自由度,分别是XYZ上3个方向位移自由度和1个压力自由度,为各向同性材料。输入材料属性时,需要输入流体的材料密度(作为DENS 输入)及流体声速(作为SONC输入),流体粘性产生的损耗效应忽略不计。FLUID29是FLUID30单元在二维上的简化,少了一个Z向的位移。SOLID45单元用于构造三维实体结构。单元通过8 个节点来定义,每个节点有 3 个沿着XYZ方向平移的自由度。PLANE42是SOLID45单元在二维上的简化。 在利用ANSYS建模分析时,流场域单元属性分为2种,由KEYOPT(2)(指定流体和结构分界面处结构是否存在) 控制,在流固耦合交界面上的单元KEYOPT(2) = 0 ,表示分界面处有结构,其他流体单元KEYOPT(2)=1,表示分界面处无结构。流体-结构分界面通过面载荷标志出来,指定FSI label可以把分界面处的结构运动和流体压力耦合起来,分界面标志在分界面处的流体单元标出。 数值分析的步骤 1) 建立流体单元的实体模型。建立流体模型,需要确定流体域的范围,可以把无限边界流体简化成流体区域的半径为固体结构半径的10倍。 2) 标记流固耦合界面。选取流体单元中流固交界面上的节点,执行FSI 命令,流固耦合交界面的处理:流体与固体是两个独立的实体,在划分单元时在两者交界面上的单元网格要划分一致,这样在交界面上的同一位置一般就有两个重合的节点,一个节点属于流体单元,一个节点属于固体单元,这两个重合节点在交界面的位移强制保持一致。 3) 建立固体结构实体模型。建立固体结构模型,定义单元属性,采用映射方式进行网格的划分。 4) 施加约束条件。由于流体区域的尺寸远大于固体结构尺寸,故可以不考虑流体液面的重力的影响,将流体边界处的单元节点上施加压力(PRES) 为零的约束。因为选择的算例为悬臂结构,在固体结构底部加全约束。 5) 选择求解算法,进行求解。定义分析类型为模态分析,设定提取频率阶数和提取模态的方法。因为耦合问题的刚度矩阵,质量矩阵都不对称,需要采用非对称矩阵法(UNSYMMETRIC)求解。 6) 查看结果。进入后处理模块,查看结构模型的频率及振型。 以半浸没与水中的桥墩模态问题为背景,并假设: 1. 桥墩为实心等截面的实体,实际桥墩模型应该是空心壳体,截面尺寸也 非常复杂,因而需要分块划分单元。

SAMCEF有限元分析 耦合分析典型实例

第18章耦合分析典型实例 18.1 压电材料与控制装置的强耦合模型 18.1.1 模型描述 如图18-1所示结构,外部环型结构为激励控制部件,其在外力作用下变形产生挤压效应。结构中间横梁为压电材料,其在控制部件的挤压下产生电势差。 18.1.2 分析目标 ●线性静力分析 ●耦合系统的模态分析 ●耦合系统谐响应分析 ●耦合系统瞬态响应分析 18.1.3 线性静力分析 我们将分析压电材料对指定常电压的反应,同时也会得到环形激励控制部件的变形。 1. 模型处理及准备 (1)打开Samcef Field软件 (2)在求解模块对话框Solver Driver Setting中进行如下设定: Domain=Piezeoelectric Analysis Solver=OOFELIE

Analysis type=Linear Static (3)点击工具菜单中File\Import geometry,在出现的Input对话框中,直接进入 samceffield安装路径下的example文件夹,选择brep\piezoelectrics\Actuator.brep,点 击导入模型。 (4)点击工具菜单中的File\Save As将文件另存到一新的工作路径中。 2. 建立分析数据 点击按钮进入分析数据定义模块。 (1)几何属性和材料属性 我们先来定义两个非压电部分的属性,对于本例来说就是外部的环和两个压电体中间的连接部分。 1)在数据树或3D显示区域选择”Actuator”,点击工具栏中的几何属性定义 按钮,在几何属性定义对话框中,几何类型(behavior)选择体,对象选择类型(place on) 选择。几何属性(type)选择柔体,在3D显示区域选择以下两个高 亮显示部分作为Flexible part(如图18-2)。 图18-2 非压电部分的选择 2)点击工具栏中的材料特性定义按钮,Place on仍旧选择上面提到的两个实体,材料 特性定义对话框中进行如下定义。 - Type : Isotropic - Young modulus : 2.1e11 Pa - Poisson ratio : 0.3 - Mass density : 7800 kg/m3 注:在OOFELIE定义材料特性时,和结构分析中定义材料有一定区别,对于一般线弹性本构模型,只需给出弹性模量,柏松比以及材料的密度。 接下来我们来定义压电材料的属性。 3)仍旧点击工具栏中的几何属性定义按钮,几何类型(Behavior)选择 导体,对象选择类型(place on)选择。导体类型(Type of orientation)选择(横向同性)。Transverse isotropic direction选

有限元热力学常见概念汇总

Film Coefficient(对流换热系数) 流体与固体表面之间的换热能力,比如说,物体表面与附近空气温差1℃,单位时间单位面积上通过对流与附近空气交换的热量。单位为W/(m^2·℃)。表面对流换热系数的数值与换热过程中流体的物理性质、换热表面的形状、部位、表面与流体之间的温差以及流体的流速等都有密切关系。物体表面附近的流体的流速愈大,其表面对流换热系数也愈大。如人处在风速较大的环境中,由于皮肤表面的对流换热系数较大,其散热(或吸热)量也较大。对流换热系数可用经验公式计算,通常用巴兹公式计算 1、详细内容 对流传热系数也称对流换热系数。对流换热系数的基本计算公式由牛顿于1701年提出,又称牛顿冷却定律。牛顿指出,流体与固体壁面之间对流传热的热流与它们的温度差成正比,即: q = h*(tw-t∞) Q = h*A*(tw-t∞)=q*A 式中: q为单位面积的固体表面与流体之间在单位时间内交换的热量,称作热流密度,单位W/m^2; tw、t∞分别为固体表面和流体的温度,单位K; A为壁面面积,单位m^2; Q为面积A上的传热热量,单位W; h称为表面对流传热系数,单位W/(m^2.K)。 2、理论发展 对流换热系数h的物理意义是:当流体与固体表面之间的温度差为1K时,1m*1m壁面面积在每秒所能传递的热量。h的大小反映对流换热的强弱。 如上所述,h与影响换热过程的诸因素有关,并且可以在很大的范围内变化,所以牛顿公式只能看作是传热系数的一个定义式。它既没有揭示影响对流换热的诸因素与h之间的内在联系,也没有给工程计算带来任何实质性的简化,只不过把问题的复杂性转移到传热系数的确定上去了。因此,在工程传热计算中,主要的任务是计算h。计算传热系数的方法主要有实验求解法、数学分析解法和数值分析解法。 影响对流传热强弱的主要因素有: 1. 对流运动成因和流动状态; 2. 流体的物理性质(随种类、温度和压力而变化); 3. 传热表面的形状、尺寸和相对位置; 4. 流体有无相变(如气态与液态之间的转化)。 3、实例应用 在不同的情况下,传热强度会发生成倍直至成千倍的变化,所以对流换热是一个受许多因素影响且其强度变化幅度又很大的复杂过程。 4、对流换热系数的大致量级: 空气自然对流 5 ~ 25 气体强制对流 20 ~ 100 水的自然对流 200 ~1000 水的强制对流 1000 ~ 15000

matlab有限元分析实例

MATLAB: MATLAB是美国MathWorks公司出品的商业数学软件,用于数据分析、无线通信、深度学习、图像处理与计算机视觉、信号处理、量化金融与风险管理、机器人,控制系统等领域。 MATLAB是matrix&laboratory两个词的组合,意为矩阵工厂(矩阵实验室),软件主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等。MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C,FORTRAN,C++,JAVA的支持。 MATLAB有限元分析与应用:

《MATLAB有限元分析与应用》是2004年4月清华大学出版社出版的图书,作者是卡坦,译者是韩来彬。 内容简介: 《MATLAB有限元分析与应用》特别强调对MATLAB的交互应用,书中的每个示例都以交互的方式求解,使读者很容易就能把MATLAB用于有限分析和应用。另外,《MATLAB有限元分析与应用》还提供了大量免费资源。 《MATLAB有限元分析与应用》采用当今在工程和工程教育方面非常流行的数学软件MATLAB来进行有限元的分析和应用。《MATLAB有限元分析与应用》由简单到复杂,循序渐进地介绍了各种有限元及其分析与应用方法。书中提供了大量取自机械工程、土木工程、航空航天工程和材料科学的示例和习题,具有很高的工程应用价值。

基于MpCCI的Abaqus和Fluent流固耦合案例1

CAE联盟论坛精品讲座系列 基于MpCCI的Abaqus和Fluent流固耦合案例 主讲人:mafuyin CAE联盟论坛总监 摘要:通过MpCCI流固耦合接口程序,对某薄壁管道流动中的传热过程进行了Abaqus和Fluent相结合的流固耦合仿真分析。信息介绍了从建模、设置到求解计算和后处理的全过程,对相关研究人员具有参考意义。 1 分析模型 用三维建模软件solidworks建立了一个管径为1m的弯管,结构尺寸如图1a所示,管的结构如图1b所示,流体的模型如图1c所示。值得注意的是,由于拓扑特征的原因,这样的管壁模型无法通过对圆环扫略直接生成,而需先通过对大圆的扫略生成实心的模型(类似于流体模型),然后进行抽壳得到管壁的模型。用同样的方法对大圆半径减去管壁厚度的圆进行扫略得到流体模型。 a. 尺寸关系 b. 管壁结构 c. 流体模型 图1. 几何模型示意图 图2. 流固耦合传热分析模型示意图 内壁面(耦合面) 速度入口 v=6m/s; T in=600K 外壁面 压力出口 P=0Pa;T out=300K

由于管壁结构和流体的热学行为不同,传热系数等都不一样,所以属于典型的流固耦合传热问题,热学模型如图2所示。即管的一端为流体速度入口,一端为压力出口,给定流体外壁面一个初始温度600K,流体入口速度为6m/s,温度为600K,出口相对大气压力为0Pa,出口温度为300K。需要求解流体和管壁的温度场分布情况。 2 流体模型 将图1c的流体模型以Step格式导入Fluent软件通常使用的前处理器Gambit中,如图3a所示。设置求解器为,然后划分体网格,网格尺寸为100mm,类型为六面体单元,一共生成4895个体单元,网格如图3b所示。 a. 导入Gambit软件中的流体模型 b. 流场的网格模型 图3. 流体模型及网格示意图 进行网格划分后,需定义边界条件,在Gambit软件中先分别定义速度入口(VELOCITY_INLET)、压力出口(PRESSURE_OUTLET)和壁面(Wall)三组边界条件,具体参数设置在Fluent软件中进行。然后定义流体属性,名称定义为air,类型为Fluid。这些定义的目的是能够在Fluent软件中识别出这些特征,具体类型和参数都可以在Fluent软件中进行设置和修改。定义完后点击【Export】,选择【Mesh】,选择路径和文件名称并进行输出。 打开Fluent6.3.26或以上的版本,选择3D求解器,点击【File】→【Read】→【Case】,然后选择Gambit中输出的msh文件,即可将网格文件读入Fluent 软件中。读入模型后,进行求解参数和条件的设置。

Ansys有限元分析实例[教学]

Ansys有限元分析实例[教学] 有限元分析案例:打点喷枪模组(用于手机平板电脑等电子元件粘接),该产品主要是使用压缩空气推动模组内的顶针作高频上下往复运动,从而将高粘度的胶水从喷嘴中打出(喷嘴尺寸,0.007”)。顶针是这个产品中的核心零件,设计使用材料是:AISI 4140 最高工作频率是160HZ(一个周期中3ms开3ms关),压缩空气压力3-8bar, 直接作用在顶针活塞面上,用Ansys仿真模拟分析零件的强度是否符合要求。 1. 零件外形设计图:

2. 简化模型特征后在Ansys14.0 中完成有限元几何模型创建:

3. 选择有限元实体单元并设定,单元类型是SOILD185,由于几何建模时使用的长度单位是mm, Ansys采用单位是长度:mm 压强: 3Mpa 密度:Ton/M。根据题目中的材料特性设置该计算模型使用的材料属性:杨氏模量 2.1E5; 泊松比:0.29; 4. 几何模型进行切割分成可以进行六面体网格划分的规则几何形状后对各个实体进行六面体网格划分,网格结果: 5. 依据使用工况条件要求对有限元单元元素施加约束和作用载荷:

说明: 约束在顶针底端球面位移全约束; 分别模拟当滑块顶断面分别以8Bar,5Bar,4Bar和3Bar时分析顶针的内应力分布,根据计算结果确定该产品允许最大工作压力范围。 6. 分析结果及讨论: 当压缩空气压力是8Bar时: 当压缩空气压力是5Bar时:

当压缩空气压力是4Bar时: 结论: 通过比较在不同压力载荷下最大内应力的变化发现,顶针工作在8Bar时最大应力达到250Mpa,考虑到零件是在160HZ高频率在做往返运动,疲劳寿命要求50百万次以上,因此采用允许其最大工作压力在5Mpa,此时内应力为156Mpa,按线性累积损伤理论[3 ]进行疲劳寿命L-N疲劳计算,进一部验证产品的设计寿命和可靠性。

血管流固耦合分析实例

Ansys14 workbench血管流固耦合实例 根据收集的一些资料,进行学习后,试着做了这个ansys14workbench的血管流固耦合模拟,感觉能够耦合上,仅是熟悉流固耦合分析过程,不一定正确,仅供参考,希望大家多讨论。谢谢! 1、先在proe5中建立血管与血液流体区的模型(两者装配起来),或者直接在workbench中建模。 图1 模型图 2、新建工程。在workbench中toolbox中选custom system,双击FSI: FluidFlow(fluent)->static structure. 图2 计算工程 3、修改engineering data,因为系统缺省材料是钢,需要构建血管材料,如图3所示。先复制steel,而后修改密度1150kg/m3,杨氏模量4.5e8Pa,泊松比0.3,重新命名,最后在主菜单中点击“update project”保存.

图3 修改工程材料 4、模型导入,进入gemetry模块,import外部模型文件。 图4 模型导入图 5、进入FLUENT网格划分。 在workbench工程视图中的Mesh上点击右键,选择Edit…,如图5所示,进入网格划分meshing界面,如图6所示。我们这里需要去掉血管部分,只保留血液几何。

图5 进入网格划分

图6 禁用血管模型 6、设置网格方法。 默认是采用ICEM CFD进行网格划分,设置方式如图7所示,截面圆弧边分为12份,纵截面的边均分为10份,网格结果如图8所示。另外在这个界面中要设置边界的几何面,如inlet、outlet、symmetry 图7 设置网格划分方式 图8 最终出网格

第二章 压电复合材料有限元分析方法

第二章压电复合材料有限元分析方法 2.1 1—3型压电复合材料常用的研究方法 第一、理论研究,包括利用细观力学和仿真软件进行数值分析的方法。人们对1-3型压电复合材料宏观等效特征参数进行研究时,从不同角度出发采用了形式多样的模型和理论,其中夹杂理论和均匀场理论具有代表性。夹杂理论的思想是,从细观力学出发,将1-3形压电复合材料的代表性体积单元(胞体)作为夹杂处理。求解过程中,使用的最著名的两个模型为:Dilute模型和Mori-Tanaka模型。夹杂理论的优点是其解析解能较好地反映材料的真实状况,解精度较高;缺点是其解题和计算过程烦琐,有时方程只能用数值方法求解。均匀场理论的思想是基于均匀场理论和混合定律,同时借助1-3型压电复合材料的细观力学模型导出其宏观等效特征参数。其基本的研究思路是:假设组成复合材料的每一相中力场和电场均匀分布,结合材料的本构方程得到1-3型压电复合材料的等效特征参数。Smith,Auld采用此理论研究了1-3型压电柱复合材料的弹性常数、电场、密度等等效特征参数。Gordon,John采用此理论研究了机电耦合系数、耗损因子、电学品质因子等等效特征参数。Bent, Hagood和Yoshikawa等基于此理论对交叉指形电极压电元件等效特征参数进行了研究。均匀场理论优点在于物理模型简单,物理概念清晰,计算也不复杂,并具有相当的精度和可靠性;不足在于其假设妨碍了两相分界面上的协调性。有限元作为一种广泛应用于解决实际问题的数值分析方法,将其引入压电复合材料研究中具有重要的意义。John,Gordon等用有限元方法分析了1-3型压电柱复合材料中压电柱为方形柱、圆形柱、二棱柱时的力电耦合系数及其波速特性,得到了压电柱在几何界面不同的情况下的等效力电耦合系数及等效波速曲线。 第二、实验研究。Helen,Gordon等对1-3型压电复合材料的宏观等效特征参数进行了理论和实验研究,结果表明两者符合良好;LVBT等运用了1-3型压电复合材料进行了声学方面的控制取得了良好的效果;John,Bent等对压电纤维复合材料的性能进行了深入的研究,结果显示压电纤维复合材料在高电场、大外载荷环境下具有优良的传感和作动性能。参数辨识研究是试验研究中重要的一种方法,基本思路是:分析1-3型压电纤维复合材料的响应特性,从中得到其等效宏观的模态和弹性波的传播特性参数。Guraja,Walter等采用的就是这种方法,他们研究了1-3型压电纤维复合材料薄板、厚板、变截面板的响应特性,得到了其相应的声波传播速度c,频率f,机械品质因素Q等参数的表达式,为1-3型压电纤维复合材料在超声波方面的应用提供了依据。 综合对比以上的研究方法,夹杂理论得出的结果比较接近实际结果,但是计算烦琐,而且对于高体积百分比的复合材料其计算结果跟实际相差较大;均匀场理论计算较为简单,但是模糊了两相材料之间的界面作用;实验研究方法是最接近实际的一种方法,但是由于实验条件、测试技术等一系列因素的制约使其不能广泛应用十实际中。由于交叉指形电极压电复合材料的复杂性,利用上面提到的夹杂理论和均匀场理论的方法,很难得到压电元件整体模型的性能状况。而数值研究有限元法,利用先进的分析软件ANSYS进行压电复合材料性能分析,可以超越目前现有的生产工艺和测试技术水平得到比较准确的分析结果,又可以减小压电元件的设计周期,减少实验制作压电元件的材料浪费和设备损耗。 2.2 有限元分析方法概述 有限元法(又称为有限单元法或有限元素法)是利用计算机进行数值模拟分析的方法。诞生于20世纪50年代初,最初只应用于力学领域中,现在广泛应用于结构、热、流体、电磁、声学等学科的设计分析及优化,有限元计算结果已成为各类工业产品设计和性能分析的

基于LSDYNA及FLUENT的板壳结构流固耦合分析

基于 LS-DYNA 及 FLUENT 的板壳结构流-固耦合分析
汪丽军 北京航空航天大学,交通科学与工程学院 100191
[摘 要]: 本文采用 ANSYS 显示动力分析模块 LS-DYNA 及流场分析模块 FLUENT,对水下的板壳 结构运动及其界面的流-固耦合现象进行了仿真分析。流场计算得到的界面压强数据以外载荷 的形式施加于结构表面,使其产生位移及变形;同时,结构的变化又进一步影响了流场的分 布。通过往复的双向耦合迭代,得到了板壳结构的动力学响应以及流场的分布情况。仿真结 果与试验结果的对比表明,此方法适用于解决兼有大位移及较大变形特征的流-固耦合问题。 [关键词]: 板壳结构 流-固耦合 有限元方法 ANSYS
Analysis of Fluid-Structure Interaction for Plate/Shell Structure Based on LS-DYNA and FLUENT
Wang Lijun School of Transportation Science & Engineering, Beihang University 100191
Abstract: In this paper,the movement of plate under water and the fluid-structure interaction(FSI) is simulated numerically by combining explicit dynamic solver LS-DYNA and computational fluid dynamics solver FLUENT in ANSYS. The pressure obtained from the calculation of flow field are applied as external loads on the surface of the plate, then the structural deformation and displacement can be calculated as well, which will affect the shape and pressure distribution of the flow field reversely. After sequential coupling iterations the dynamic response of the structure and flow field distribution are obtained consequently. By comparing numerical and experimental results it is proved that this proposed coupling method is suitable for solving such a kind of FSI problems considering both large displacement and comparatively large deformation. Keyword: Plate/shell structure, Fluid-Structure Interaction, Finite element method,ANSYS
1
前言
在自然界中,流-固耦合现象广泛存在于航空、航天、汽车、水利、石油、化工、海洋 以及生物等领域。很多实际问题中流体载荷对于结构的影响不可忽略;同时,结构的位移 和变形也会对流场的分布产生重要影响。例如各种水下运动机构都需要考虑这种现象。

有限元分析案例

有限元分析案例 图1 钢铸件及其砂模的横截面尺寸 砂模的热物理性能如下表所示: 铸钢的热物理性能如下表所示: 一、初始条件:铸钢的温度为2875o F,砂模的温度为80o F;砂模外边界的对流边界条件:对流系数0.014Btu/hr.in2.o F,空气温度80o F;求3个小时后铸钢及砂模的温度分布。 二、菜单操作: 1.Utility Menu>File>Change Title, 输入Casting Solidification; 2.定义单元类型:Main Menu>Preprocessor>Element Type>Add/Edit/Delete, Add, Quad 4node 55; 3.定义砂模热性能:Main Menu>Preprocessor>Material Props>Isotropic,默认材料编号1, 在Density(DENS)框中输入0.054,在Thermal conductivity (KXX)框中输入0.025,在S pecific heat(C)框中输入0.28; 4.定义铸钢热性能温度表:Main Menu>Preprocessor>Material Props>-Temp Dependent->Temp Table,输入T1=0,T2=2643, T3=2750, T4=2875; 5.定义铸钢热性能:Main Menu>Preprocessor>Material Props>-Temp Dependent ->Prop Table, 选择Th Conductivity,选择KXX, 输入材料编号2,输入C1=1.44, C2=1.54, C3=1.22, C4=1.22,选择Apply,选择Enthalpy,输入C1=0, C2=128.1, C3=163.8, C4=174.2; 6.创建关键点:Main Menu>Preprocessor>-Modeling->Create>Keypoints>In Active

板结构有限元分析实例详解

板结构有限元分析实例详解1:带孔平板结构静力分析本节介绍带孔平板结构静力分析问题,同时介绍布尔操作的基本用法。 8.3.1 问题描述与分析 有孔的矩形平板,左侧边缘固定,长400mm,宽200 mm,厚度为10 mm,圆孔在板的正中心,半径为40 mm,左侧全约束,右侧边缘均布应力1MPa,如图8.7所示。求板的变形、位移及应力变化情况。(材料的材料属性为:弹性模量为300000 MPa,剪切模量为0.31。) 图8.7 带孔的矩形平板 由于小孔处边缘不规则,本文采用PLANE82高阶平面单元进行分析。 8.3.2 求解过程 8.3.2.1 定义工作目录及文件名 启动ANSYS Mechanical APDL Product Launcher窗口,如图8.8所示。在License下 拉选框中选择ANSYS Multiphysics产品,在Working Directory输入栏中输入工作目 录:C:\ANSYS12.0 Structural Finite Elements Analysis and Practice\Chapter 8\8-1,在Job Name一栏中输入工作文件名:Chapter8-1。以上参数设置完毕后,单 击Run按钮运行ANSYS。

图8.8 ANSYS设置窗口菜单 可以先在目标文件位置建立工作目录,然后单击Browse按钮选择工作目录;也 可以通过单击Browse按钮选择工作文件名。 8.3.2.2 定义单元类型和材料属性 选择Main Menu>Preferences命令,出现Preferences for GUI Filtering对话框, 如图8.9所示,在Individual discipline(s) to show in the GUI中勾选Structural,过滤掉ANSYS GUI菜单中与结构分析无关的选项,单击OK按钮关闭该对话框。 图8.9 Preferences for GUI Filtering对话框

abaqus与fluent流固耦合

基于MPCCI的流固耦合成功案例 基于MPCCI的流固耦合成功案例 (一)机翼气动弹性分析 1 问题陈述 机翼绕流问题是流固耦合中的经典问题。以前由于缺乏考虑流固耦合的软件,传统的分析方法是将机翼视为刚体,不考虑其弹性变形,通过CFD软件来计算机翼附近的流场。这个强硬的假设很难准确的描述流场的实际情况。更无法预测机翼的振动。MPCCI是基于代码耦合的并行计算接口,它可以同时调用结构和流体的软件来实现流固耦合。我们通过MPCCI,能很好的预测真实情况下的机翼绕流问题。采用ABAQUS结构分析软件来求解结构在流畅作用下的变形和应力分布,通过Fluent软件来计算由于固体运动和变形对整个流场的影响。 2 模拟过程分析顺序 MpCCI的图形用户界面可以方便的读入结构和流体的输入文件。后台调用ABAQUS和FLUENT。在MPCCI耦合面板中选择耦合面,然后选择在相应耦合面上流体和固体需要交换的量。启动MpCCI进行耦合。 3 边界条件设置

图1 无人机模型和流体计算模型 结构部分单个机翼跨度在1.5m左右,厚度为0.1m左右。边界条件为机翼端部的固定,三个方向的位移完全固定,另一端完全自由。在固体中除了固定端的面外,其他三个面为耦合面。流体部分采用四面体网格,采用理想气体作为密度模型。流体的入口和出口以及对称性边界条件如下图所示。 图2 固体有限元模型 4 计算方法的选择 通过结合ABAQUS和FLUENT,使用MPCCI计算流固耦合。在本例中,固体在流场作用下产生很大的变形和运动。在耦合区域,固体结构部分计算耦合面上的节点位移,通过MPCCI传输给FLUENT的耦合界面,FLUENT 计算出耦合区域上的节点力载荷,然后通过MPCCI传给结构软件ABAQUS。在MPCCI的耦合面板中选择的耦合面如图所示,交换量为:节点位移、相对受力。采用ABAQUS中的STANDARD算法,时间增量步长为0.1毫秒。 5 计算结论 通过MPCCI结合ABAQUS和FLUENT,成功地计算在几何非线性条件下的气动弹性问题,得到了整个流体区域的流场分布以及结构的动态响应历程。

基于NX有限元分析实验报告

有限元分析及应用 专业:机械 姓名:你喝 学号:2 0 1 3 X X 指导老师:没意义 工字梁热力学与结构学耦合分析

有限元分析(FEA,Finite Element Analysis)将物体划分成有限个单元,这些单元之间通过有限个节点相互连接,单元看作是不可变形的刚体,单元之间的力通过节点传递,然后利用能量原理建立各单元矩阵;在输入材料特性、载荷和约束等边界条件后,利用计算机进行物体变形、应力和温度场等力学特性的计算,最后对计算结果进行分析,显示变形后物体的形状及应力分布图。有限元分析的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的(较简单的)近似解,然后推导求解这个域总的满足条件(如结构的平衡条件),从而得到问题的解。这个解不是准确解,而是近似解,因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。 热——结构耦合分析是指求解温度场对结构中应力、应变和位移等物理量的影响,热——结构耦合问题是结构分析中较常见的一类耦合分析问题。由于结构温度场的分布不均会引起结构的热应力,或者是结构件在高温环境中工作,材料受到温度的影响会发生性能的改变,这些都是进行结构分析时需要考虑的因素。为此需要先进行相应的热分析,然后再进行结构分析。在NX环境中进行热——结构耦合分析,首先进行热分析求得结构的温度场,然后再进行结构分析,并将前面得到的温度场作为体载荷加到结构中,求解结构的应力分布。 1.模型建立 2.热分析 新建FEM和仿真 点击开始按钮,选择“高级仿真”,激活高级仿真模块。在仿真导航器中选择“新建FEM 和仿真” 解算方案 网格收集器 添加材料属性,从材料清单中选择“Steel”,单击“确定” 划分网格 添加约束(进入仿真环境) 所有外表面添加对流约束,环境温度为45,对流系数为100W/m^2-C 添加热约束 在工字梁顶端设置65恒温 解算方案求解

ANSYS有限元分析实例

有限元分析 一个厚度为20mm的带孔矩形板受平面内张力,如下图所示。左边固定,右边受载荷p=20N/mm作用,求其变形情况 P 一个典型的ANSYS分析过程可分为以下6个步骤: ①定义参数 ②创建几何模型 ③划分网格 ④加载数据 ⑤求解 ⑥结果分析 1定义参数 1.1指定工程名和分析标题 (1)启动ANSYS软件,选择File→Change Jobname命令,弹出如图所示的[Change Jobname]对话框。 (2)在[Enter new jobname]文本框中输入“plane”,同时把[New log and error files]中的复选框选为Yes,单击确定 (3)选择File→Change Title菜单命令,弹出如图所示的[Change Title]对话框。 (4)在[Enter new title]文本框中输入“2D Plane Stress Bracket”,单击确定。 1.2定义单位

在ANSYS软件操作主界面的输入窗口中输入“/UNIT,SI” 1.3定义单元类型 (1)选择Main Menu→Preprocessor→Element Type→Add/Edit/Delete命令,弹出如图所示[Element Types]对话框。 (2)单击[Element Types]对话框中的[Add]按钮,在弹出的如下所示[Library of Element Types]对话框。 (3)选择左边文本框中的[Solid]选项,右边文本框中的[8node 82]选项,单击确定,。 (4)返回[Element Types]对话框,如下所示 (5)单击[Options]按钮,弹出如下所示[PLANE82 element type options]对话框。

几个ansys流固耦合的例子

一般说来,ANSYS的流固耦合主要有4种方式: 1,sequential 这需要用户进行APDL编程进行流固耦合 sequentia指的是顺序耦合 以采用MpCCI为例,你可以利用ANSYS和一个第三方CFD产品执行流固耦合分析。在这个方法中,基于网格的平行代码耦合界面(MpCCI) 将ANSYS和CFD程序耦合起来。即使网格上存在差别,MpCCI也能够实现流固界面的数据转换。ANSYS CD中包含有MpCCI库和一个相关实例。关于该方法的详细信息,参见ANSYS Coupled-Field Analysis Guide中的Sequential Couplin 2,FSI solver 流固耦合的设置过程非常简单,推荐你使用这种方式 3,multi-field solver 这是FSI solver的扩展,你可以使用它实现流体,结构,热,电磁等的耦合 4,直接采用特殊的单元进行直接耦合,耦合计算直接发生在单元刚度矩阵 一个流固耦合的例子 length=2 width=3 height=2 /prep7 et,1,63 et,2,30 !选用FLUID30单元,用于流固耦合问题 r,1,0.01 mp,ex,1,2e11 mp,nuxy,1,0.3 mp,dens,1,7800 mp,dens,2,1000 !定义Acoustics材料来描述流体材料-水 mp,sonc,2,1400 mp,mu,0, ! block,,length,,width,,height esize,0.5 mshkey,1 ! type,1 mat,1 real,1 asel,u,loc,y,width amesh,all alls ! type,2 mat,2 vmesh,all

有限元法的应用现状研究

有限元法的应用现状研究3 于亚婷,杜平安,王振伟 (电子科技大学机械电子工程学院,四川成都 610054) 摘要:有限元法(FEM:Finite Element Method)作为一种最有效的数值方法,在工程实际中得到了广泛、深入的应用。以应用为主线,首先回顾了FEM的发展历程,然后从FEM的应用过程和应用领域两个方面详细地论述了FEM应用的有关问题,并例举了相关的应用实例。最后总结了FEM的国内外应用现状及研究热点问题。 关键词:有限元法;应用过程;应用领域;现状 中图分类号:TP391.7 文献标识码:A 文章编号:1001-2354(2005)03-0006-04 1 F EM的发展历程 FEM作为求解数学物理问题的一种数值方法,已经历了50余年的发展。20世纪50年代,它作为处理固体力学问题的方法出现。1943年,Courant第一次提出单元概念[1]。1945~1955年,Argyris等人在结构矩阵分析方面取得了很大进展[1]。1956年,Turner、Clough等人把刚架位移法的思路推广应用于弹性力学平面问题[1]。1960年,Clough首先把解决弹性力学平面问题的方法称为“有限元法”[1],并描绘为“有限元法= Rayleigh Ritz法+分片函数”。几乎与此同时,我国数学家冯康也独立提出了类似方法。 FEM理论研究的重大进展,引起了数学界的高度重视。自20世纪60年代以来,人们加强了对FEM数学基础的研究。如大型线性方程组和特征值问题的数值方法、离散误差分析、解的收敛性和稳定性等。 FEM理论研究成果为其应用奠定了基础,计算机技术的发展为其提供了条件。20世纪70年代以来,相继出现了一些通用的有限元分析(FEA:Finite Element Analysis)系统,如SA P、ASKA、NASTRAN等,这些FEA系统可进行航空航天领域的结构强度、刚度分析,从而推动了FEM在工程中的实际应用。 20世纪80年代以来,随着工程工作站的出现和广泛应用,原来运行于大中型机上的FEA系统得以在其上运行,同时也出现了一批通用的FEA系统,如ANSYS-PC、N ISA,SU PER SA P等[3]。20世纪90年代以来,随着微机性能的显著提高,大批FEA系统纷纷向微机移植,出现了基于Windows的微机版FEA系统。 经过半个多世纪的发展,FEM已从弹性力学平面问题扩展到空间问题、板壳问题;从静力问题扩展到动力问题、稳定问题和波动问题;从线性问题扩展到非线性问题;从固体力学领域扩展到流体力学、传热学、电磁学等其他连续介质领域;从单一物理场计算扩展到多物理场的耦合计算。它经历了从低级到高级、从简单到复杂的发展过程,目前已成为工程计算最有效的方法之一。 2 F EM应用的有关问题 2.1 FEM的应用过程 FEM应用于实际问题须经历以下过程,如图1所示 。 图1 FEM的应用过程 (1)问题的数学描述。对问题客观规律的数学描述(通常是微分方程及边界条件)是建立有限元方程的前提。单元特性矩阵和整体有限元方程都是基于数学模型建立的。常见的弹性力学基本方程、运动方程、热传导方程等都是对客观现象的数学描述。 (2)有限元方程的建立。利用变分原理,通过离散、单元分析、整体分析等过程,建立数学模型的有限元方程,它通常是一组易于用数值方法求解的代数方程。 (3)算法研究。有限元方程的计算量庞大,须有有效的算法来保证计算效率和精度,同时考虑对计算条件的要求。如求解大型线性方程组的带宽法、波前法,求解大型特征值问题的分块Lanczos法等。 第22卷第3期2005年3月 机 械 设 计 J OU RNAL OF MACHIN E DESIGN Vol.22 No.3 Mar. 2005 3收稿日期:2004-07-22;修订日期:2004-09-02 基金项目:四川省学术与技术带头人培养基金资助项目(2200104) 作者简介:于亚婷(1979-),女,陕西人,电子科技大学博士研究生,研究方向:CAD/CAM/CA E有限元法应用等。

ansys有限元分析作业经典案例

有 限 元 分 析 作 业 作业名称 输气管道有限元建模分析 姓 名 陈腾飞 学 号 3070611062 班 级 07机制(2)班 宁波理工学院

题目描述: 输气管道的有限元建模与分析 计算分析模型如图1所示 承受内压:1.0e8 Pa R1=0.3 R2=0.5 管道材料参数:弹性模量E=200Gpa;泊松比v=0.26。 图1受均匀内压的输气管道计算分析模型(截面图) 题目分析: 由于管道沿长度方向的尺寸远远大于管道的直径,在计算过程中忽略管道的断面效应,认为在其方向上无应变产生。然后根据结构的对称性,只要分析其中1/4即可。此外,需注意分析过程中的单位统一。 操作步骤 1.定义工作文件名和工作标题 1.定义工作文件名。执行Utility Menu-File→Chang Jobname-3070611062,单击OK按钮。 2.定义工作标题。执行Utility Menu-File→Change Tile-chentengfei3070611062,单击OK按钮。 3.更改目录。执行Utility Menu-File→change the working directory –D/chen 2.定义单元类型和材料属性 1.设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK

2.选择单元类型。执行ANSYS Main Menu→Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 8node 82 →apply Add/Edit/Delete →Add →select Solid Brick 8node 185 →OK Options…→select K3: Plane strain →OK→Close如图2所示,选择OK接受单元类型并关闭对话框。 图2 3.设置材料属性。执行Main Menu→Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic,在EX框中输入2e11,在PRXY框中输入0.26,如图3所示,选择OK并关闭对话框。 图3 3.创建几何模型 1. 选择ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标:input:1(0.3,0),2(0.5,0),3(0,0.5),4(0,0.3) →OK

相关文档
最新文档