有限元、边界元、有限差分法的区别

有限元、边界元、有限差分法的区别
有限元、边界元、有限差分法的区别

有限元法、边界元法、有限差分法的区别和各自的优点

请问:有限元法、边界元法、有限差分法等方法有哪些区别和各自的优点?尤其是在声学方面。

谢谢!

网格的跑分上不同,差分要求模型规则,有限元可以是任意不规则模型,

FEM: irregular grid-> easy to describe complex shape, hard in mesh generation \.a4hj

FDM: regular mesh -> easy in grid generation, hard to describe complex shape=> less accurate than FEM

BEM: irregular mesh in boundary -> mesh generation much easier than that of FEM. need much less computation resource than the above two. BUT need basic solution (Green function) at the boundary.

对于这个基础问题一定要搞清楚,不然有限元就无从谈起。

有限元法的优点是适应性强,自由边界条件自动满足,但是不适合计算大尺度,对于透射边界需单独处理,单

元太多的模型,计算速度慢i7g c1T `w5v

边界元法的优点是域内二维问题化成了边界一维问题来处理,自动满足透射边界,但是构造G函数非常麻烦有限差分法适合大尺度(如地震波),方法简单,计算速度快,但是边界处理太麻烦.

:) :( :D :'(

[quote]原帖由[i]jonewore[/i] 于2007-10-1 20:31 发表

[url=https://www.360docs.net/doc/c517984200.html,/forum/redirect.php?goto=findpost&pid=1152036&ptid=7785 04][img]https://www.360docs.net/doc/c517984200.html,/forum/images/common/back.gif[/img][/url]

有限元法的优点是适应性强,自由边界条件自动满足,但是不适合计算大尺度,对于透射边界需单独处理,单元太多的模型,计算速度慢

边界元法的优点是域内二维问题化成了边界一维问题来处理,自动满足透射边界,但是构造 ... [/quote] -O

你说自动满足透射边界是什么意思?是说边界的反射波可以完全吸收吗(不用再使用人工边界?)?能不能详细说一下呢。。。

求解微分方程时,

有限元从微分方程的等效积分形式出发,个人理解其实就是把微分方程转化成一个泛函问题8k%h K p+M"J

有限差分是直接在网格结点上采用差分方程近似微分方程

有限元法适用于拉格朗日坐标系下建立的微分方程

有限差分法适用于欧拉坐标系下建立的微分方程

不过,据个人经验,等效积分方法与差分法混合使用威力也很强大

有限元法适用于拉格朗日坐标系下建立的微分方程,也就意味着有限元法适应性很强

此外,对于求解微分方程,个人认为应该顺便浏览一下“坐标函数(或者函数空间)”的概念,函数空间意味着把原函数在一个函数形成的空间中分解,即代数几何化,对于对空间(或几何)有一定感觉而代数能力稍弱的人来说,应该会获益匪浅。当然,这仅是个人观点

离散的对象不同,差分法直接离散方程,有限元将物理问题转化成求势能泛函的最小值,对泛函的积分区域进行离散,边界元法将方程的解表示成边界积分的形式,然后离散边界区域。

基本原理相同,但应用范围不一样。

有限元法、边界元法、有限差分法等方法有哪些区别和各自的优点?尤其是在声学方面。

答:首先三者的算法不同。从应用角度看,有限元法和有限差分法适合用于有限区域,而边界元法适用于无限域或很大区域。

从划分网格方面看,有限差分法对网格质量要求高于有限元法。

但是在求解非线性程度很高的问题时,有时用有限元法迭代求解如果失败则得不到任何结果。而有限差分法只要网格够细密划分单元形状质量够好,即使得不到最终结果,还是可以得到部分结果。m kH w$o(~

声学问题好像多用有限差分方法以及边界元法。

UG有限元分析

UG有限元分析 第1章有限元分析方法及NX Nastran的由来 1.1 有限元分析方法介绍 计算机软硬件技术的迅猛发展,给工程分析、科学研究以至人类社会带来急剧的革命性变化,数值模拟即为这一技术革命在工程分析、设计和科学研究中的具体表现。数值模拟技术通过汲取当今计算数学、力学、计算机图形学和计算机硬件发展的最新成果,根据不同行业的需求,不断扩充、更新和完善。 1.1.1 有限单元法的形成 近三十年来,计算机计算能力的飞速提高和数值计算技术的长足进步,诞生了商业化的有限元数值分析软件,并发展成为一门专门的学科——计算机辅助工程CAE(Computer Aided Engineering)。这些商品化的CAE软件具有越来越人性化的操作界面和易用性,使得这一工具的使用者由学校或研究所的专业人员逐步扩展到企业的产品设计人员或分析人员,CAE在各个工业领域的应用也得到不断普及并逐步向纵深发展,CAE工程仿真在工业设计中的作用变得日益重要。许多行业中已经将CAE分析方法和计算要求设置在产品研发流程中,作为产品上市前必不可少的环节。CAE仿真在产品开发、研制与设计及科学研究中已显示出明显的优越性: ?CAE仿真可有效缩短新产品的开发研究周期。 ?虚拟样机的引入减少了实物样机的试验次数。 ?大幅度地降低产品研发成本。 ?在精确的分析结果指导下制造出高质量的产品。 ?能够快速对设计变更作出反应。 ?能充分和CAD模型相结合并对不同类型的问题进行分析。 ?能够精确预测出产品的性能。 ?增加产品和工程的可靠性。 ?采用优化设计,降低材料的消耗或成本。 ?在产品制造或工程施工前预先发现潜在的问题。 ?模拟各种试验方案,减少试验时间和经费。

有限差分法、有限单元和有限体积法简介

有限差分法、有限单元法和有限体积法的简介 1.有限差分方法 有限差分方法(Finite Difference Method,FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。 对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑,可分为中心格式和逆风格式。考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。 2.有限元方法 有限元方法(Finite Element Method,FEM)的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。 有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。 在数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的

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

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

排气系统的声学模型

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

电磁仿真算中的有限元法

1电磁仿真算法中的有限元法 1.1常规的电磁计算方法简介 从上世纪50年代以来,伴随着计算机技术的进步,电磁仿真算法也蓬勃发展起来,这其中主要包括:单矩法、矩量法和有限元法等属于频域技术的算法; 传输线矩阵法、时域积分方程法以及时域有限差分法等属于时域技术的算法。除了这些以外, 还有属于高频技术的集合衍射理论等。本文根据国内外计算电磁学的发展状况,对日常生活中比较常用的电磁计算方法做了介绍,并对有限元法做了重点说明。 ⑴矩量法 矩量法属于电磁场的数值计算方法中频域技术的一种, 它的基本原理是利用把待解的微积分方程转化成的算子方程, 然后将由一组线性组合表示的待求函数代入第一步中的算子方程, 然后将算子方程转化成矩阵方程, 最后再通过计算机进行大量的数值计算从而得到数值结果。该方法在求解非均勻和不规则形状对象时,面很广,但会生成病态矩阵,所以会在一定程度上受到限制。矩量法的特点就是适用于求解微积分方程, 并且求解方法统一简单。但缺点就是会占用大量计算机内存,影响计算速度。 (2)单矩法 单矩法是一种解析方法和数值方法相结合的混合数值算法法,该方法的关键在于,如何合理的选择一个球面最小的半径,使得能够将分析对象的结构全部包含在内,以便将内外场进行隔离。外边的散射场单独使用其他函数表示,而包围的内部区域使用有限元法亥姆赫兹(Helmholtz)方程。此方法对于计算复杂形体乃至复杂埋入体内的电磁散射是种极为有效的手段。 (3)时域有限差分法 时域有限差分法(FDTD)近几年来越来越受到各方的重视, 因为一方面它处理庞大的电磁福射系统方面和复杂结构的散射体时很突出,另外一方面则在于它不是传统的频域算法, 它是种时域算法, 直接依靠时间变量求解麦克斯韦方程组,可以在有限的时间和体积内对场进行数据抽样, 这样同时也能够保证介质边界

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

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

有限元边界条件和载荷

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.牵引力(或斜压力) 牵引力是作用在一块区域上任意方向而不仅仅是垂直于此区域的力。垂直于此区域的力称为压力。

有限差分和有限体积的 有限元等

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

有限差分、有限元区别

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

边界元与有限元

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

内燃机零部件有限元计算中边界条件处理的研究

内燃机零部件有限元计算中边界条件处理的研究 * 孙 军 汪景峰 桂长林 (合肥工业大学机械与汽车工程学院 合肥 230009) 摘 要:有限元方法已经成为内燃机零部件应力和变形计算的主要手段,但是目前在内燃机零部件有限元分析中采用的边界条件是否合理,有无必要采用更符合实际的边界条件?本文以曲轴为例,模拟实际 状况,采用不同的边界条件进行了有限元计算。计算结果表明,边界条件处理对曲轴有限元分析结果影响很大。因此,为了提高内燃机零部件有限元计算结果的精度,非常有必要根据实际情况确定边界条件。 关键词:边界条件 有限元 内燃机中图分类号:TK412.4 文献标识码:A 文章编号:1671-0630(2005)03-0006-03 Study on Boundary Condition in Finite Ele ment Calculation for Parts of Internal Co mbustion Engi ne Sun Jun ,W ang Jingfeng ,Gui Changlin H efeiUn i v ersity of Techno l o gy (H efei 230009) Abst ract :The fi n ite ele m ent m et h od has beco m e the m a i n m eans to calcu late t h e stress and de f o r m ation o f parts for inter na l co m bustion engine .Bu,t whether the boundary conditi o ns used i n FE ana l y sis on parts o f i n -ter nal co m busti o n eng ine are reasonable ?Is it necessary to use the boundary condition ,wh ich ism ore adapta -b le to the facts ?As an exa m p le ,the crankshaft is ca lculated by FE usi n g d ifferent boundary conditi o ns that si m ulate factual conditi o ns .The resu lts sho w t h at the boundary conditi o ns have i m portant effects on the results of FE analysis o f crankshaf.t Therefo re ,it is necessary to choose boundary cond itions acco r d i n g to factua l con -d iti o n i n o r der to i m prove the prec isi o n of calcu l a ti n g resu lts for parts o f i n ternal co m bustion eng i n e .K eyw ords :Boundary conditi o n ,F i n ite ele m en,t I C eng i n e 前言 随着有限元计算技术的进步,有限元方法目前已 经成为内燃机零部件应力和变形计算的主要手段。内燃机零部件的有限元分析,类似于其他问题的有限元分析,边界条件的处理是否合理直接影响计算结果的精确性。本文以曲轴为例,分析目前采用的边界条件是否合理,有无必要采用更符合实际的边界条件。 目前在曲轴有限元计算中,载荷边界条件的处理(重点是作用在轴颈表面的力处理)基本采用的是定 型模式,其假设作用在轴颈上的载荷(其与曲轴轴承油膜压力对应)为分布载荷,沿轴线方向均布或呈抛物线分布,沿圆周方向呈余弦分布 [1~4] 。这种处理方 法简单易行,但其属于较理想的状况,因为实际曲轴轴承的油膜压力分布规律复杂,且随时间变化。沿轴向抛物线型的油膜压力分布规律仅适合于无限短且轴颈轴线与轴承孔中心线平行的滑动轴承,实际的曲轴轴承为有限长轴承,且由于受到诸多因素的影响,如载荷作用下轴的变形、轴承的制造与装配误差和轴的热变形 * 基金项目:国家自然科学基金资助项目(50175023) 作者简介:孙军(1960-),男,硕士,研究方向,内燃机现代设计理论与方法。 第34卷 第3期2005年6月小型内燃机与摩托车 S MALL I N TERNAL COM B UST I O N ENG I N E AND MOTORCYCLE Vo.l 34No .3 June .2005

有限差分,有限元,有限体积等的区别介绍

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

IDESA有限元分析_第6篇第26章 基于几何施加边界条件

第26章MasterFEM 教程:定义边界条件 前面的教程简单介绍了仿真分析的流程。本篇将介绍更多高级定义边界条件的内容(载荷和约束)。 用户将学会: ?创建约束和约束集。 ?创建载荷和载荷集。 ?创建边界条件集。 ?解算定义以上边界条件的模型。 ?创建均布载荷。 ?解算定义以上边界条件的模型。 ?比较不同工况下的结果。 开始前必备知识: 熟悉MasterFEM界面和创建零件。 熟悉在模型文件中管理零件。 熟悉拉伸特征和旋转特征的布尔运算。 熟悉仿真分析流程。 熟悉自由网格划分。 设置1/3 如果还没有运行一个新的模型文件,创建一个新文件并命名。 ·1·

·2· File Open 打开模型文件菜单 确信用户是在以下工作状态和任务当中 : 设置工作单位为毫米(mm) Options Units 设置2/3 工作内容:按照以下尺寸草绘封闭形状的图形。 提示 : 为什么:这个零件代表了典型机构连杆的应力集中部位。

工作内容: 命名零件 提示: 命名菜单 设置3/3 工作内容:创建一个和零件关联的有限元模型(FEM1)。 提示 保存模型文件。 File Save 警告! 如果软件提示用户保存模型文件,用户应选择:No 记住:只有教程中提示保存模型文件,而不是软件提示保存的时候,用户才可以执行保存文件操作。 为什么: 在上一次保存以后的错误操作不能撤销恢复,用户可以选择重新打开文件,恢复到上一次保存时的状态。 提示: ·3·

重新打开模型文件的快捷键:按Control-Z。 创建约束和约束集1/3 工作内容:全约束以下高亮表面。 怎样做: 表面上定义约束的菜单 OK 创建约束和约束集2/3 注意事项: 会产生约束符号。 在几何边缘、表面、顶点的约束用不同的颜色和符号表示。 ·4·

有限差分,有限元,有限体积等离散方法的区别介绍

https://www.360docs.net/doc/c517984200.html,/s/blog_501a61220100f9rs.html 有限差分,有限元,有限体积等离散方法的区别介绍 (2009-10-25 22:07:18) 转载 以下介绍是本人从网络上搜集的,供计算数学虫子参考。也许小木虫论坛有,我没搜索到。欢迎大家补充内容。 转自https://www.360docs.net/doc/c517984200.html,/bbs/viewthread.php?tid=1618917&pid=16196206&page=1#pid16196206 1 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑,可分为中心格式和逆风格式。考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。 2 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。 在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。 根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分

有限差分及有限单元法的区别

1 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑,可分为中心格式和逆风格式。考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。 2 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。 在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。 根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。从

对有限元法 有限差分法 边界元法和模拟电荷法的粗略总结

对有限元法、有限差分法、边界元法和模拟电荷法的粗略总结: 有限元法(finite element method):将连续的求解域离散为一组单元的组合体,用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数,近似函数通常由未知场函数及其导数在单元各节点的数值插值函数来表达。从而使一个连续的无限自由度问题变成离散的有限自由度问题。缺点是有限元必须同时对所有域内节点和边界节点联立求解,待求未知数多,要求解的方程规模大,导致输入数据多,计算的准备工作量大。 有限差分法(finite difference method):直接从微分方程出发,将求解区域划分为网格,近似地用差分、差商代替微分、微商,于是无限度的问题化成有限自由度的问题。这种方法在解决规则边界的问题时极为方便,但是正是由于这种限制而增加了它的局限性,即对于非规则边界的问题适用性较差。 边界元法(boundary element method):边界元法是在有限元法之后发展起来的一种较精确有效的工程数值分析方法。它以定义在边界上的边界积分方程为控制方程,通过对边界分元插值离散,化为代数方程组求解。它与基于偏微分方程的区域解法相比,由于降低了问题的维数,而显著降低了自由度数,边界的离散也比区域的离散方便得多,可用较简单的单元准确地模拟边界形状,最终得到阶数较低的线性代数方程组。又由于它利用微分算子的解析的基本解作为边界积分方程的核函数,而具有解析与数值相结合的特点,通常具有较高的精度。特别是对于边界变量变化梯度较大的问题,如应力集中问题,或边界变量出现奇异性的裂纹问题,边界元法被公认为比有限元法更加精确高效。由于边界元法所利用的微分算子基本解能自动满足无限远处的条件,因而边界元法特别便于处理无限域以及半无限域问题。边界元法的主要缺点是它的应用范围以存在相应微分算子的基本解为前提,对于非均匀介质等问题难以应用,故其适用范围远不如有限元法广泛,而且通常由它建立的求解代数方程组的系数阵是非对称满阵,对解题规模产生较大限制。对一般的非线性问题,由于在方程中会出现域内积分项,从而部分抵消了边界元法只要离散边界的优点。 模拟电荷法(charge simulation method):在实际工程计算中,电极表面上连续分布的束缚电荷的分布情况是未知的,不能直接由给定的边界条件解出。如果在计算场域之外设置n个被称为模拟电荷的离散电荷来等效代替这些待求的连续电荷分布,则根据等值替代前后条件不变的前提条件,即可求得各模拟电荷的量值,从而使场域内任意一点的电位与场强便可由各模拟电荷所产生的场量叠加而获得,以此作为原场的逼近解。相比较于有限元法和有限差分法,模拟电荷法的优点是无需封边、使计算问题的维数降低一维、能直接求解出场域内的任意点的场强、计算精度高。

计算流体力学中有限差分法、有限体积法和有限元法的区别

有限元法,有限差分法和有限体积法的区别 1. FDM 1.1 概念 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。 1.2 差分格式 (1)从格式的精度来划分,有一阶格式、二阶格式和高阶格式。 (2)从差分的空间形式来考虑,可分为中心格式和逆风格式。 (3)考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 1.3 构造差分的方法 构造差分的方法有多种形式,目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。 2. FEM 2.1 概述 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。 2.2 原理 有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。 根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。(1)从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法;(2)从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格;

相关文档
最新文档