结构的多机并行分析I——结构的并行有限元方法

本文详细分析和讨论了结构分析的并行有限元方法—并行预处理共轭梯度法(以下简称PPCG法)。着重讨论了基于自带存储器的多处理机系统的并行预处理算法问题,并由此提出了两种PPCG法:PPCG1和PPCG2法。这两种方法适用于以单道剖分(one-way dissection)的子结构法为基础的并行分析。由于这种剖分法产生的结构刚度矩阵具有箭头形状,可独立地消除各子结构的内部自由度,并且不会在刚度矩阵中产生新的非零元素,因此很适合具有较多处理机的并行机系统对复杂结构进行的并行分析。     

陀螺仪漂移特性的小波分析

在简述小波分析中Mallat算法的基础上,用其实现陀螺仪漂移特性的小波分析,并对有关的工程问题进行了讨论,给出了某型陀螺实测数据的应用实例,证明了此方法的有效性。对今后的应用提出了一些看法。     

GPS／捷联惯性组合导航系统的一种半实时算法

本文讨论了ＧＰＳ／捷联惯性组合导航系统的原理和特点，针对其理论仿真和工程实现中的衔接问题，提出了一种半实时算法。该算法既继承了理论仿真方便、直接的优点；又充分利用了实测数据。系统调试成功后，可以直接将全套软件转换到导航计算机中。文中给出了半实时算法的原理和方案，讨论了组合卡尔曼滤波器的实现和ＧＰＳ仿真器、航迹仿真器、惯性器件等的设置。最后给出了有关结论。     

无网格径向点插值法在拓扑优化中的应用

论文基于无网格径向点插值法(RPIM)对连续体结构进行拓扑优化设计.以高斯点的相对密度为设计变量,以结构的柔度最小化为目标函数,利用带惩罚的各向同性固体材料模型(SIMP)和优化准则法,建立了设计变量的迭代格式.利用灵敏度过滤技术有效地消除了点状棋盘格现象.给出了相应的计算流程,并用Fortran程序语言实现其算法.算例证明,应用无网格径向点插值法能够有效地对连续体结构进行拓扑优化设计.     

结构非概率可靠性指标的求解方法

在结构的非概率可靠性方法中,结构的可靠性是用非概率可靠性指标来度量的。本文研究了区间变量的一般运算规则,给出了非概率可靠性指标的三种求解方法：定义法、转换法和优化法。给出了定义法的具体描述和实现方法,并给出了转换法的一般求解步骤和一种可通用的区间优化算法的实现过程。三种方法的给出,基本上解决了结构可靠性分析中非概率可靠性指标的求解问题。实例计算表明所提方法是有效和可行的。     

以u—U形式方程为基础的结构—饱和土壤相互作用的固结有限元分析

提出用u-U形式方程对饱和土固结问题进行分析,对于结构则采取多重子结构技术实现向饱和土交界面上的凝聚,从而建立了结构与饱和土相互作用的固结控制方程。文中给出了方程的求解方案。文末给出了数值算例。本文的工作已在DIASS程序系统中实现。     

拟协调分层薄壳单元在圆柱薄壳结构弹-塑性分析中的应用

本文基于拟协调元的基本原理,引入薄壳分层子单元的概念,提出了用于圆柱薄壳结构弹-塑性分析的三角形和矩形拟协调分层圆柱薄壳单元的模型,给出了第k层薄壳分层子单元的弹-塑性矩阵和物理矩阵[D~(k)]的表达式,并最后给出了整个薄壳单元的物理矩阵和刚度矩阵的表达式。 根据上述原理和方法用Fortran语言编制了圆柱薄壳结构弹-塑性有限元分析程序,实例计算与实验结果较为吻合,从而说明用本方法进行圆柱薄壳结构的弹-塑性分析是可靠而行之有效的。     

推进－迭代算法计算效率的研究

本文对求解三维定常超音速动性流场的一次空间推进，在每一个推进站沿伪时间层局部迭代的推进－迭代算法作了进一步的研究．在每一推进站（侧向平面）沿伪时间层局部迭代时，给出了四种不同的隐式迭代方法，即沿侧面两个方向（法向和周向）全用隐式；法向隐式而周向采用Ｇａｕｓｓ－Ｓｉｌｄｌｅ来回扫描迭代；法向隐式而周向显式及以系数矩阵谱半径代替系数矩阵的简化标量隐式算法．用这四种算法模拟了三维球锥黏性绕流，给出了四种不同算法的计算效率和收敛特性比较．     

树形多体系统动力学的隐式数值算法

研究了树形多体系统动力学的隐式算法．用矩阵形式给出了多体系统的正则方程及其右端函数的Ｊａｃｏｂｉ矩阵，并给出该矩阵的分块算法和对角隐式Ｒｕｎｇｅ－Ｋｕｔａ法（ＤＩＲＫＭ）以及隐式辛Ｒｕｎｇｅ－Ｋｕｔａ法（ＩＳＲＫＭ）．该算法便于编程计算，能提高计算效率，保持长期计算的稳定性．并用算例说明该算法的有效性     

轴对称载荷作用下旋转壳的近似屈服条件

文献[1—3]已给出了轴对称旋转壳以四个广义内力n_1,n_2,m_1,m_2表达的精确屈服条件.本文对服从最大切应力准则的壳,给出了以膜力强度和力矩强度表达的静力屈服面n~2+m=1和机动屈服面2(n-1/2)~2十m=1,机动屈服面不超过静力屈服面的1.31倍.本文还应用这两个屈服面计算了受均布法向载荷作用下的简支和固定边球壳的极限载荷,改进了文献[4]给出的上下界.从给果中发现了球壳的极限载荷基本上只与一个壳体参数有关.最后还给出了周边简支可移的扁球壳极限载荷的上下界.     

位移法解受压固支正交异性矩形分层的屈曲和后屈曲

在足够大的压力作用下,层合板中的分层会产生屈曲和后屈曲。本文用位移法求解了受压固支正交异性分层的屈曲和后屈曲问题,并给出了有关的数值结果。     

固体力学应力分析中的有限域法

本文给出了用有限域法进行了和学应力分析的一般原理，并与目前广泛采用的有限元不做了比较。文中的研究表明：有限域法与有限元法极为相似，前者使用单位位移加权，后者从虚位移原理出发；两种方法的实施过程也有共同点，它们都进行网络部分及逼近。     

基于开放式结构有限元系统SiPESC.FEMS的动力时程分析通用算法构架

基于开放式结构有限元分析软件系统SiPESC.FEMS,针对动力时程分析算法共性,采用C++面向对象程序设计方法和软件设计模式,研发了一种结构动力分析通用算法构架。构架的核心思想是算法与数据模型相分离,从而实现算法通用性,整个构架由四个基本类及子类构成。本文重点阐述了基本类的抽象过程和通用接口的设计思想,给出了利用插件技术实现算法构架的步骤。利用该构架已实现了Newmark法、Wilson-θ法、中心差分法及改进中心差分法,并进行数值验证。研究工作表明,算法构架适用于通用时程积分算法,可方便地进行动态扩展,具备良好的开放性和重用性。     

基于有限元分析的特征值反问题求解的逆摄动方法

本文研究特征值反问题的求解方法,根据广义特征值反问题理论和有限元法的特点,以转子系统平面梁单元有限元模型结构分析的特征值反问题求解为例,给出一种新的逆摄动方法,给出了本逆摄动法较完整的理论基础,给出了其逆摄动参数的显式计算公式及相应的取值方法,本逆摄动法也可推广到其他单元类型的有限元模型特征值反问题的求解。     

基于ICCP算法及其推广的重力定位

基于ICCP算法的重力匹配定位可以用于限制推算定位随航行时间增长的位置误差。给出了ICCP算法的设计思想,同时针对算法的假设前提进行了推广,使算法能够在考虑重力传感器测量数据存在误差的情况下,实现推位航法的误差校正。仿真结果证明这种推广具有较好的定位精度,能够满足AUV的导航要求,对实现AUV的自主无源导航有重要意义。     

关于《再论约束最小二乘法》的注记

约束最小二乘法在计算力学和工程问题的理论研究和实际应用中都起到重要的作用。本文借助矩阵广义奇异值分解的方法，研究了一般矩阵的约束最小二乘问题，给出了一般矩阵约束最小二乘问题的算法，推广了文献[1]的结果。     

新型陀螺仪数字式信号提取图形的研究

介绍了用在新型陀螺仪上一种赤道三角图形光电信号读取方案，给出了三种三角图形的实现方法，就镀层金属法和高精密机床刻制技术进行了研究和试验，分析了由于边缘效应和伪信号带来的测角误差，就两种方法实现的图形进行了比较，分析了各自对信号检测带来的影响。建立了刻线误差数学模型且进行了仿真。为了消除边缘效应和伪信号影响，提出了一种新的信号提取图形法一组合图形法，并给出了它的实现模式和试验结果。最后分析了由于图形的存在而使得转子表面不规则对陀螺性能产生的负影响，可以用这些参数来预测陀螺漂移性能。     

双参考光全息—散斑干涉术的原理和实验

本文用统计光学方法对双参考光全息——散斑干涉术进行了详尽的理论分析,给出了全息和散斑干涉场区域平均光强分布与光学系统主要参数及三维变形场各分量之间关系的解析表达式,进而讨论了离面位移和面内位移测量的上、下限,最后还给出了有关的实验结果。     

地磁导航系统中等值线的提取及应用

地磁导航系统采用等值线最近点迭代算法达到最终匹配定位的目的,等值线最近点迭代算法是在等值线的基础上实现的,因此,等值线的提取是实现地磁导航的“必经之路”。文中介绍了常用的两种等值线提取方法,即线性插值方法和等值线延伸寻找最近点法,其中前者的计算量大,实时匹配效果差,而后者能在地磁匹配过程中快速提取匹配需要的等值线。此外,文中还给出了地磁等值线图层在基于电子海图显示信息系统的地磁导航系统中的生成和添加方法。     

某型双立尾战斗机的网格生成及流场解

以某型双立尾战斗机为例，给出了实际复杂型号飞行器结构网格生成的一种分块方案及实现过程。用一种矢性三次多项式插值的方法生成相邻块公共交界面的网格，用求解椭圆型方程的方法生成块内空间网格。在所生成的双立尾战斗机分块结构网格中分区求解Ｅｕｌｅｒ方程，获得了合理的全机气动力系数、翼面压力分布和绕流图谱结果