结构局部修改重分析算法的设计、优化与并行

来源 :北京大学 | 被引量 : 0次 | 上传用户:vickyfucandy
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
有限元分析是工程设计及科学研究中最重要和最成功的数值方法,以其为核心的CAD/CAE软件对国家的经济建设和国家安全起到了非常关键的作用,而自主研发此类软件并提高有限元分析效率一直以来都是工程师和科研人员关注的焦点问题之一。  实际工程中经常遇到结构修改的情形,如材料非线性分析、结构拓扑优化设计、施工模拟、抗连续倒塌设计等等。在有限元分析中,我们可以将上述的多步计算要求描述为:在前步结构上进行修改然后重新进行分析计算。结构修改也可以直接用矩阵变化的形式提出并转化为数学问题:一个线性方程组已经完成求解,当系数矩阵中元素的值发生改变或者增加、删除了某些行(列)时,如何快速得到修改后方程组的解?  结构修改一直是计算力学研究的热点之一,然而目前已有的重分析算法尚不能完全满足实际工程的需求,对于大规模问题的大量修改或者拓扑修改往往力不从心。针对这个问题,本文在目前通用的稀疏矩阵求解方案基础上,借助图剖分填充元优化的结果,提出了一种新的结构局部修改算法,并从软件和硬件两个方面进行了大量优化,进一步提高计算效率。  本文的主要工作如下:  首先,本文提出了结构修改算法的新思路:更新刚度矩阵的三角分解因子;当分解修改刚度矩阵时,只需要在原三角分解因子的基础上,沿着剖分二叉树从修改节点回溯到根节点的路径进行修改即可,其余部分则保持不变。基于该思路,从非拓扑修改入手,共提出了三种不同的修改算法。算法Ⅰ直接跳过未修改行的重分解,只计算修改行;算法Ⅱ改为计算修改行的变化量;算法Ⅲ则结合了前两种算法的优点,依计算量对各行分别选择算法Ⅰ或算法Ⅱ进行计算,效率最高。对于修改量较少的局部修改问题,算法Ⅲ能够减少大约70%到90%的计算量。  其次,本文将非拓扑修改算法扩展到整体构架不变的拓扑修改情形,并设计出统一的计算流程,降低了拓扑修改算法的应用难度。其核心思想是:将有限元方程系数矩阵中添加或删除的行(列)视为一种特殊的修改。算法Ⅰ只需重新计算修改行即可;算法Ⅲ则建立了精确的行间关系和计算式。拓扑修改与非拓扑修改的统一大大地拓宽了修改算法的应用范围。  为了充分利用多核CPU的硬件资源,本文实现了稀疏矩阵求解与结构修改算法的共享内存式并行。采用的方法是流水线并行模式:将分解一行拆分成行间消元与自消元,形成多个小任务,当一个小任务所需数据准备好时即开始计算,而不必等前面所有行都分解完才开始。该并行算法在常见双核或四核CPU上的加速比通常可以达到3左右。  最后,本文介绍了结构修改算法的几个应用实例,提出了适用性强的计算流程方案,为本算法在实际工程中的应用与推广提供了参考与便利。  综合软件和硬件两方面的优化成果,本文结构局部修改并行算法的运算速度将达到重新完整计算的5到30倍,计算效率的提升非常显著。与现有的结构修改直接解法相比,本文算法还具有以下优点:更适合处理大规模问题的大量修改,额外占用的内存极少,适合依次连续修改。更为重要的是,本文算法与现有的稀疏矩阵有限元方程直接求解方案完全兼容,实现简单,因此有望成为结构局部修改的标准算法。
其他文献
本文的研究内容是太原市科技局大学生创新创业项目“刮板输送机驱动系统”(项目编号:110148032)的重要组成部分,是为了进一步提高开关磁阻电机的动态性能以及带载能力,增加电机的制动功能而提出的。本文结合SRD系统的实际需求,研究了电机在CCC控制方式下的单双相励磁运行模式,提出了电压PWM控制方式下的开通角、关断角优化方法,增加了电机的制动功能。主要研究内容如下:在查阅大量相关资料并进行多次实地
学位
随着电力系统的飞速发展和电力市场的逐步实现,大电网的稳定问题(尤其是暂态稳定问题)变得更加复杂。这些变化使得在电力系统中应用安全稳定控制装置(系统)已经成为一种趋势,因此,安全稳定控制装置(系统)的研究与开发也得到了长足的发展。在安全稳定控制装置(系统)开发的过程中,不可避免地遗留下一些错误和缺陷,如果不及时将其找出与改正,这个产品就可能会在现场中工作失常,甚至会给电力系统带来隐患。这使得对安全稳
涉及多相流的工程技术已广泛渗入我国国民经济各部门。以均速管流量计为代表的压差式计量仪器,结构简单、计量精度较高,已被应用到两相流动的测量中。但目前与均速管流量计相应的多相流动的理论、实验及数值研究工作,还相对较为薄弱。实际应用过程中,均速管流量计在多相流的测量中,堵塞还是影响其高效广泛应用的重大问题。本课题组设计的新型双体均速管流量计,在威力巴原型的基础上,引入局部文丘里效应,使有效输出压差上有了
随着全国联网的不断进展,实时功角的测量和应用越来越受到关注。本文的工作是研发一种基于GPS的新型结构的广域动态监测系统,主要侧重于后台监测软件的开发工作。通过提出系统的功能需求,制定系统的整体设计方案,选择工业控制软件平台实现后台监测和数据处理,从而构成了一个架构灵活配置、易于扩展、分层分布式的电网广域动态监测系统。并利用具有国际水准的RTDS数模仿真系统对整个监测系统进行了双端模型和三端环网模型
学位
学位
由于超声波电机原理和结构的新颖性及拥有的直接驱动、结构多样性、无电磁干扰及不受电磁干扰影响的特性,虽然目前还不能完全替代传统的电磁电机,但在一些特定的应用场合取代电磁电机是未来微特电机的发展趋势。在电机的研究过程中,必要的数学模型和分析方法是研究系统控制特性和性能分析的基础,是电机结构参数优化设计的有效工具。本文对环形行波型超声波电机的工作原理进行数学推导,将电机的定转子接触简化为二维弹性接触,推
德国媒体近日刊登了一篇题为的报道,该报道援引美国航空航天专家的观点,指出亚洲近期接连发生的海底大地震是由天王星的异常而引起的.
学位