计算力学的数值模拟
发布网友
发布时间：2022-04-29 18:46
共1个回答
热心网友
时间：2022-06-19 17:24
力学现象的数学模拟，常常归结为求解常微分方程、偏微分方程、积分方程、或代数方程。求解这些方程的方法有两类：一类是求分析解，即以公式表示的解；另一类是求数值解，即以成批数字表示的解。很多力学问题相当复杂，特别是复杂的偏微分方程组，一般难以得出它们的分析解，而用数值方法求解则运算步骤繁复，耗用人力很多，因此在电子计算机出现以前，非不得已不用。20世纪50年代以来，出现了配有现代程序设计语言的通用数字计算机。计算机的快速运算和大存贮量，使解复杂的力学问题成为可能。三十多年来，随着计算机的改进，数值方法得到广泛的应用和很大的发展；主要是考虑算得更快、更准、省钱，并为原先不能算的问题构造算法。数值方法很多，求解偏微分方程数值解，以有限差分方法和有限元法使用最广；此外，还有变分方法、直线法、特征线法和谱方法，等等。这些方法的实质绝大多数是将偏微分方程问题化成代数问题，然后再用计算机求未知函数的数值解。 有简单、灵活和通用性强等特点。用差分方法求数值解时，须先将自变量的定义域“离散化”，即只企图算自变量定义域中有限个点的未知函数的近似值。如果自变量只有一个，则可把要计算的区间离散成个线段。如果自变量有两个，而计算区域是图1[二变量区域的离散化]所示的矩形，则最简单的离散方式是把区域分成乘个小矩形。小矩形的长 和宽分别叫作方向和方向的步长。微分方程中出现的偏导数(，)， 在微积分中是差商的极限，在有限差分方法中则代以差商。如图1[二变量区域的离散化]中点的有的情形可代以差商(()-())/2，有的情形可代以(()-())/，如果有二阶偏导数，常常可代以二阶差商(()-2()+())/2，其中()、()和()分别表示相应点的值。 如以适当的差商来代替微分方程每一个导数，就得到对应于
原微分方程的差分方程怎样选差商至关重要。此外，偏微分方程总还要附加边界或初始条件，这些条件也要用差分形式表示。这样，对于每个网格点的未知函数值作出未知量的代数方程组。如果网格分得较密，即步长和都比较小，或与 的数值都比较大，则所得代数方程组的未知量的数目将很大，但借助计算机，还是可以很快求出解来。由于步长无法取为零，因此用差分方法只能求得原微分方程的近似解。但只要选择合理的差商和步长，计算结果仍能令人满意，有时还能得到精度很高的解。有限元法
这种方法是把计算区域剖分成大小不等的三角形（或其他形状的）单元，然后在各单元上用适当的插值函数来代替未知函数。根据变分原理，可将偏微分方程化成代数方程来求解。这种方法具有很广泛的适应性，特别适于求解具有复杂边界形状和物理条件的问题，而且很容易在计算机上实现。1970年以来已研究出一些适用于广泛的线性问题的有限元通用程序，对工程设计起很大作用。按照有限元法剖分的思想，把汽车外壳剖分成大小不等的许多三角形单元，而对弯曲边界只须裁弯取直即可。在应力变化剧烈和要求精确计算的地方，须把单元取得小些；在变化不剧烈的地方则可取得大些。用这种方法不仅可以适应复杂的区域，还可以尽量减少总的单元数目，从而减少未知量的数目。如果在有限差分方法中用矩形网格，则较难处理如此复杂的区域。
