弹性波动方程数值解的有限元并行算法

弹性波正演模型的快速有限元算法
钱勇先
【期刊名称】《石油地球物理勘探》
【年(卷),期】1990(000)003
【摘要】用有限元法做弹性波正演模型计算已进入实用阶段,但有限元法所要求的大存储量和计算速度慢的缺点依然存在。

本文在深入分析有限元算法的基础上,注意到刚度矩阵,耦合矩阵在特殊单元中的对称性、耦合矩阵具有的转置性及某些特殊特征,采用正方形单元并将与节点相乘的元素归并到时间离散中运算,使计算速度提高一倍,使整体矩阵的存放量减少一半。

文中最后给出了一个计算实例,证实该快速算法的正确性。

【总页数】8页(P260-266,378)
【作者】钱勇先
【作者单位】江汉石油学院物探系
【正文语种】中文
【中图分类】TE13
【相关文献】
1.弹性波动方程数值解的有限元并行算法 [J], 王月英;孙成禹
2.一种基于矢量有限元与多层快速多极子技术的电磁散射快速并行算法 [J], 袁军;刘其中;郭景丽
3.基于能量最小化原理的弹性波CT成像频域有限元反演算法 [J], 薛龙;刘天云;张建民
4.拟地震法电磁波场有限元法正演模型 [J], 沈飚;何继善
5.基于能量最小化原理的弹性波CT成像频域有限元反演算法 [J], 薛龙;刘天云;张建民;
因版权原因，仅展示原文概要，查看原文内容请购买。

205第八章 弹性连续体振动的有限元解上一章讨论的弹性连续体振动问题都只是在简单的特殊边界情况下才能得到精确解，而对于复杂弹性连续体的振动，通常无法得到精确解。

因此，只能采用近似解，近似解方法很多，其要旨在于将无限自由度系统（连续体）变换成为有限多自由度系统（离散系统）来处理。

下面仅介绍现代工程中常用的有限元分析法（简称有限元）。

有限元的基本思想是将一个复杂结构（连续系统）看成是有限个基本元素（单元）在有限个结点彼此相联结的组合结构。

每个单元都是一个弹性体。

有限元法通常是采用位移法，即以结点处的位移作为基本未知量，单元的位移是用结点位移的插值函数表示，单元以至整个结构的一切参数包括位移、应变、应力等都通过结点位移表示出来。

从振动问题来看，最后是将一个连续体的振动问题变成了一个以有限个结点位移为广义坐标的多自由度系统的振动问题。

由于平面刚架所采用的单元综合有受轴向力杆及横向弯曲的两种作用，因此本节着重通过平面刚架（杆系结构）的振动为例来阐明有限元的基本概念与分析步骤。

有限单元法分析过程基本上可分为结构离散化、单元分析、整体分析三个步骤。

一. 结构的离散化如图8-1所示平面刚架，它由三根杆件组成，杆件之间通过刚结点2、3相互连接，若在某些结点上受有作用在刚架平面内的激扰力，这就是一个平面刚架的振动问题。

结构的离散化就是将结构分成若干个单元。

在此例中，取三个杆各为一个单元。

共分成三个单元，在图中分别以①、②、③标明各单元的标号。

各单元之间以及它们与基础的联结点称为结点。

结点编号分别为1、2、3、4。

结点1、4是杆件和基础固定联结的结点，它们不能运动，称为固定结点。

2、3点是杆件间的联结点，是可以运动的，称为可动结点。

可动结点的数目乘以每个结点的基本未知量的数目就是整个结构的自由度数。

对于工程结构来说，结构划分的单元愈多，则结构的自由度数愈多，计算工作量就愈大，但计算结果愈有可能趋于精确解，因此在结构离散化过程中，划分单元的多少，应综合考虑精确度与计算工作量的因素。

第二章 弹性力学问题有限元方法的一般原理和表达格式 2.1 引言本章将讨论通过弹性力学变分原理建立弹性力学问题有限元法列式的基本步骤。

最小位能原理的未知场变量是位移，以结点位移为基本未知量，并以最小位能为基础建立的有限单元位移元。

它是有限元方法中应用最普遍的单元。

对于一个力学或物理问题，在建立其数学模型以后，用有限元方法对它进行分析的首要步骤是选择单元形式。

平面问题三结点三角形单元是有限元方法最早采用，而且至今仍经常采用的单元形式。

我们将以此作为典型，讨论如何应用广义坐标建立单元位移模式与位移插值函数，以及如何根据最小位能原理建立有限元求解方程的原理、方法与步骤，并进而导出弹性力学问题有限元方法的一般列式。

2.2 弹性力学平面问题的有限元列式2.2.1 单元位移模式及插值函数典型的三结点三角形单元结点编码为i,j,m 。

每个结点有两个位移分量，如图2.2所示。

每个结点的位移可用位移矢量i α表示，即⎥⎦⎤⎢⎣⎡=i i i v u α ),,(m j i每个单元有6个结点位移分量（称为6个自由度），于是单元结点的位移向量可表示为[]Tm m j j i im j i e v u v u v u =⎥⎥⎦⎤⎢⎢⎣⎡=ααααe α为单元结点位移列阵。

1．单元的位移模式和广义坐标在有限元方法中单元的位移模式，是指在单元内位移的插值函数，其一般形式采用多项式作为近似函数，因为多项式运算简单，并且随着项数的增多，可以逼近任何一段光滑的函数曲线。

假设3结点三角形单元位移模式选取一次多项式y x u 321βββ++=y x v 654βββ++= （2.2.1）它的矩阵形式是φβ=u （2.2.2）其中⎥⎦⎤⎢⎣⎡=v u u ，⎥⎦⎤⎢⎣⎡=ϕϕφ00 []y x 1=ϕ[]T 654321βββββββ=由于三个结点也在单元内，满足位移模式，于是得i i i y x u 321βββ++=j j j y x u 321βββ++= （2.2.3） m m m y x u 321βββ++=上式是关于321,,βββ的线性方程组。

弹性体力学数值模拟算法弹性体力学数值模拟算法是一种用于模拟弹性材料力学行为的数值计算方法。

该算法可以通过计算机模拟得出弹性体的应力分布、变形情况以及其他相关参数，为工程设计、科学研究以及现实世界中的弹性材料应用提供了有力支持。

一、弹性体力学概述弹性体力学是研究物体由于外力而发生形变时的力学行为，主要关注物体的应力、应变以及相应的力学性质。

弹性体力学数值模拟算法正是基于这些力学特性进行计算。

在弹性体力学中，弹性模量、剪切模量以及泊松比等参数是十分关键的。

弹性体力学数值模拟算法通过考虑这些参数，结合物体的几何形状以及外力条件，利用计算方法来模拟物体的力学行为。

二、有限元法数值模拟算法有限元法是弹性体力学数值模拟算法中常用的一种方法。

该方法将弹性体划分为许多小的有限单元，在每个单元内进行数学计算，再通过组合这些小单元的结果，得出整个弹性体的力学行为。

具体而言，有限元法首先进行网格划分，将弹性体分割成许多小单元，然后在每个单元内利用数学公式推导出该单元的力学行为。

这些小单元之间通过边界条件相互连接，构成整个体系。

最后，通过求解线性方程组，得到弹性体的应力分布和变形情况。

有限元法数值模拟算法的优点是适用于各种复杂形状的弹性体，并且精度较高。

然而，该方法在处理大规模模拟问题时，计算量较大，需要较高的计算资源。

三、有限差分法数值模拟算法有限差分法是另一种常用的弹性体力学数值模拟算法。

该方法将弹性体通过网格划分为一系列小的节点，在每个节点处利用差分近似推导出该节点的力学行为。

有限差分法通过离散化的方式，将偏微分方程转化为差分方程，再通过求解这些差分方程得到相应的力学参数。

这个过程类似于对输出信号进行采样，将连续的物理问题离散化为一个个节点上的差分问题。

有限差分法数值模拟算法的优点是计算较为简单，容易实现。

但是对于复杂形状的物体，可能需要更多的节点和更小的步长来保证计算精度。

四、其他数值模拟算法除了有限元法和有限差分法，还有其他一些弹性体力学数值模拟算法，如边界元法、网格无关法等。

§5 弹性理论与有限元解法有限元法用来解弹性体的小变形问题是很有效的.这是由于变形能与外力势能或总势能可以表示为形式划一的二次泛函.由于弹性力学问题提法,特别是边界条件的复杂性,要使有限元解法灵活运用,还需要对各种问题作些分析与讨论.因此本节除了介绍各类问题与§1中B 、D 等有关的基本关系式外，还着重分析边界条件在有关变分问题中所起的作用，并讨论各类变分问题与微分方程定解问题的等价性。

一、 三维的弹性问题本节在直角坐标系(x ,y ,z )中讨论弹性体内部受力与形变的情况.为了简化论述,假设弹性体是单连通、均匀并各向同性的,必要时还用通常的矢量记号以突出其力学意义.[应力与平衡方程] 弹性体在外力、温差等作用下，各部分之间将产生内力，表示内力的大小是作用面上受力的强度或单位面积所受的力,即所谓应力.对弹性体内的一点P ,可在其邻近作一微元六面体，其棱边平行于坐标轴。

六面体的六个面中有三个面其外法线方向分别与x 轴、y 轴或z 轴同向（其余三个则与坐标轴反向）,可分别称为x ,y ,z 坐标面(或第1,2,3坐标面).设321,,σσσ分别表示作用于x ,y ,z 坐标面的应力,并记i σ的三个分量为321,,i i i σσσ(对其余三个面的应力取号相反.例如,图19.18上所示的分量全是正的),这九个分量ij σ(i ,j =1,2,3) 构成一个张量,称为应力张量.从图可看出,分量ii σ表示第i 坐标面上的正应力(受拉取正,受压为负);而分量ij σ则表示沿第i 坐标面剪应力的二分量(使扭角变成锐角的为正).由微元六面体力矩的平衡可得剪应力互等定律,即ij σ =ji σ (j i ≠)因此应力张量是对称的,其分量只有六个是独立的.在有限元解法中,一般是把这六个分量按如下次序排成列矢量,并记作τσσσσσσσ),,,,,(}{123123332211=设过点P 作一任意的斜面,其法线n 的方向余弦为(321,,n n n ),则利用与三坐标面围成的四面体的平衡条件可得作用于该斜面的应力n σ的三个分量∑=++=31332211i j j j iij n n n n σσσσ(j =1,2,3)于是n σ在方向余弦为(321,,m m m )的直线上的投影,即其分量为∑==31,j i j i ij nm m n σσ这表明,一个点的应力状态完全由该点的六个应力分量{σ}所决定.例如,作用于上述斜面的正应力为∑==31,j i j i ij nn n n σσ而沿该斜面的剪应力方向取n n ⨯⨯n σ,大小等于∑∑∑===-=-=313131,2222)()(j i j i j i ij iij nnn n n n n σσσστ反过来,对于0=n τ的斜面称为P 点的应力主面,相应的法线称为P 点的应力主轴,而其正应力nn σ称为P 点的主应力.可以证明,在弹性体内任意一点,一定存在三个互相正交的主应力,而且其中最大(小)的一个就是该点的极大(小)正应力.三个正应力之和332211σσσΘ++= 称为体积应力,它在坐标变换下是个不变量,因而等于三个主应力之和.设τ),,(}{321f f f f ==f ,表示作用于P 点的体力(单位体积的外力),对体积元素作积分或均值分析都可以推出力的平衡方程0321=+∂∂+∂∂+∂∂i ii i f zy x σσσ (i =1,2,3) (26) [应变与几何方程] 弹性体内任一点P (x ,y ,z )在小变形后移到),,(z y x P ''''其位移函数为τ),,(w v u =u式中z z w y y v x x u -'=-'=-'=,,它们是(x ,y ,z )的微量函数,假定有一微小线段PN =dr 其方向余弦为(321,,n n n ),经过小变形变为r d N P '='',则沿该方向的正应变n ε定义为单位伸长,即drdr r d n -'=ε从变形前后的r d '与dr 关于位移的表达式不难得出1}])1([])1([])1({[1232123212321-+++++++++++=n w w n w n v n n v v n u n u n u n z y x z y x z y x n ε展开右端的根式并略去高阶无穷小量(即位移导数的高次项),就得到)()()(211332232221y x x z z y z y x n u v n n w u n n v w n n w n v n u n ++++++++=ε设另一线段PM 的方向余弦为(321,,m m m ),变形前的夹角θ=∠MPN ,则33221131cos m n m n m n m n i i i ++==∑=θ (27)设变形后的夹角θ'='''∠N P M (图19.19),则∑=''='31cos i i i m n θ根据变形前后两线段大小、方向的变化,不难得出))(())(()(2cos )1(cos 31132332332211x x z y z y x m n w u m n m n v w m n m n w m n v m n u m n +++++++++--='θεεθ))((1221y x u v m n m n +++ (28)对照(27),(28)可知,只要在P 点给定如下六个导数值:yx x z xy x y x u v w u v w wv u +==+==+==⎪⎩⎪⎨⎧===211213313223332211εεεεεεεεε (29)就可以完全确定P 点邻近的变形状态.ii ε表示沿第i 坐标轴的正应变,ij ε表示经过小变形第i ,j 坐标方向之间的直角改变量即所谓剪应变(12ε如图19.20所示,变形后成锐角为正,成钝角为负).这六个量称为应变分量,记作τεεεεεεε),,,,,(}{123123332211=同样可以证明,在弹性体内任意一点,一定存在互相正交的应变主轴,变形后三轴交角仍然保持直角,即剪应变为零;三主轴的正应变称为主应变,而且其中最大(小)的一个就是该点的极大(小)的正应变.三个正应变之和332211εεε++=e称为体积应变,也是个不变量,而且表示微元中每单位体积的改变量.对于各向同性体来说,应力主轴与应变主轴的方向还是一致的.§1中所述的变形能密度}{}{21}{}{21d εσσεττ==V正是把坐标系变换到共同的主轴方向并应用虎克定律而推导出来的.关系式(29)称为几何方程,其矩阵形式为)92(000000000}{123123332211'≡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=u B w v u xy x z y z z y xεεεεεεε[物理方程与弹性系数] 对杆件作简单拉、压的小变形实验证明,单独的轴向(取为x 轴)力不会引起剪应变,其正压力x σ与正应变x ε之间具有如下的线性关系x x E εσ=或1111εσE =如果考虑到沿x 方向的伸长还伴有侧向收缩,则产生沿y ,z 方向的正应变x z y νεεε-==或113322νεεε-==这就是虎克定律,式中系数E ,v 分别称为弹性模数和泊松比.从胡克定律就可以推导出应变与应力之间的一般关系式,即所谓物理方程)30(1)1(2)]([11)1(2)]([11)1(2)]([1121212221133333131313311222223232333221111⎪⎪⎪⎩⎪⎪⎪⎨⎧=+=+-==+=+-==+=+-=σσεσσνσεσσεσσνσεσσεσσνσεGE v E G E v E G E v E式中)1(2v E G +=称为剪切弹性模数. 从(30)求逆得出应力与应变之间的关系式)(),3,2,1(2j i G i G e ijij ii ii ≠==+=εσελσ或写成矩阵形式:)31(}{}{000000000000000000200020002}{εελλλλλλλλλσD G G G G G G≡⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+++=式中)21)(1(v v Ev -+=λ与G 称为拉梅系数,e 就是体积应变.它与体积应力Θ成正比:ΘE v e 21-= 或 e vE 21-=Θ式中比例常数vE21-称为体积弹性模数.[边界条件] 弹性体Ω的边界Ω'承受面力τ),,(}{321q q q q ==q有三种方式:固定支承,荷载支承和弹性支承.假定Ω'接受这三种支承的部分分别记作2,ΩΩ''1和3Ω',则其边界条件可表示为 1° 几何约束条件: 在1Ω'上给定位移,即w w v v u u ===',,1Ω2° 面力平衡条件: 在2Ω'上给定荷载即面力q ,(321,,n n n )表示2Ω'上任一面积元素的外法线方向余弦.由于应力要与面力平衡,从(27)可得其条件为)3,2,1(∑312=='=j q n i ji ijσΩ3° 耦合平衡条件: 在3Ω'上弹性体与另一弹性结构耦合,这些耦合边界上的位移既不受约束也不完全自由,而是接受与其位移偏差(相对于某一给定的位移值w v u ,,)成正比的弹性反力.每单位面积上它的三个分量可表示为)3,2,1(][)]()()([321321=+++-=-+-+--j q w c v c u c w w c v v c u u c jj j j j j j这里弹性支承系数矩阵C =(ji c )是正定的,而w c v c u c q j j j j 321++=可看作给定的面力.同样,这反力应由3Ω'上的应力来平衡.于是其条件可写成)3,2,1()(-∑313213=+++='=j q w c v c u c n i jj j j i ijσΩ[外力势能及其计算公式] 弹性体内及其边界上,凡给定的外力因变形而作功的部分都要累加起来,再补上负号可看作弹性体相对于外力系统的势能,即所谓外力势能-F (u ).例如在弹性体内部除体力f 外,也可能有部分施加集中的面力、线力和点力，这些力与该部分的位移的内积就是它们所作的功，因此都要算进去。

求解弹性波动方程弹性波动方程是波动力学中的重要概念，描述了固体、液体等弹性体中的波动传播情况。

求解弹性波动方程可以帮助我们更好地了解地震、地质勘探、材料力学等领域中的问题。

本文将从弹性波动方程的概念、求解方法及应用方面进行探讨。

一、弹性波动方程的概念弹性波动是指在弹性介质中，由于外力的作用导致的能量传递方式，其传播速度与介质的物理性质有关。

弹性波动方程则是描述弹性波传播过程的数学模型，通常也称为弹性振动方程。

弹性波动方程通常包含两个基本变量：质点的位移和质点的速度。

在三维坐标系中，弹性波动方程可以写成如下形式：ρ∂²u/∂t² = ∂iτij/∂xj其中，ρ为介质密度，u为波动的位移，τij为应力张量，xj表示坐标轴方向，i、j为坐标轴的编号。

二、求解弹性波动方程的方法求解弹性波动方程可以采用有限元法、有限差分法、声波分解法等多种方法。

其中，有限元法和有限差分法是比较常用的求解方法。

1. 有限元法有限元法是一种比较通用的求解偏微分方程的方法，其主要思想是将求解区域离散化成若干个单元，然后对每个单元内的位移和应力进行求解。

最后将所有单元的解合并起来，就可以得到整个区域内的位移和应力分布情况。

2. 有限差分法有限差分法是一种利用差分近似方法来求解微分方程的数值方法。

其思想是将求解区域划分为网格，然后在每个网格节点处用差分公式来逼近微分方程的解。

最后通过数值计算来得到整个区域内的解。

三、弹性波动方程的应用弹性波动方程在地质勘探、地震学、材料力学等领域中有广泛的应用。

1. 地质勘探地质勘探中需要了解地下结构，以便确定油气、煤炭等矿产资源的位置和规模。

在勘探中，我们可以通过广播一定频率的声波，然后通过测量接收到的反射波、折射波等来推导地下结构的性质。

弹性波动方程能够有效地描述声波在地下结构中的传播情况，从而帮助勘探人员更准确地掌握地下资源的分布情况。

2. 地震学地震学是研究地震、地震波等现象的学科，也是弹性波动方程的重要应用领域。

波动方程有限元—差分法数值解有限元差分法是新兴的数值解法，可以用于计算复杂的波动方程。

本文研究了有限元差分法在解决波动方程的数值方面的应用，主要从以下几方面来探讨。

首先，介绍了波动方程的定义及其特性，简要介绍了波动方程模型的建立和求解，以及有限元差分法在求解波动方程中的应用。

其次，提出了三种典型的有限元差分法求解波动方程的算法：时间歉分差分法、时间步长变化差分法和双时间空间差分法。

之后给出了三种算法的实现过程，通过计算了实例，验证了有限元差分法的有效性。

最后，归纳总结了有限元差分法的特点和优点，指出除了有限元差分法外，还可以使用基于全时空的解析方法求解波动方程的有效性。

波动方程是用于描述物理系统的变化的常微分方程，例如用于描述热传导问题的拉普拉斯方程，用于描述电磁学场的马尔科夫方程，以及用于描述流体动力学问题的贝尔加拉里散度方程等。

由于波动方程的解一般是复杂的，有限元差分法可以利用有限个均匀分布的有限元素网格来模拟波动方程，然后采用差分法进行数值求解。

有限元差分法的使用可以使波动方程的解变得更加精确，同时由于其求解的过程是近似的，因此也可以节省计算时间，节约资源。

此外，有限元差分法还可以应用于并行计算，以提高计算速度。

由于近年来随着计算机性能的不断提高，有限元差分法在解决波动方程的数值解中逐渐被认可，逐渐被广泛应用。

综上所述，有限元差分法是一种新兴的数值解法，能够比较有效地求解复杂的波动方程，有助于深入研究物理系统的变化，在实际应用中有重要的意义。

因此，有关有限元差分法在波动方程求解中的应用，以及这种求解方法的特点和优点，应予以深入研究，以挖掘这种数值方法的潜力，为解决实际物理问题提供有效的数值支持。

1.1控制方程在经典的线弹性理论当中，线弹性均匀介质的运动方程为⋅⋅=+∂∂i i jij u f x ρρσ (3,2,1=i ) (1-1)其中，应力张量ij σ和应变张量kl ε满足Hooke 定律，即满足本构方程:kl ijkl ij c εσ= (1-2)应变张量kl ε可通过位移矢量i u 确定，从而给出几何方程:)(21,,k l l k kl u u +=ε (k,l=1,2,3) (1-3)各向同性线弹性材料的弹性特征可以用下列材料常数进行描述:弹性模量E,剪切 模量μ、泊松比v 、体积模量K 和Lame 常数λ等。

这些常数中任意两个独立，而其它 几个常数能够用它们间接进行表示。

在各向异性的情形下，则需要采用21个相互独立 的弹性常数ijkl C ，进行表示。

介质表面S 上的边界条件为u i i S x t x u t x u ∈=),,(),( (1-4)σσS x t x p n t x i j ij ∈=),,(),( (1-5)由方程(2-1)-(2-3)这三个线弹性动力学的控制方程所描述的线弹性动力问题的初始 条件为V x x u t x u i i ∈=),(),(00 (1-6)V x x v t x u i i ∈=⋅),(),(00 (1-7)其中，)(x u i ，)(0x v i ，),(t x u i 和),(t x p i 是已知量。

由弹性动力学问题解的唯一性定理可知:若弹性体(体积为Y ，表面为S)的解能 满足方程(1-1)-式(1-7)，则其位移场、应力场和应变场的解答是唯一的。

将三个控制方程合并，可以得到用位移表示的运动方程:⋅⋅=+i i lj k ijkl u f u c ρρ, (3,2,1=i ) (1-8)方程(1-8)所示的就是著名的Navie-Cauchy 运动方程。

而在各向同性线弹性的情形 下，则有)(jk il jl ik kl ij ijkl c δδδδμδλδ++= (1-9)因此，各向同性材料的本构关系也为表示为⎩⎨⎧=≠=+=)(,1)(,0,2j i j i ij ij kk ij ij δμεελδσ (1-10)1.2有限元方程现在我们考虑一个任意形状的封闭区域，并将其用一定形式的网格划分为相应的有 限单元，各单元的单元矩阵可由如下所示的虚功方程求得。

弹性力学与有限元法分析弹性力学是固体力学的一个重要分支，是研究弹性固体在受外力作用、温度改变、边界约束或其他外界因素作用下而发生的应力、形变和位移状态的科学。

有限单元法是力学、数学、物理学、计算方法、计算机技术等多种学科综合发展和结合的产物，是随着计算机技术的广泛应用而迅速发展起来的一种数值分析方法。

有限元法的基本思想就是化整为零，分散分析，再集零为整。

即用结构力学方法求解弹性力学问题，实质是将复杂的连续体划分为有限多个简单的单元体，单元体之间仅仅通过结点相连，实现化无限自由度问题为有限稀有度问题，将连续场函数的（偏）微分方程的求解问题转化为有限个参数的代数方程组的求解问题。

有限元方法经过近半个世纪的发展，目前已经成为各种工程问题特别是结构分析问题的标准分析方法，而有限元软件也已成为现代结构设计中不可缺少的工具。

有限元软件是有限元理论通向实际工程应用的桥梁，它的应用极大地提高了力学学科解决自然科学和工程实际问题的能力，进一步促进了有限元方法的发展。

ANSYS 软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，广泛用于机械制造、石油化工、航空航天、汽车交通、土木工程、造船、水利等一般工业及科学研究。

ANSYS 软件的组成：（一）前处理模块该模块为用户提供了一个强大的实体建模及网格划分工具，可以方便的构造有限元模型，软件提高了100种以上的单元类型，用来模拟工程中的各种结构和材料。

包括：1.实体建模：参数化建模，布尔运算及体素库，拖拉、旋转、拷贝、蒙皮、倒角等。

2.自动网格划分，自动进行单元形态、求解精度检查及修正。

3.在集合模型上加载：点加载、分布载荷、体载荷、函数载荷。

4.可扩展的标准梁截面形状库。

（二）分析计算模块该模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力。