基于减基概念的凸集结构可靠性分析方法

针对椭球凸集参数域结构的可靠性分析问题,提出了一种基于减基概念的快速求解方法。首先,将椭球参数域进行坐标正交变换,获得标准的椭球域及其相应的矩形域,在矩形域采样且通过坐标逆向变换获得原椭球参数域的样本参数点集,并以此构建结构的减基空间及其相应的减基算法;随后,在标准椭球域产生均匀的等概率抽样点,并通过坐标逆向变换和相应的减基算法进行蒙特卡洛减基模拟来分析结构的可靠度及其可靠域。由于是在低维的逼近空间中进行椭球参数域结构的位移向量解分析,故而较之有限元法能够获得较高的计算效率。算例测试验证了本文方法的有效性。     

一种动力响应混合快速分析方法及其在功能梯度材料板瞬态响应分析中的应用

提出了一种把减基法嵌套在Newmark法中快速求解动力问题的混合算法.这种算法把Newmark法作为动力问题的基本求解器.而在嵌套减基法时,把参数样本空间对应的解集做成一个训练空间,采用贪婪算法进行减基空间基向量的自适应选取,进而由减基空间得到减缩系统.通过求解减缩系统对应的Newmark增量式来达到快速求解动力问题的目的.最后采用了一个无限大的功能梯度材料板的瞬态响应分析来验证文中方法的高效性和精确性.     

一种求解大型结构动力响应的快速计算方法

针对大型工程结构动力响应求解效率较低的问题,提出了一种基于减基法的快速求解方法。该方法对动力学方程进行时间域积分构建减基空间,利用Galerkin映射向减基空间进行投影得到减缩方程,利用减缩方程快速求得原系统的逼近解,从而大大提高了动力学方程的求解效率。该方法还考虑了Galerkin映射下的奇异性计算,文中算例验证了该...     

一种分析区间参数结构可靠性的减基快速计算方法

将减基法与蒙特卡洛模拟结合,提出了一种快速计算区间不确定结构可靠性的方法。该方法分为离线和在线计算两个阶段,离线阶段利用减基法建立减基空间,进而形成减缩模型;而在线阶段则将减缩模型融入蒙特卡洛方法,进而快速求解区间失效概率及相应的区间参数失效域。该方法在减基空间而非有限元空间中分析区间结构的可靠性,减少了求解时间,提高了计算效率,通过算例验证了该方法的可行性和有效性。     

GPS/SINS系统空中对准姿态角误差可观测性研究

提出一种新的判断系统可观测性和可观测度的方法,详细分析了机体各种运动对系统姿态角误差可观测性和可观测度的影响,并把该方法应用于组合导航系统的可观测性和可观测度的研究中。该方法利用了误差状态的最小二乘估计均方误差阵的特征值和特征向量,能判断系统的可观测度,避免进行卡尔曼滤波计算求协方差阵。仿真结果表明该方法简单、快速、有效。     

计算特征向量摄动量的混合基展开法

在结构修改和模型校正中，模态展开法是计算特征向量摄动量的常用方法之一，但当高阶模态被截断时，它会带来很大的截断误差。本文利用已知的有限阶模态，构造了Ｎ维欧氏空间的一个新基－混合基，并将特征向量的摄动量在新基上展开来计算特征向量的一、二阶 摄动量。该方法使得不管截模态个数的多少，其精度总与全模态展开法相同，且计算量都远少于全模记展开法；与改进的部分民开法相比，本方法不要求所截留的模态边连续的低阶模态     

智能结构密频系统振动控制及其摄动分析

提出一种智能结构密频控制规律的设计方法。其主要工作为：将密频子空间转化为重频子空间,把密频系统作为重频系统上的摄动;为了解决重频密频子空间相对应的特征向量选取的敏感性问题,通过摄动分析得到与重频密频子空间相对应特征向量的线性组合并利用闭环系统特征值极点配置的方法得到重频密频系统的振动规律;把所设计的振动控制规律作用于原系统和摄动系统上,讨论了结构参数改变后对系统的动态特性的影响;最后通过算例证明该方法的有效性。     

基于Lanczos算法的模态重分析法及其在车身结构设计中的应用

模态重分析是指在结构修改之后不需要重新求解广义特征值方程,仅需要根据初始计算结果对修改后的问题进行求解,并能够在保证精度的前提下,提高计算速度。随着结构复杂度和修正量的增加,传统重分析方法的求解精度和稳定性随之下降。为此,利用初始结构模态分析结果,结合Lanczos算法和投影技术,采用缩减基方法求解修改结构的特征值和特征向量,使其同时具备了Lanczos向量快速收敛的优点和基于全局近似的缩减基向量的高精度。为了验证该方法的性能和准确性,对本文方法基于扩展基向量和瑞利-里兹分析的模态重分析法以及改进的单步摄动瑞利商逆迭代法进行了测试。测试结果表明,该方法具有最高的计算精度。同时,将该方法成功用于车架和车门的前期设计中,计算结果表明,该方法具备处理计算规模大、拓扑修改变化量大的结构分析问题的潜力。     

基于改进残余力向量法的桁架结构损伤诊断

提出一种基于改进残余力向量法的桁架结构损伤诊断方法. 先由灵敏度分析, 求出结构刚度联系矩阵,再由刚度联系矩阵将损伤后的刚度摄动矩阵展开成对角矩阵,代入 残余力向量方程,得到由刚度联系矩阵表示的新的残余力向量方程,此方程可以直接求解, 即可诊断出桁架结构的损伤杆件及其损伤程度. 对于实测中难以获得完备振型的情况,采用 模态扩阶的方法来获得完备的测试振型. 最后以一桁架结构进行数值仿真分析,证实了该方 法的有效性.     

基于双向重投影的双目视觉里程计

针对视觉里程计中的特征点漂移与累积误差问题,提出一种基于双向重投影的双目视觉里程计方法。为降低特征点匹配的高错误率及改善匹配效率,提出一种改进的随机采样一致(RANSAC)算法。该算法依据特征点寿命长短、环形匹配、优质匹配等方法选取优质特征点以降低特征点样本容量,从而改善RANSAC算法的效率;为抑制视觉里程计累积误差、提升位姿估计精度,提出一种双向重投影的位姿估计算法,该方法通过解算下一帧图像中的最小化重投影误差得到相机位姿以更新空间点集,更新后的空间点重新投影至当前帧图像以提升位姿解算精度。最后利用KITTI数据集进行了仿真。结果表明该算法的150 m内定位误差在X方向小于1 m,在Z方向小于0.9 m;相对于ORB-SLAM2算法,定位精度提升了约20%,系统的稳定性更好,有效削弱了累积误差,验证了该算法的正确性与可行性。     

基于减基法的层合板瞬态响应快速分析方法

用混合数值法分析复合材料层合板结构的瞬态响应时,在波数域中引入减基法、构造减基空间,把原大型特征值问题映射到减基空间进行降阶,在保证精度的前提下快速近似求解原大型特征值问题,从而快速得到结构的时域响应. 算例表明,引入减基法后的混合数值法能更有效地解决结构瞬态响应问题. 该方法可以推广用于快速解决大型结构动力优化设计和结构动力反问题等结构动力分析问题.      

粘弹性阻尼结构系统振动特性的迭代分析

含粘弹性阻尼材料结构系统的振动特性分析是相当复杂的,求解比较困难,本文给出了一种有效的迭代解法。文中首先定义了广义Rayleigh商,指出求广义Rayleigh商的驻值问题等价于原特征问题的求解;由近似估计得到的特征向量,广义Rayleigh商给出一个较好的特征值近似估计。据此构造了迭代方法,并成功地计算了一些例子。     

大规模边界元模态分析的高效数值方法

随着大规模快速边界元计算技术的发展,在复杂结构的动态设计、振动与噪声分析中愈来愈多地采用边界元法,因此求解大规模边界元特征值问题、进行复杂结构和声场模态分析,成为工程应用中一个十分重要,但却极具挑战性的课题,目前国际上还没有十分有效的数值方法.本文针对边界元法中典型的非线性特征值问题,提出了一种通用、高效的数值解法,称为基于预解矩阵采样的Rayleigh-Ritz投影法,记为RSRR.首先,通过求解一系列频域边界元问题来构造特征向量搜索空间,进而可以采用Rayleigh-Ritz投影,将原问题转化为一个可以采用现有方法求解的小规模缩减特征值问题;其次,为了降低Rayleigh-Ritz投影过程的计算量,基于解析函数的Cauchy积分公式,构造了边界元系数矩阵的插值近似方法,以及缩减特征值问题系数矩阵的快速计算方法,给出了插值项数的估计策略;最后,将RSRR与声学快速边界元法结合,应用于大规模吸声结构的复模态分析.数值算例表明,RSRR方法能够可靠地求出给定频段内的全部特征值和特征向量,具有计算效率高、精度高、通用等优点.     

Computational study of the axial instability of rimming flow using Arnoldi method

Axial instability of rimming flow has been investigated by solving a linear generalized eigenvalue problem that governs the evolution of perturbations of two‐dimensional base flow. Using the Galerkin finite element method, full Navier–Stokes equations were solved to calculate base flow and this base flow was perturbed with three‐dimensional disturbances. The generalized eigenproblem formulated from these equations was solved by the implicitly restarted Arnoldi method using shift‐invert technique. This study presents instability curves to identify the critical wavelength of the neutral mode and the critical β, which measures the importance of gravity relative to viscosity. The axial instability of rimming flow is examined and three‐dimensional flow was reconstructed by using eigenvector and growth rate at a critical wave number. The critical β value in the axial instability analysis was observed to be comparable to the onset β value of the transition between the bump and the homogeneous film state in 2‐D base flow calculations. Copyright © 2006 John Wiley & Sons, Ltd.     

基于交叉降阶取样的改进响应面方法

针对经典响应面法和加权响应面法存在的计算精度不高和稳定性不强的问题,提出一种新的样本点选择策略,即交叉降阶思想,并基于此发展出一种新的改进响应面方法。该方法以均值点加上0.5倍的标准差作为初始迭代点。在Bucher试验设计的基础之上,将n维坐标系降阶为n个一维坐标系。然后,依据样本点的重要性,筛选出n+1个靠近真实失效面的优秀样本点,并用这些样本点去拟合线性响应面函数。最后依据样本中心间的距离与上一次迭代得到的设计点的范数之商作为其收敛准则。通过5个数值算例与工程算例,验证了交叉降阶响应面法具有较高的精度、效率、稳定性和一定的工程适用性。     

考虑尺寸效应的双材料轴对称界面端应力奇异性

提出了一种分析双材料轴对称界面端的应力奇异行为的特征值法.基于弹性力学空间轴对称问题的基本方程和一阶近似假设,利用分离变量形式的位移函数和无网格算法,导出了关于应力奇异性指数的离散形式的奇异性特征方程.由奇异性特征方程的特征值和特征向量,即可确定应力奇异性指数、位移角函数和应力角函数.数值求解了纤维/基体轴对称界面端模型的奇异性特征方程, 结果表明:尺寸效应参数δ(奇异点与轴对称轴的距离和应力奇异性支配区域大小的比值)影响着应力奇异性的强弱与阶次, 准一阶近似解析解只是δ>>1时的一个特例.