计算力学及算法研究.pptVIP

计算力学的算法研究 　　算法是计算力学研究的重点。算法的研究沿着两个途径进行：一是不断从计算数学中学习和移植有用的算法，如求解线性代数方程组和特征值问题的各种算法;二是不断提出新的高效算法，如分片插值原理在有限元法和边界元法中的应用，这一方法已被广泛应用于温度场、电磁场和其它物理问题的分析中。 　　有限元法的数学理论研究表明，对于适当构造的单元，当有限元网格无限加密时有限元解应收敛于真实解，但事实上对于一个特定的网格无法对有限元解的精度做出可靠的估计，自适应有限元法较好地解决了这一问题。自适应有限元法的做法是在求得给定网格下的有限元位移和应力之后采用应力恢复方法得到有限元结点应力的改进值，然后借助于形函数插值求得应力分布的改进估计，称为最佳猜测应力。理论研究表示，改进前的应力解与最佳猜测应力解的差可以用来估计有限元解的误差。当按上述方法求得的误差估计大于规定的误差允许值时，可以根据误差预测出结构中各个区域中有限单元的合理尺寸，然后直接采用有限元网格自动生成器重新划分单元，再次对问题求解。这种应力恢复、误差估计和网格重划分过程可以重复多次,直到得出满意的解答为止。可以看出，自适应有限元法的主要特点是依靠有限元分析结果的误差估计，计算出全局或局部单元的误差，依此误差预示下一步有限元分析所需网格尺寸。在结构优化、自由边界等问题中，区域的形状在不断修改，采用不同的网格离散形状不同的区域会得出不同的解答，自适应有限元分析方法有助于得到更为可靠的算法。 自适应有限元法的研究 解析法与数值法的结合 　　计算机软硬件的飞速发展为计算力学的发展提供了有力的工具，但是研究计算效率高、存储占用低的有效算法仍是计算力学研究的重点，解析方法与数值方法的结合提供了一系列高效算法。有限条法就是这样一种半解析方法。早期的有限条法成功的用来解决板壳和箱形梁柱这样一类薄壁结构问题，样条有限条法更易于适应任意形状和支撑条件，能够精确的处理集中荷载。采用群论方法求解具有更为复杂条件的结构取得了显著成效。早期的康托洛维奇法又得到新的发展。康托洛维奇法是将求解偏微分方程的问题化为求解常微分方程问题的一种技术，将未知函数代入泛函求得一组欧拉方程及相应的边界条件，得到未知函数的一级近似解。改进的方法有二：其一是求未知函数的二级近似解或三级近似解;其二是采用变量轮换法。线法也可认为是一种半解析半数值方法。如今，有上百种有效的解析与数值结合的方法在结构计算中得到了广泛的应用。 并行计算力学 　　并行计算方法是在工程计算规模的大幅度增加与计算机的计算能力受到限制的矛盾日益突出的情况下产生与发展的。并行计算力学的研究要解决两个方面的问题：一是并行计算工具;二是并行计算方法。并行计算工具的研究有两个方向：一种是传统的并行计算机系统;另一种是利用局域网或互联网在并行程序设计环境支持下的并行处理系统。并行计算力学的研究开始于子结构并行分析方法，逐渐推广到结构特征值问题、结构动力方程、结构弹塑性分析等方面的并行方法研究。目前的研究热点集中在扩展并行计算的应用领域、高效并行计算方法的研究和并行程序生成器的实现。 结构优化 　　结构优化致力于研究系统而高效率的算法和设计工具，使得设计在满足各种约束条件下实现某种目标的最优化结果。当约束函数和目标函数都是设计变量的线性函数时，响应的优化问题是线性规划问题，这类问题已经得到了很好的解决。而当约束函数或目标函数是设计变量的非线性函数时，相应的问题是非线性规划问题，非线性规划比线性规划的求解要困难的多。工程中较多的是采用准则法和序列线性规划法。在序列线性规划中适当的选取中间变量，在初始设计点附近将原来的非线性规划用一个线性规划来近似，从而求出原问题的近似解。从这个近似解出发又可构造出新的近似线性规划，求解后得到新的近似解，如此迭代下去直到满足足够的精度为止。广泛采用的序列二次规划是收敛速度更高的序列近似规划，目前已研究出了一套十分有效的结构优化算法和相应的结构优化设计软件。优化方法的研究也带动了结构分析方法的发展，利用数学规划工具可以有效的求解能够归结为非线性规划的许多非线性问题。 Thank you