边界条件的使用方法

边界条件的使用方法
边界条件的使用方法

FLUENT中各种边界条件的适用范围

速度入口边界条件:用于定义流动入口边界的速度和标量。

压力入口边界条件:用来定义流动入口边界的总压和其它标量。

质量流动入口边界条件:用于已知入口质量流速的可压缩流动。在不可压缩流动中不必指定入口的质量流,因为当密度是常数时,速度入口边界条件就确定了质量流条件。

压力出口边界条件:用于定义流动出口的静压(在回流中还包括其它的标量)。当出现回流时,使用压力出口边界条件来代替质量出口条件常常有更好的收敛速度。

压力远场边界条件:用于模拟无穷远处的自由可压缩流动,该流动的自由流马赫数以及静态条件已知。这一边界类型只用于可压缩流。

质量出口边界条件:用于在解决流动问题之前,所模拟的流动出口的流速和压力的详细情况还未知的情况。在流动出口是完全发展的时候这一条件是适合的,这是因为质量出口边界条件假定出了压力之外的所有流动变量正法向梯度为零。不适合于可压缩流动。进风口边界条件:用于模拟具有指定的损失系数、流动方向以及周围(入口)环境总压和总温的进风口。

进气扇边界条件:用于模拟外部进气扇,它具有指定的压力跳跃、流动方向以及周围(进口)总压和总温。

通风口边界条件:用于模拟通风口,它具有指定的损失系数以及周围环境(排放处)的静压和静温。

排气扇边界条件:用于模拟外部排气扇,它具有指定的压力跳跃以及周围环境(排放处)的静压。

速度入口边界条件:速度入口边界条件用于定义流动速度以及流动入口的流动属性相关标量。这一边界条件适用于不可压缩流,如果用于可压缩流它会导致非物理结果,这是因为它允许驻点条件浮动。应该注意不要让速度入口靠近固体妨碍物,因为这会导致流动入口驻点属性具有太高的非一致性。

压力入口边界条件:压力入口边界条件用于定义流动入口的压力以及其它标量属性。它即可以适用于可压缩流,也可以用于不可压缩流。压力入口边界条件可用于压力已知但是流动速度和/或速率未知的情况。这一情况可用于很多实际问题,比如浮力驱动的流动。压力入口边界条件也可用来定义外部或无约束流的自由边界。

质量流动入口边界条件:用于已知入口质量流速的可压缩流动。在不可压缩流动中不必指定入口的质量流,因为当密度是常数时,速度入口边界条件就确定了质量流条件。当要求达到的是质量和能量流速而不是流入的总压时,通常就会使用质量入口边界条件。调节入口总压可能会导致解的收敛速度较慢,所以如果压力入口边界条件和质量入口条件都可

以接受,应该选择压力入口边界条件。

压力出口边界条件:压力出口边界条件需要在出口边界处指定静(gauge)压。静压值的指定只用于亚声速流动。如果当地流动变为超声速,就不再使用指定压力了,此时压力要从内部流动中推断。所有其它的流动属性都从内部推出。在解算过程中,如果压力出口边界处的流动是反向的,回流条件也需要指定。如果对于回流问题指定了比较符合实际的值,收敛性困难就会被减到最小。

压力远场边界条件:FLUENT中使用的压力远场条件用于模拟无穷远处的自由流条件,其中自由流马赫数和静态条件被指定了。压力远场边界条件通常被称为典型边界条件,这是因为它使用典型的信息(黎曼不变量)来确定边界处的流动变量。这一边界条件只应用于密度是用理想气体定律计算出来的情况,不可以适用于其它情况要有效地近似无限远处的条件,必须建立的这个远场放到所关心的计算物体的足够远处。例如,在机翼升力计算中远场边界一般都要设到20倍弦长的圆周之外。

质量出口边界条件:当流动出口的速度和压力在解决流动问题之前是未知时,FLUENT会使用质量出口边界条件来模拟流动。你不需要定义流动出口边界的任何条件(除非你模拟辐射热传导、粒子的离散相或者分离质量流):FLUENT会从内部推导所需要的信息。然而,重要的是要知道这一边界类型的限制。

注意:下面的几种情况不能使用质量出口边界条件: 1.如果包含压力出口,请使用压力出口边界条件 2.如果模拟可压缩流 3.如果模拟变密度的非定常流,即使流动是不可压的也不行。

各类边界条件fluent

Fluent技巧 边界条件 定义边界条件概述 边界条件包括流动变量和热变量在边界处的值。它是FLUENT分析得很关键的一部分,设定边界条件必须小心谨慎。 边界条件的分类:进出口边界条件:压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件、质量出口、通风口、排气扇;壁面、repeating, and pole boundaries:壁面,对称,周期,轴;内部单元区域:流体、固体(多孔是一种流动区域类型) ;内部表面边界:风扇、散热器、多孔跳跃、壁面、内部。(内部表面边界条件定义在单元表面,这意味着它们没有有限厚度,并提供了流场性质的每一步的变化。这些边界条件用来补充描述排气扇、细孔薄膜以及散热器的物理模型。内部表面区域的内部类型不需要你输入任何东西。) 下面一节将详细介绍上面所叙述边界条件,并详细介绍了它们的设定方法以及设定的具体合适条件。周期性边界条件在本章中介绍,模拟完全发展的周期性流动将在周期性流动和热传导一章中介绍。 使用边界条件面板 边界条件(Figure 1)对于特定边界允许你改变边界条件区域类型,并且打开其他的面板以设定每一区域的边界条件参数 菜单:Define/Boundary Conditions... Figure 1: 边界条件面板 改变边界区域类型 设定任何边界条件之前,必须检查所有边界区域的区域类型,如有必要就作适当的修改。比方说:如果你的网格是压力入口,但是你想要使用速度入口,你就要把压力入口改为速度入口之后再设定。 改变类型的步骤如下:: 1.在区域下拉列表中选定所要修改的区域 2.在类型列表中选择正确的区域类型 3.当问题提示菜单出现时,点击确认 确认改变之后,区域类型将会改变,名字也将自动改变 (如果初始名字时缺省的请参阅边界条件区域名字一节),设定区域边界条件的面板也将自动打开。 !注意:这个方法不能用于改变周期性类型,因为该边界类型已经存在了附加限制。创建边界条件一节解释了如何创建和分开周期性区域。需要注意的是,只能在图一中每一个类别中改变边界类型(注意:双边区域表面是分离的不同单元区域.) Figure 1: 区域类型的分类列表 设定边界条件 在FLUENT中,边界条件和区域有关而与个别表面或者单元无关。如果要结合具有相同边界条件的两个或更多区域请参阅合并区域一节。 设定每一特定区域的边界条件,请遵循下面的步骤: 1.在边界条件区域的下拉列表中选择区域。 2. 点击Set...按钮。或者,1.在区域下拉列表中选择区域。 2.在类型列表中点击所要选择的类型。或者在区域列表中双击所需区域.,选择边界条件区域将会打开,并且你可以指定适当的边界条件

FLUENT中各种边界条件的适用范围

FLUENT中各种边界条件的适用范围 速度入口边界条件:用于定义流动入口边界的速度和标量。 压力入口边界条件:用来定义流动入口边界的总压和其它标量。 质量流动入口边界条件:用于已知入口质量流速的可压缩流动。在不可压缩流动中不必指定入口的质量流,因为当密度是常数时,速度入口边界条件就确定了质量流条件。压力出口边界条件:用于定义流动出口的静压(在回流中还包括其它的标量)。当出现回流时,使用压力出口边界条件来代替质量出口条件常常有更好的收敛速度。 压力远场边界条件:用于模拟无穷远处的自由可压缩流动,该流动的自由流马赫数以及静态条件已知。这一边界类型只用于可压缩流。 质量出口边界条件:用于在解决流动问题之前,所模拟的流动出口的流速和压力的详细情况还未知的情况。在流动出口是完全发展的时候这一条件是适合的,这是因为质量出口边界条件假定出了压力之外的所有流动变量正法向梯度为零。不适合于可压缩流动。 进风口边界条件:用于模拟具有指定的损失系数、流动方向以及周围(入口)环境总压和总温的进风口。 进气扇边界条件:用于模拟外部进气扇,它具有指定的压力跳跃、流动方向以及周围(进口)总压和总温。 通风口边界条件:用于模拟通风口,它具有指定的损失系数以及周围环境(排放处)的静压和静温。 排气扇边界条件:用于模拟外部排气扇,它具有指定的压力跳跃以及周围环境(排放处)的静压。 速度入口边界条件:速度入口边界条件用于定义流动速度以及流动入口的流动属性相关标量。这一边界条件适用于不可压缩流,如果用于可压缩流它会导致非物理结果,这是因为它允许驻点条件浮动。应该注意不要让速度入口靠近固体妨碍物,因为这会导致流动入口驻点属性具有太高的非一致性。 压力入口边界条件:压力入口边界条件用于定义流动入口的压力以及其它标量属性。它即可以适用于可压缩流,也可以用于不可压缩流。压力入口边界条件可用于压力已知但是流动速度和/或速率未知的情况。这一情况可用于很多实际问题,比如浮力驱动的流动。压力入口边界条件也可用来定义外部或无约束流的自由边界。 质量流动入口边界条件:用于已知入口质量流速的可压缩流动。在不可压缩流动中不必指定入口的质量流,因为当密度是常数时,速度入口边界条件就确定了质量流条件。当要求达到的是质量和能量流速而不是流入的总压时,通常就会使用质量入口边界条件。调节入口总压可能会导致解的收敛速度较慢,所以如果压力入口边界条件和质量入口条件都可以接受,应该选择压力入口边界条件。 压力出口边界条件:压力出口边界条件需要在出口边界处指定静(gauge)压。静压值的指定只用于亚声速流动。如果当地流动变为超声速,就不再使用指定压力了,此时压力要从内部流动中推断。所有其它的流

COMSOL周期性边界条件的应用

COMSOL周期性边界条件的应用 在将真实的物理问题转化为仿真模型时,为了通过有限的计算资源获得尽可能高的计算精度,模型简化是必要的。模型简化的前提是所模拟的物理问题具有结构、材料属性及边界条件的对称性或均匀性,以此为基础,可通过特定的方程及边界条件建立模型,例如降维方程,镜像/周期性/旋转对称边界条件,或根据工程经验将某些计算域简化为边界等等。 当处理空间或时间上具有周期性的物理问题时,采用周期性边界条件(Periodic/Cyclic Condition),可将复杂结构的模拟简化为周期单元,在不失精确度的前提下,大大降低计算量。 COMSOL提供的周期性边界条件包括四种类型: ?连续性周期边界(Continuity),指在源和目标边界上的场值相等; ?反对称周期边界(Antiperiodicity),源和目标边界上场值符号相反; ?弗洛奎特周期性边界(Floquet periodicity),源和目标边界上场值相差一个位相因子,位相因子由波矢和边界相对距离确定。Continuity和Antiperiodicity边界可以认 为是Floquet periodicity边界在位相分别为0和π情况下的两个特例。 ?循环对称性边界(Cyclic Symmetry),源和目标边界上场值相差一个位相因子,位相因子由计算域所对应的扇形角和角向模式数决定。 以下是几个典型应用: 1.微纳光学领域内的光子晶体(Photonic Crystal)、表面等离子体激元(Surface Plasmon) 阵列结构及超材料(Metamaterial),这几种结构均由空间上周期性重复的散射体构 成,当计算透射率及能带结构时,常常可采用Floquet perioidcity边界将结构简化。 超材料能带分析 Metamaterial.mph 2.作为压电传感器件的声表面波器件(Surface Acoustic Wave, SAW)的本征频率问题 计算。

基于有限元和边界元的噪声分析

half 重登录 隐身 控制面板 搜索 状态 展区 振动博客 论坛服务 退出 振动论坛 → 专题讨论区→ 噪声分析及控制→声学基础理论→[转帖]基于有限元和边界元的噪声 分析 复制本页地址 粘贴我的收件箱 (0) 您 是本帖的第42个阅读者 标题:[转帖]基于有限元和边界元的噪声 分析树形 打印 收藏 推荐 提交网摘 等级:本科生 威望:18 现金:308 经验:1107 魅力:627 文章:109 注册: 2005-07-24 活跃度: 活跃等级:①年迈乌龟 在线等级: van321 ▼楼主 物体受到激励后,必将会产生振动,由物体的振动而引起与之相接触的流体的振动(如空气),从而在流体中产生噪声。对流体的噪声分析可以在频率域内或者时间域内进行,可以采用流体与结构耦合的形式进行分析,也可以只采用流体的形式进行计算分析,可以计算内声场也可以计算外声场,例如对于汽车而言,可以计算内声场,也可以计算外声场。在低频范围内采用边界元或者有限元的方法,在高频内采用统计能量的方法,计算结果包括声场中任意一点处的声学响应,如声压、声强、声功率,还可以是某点处的响应函数,如声压函数、模态贡献量函数,还可以进行一些特殊的分析,如声学传递矢量分析、面板贡 献量分析和灵敏*分析,以及高频域内的统计 能量分析。 如图所示是某轿车的排气系统的有限元声学模型,图所示是该排气系统中消声器的声学 模型。 [转帖]基于有限元和边界元的噪声分析

排气系统的声学模型

消声器的声学模型 ?声学模态分析 声学模态类似于结构模态,声波在流体团中传播时,会引发流体的振荡,流体的振荡也是有一定的固有频率和振动样式(振型),通过声学模态计算可以计算出流体的声学共振频率,防止流体和流体周围的结构产生共振而引发共鸣。 图所示是排气系统的声学模态云纹图。

fluent边界条件(二)

周期性边界条件 周期性边界条件用来解决,物理模型和所期待的流动的流动/热解具有周期性重复的特点。FLUENT提供了两种类型的周期性边界条件。第一种类型不允许通过周期性平面具有压降(对于FLUENT4用户来说:这一类型的周期性边界是指FLUENT4中的圆柱形边界)。第二种类型允许通过平移周期性边界具有压降,它是你能够模拟完全发展的周期性流动(在FLUENT4中是周期性边界)。 本节讨论了无压降的周期性边界条件。在周期性流动和热传导一节中,完全发展的周期性模拟能力得到了详尽的描述。 周期性边界的例子 周期性边界条件用于模拟通过计算模型内的两个相反平面的流动是相同的情况。下图是周期性边界条件的典型应用。在这些例子中,通过周期性平面进入计算模型的流动和通过相反的周期性平面流出流场的流动是相同的。正如这些例子所示,周期性平面通常是成对使用的。 Figure 1: 在圆柱容器中使用周期性边界定义涡流 周期性边界的输入 对于没有任何压降的周期性边界,你只需要输入一个东西,那就是你的所模拟的几何外形是旋转性周期还是平移性周期。(对于有周期性压降的周期流还要输入其它的东西,请参阅周期性流动和热传导一节。) 旋转性周期边界是指关于旋转对称几何外形中线形成了一个包括的角度。本节中的图一就是旋转性周期。平移性周期边界是指在直线几何外形内形成周期性边界。下面两图是平移性周期边界:

Figure 1: 物理区域 Figure 2: 所模拟的区域 对于周期性边界,你需要在周期性面板(下图)中指定平移性边界还是旋转性边界,该面板是从设定边界条件菜单中打开的。 Figure 3: 周期性面板 (对于耦合解算器,周期性面板中将会有附加的选项,这一选项允许你指定压力跳跃,详细内容请参阅周期性流动和热传导一节。) 如果区域是旋转性区域,请选择旋转性区域类型。如果是平移性就选择平移性区域类型。对

FLUENT进行流体动力学分析时,分析边界条件的种类及应用要点

FLUENT进行流体动力学分析时,分析边界条件的种类及应用要点。答:FLUENT 软件提供了十余种类型的进、出口边界条件,分别如下: (1) 速度入口(velocity-inlet):给出入口边界上的速度。 给定入口边界上的速度及其他相关标量值。该边界条件适用于不可压速流动问题,对可压缩问题不适合,否则该入口边界条件会使入口处的总温或总压有一定的波动。 (2) 压力入口(pressure-inlet):给出入口边界上的总压。 压力入口边界条件通常用于流体在入口处的压力为已知的情形,对计算可压和不可压问题都适合。压力进口边界条件通常用于进口流量或流动速度为未知的流动。压力入口条件还可以用于处理自由边界问题。 (3) 质量入口(mess-flow-inlet):给出入口边界上的质量流量。 质量入口边界条件主要用于可压缩流动;对于不可压缩流动,由于密度是常数,可以用速度入口条件。质量入口条件包括两种:质量流量和质量通量。质量流量是单位时间内通过进口总面积的质量。质量通量是单位时间单位面积内通过的质量。如果是二维轴对称问题,质量流量是单位时间内通过2π弧度的质量,而质量通量是通过单位时间内通过1 弧度的质量。 (4) 压力出口(pressure-outlet):给定流动出口边界上的静压。 对于有回流的出口,该边界条件比outflow 边界条件更容易收敛。给定出口边界 上的静压强(表压强)。该边界条件只能用于模拟亚音速流动。如果当地速度已经超过音速,该压力在计算过程中就不采用了。压力根据内部流动计算结果给定。其他量都是根据内部流动外推出边界条件。该边界条件可以处理出口有回流问题,合理的给定出口回流条件,有利于解决有回流出口问题的收敛困难问题。(5) 无穷远压力边界 (pressure-far-field):该边界条件用于可压缩流动。 如果知道来流的静压和马赫数,FLUENT 提供了无穷远压力边界条件来模拟该类问题。该边界条件适用于用理想气体定律计算密度的问题。为了满足无穷远压力边界条件,需要把边界放到我们关心区域足够远的地方。

[考试]在fluent中修改周期性边界条件

[考试]在fluent中修改周期性边界条件中国振动联盟 标题: 在fluent中修改周期性边界条件,怎么不对啊 [打印本页] 作者: skgk-qqq 时间: 2012-2-26 09:39 标题: 在fluent中修改周期性边界 条件,怎么不对啊 我是在fluent主界面输入命令:grid mod check,然后回车,得到periodic zone[()],我再输入3,回车,shadow zonezone[()],我再输入10,回车,得到Rottional periodic,(if no,translational)[yes],然后回车,得到Create periodic zones?[yes],然后回车,得到zone 3;matched 0 out of 10854 faces. zone 10:matched 0 out of 10854 faces. Error: Failed to make zones periodic.ERROE:object:#f.请教各位了,着急啊~~~ 作者: skgk-qqq 时间: 2012-2-26 09:51 回复 1 # skgk-qqq 的帖子 各位大哥,帮帮忙啊,着急啊 作者: Seventy721 时间: 2012-2-26 11:01 大概是因为你的两个periodic面上的网格不完全一致,导致不能match。这两 个面的几何尺寸和网格划分必须完全一致。建议划分网格之前在两个面上建立 hard link,这样网格就会完全一样了。如果还不行就调整判断网格差异的tolerance,我记得用户手册里有说明,你找找看。 作者: skgk-qqq 时间: 2012-2-26 16:15 回复 3 # Seventy721 的帖子 我已经建立了link了啊,经过网格检查,网格单元数量也是一致的,而且输 出meh文件也正确,请问怎么调整tolerance啊,着急啊

浅话边界条件与初始条件

浅话边界条件与初始条件 边界条件 在说边界条件之前,先谈谈初值问题和边值问题。 初值和边值问题: 对一般的微分方程,求其定解,必须引入条件,这个条件大概分两类---初始条件和边界条件,如果方程要求未知量y(x)及其导数y′(x)在自变量的同一点x=x0取给定的值,即y(x0 )=y0,y′(x0)= y0′,则这种条件就称为初始条件,由方程和初始条件构成的问题就称为初值问题; 而在许多实际问题中,往往要求微分方程的解在在某个给定的区间a ≤ x≤b 的端点满足一定的条件,如y(a) = A , y(b) = B则给出的在端点(边界点)的值的条件,称为边界条件,微分方程和边界条件构成数学模型就称为边值问题。 三类边界条件: 边值问题中的边界条件的形式多种多样,在端点处大体上可以写成这样的形式,Ay+By'=C,若B=0,A≠0,则称为第一类边界条件或狄里克莱(Dirichlet)条件;B≠0,A=0,称为第二类边界条件或诺依曼(Neumann)条件;A≠0,B≠0,则称为第三类边界条件或洛平(Robin)条件。 总体来说, 第一类边界条件:给出未知函数在边界上的数值; 第二类边界条件:给出未知函数在边界外法线的方向导数; 第三类边界条件:给出未知函数在边界上的函数值和外法向导数的线性组合。 对应于comsol,只有两种边界条件: Dirichlet boundary(第一类边界条件)—在端点,待求变量的值被指定。

Neumann boundary(第二类边界条件)—待求变量边界外法线的方向导数被指定。 再补充点初始条件: 初始条件,是指过程发生的初始状态,也就是未知函数及其对时间的各阶偏导数在初始时刻t=0的值.在有限元中,好多初始条件要预先给定的。不同的场方程对应不同的初始条件。 总之,为了确定泛定方程的解,就必须提供足够的初始条件和边界条件.边界条件与初始条件是控制方程有确定解的前提。边界条件是在求解区域的边界上所求解的变量或其导数随时间和地点的变化规律。对于任何问题,都需要给定边界条件。初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况,对于瞬态问题,必须给定初始条件,稳态问题,则不用给定。 对于边界条件与初始条件的处理,直接影响计算结果的精度。 在CFD模拟时,基本边界条件有: 1流动进口边界 包括速度进口边界,压力进口边界,质量进口边界(可压流动)。 在使用流动进口边界时,需要涉及到某些流动参数,如绝对压力,湍动能及耗散率,这些参数要做特殊考虑。关于参考压力,在流场数值计算中,压力总是按相对值表示的,实际求解的压力并不是绝对值,而是相对于进口压力而言的。 在有些情况下,可以通过设定进口压力为0,求解其他点的压力。还有时,为了减小数字截断误差,往往故意抬高或降低参考压力场的值,可使其余各处的计算压力场与整体数值计算的量级相吻合。 2流动出口边界 一般选在离几何扰动足够远的地方来施加。在这样的位置,流动是充分发展的,沿流动方向没有变化。该边界只有在进入计算域的流动是以进口边界条件给定时才使用,而且在只有一个出口的计算域中使用。

有限元边界条件和载荷

X边界条件和载荷 10.1边界条件 施加的力和/或者约束叫做边界条件。在HyperMesh中,边界条件存放在叫做load collectors的载荷集中。Load collectors可以通过在模型浏览器中点击右键来创建(Create > Load Collector)。 经常(尤其是刚开始)需要一个load collector来存放约束(也叫做spc-单点约束),另外一个用来存放力或者压力。记住,你可以把任何约束(比如节点约束自由度1和自由度123)放在一个load collector中。这个规则同样适用于力和压力,它们可以放在同一个load collector中而不管方向和大小。 下面是将力施加到结构的一些基本规则。 1.集中载荷(作用在一个点或节点上) 将力施加到单个节点上往往会出现不如人意的结果,特别是在查看此区域的应力时。通常集中载荷(比如施加到节点的点力)容易产生高的应力梯度。即使高应力是正确的(比如力施加在无限小的区域),你应该检查下这种载荷是不是合乎常理?换句话说,模型中的载荷代表了哪种真实加载的情形? 因此,力常常使用分布载荷施加,也就是说线载荷,面载荷更贴近于真实情况。 2.在线或边上的力 上图中,平板受到10N的力。力被平均分配到边的11个节点上。注意角上的力只作用在半个单元的边上。

上图是位移的云图。注意位于板的角上的红色“热点”。局部最大位移是由边界效应引起的(例如角上的力只作用在半个单元的边上),我们应该在板的边线上添加均匀载荷。 上述例子中,平板依然承受10N的力。但这次角上节点的受力减少为其他节点受力的一半大小。 上图显示了由plate_distributed.hm文件计算得到的平板位移的云图分布。位移分布更加均匀。 3.牵引力(或斜压力) 牵引力是作用在一块区域上任意方向而不仅仅是垂直于此区域的力。垂直于此区域的力称为压力。

周期边界条件

周期边界条件 aresaran (答网友问) (1)、究竟什么是"周期性边界条件"?如何去定义它的,为什么要引入这样一个定义。 周期边界条件源于这样的问题:宏观结构的信息不足以描述问题的细节,所以引入微观结构的信息来统计物质的宏观性质。周期边界条件广泛用于molecular dynamics & micromechanics. Fig1.细观力学的RVE 代表单元 尽管目前计算机的运算速度极大提高,但是仍然不能够用于进行大规模的宏微观联合计算。 因此引入了代表单元的概念,代表单元RVE 就如同是一个打开微观世界的一个窗口,看到的只是窗户里面的东西,我们假设整个微观世界是统计均匀的,因此无限量的复制了这个窗口,就可以得到所有微观信息。当然这个代表单元有要求,如上图,宏观结构尺寸远远尺寸,但是这个达标单元的尺寸又要能 足够多的包含微观颗粒的信息,有代表性,所以要求l L >>l A <<这是个一般性定义。 (2)、"周期性边界条件" 是不是只是在处理复合材料问题时才用,而且从众位大侠的讨论中似乎让我觉得这有点像"子结构"? Fig2. 2D or 3 D RVE

子结构和代表单元根本不在一个层次上,RVE 的建模与普通建模没什么区别,当然你想得到随机的微观结构,就需要用外部程序比如matlab 书写相应的inp 文件。 Fig3. Ref. Frederic Feyel. Multiscale elastoviscoplastic analysis of composite structures. Computational Materials Science,1999,16: 344~354 2FE 子结构模型适合多尺度计算。如图三,是一个发动机叶片,局部区域希望能够用细观微结构描述,其余结构希望是均匀材料。 这个问题的模型就可以将复合材料区域SiC/Ti 用子模型/子结构实现代表单元,子结构传递边界条件给代表单元, 实现微观和宏观的关联。 (3)、"周期性边条"是不是"旋转周期结构"里所需施加的边界条件? 对于复合材料层合壳体结构的旋转周期结构,相当于直角坐标周期结构的球坐标变换,物理意义等同。 (4)、为什么有些"轴对称单元"也在用这个? 因该是指对称性条件和周期性条件的关系,下面的例子会给出解释。 【1】周期边界条件的推导实例: ij 是边界上施加的的宏观应变条件 Displacement BC. j ij i i l x u y u ε+=)()( Traction BC. )()()()(x n x y n y j ij j ij σσ?=

定解条件和定解问题

定解条件和定解问题 含有未知函数的偏导数的方程叫偏微分方程,常微分方程可以看成是特殊的偏微分方程。方程的分数是1的称为方程式,个数多于1的叫做方程组。方程(组)中出现的未知函数的最高阶偏导数的阶数称为方程(组)的阶数。如果方程(组)中的项关于未知函数及其各阶偏导数的整体来讲是线性的,就称方程(组)为线性的,否则就称为非线性的。非线性又分为半线性、拟线性和完全非线性。 一、定解条件 给定一个常微分方程,有通解和特解的概念。通解只要求满足方程,即满足某种物理定律,而不能完全确定一个物理状态。特解除了要求满足方程还要满足给定的外加(特殊)条件。对偏微分方程也是如此,换句话说,只有偏微分方程还不足以确定一个物理量随空间和时间的变化规律,因为在特定情况下这个物理量还与它的初始状态和它在边界受到的约束有关。描述初始时刻的物理状态和边界的约束情况,在数学上分别称为初始条件(或初值条件)和边界条件(或边值条件),他们统称为定解条件。 初始条件:能够用来说明某一具体物理现象初始状态的条件,即描述物理过程初始状态的数学条件。 边界条件:能够用来说明某一具体物理现象边界上的约束情况的条件,即描述物理过程边界状态的数学条件。 定解条件:初始条件和边界条件的统称。 非稳态问题:定解条件包括初始条件和边界条件。

稳态问题:定解条件为边界条件。 1、弦振动方程 ( 2(,),0,0tt xx u a u f x t x l t -=<<>) 初始条件是指初始时刻(0t =)弦的位移和速度。若以()x ?, ()x ψ分别表示弦上任意点x 的初始位移和初始速度,则初始条件为: 边界条件是指弦在两端点的约束情况,一般有三种类型。 (1)第一类边界条件(狄利克雷(Dirichlet )边界条件):已知端点()x a a o a l ===或处弦的位移是()a g t ,则边界条件为: (0,)(0,)u t g t = 或 (,)(,)u l t g l t = 当0()0()0l g t g t ≡≡或时,表示在该点处弦是固定的。 (2)第二类边界条件(诺伊曼(Neumann )边界条件):已知端点0x x l ==或处弦所受的垂直于弦线的外力0()g t 或()l g t ,则边界条件为: 0(0,)()x Tu t g t -= 或 (,)()x l Tu l x g t = 当00()0l g g t ≡≡或时,表示弦在端点0x x l ==或处自由滑动。 (3)第三类边界条件(混合边界条件或罗宾(Robin )边界条件:已知端点处弦的位移和所受的垂直于弦线的外力的和: 000(0,)(0,)g (t),0,x Tu t k u t k -+=> 或 (,)(,)(),0x l l l Tu l t k u l t g t k +=>, (,0)(),0(,0)(), t u x x x l u x x ?ψ=?<

fluent边界条件(一)

边界条件 定义边界条件概述 边界条件包括流动变量和热变量在边界处的值。它是FLUENT分析得很关键的一部分,设定边界条件必须小心谨慎。 边界条件的分类:进出口边界条件:压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件、质量出口、通风口、排气扇;壁面、repeating, and pole boundaries:壁面,对称,周期,轴;内部单元区域:流体、固体(多孔是一种流动区域类型) ;内部表面边界:风扇、散热器、多孔跳跃、壁面、内部。(内部表面边界条件定义在单元表面,这意味着它们没有有限厚度,并提供了流场性质的每一步的变化。这些边界条件用来补充描述排气扇、细孔薄膜以及散热器的物理模型。内部表面区域的内部类型不需要你输入任何东西。) 下面一节将详细介绍上面所叙述边界条件,并详细介绍了它们的设定方法以及设定的具体合适条件。周期性边界条件在本章中介绍,模拟完全发展的周期性流动将在周期性流动和热传导一章中介绍。 使用边界条件面板 边界条件(Figure 1)对于特定边界允许你改变边界条件区域类型,并且打开其他的面板以设定每一区域的边界条件参数 菜单:Define/Boundary Conditions... Figure 1: 边界条件面板 改变边界区域类型 设定任何边界条件之前,必须检查所有边界区域的区域类型,如有必要就作适当的修改。比方说:如果你的网格是压力入口,但是你想要使用速度入口,你就要把压力入口改为速度入口之后再设定。 改变类型的步骤如下:: 1.在区域下拉列表中选定所要修改的区域

2.在类型列表中选择正确的区域类型 3.当问题提示菜单出现时,点击确认 确认改变之后,区域类型将会改变,名字也将自动改变(如果初始名字时缺省的请参阅边界条件区域名字一节),设定区域边界条件的面板也将自动打开。 !注意:这个方法不能用于改变周期性类型,因为该边界类型已经存在了附加限制。创建边界条件一节解释了如何创建和分开周期性区域。需要注意的是,只能在图一中每一个类别中改变边界类型(注意:双边区域表面是分离的不同单元区域.) Figure 1: 区域类型的分类列表 设定边界条件 在FLUENT中,边界条件和区域有关而与个别表面或者单元无关。如果要结合具有相同边界条件的两个或更多区域请参阅合并区域一节。 设定每一特定区域的边界条件,请遵循下面的步骤: 1.在边界条件区域的下拉列表中选择区域。 2. 点击Set...按钮。或者,1.在区域下拉列表中选择区域。 2.在类型列表中点击所要选择的类型。或者在区域列表中双击所需区域.,选择边界条件区域将会打开,并且你可以指定适当的边界条件 在图像显示方面选择边界区域 在边界条件中不论你合适需要选择区域,你都能用鼠标在图形窗口选择适当的区域。如果你是第一次设定问题这一功能尤其有用,如果你有两个或者更多的具有相同类型的区域而且你想要确定区域的标号(也就是画出哪一区域是哪个)这一功能也很有用。要使用该功能请按下述步骤做: 1.用网格显示面板显示网格。 2.用鼠标指针(默认是鼠标右键——参阅控制鼠标键函数以改变鼠标键的功能)在图形窗口中点击边界区域。在图形显示中选择的区域将会自动被选入在边界条件面板中的区域列表中,它的名字和编号也会自动在控制窗口中显示改变边界条件名字 每一边界的名字是它的类型加标号数(比如pressure-inlet-7)。在某些情况下你可能想要对边界区域分配更多的描述名。如果你有两个压力入口区域,比方说,你可能想重名名它们

晶格振动、金属电子论、能带理论三个地方都用到了周期性边界条件

1.固体物理教材在晶格振动、金属电子论、能带理论三个地方都用到了周期性 边界条件,试比较其异同并阐述你的理解。 周期性边界条件是边界条件的一种,反映的是如何利用边界条件替代所选部分(系统)受到周边(环境)的影响。可以看作是如果去掉周边环境,保持该系统不变应该附加的条件,也可以看作是由部分的性质来推广表达全局的性质。 周期性边界条件的引入有两个目的:在粒子的运动过程中,若有一个或几个粒子跑出模型,则必有一个或几个粒子从相反的界面回到模型中,从而保证该模拟系统的粒子数恒定;计算原子间作用力的时候采取最近镜像方法,这样模型中处于边界处的原子受力就比较全面,从而消除了边界效应。这种方法在计算机分子动力学模拟中使用非常广泛。 由此,在讨论晶格振动、金属电子论、能带理论的周期性边界条件时只是在不同的范围中周期性边界条件具体的定义、应用以及意义。 晶格振动的周期性边界条件:由N个原子组成一个模型——原子数目有限,但各原子完全等价。第j个原子的运动与第 mN+j个原子的运动情况完全一样。对于原子的自由运动,边界上的原子与其它原子一样,无时无刻不在运动,对于有N个原子原子链,硬性设定u1=0,uN=0的边界条件是不符合事实的。其实不论什么边界条件都与事实不符合,但为了求近似解,必须选取一个边界条件,晶格振动谱的实验测定是对晶格振动理论的最有力验证,周期性边界条件是晶格振动理论的前提条件。 金属电子论的周期性边界条:.金属中自由电子气应该服从量子力学规律,在保留独立电子近似和自由电子近似基础上应通过求解薛定愕方程给出电子本征态和本征能量,从而来解释金属性质。我们把自由电子气等效为在温度 T=0K,V =L3的立方体内运动的 N个自由电子。独立电子近似使我们可以把 N个电子问题转换为单电子问题处理。要计算一系列想关函数都与波矢 k有关。波矢 k 的取值要由边界条件决定,边界条件的选取既要反映出电子是在有限体积中运动的特点,又要在数学上便于操作,因此,类似于晶格振动是的情况,周期性边界条件(Born-Karman边界条件)是人们通常采用的最适合的方法。 能带理论的周期性边界条件:能带论的基本出发点是认为固体中的电子不再是完全被束缚在某个原子周围,而是可以在整个固体中运动的,称之为共有化电子。但电子在运动过程中并也不像自由电子那样,完全不受任何力的作用,电子在运动过程中受到晶格原子势场和其它电子的相互作用。能带理论是基于三个基本(近似)假设:1)Born-Oppenheimer 绝热近似:离子的波函数与电子的位置及状态无关:多粒子问题→多电子问题2)Hatree-Fock平均场近似:忽略电子与电子间的相互作用,用平均场代替电子与电子间的相互作用:多电子问题→单电子问题。3)周期场近似:单电子问题→单电子在周期场中运动问题。由于这三个基本假设,每个电子都处在完全相同的严格周期性势场中运动,因此每个电子的运动都可以单独考虑。在计算电子运动的薛定谔方程时,由于势场的周期性反映了晶格的平移对称性,可定义一个平移算符,为了确定平移算符的本征值,引入周期性边界条件。

(完整版)fluent边界条件设置

边界条件设置问题 1、速度入口边界条件(velocity-inlet):给出进口速度及需要计算的所有标量值。该边界条件适用于不可压缩流动问题。 Momentum 动量?thermal 温度radiation 辐射species 种类 DPM DPM模型(可用于模拟颗粒轨迹)multipahse 多项流 UDS(User define scalar 是使用fluent求解额外变量的方法) Velocity specification method 速度规范方法:magnitude,normal to boundary 速度大小,速度垂直于边界;magnitude and direction 大小和方向;components 速度组成?Reference frame 参考系:absolute绝对的;Relative to adjacent cell zone 相对于邻近的单元区 Velocity magnitude 速度的大小 Turbulence 湍流 Specification method 规范方法 k and epsilon K-E方程:1 Turbulent kinetic energy湍流动能;2 turbulent dissipation rate 湍流耗散率 Intensity and length scale 强度和尺寸:1湍流强度 2 湍流尺度=0.07L(L为水力半径)intensity and viscosity rate强度和粘度率:1湍流强度2湍流年度率 intensity and hydraulic diameter强度与水力直径:1湍流强度;2水力直径

边界元与有限元

边界元与有限元 边界元法boundary element method 定义:将力学中的微分方程的定解问题化为边界积分方程的定解问题,再通过边界的离散化与待定函数的分片插值求解的数值方法。 所属学科:水利科技(一级学科) ;工程力学、工程结构、建筑材料(二级学科) ;工程力学(水利)(三级学科) 边界元法(boundary element method)是一种继有限元法之后发展起来的一种新数值方法,与有限元法在连续体域内划分单元的基本思想不同,边界元法是只在定义域的边界上划分单元,用满足控制方程的函数去逼近边界条件。所以边界元法与有限元相比,具有单元个数少,数据准备简单等优点.但用边界元法解非线性问题时,遇到同非线性项相对应的区域积分,这种积分在奇异点附近有强烈的奇异性,使求解遇到困难。 简介 边界元法是在有限元法之后发展起来的一种较精确有效的工程数值分析方法。又称边界积分方程-边界元法。它以定义在边界上的边界积分方程为控制方程,通过对边界分元插值离散,化为代数方程组求解。它与基于偏微分方程的区域解法相比,由于降低了问题的维数,而显著降低了自由度数,边界的离散也比区域的离散方便得多,可用较简单的单元准确地模拟边界形状,最终得到阶数较低的线性代数方程组。又由于它利用微分算子的解析的基本解作为边界积分方程的核函数,而具有解析与数值相结合的特点,通常具有较高的精度。特别是对于边界变量变化梯度较大的问题,如应力集中问题,或边界变量出现奇异性的裂纹问题,边界元法被公认为比有限元法更加精确高效。由于边界元法所利用的微分算子基本解能自动满足无限远处的条件,因而边界元法特别便于处理无限域以及半无限域问题。边界元法的主要缺点是它的应用范围以存在相应微分算子的基本解为前提,对于非均匀介质等问题难以应用,故其适用范围远不如有限元法广泛,

ABAQUS旋转周期对称边界条件的设置

ABAQUS旋转周期对称边界条件的设置 旋转周期对称设置包括:旋转周期对称设置,外加主面上的对称面约束,两者一起构成旋转对称的边界条件。下面所述的两种方法是仅针对旋转周期对称的设置。 两种方法: 1)修改inp文件: 找到*End Assembly,将之替换为 *TIE,CYCLIC SYMMETRY,NAME=TIE-CYCLIC Surf-Cylic-SLAVE,Surf-Cylic-MASTER ** *End Assembly ** *CYCLIC SYMMETRY MODEL,N=60 0,0,0,0,0,1 --------------------------- 上面设置中包括:主面的设置,从面的设置,模型周期的数目,以及旋转轴。因此需要建立这两个面的集合:Surf-Cylic-MASTER,Surf-Cylic-SLAVE。N=60表示有60个。0,0,0为旋转轴的起点,0,0,1为旋转轴的终点。 2)直接在前处理cae中设置 首先,建立主面和从面的集合,便于选取; 其次,为旋转轴的起点和终点建立参考点(RP),旋转轴一定要设在整个模型的旋转中心上;参考点可通过输入坐标的方式建立。注意:其他方式建立点都不可行,以下详述。 最后,输入周期的数目,本模型为整体模型的多少分之一,即输入倒数即可。 以上步骤参见下图。 【旋转轴起点和终点的建立】 1)除参考点以外其他的建点的方式不行,比如建立datum point,无法在viewport中直接选中,同样建立集合时也选不中datum point。 2)使用attachment point建立的点虽然可以直接在viewport中选中,建立集合时也可选中,但无法写入inp文件,当write inp 文件时就造成cae崩溃直接退出软件! 总之,旋转轴的设置,直接在前处理cae界面中设置,不如直接在inp文件中修改方便!因为修改inp旋转轴只要直接给定起点和终点坐标就OK,省去先建立RP点的步骤。 【主面上设置对称面】 在边界条件中选对称面设置即可。先要建立一个柱坐标系为好。将柱坐标系的Z轴建在旋转中心上,R轴在模型两对称侧面的平分线上,T轴即自动建好为切线方向。对称边界设置时,选取之前建立的主面,方向为U2=UR1=UR3=0,此即为T轴为对称面的法线方向。 【补充说明】 对于一个具体的部件,除上述约束外,根据实际情况还需加上其他约束条件避免存在任何刚体位移的出现。如Z向(轴向)上避免刚体位移,径向上避免刚体位移。 下文算例中的详情看文末的总结。

有限元法与有限差分法的主要区别

有限元法与有限差分法的主要区别 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑,可分为中心格式和逆风格式。考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的平方误差最小;在配置法中,先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程,即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身,还要求它的导数值在插值点取已知值,称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标,有对称和不对称等。常采用的无因次坐标是一种局部坐标系,它的定义取决于单元的几何形状,一维看作长度比,二维看作面积比,三维看作体积比。在二维有限元中,三角形单元应用的最早,近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元,常采用的插值函数为有La grange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。对于有限元方法,其基本思路和解题步骤可归纳为(1)建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的出发点。(2)区域单元剖分,根据求解区域的形状及实际问题的物理特点,将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作,这部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标,同时还需要列出自然边界和本质边界的节点序号和相应的边界值。(3)确定单元基函数,根据单元中节点数目及对近似解精度的要求,选择满足一定插

相关文档
最新文档