设计师品牌 网站,python自学要多久,安卓手机优化,如何建设一个手机网站有限元法(Finite Element Method)是基于近代计算机的快速发展而发展起来的一种近似数值方法#xff0c;用来解决力学#xff0c;数学中的带有特定边界条件的偏微分方程问题(PDE)。而这些偏微分方程是工程实践中常见的固体力学和流体力学问题的基础。有限元和计算机发展共同构…有限元法(Finite Element Method)是基于近代计算机的快速发展而发展起来的一种近似数值方法用来解决力学数学中的带有特定边界条件的偏微分方程问题(PDE)。而这些偏微分方程是工程实践中常见的固体力学和流体力学问题的基础。有限元和计算机发展共同构成了现代计算力学 (Computational Mechanics)的基础。有限元法的核心思想是“数值近似”和“离散化” 所以它在历史上的发展也是围绕着这两个点进行的。“数值近似”由于在有限元法被发明之前所有的力学问题和工程问题中出现的偏微分方程只能依靠单纯的解析解(Analytical Solution)得到解答。这种方法对数学要求很高而且非常依赖于一些理想化的假定(Assumption)。比如在土木工程中梁柱计算中出现的平截面假定小应变假定理想塑性假定。这些假定其实是和实际工程问题有很大偏差的而且一旦工程问题稍微复杂一些我们就不能直接得到解析解或者解析解的答案误差过大。而有限元法把复杂的整体结构离散到有限个单元(Finite Element)再把这种理想化的假定和力学控制方程施加于结构内部的每一个单元然后通过单元分析组装得到结构总刚度方程再通过边界条件和其他约束解得结构总反应。总结构内部每个单元的反应可以随后通过总反应的一一映射得到这样就可以避免直接建立复杂结构的力学和数学模型了。其总过程可以描述为总结构离散化 — 单元力学分析 — 单元组装 — 总结构分析 — 施加边界条件 — 得到结构总反应 — 结构内部某单元的反应分析在进行单元分析和单元内部反应分析的时候形函数插值(shape function interpolation)和 高斯数值积分(Gaussian Quadrature)被用来近似表达单元内部任意一点的反应这就是有限元数值近似的重要体现。一般来说形函数阶数越高近似精度也就越高但其要求的单元控制点数量和高斯积分点数量也更多。另外单元划分的越精细其近似结果也更加精确。但是以上两种提高有限元精度的代价就是计算量几何倍数增加。为了提高数值近似精度同时尽量较少地提高计算量有限元法经历了很多发展和改良。下图就是一典型的有限元问题因为模型中间空洞部分几何不规则性结构用有限三角单元划分。由于在靠外区域结构反应变化程度不是很大因此划分的单元比较大和粗糙而在内部应力变化比较大划分也比较精细。而在左边单元划分最密区域有应力集中现象(如裂纹问题的奇异解现象)所以又有相应的高级理论(比如non-local theory)来指导这部分的单元应力应变计算。结构被选择性地离散和高级理论构成了有限元发展的主要研究方向。“离散化”离散化和相应单元特性和收敛研究也是有限元中一个重要研究领域总的来说有限单元和他们组装成的总体结构主要分为1-D 单元 (1-D element)杆单元 (bar element) ------ 桁架 (truss)梁单元 (beam element) ------ 框架 (frame)板单元 (plate element) ------ 壳体 (shell)2-D单元 (2-D element) ------ 平面应力体 (plain stress) 和 平面应变体 (plain strain)三角单元 (triangle element)四边形单元 (quadrilateral element)多边形单元 (polygonal element)3-D 单元 (3-D element) ----- 立体结构 (3-D problem)三角体 (tetrahedrons element)立方体单元 (hexahedrons element)多边体单元 (polyhedrons element)具体的分类和单元形状见下图可以看到每种单元又可以提高形函数的阶数(控制点 node 数量)来提高精度。很多有限元研究也集中在这个领域。比如研究新的单元引用于结构动力反应以减小数值震荡比如用3-D单元去模拟梁单元等等。其实理论上来说这个领域可以有无限可能因为对精度和数值稳定的追求可以是无限的。“光滑边界” 和 与CAD的交互问题其实这个算不上有限元的核心思想不过是现在有限元研究热的不能再热的领域了就是Hughes提出的“NURBS”有限元法它的原理是用空间样条曲线来划分单元。如第一幅图所示传统的有限元在处理不规则边界的时候一般都是较多的单元和用三角单元多边形单元来解决而且单元控制点都是和单元在一个平面上。而NURBS 单元的控制点脱离了单元本身并且利用B-spline理论上可以把单元的光滑程度(continuity)提高到无限而且不会显著提高计算量。发展NURBS的另外一个好处是在建模中常用的CAD软件是用B-spline来进行模型建立基础的而NURBS 正好也是用用B-spline作为basis。所以CAD和NURBS的交互可以非常简单和高效的甚至可以说是无缝连接。因此在工业界中十分复杂的模型都可以用CAD进行建模再用NURBS进行有限元计算如下图。现在成吨的有限元paper都来自这个领域因为有限元的基本理论基本已经成熟和robust利用高性能计算机进行大尺度(large-scale)和高复杂结构模拟也是有限元发展的一个主要方向。P.S. 需要提到一点的是没有高性能计算机技术的大力发展就不可能有有限元的发展。有限元的理论最早是出现在1960年代直到1970之后才随着计算机的发展而迅速发展。而现在发展迅速的计算力学也是得益于高性能计算机的发展。可能当某一天计算机处理速度可以强大到我们可以用最复杂最密集的单元完美快速地模拟任意结构我们也不用再操心精度问题了。所以我觉得有限元的核心还需要加上计算机技术的发展吧。作者菲兹-END-
