搜档网
当前位置:搜档网 › 国内有限元法的发展之路

国内有限元法的发展之路

国内有限元法的发展之路
国内有限元法的发展之路

国内有限元法的发展之路

我国的力学工作者为有限元方法的初期发展做出了许多贡献,其中比较著名的有:陈伯屏(结构矩阵方法),钱令希(余能原理),钱伟长(广义变分原理),胡海昌(广义变分原理),冯康(有限单元法理论)。遗憾的是由于当时环境所致,我国有限元方法的研究工作受到阻碍,有限元理论的发展也逐渐与国外拉开了距离。

20世纪60年代初期,我国的老一辈计算科学家较早地将计算机应用于土木、建筑和机械工程领域。当时黄玉珊教授就提出了“小展弦比机翼薄壁结构的直接设计法”和“力法-应力设计法”;而在70年代初期,钱令希教授提出了“结构力学中的最优化设计理论与方法的近代发展”。这些理论和方法都为国内的有限元技术指明了方向。

1964年初崔俊芝院士研制出国内第一个平面问题通用有限元程序,解决了刘家峡大坝的复杂应力分析问题。20世纪60年代到70年代,国内的有限元方法及有限元软件诞生之后,曾计算过数十个大型工程,应用于水利、电力、机械、航空、建筑等多个领域。

20世纪70年代中期,大连理工大学研制出了JEFIX有限元软件,航空工业部研制了HAJIF系列程序。80年代中期,北京大学的袁明武教授通过对国外SAP软件的移植和重大改造,研制出了SAP-84;北京农业大学的李明瑞教授研发了FEM软件;建筑科学研究院在国家“六五”攻关项目支持下,研制完成了“BDP-建筑工程设计软件包”;中国科学院开发了FEPS、SEFEM;航空工业总公司飞机结构多约束优化设计系统YIDOYU等一批自主程序。

上世纪90年代以来,大批国外CAE软件涌入国内市场,遍及国内的各个领域,国外的专家则深入到大学、院所、企业与工厂,展示他们的CAE技术、系统功能及使用技巧,因此使得国内自主研发CAE软件受到强烈打压。同时,有关管理部门在对直接为先进装备制造业服务的CAE软件核心技术的认识上产生了偏差:CAE既不属于基础科学,又不属于科技攻关,故而失去了必要的支持,使其发展举步维艰,以至于在上世纪的最后十几年国内CAE自主创新的步伐已经非常缓慢,也逐渐的拉开了与国外CAE软件的距离。

进入21世纪后,虽然国外CAE软件占据市场主流的现状短时间内已经无法撼动,但国内自主知识产权CAE软件逐渐市场化,获得了一定的发展:元计算科技发展有限公司推出的FELAC、郑州机械研究所推出的紫瑞CAE、湖南大学与吉林大学开发了针对汽车结构的KMAS分析系统;清华大学、上海交大在注塑成型CAE领域也推出了相应的分析软件。

虽然国内CAE自主研发之路历经艰辛,但是广大专家学者用锲而不舍的战斗精神顽强地生存下来。尤其是在近几年,数字化产品设计的概念逐渐深入人心,国内高校技术研究和应用水平不断提高,有限元技术已经为广大企业所认可。随着国家对发展自主CAE平台已经愈发重视,国内CAE的研究已经逐渐走出低迷状态,获得了一定的发展,而且值得注意的是,有限元技术不再仅仅停留在高校中,而是更多的走向了企业。同时,更多使用方便、操作简单的专用分析软件也得到了广泛应用。

《结构分析中的有限元法》2015-有限元习题-参考答案

本科有限元习题参考答案

2015年3月10日作业 1、简述力学课程中介绍的各种力学模型的简化条件、基本假设和适用范围(包括有拉压杆模型、弯曲梁模型、平面应力和平面应变模型、轴对称模型、板模型、壳模型等) 2、给出弹性力学问题中平衡方程、几何方程、物理方程的表达式及其意义。 (1)平衡方程:

zy yz xz zx yx xy z yz xz z y xy zy y x zx yx x f y x z f x z y f z y x ττττττττσττσττσ====+??+??+??=+??+??+??=+??+??+??,000, 物理意义:应力分量与体力分量之间的关系。 (2)几何方程: z u x w y w z v x v y u z w y v x u zx yz xy z y x ??+??=??+??=??+??=??=??=??=γγγεεε,,,, 物理意义:应变分量与位移分量之间的关系。 (3)物理方程: [] [] [] zx zx yz yz xy xy y x z z z x y y z y x x G G G E E E τγτγτγσσμσεσσμσεσσμσε1,1,1) (1 ) (1 )(1 ===+-=+-=+-= 物理意义:应变分量与应力分量之间的关系。 3、简述最小势能原理的主要内容和主要公式。 根据虚功原理得到:??=-Γ T Ω T T 0Td Γδu d Ω)F δu -σδε(,由 )(21εδσεδδεU T T =?? ? ??=则0)21((=Γ-Ω-=∏??ΩΓ)Td u d F u T T T p σεδδ 其中,??ΩΓ Γ-Ω-=∏Td u d F u T T T p )21 (σε即为系统的总势能,它是弹性体变 形势能和外力势能之和。上面变分为零式表明:在所有区域内满足几何关系,在边界上满足给定位移条件的可能位移中,真实位移使系统的总势能取驻值(可证

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)的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。 有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。 在数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的

有限元素法有限体积法有限差分法有限容积法的区别

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

电磁仿真算中的有限元法

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

结构分析及有限元分析基础知识

第一章结构分析及有限元分析基础知识 注:摘自《NX知识工程应用技术——CAD/CAE篇》 洪如瑾编译 清华大学出版社 [目标] 本章将简述结构分析及有限元分析的基础知识,为学习与应用结构分析做好准备,包括: ※ 结构与结构分析定义 ※ 结构的线性静态分析 ※ 材料行为与故障 ※ 有限元分析的基本概念 ※ 有限元模型 1.1结构分析基础知识 1.1.1结构基本概念 1.结构定义 结构可以定义为一个正承受作用的载荷处于平衡中的系统。平衡条件意味着结构是不移动的。一个自由的支架不是一个结构,它未被连接到任一物体上并无载荷作用与它。仅当它附着到外部世界,并且有作用力、压力或力矩时,支架成为一个结构。 例如横跨江面的大桥就是一个普通的结构,一个支架通过它的支撑连接到地面上,桥的重量是在结构上的一种载荷(力)。当汽车通过桥时,附加的力作用于桥的不同位置。 一个好的结构必须满足以下标准: (1) 当预期的载荷作用时,结构必须不出现故障。这个似乎是显而易见的,并意味着结构必须是“强度足够的”。故障意味着结构破裂、分离、弯曲,以及支撑作用载荷失败。 注意:考虑到意外的载荷,通常在设计中提供安全余量。余量常常利用安全因素来描述。例如,如果在结构上期待载荷是10 000磅,规定安全因素是2.0,则结构将设计成能经受住20 000磅载荷。 (2) 当载荷作用时,结构必须不产生过分变形。这意味着结构必须“刚度足够”。 变形可接受的极限(弯曲度、挠度、拉伸等)取决于特定情况。例如,在通常住宅中的地板由足够的吊带支撑,以防止当人在地板岸上行走时有“柔软”的感觉。 (3) 在它的服务生命周期,结构的行为应不会恶化。这意味着结构必须“足够耐用”,必须考虑环境影响和“磨损与破裂”。如果一座桥假定维持50年,则桥的设计必须提供整个50年寿命的结构完整性与充分的安全余量。2.结构分析 结构分析是用于决定一个结构是否将正确完成任务的工程分析过程。结构将在某些方式中进行模拟和求解描述它的行为的数学方程。分析可以人工方法或用计算机方法来完成。 结构分析的结果(答案)用于评估性能,摘要如下: (1)“强度足够吗?”:应力必须是在一可接受的范围内。 (2)“刚度足够吗?”:位移必须是在一可接受的范围内。 (3)“耐用度足够?”:对一个长的疲劳周期应力必须足够低。

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

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

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

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

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

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

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

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

有限元、有限差分法

有限元法原理 将连续的求解域离散为一组单元的组合体,用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数,近似函数通常由未知场函数及其导数在单元各节点的数值插值函数来表达。从而使一个连续的无限自由度问题变成离散的有限自由度问题。 运用步骤 步骤1:剖分: 将待解区域进行分割,离散成有限个元素的集合.元素(单元)的形状原则上是任意的.二维问题一般采用三角形单元或矩形单元,三维空间可采用四面体或多面体等.每个单元的顶点称为节点(或结点). 步骤2:单元分析: 进行分片插值,即将分割单元中任意点的未知函数用该分割单元中形状函数及离散网格点上的函数值展开,即建立一个线性插值函数 步骤3:求解近似变分方程 用有限个单元将连续体离散化,通过对有限个单元作分片插值求解各种力学、物理问题的一种数值方法。有限元法把连续体离散成有限个单元:杆系结构的单元是每一个杆件;连续体的单元是各种形状(如三角形、四边形、六面体等)的单元体。每个单元的场函数是只包含有限个待定节点参量的简单场函数,这些单元场函数的集合就能近似代表整个连续体的场函数。根据能量方程或加权残量方程可建立有限个待定参量的代数方程组,求解此离散方程组就得到有限元法的数值解。有限元法已被用于求解线性和非线性问题,并建立了各种有限元模型,如协调、不协调、混合、杂交、拟协调元等。有限元法十分有效、通用性强、应用广泛,已有许多大型或专用程序系统供工程设计使用。结合计算机辅助设计技术,有限元法也被用于计算机辅助制造中。

有限差分法the Finite Difference Method 微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组,解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 在采用数值计算方法求解偏微分方程时,若将每一处导数由有限差分近似公式替代,从而把求解偏微分方程的问题转换成求解代数方程的问题,即所谓的有限差分法。有限差分法求解偏微分方程的步骤如下: 1、区域离散化,即把所给偏微分方程的求解区域细分成由有限个格点组成的网格; 2、近似替代,即采用有限差分公式替代每一个格点的导数; 3、逼近求解。换而言之,这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程(Leon,Lapidus,George F.Pinder,1985)

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

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

结构分析中有限元法课程建设的问题和方式

结构分析中有限元法课程建设的问题和 方式 本文从网络收集而来,上传到平台为了帮到更多的人,如果您需要使用本文档,请点击下载按钮下载本文档(有偿下载),另外祝您生活愉快,工作顺利,万事如意! 随着复杂工程结构及巨型结构的大量出现,在土木、水利等工程专业本科生中加强结构分析能力(包括电算能力)的培养显得尤为迫切。笔者将结构力学中的“矩阵位移法”和弹性力学中的“有限元法”有机结合,按照培养创新型高级专门人才的要求建设了结构分析中的有限元法课程,精心编写并出版了该课程所用的教材和电算程序。根据该课程的特点,对课堂教学、上机实习和考试等教学环节进行了改革,充分发挥教师的主导作用和学生的主体作用,取得了良好的教学效果。 关键词:结构分析;有限元法;课程建设 中图分类号:G6423文献标志码:A文章编号: 10052909(2015)02005304 一、加强结构分析能力的培养是课程建设的初衷 在土木、水利等本科工程专业的教学过程中,通过前期基础课程和专业基础课程的学习,学生初步具备了对计算工作量不大的简单结构进行结构分析的能

力。如通过结构力学课程经典理论的学习,学生可以对静定平面杆系结构进行分析,也可用力法或位移法等方法分析未知量较少的超静定平面杆系结构,计算其内力和位移。通过弹性力学基本理论的学习,学生可求出几何形状规则(如矩形或圆形)、边界条件简单(如四边固支或四边铰支)的结构在单一荷载(如均布荷载)作用下的内力解析值。 随着各行各业现代化建设的深入开展和城镇化建设的加速推进,房地产和土木工程建筑行业已成为国家的重要支柱产业,与之相伴的是大跨度结构、高层高耸结构等各种复杂结构和巨型结构的出现。这些结构中,有些是形状、边界、荷载等较复杂的连续体板壳结构或实体结构,有些虽是杆系结构,但却是空间杆系结构或计算工作量庞大的平面杆系结构,还有些是杆系结构与连续体结构的组合体。无论对以上哪种结构进行分析,都必须利用数值分析法才能进行,学生既有知识已明显不足。 在一般本科院校开设的结构力学课程体系中,通常都要介绍适合数值分析的矩阵位移法。但有的仅讲述了矩阵位移法基本原理而未涉及程序使用,有的虽然让学生使用了部分程序,却较少或几乎不涉及对空间杆系结构的分析。实际工程中较为复杂的连续体结

有限差分法

有限差分法有限差分法 finite difference method 微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组,解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 有限差分法的主要内容包括:如何根据问题的特点将定解区域作网格剖分;如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确,还需从理论上分析差分方程组的性态,包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式,为了有实用意义,一个基本要求是它们能够任意逼近微分方程,这就是相容性要求。另外,一个差分格式是否有用,最终要看差分方程的精确解能否任意逼近微分方程的解,这就是收敛性的概念。此外,还有一个重要的概念必须考虑,即差分格式的稳定性。因为差分格式的计算过程是逐层推进的,在计算第n+1层的近似值时要用到第n层的近似值,直到与初始值有关。前面各层若有舍入误差,必然影响到后面各层的值,如果误差的影响越来越大,以致差分格式的精确解的面貌完全被掩盖,这种格式是不稳定的,相反如果误差的传播是可以控制的,就认为格式是稳定的。只有在这种情形,差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有以下3种方法。最常用的方法是数值微分法,比如用差商代替微商等。另一方法叫积分插值法,因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理,一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛

相关主题