带位移偶次方项非线性随机振子在位移参数激励下的概率解

论文用指数多项式闭合或EPC(Exponential Polynomial Closure)法分析了具非零均值响应的带位移偶次方项非线性随机振子在参数激励下响应的概率密度函数解.给出了求解过程并通过算例分析验证了指数多项式闭合法在此情况下的有效性.数值结果显示,指数多项式闭合法得到的响应概率密度结果与蒙特卡洛模拟的结果符合较好,尤其是在对系统可靠性分析起主要作用的概率密度函数尾部区域符合很好.     

非高斯激励下带附加质量的悬臂梁响应研究

大多数悬臂梁的研究都是基于确定性振动分析或将外部激励假定为高斯过程，鲜少涉及非高斯激励。针对非高斯激励下带附加质量的悬臂梁的响应特征进行研究。首先，通过Lagrange方程建立考虑曲率非线性的悬臂梁振动方程；然后，通过指数多项式闭合法求解结构响应的概率密度函数，并分析了平均到达速率、激励强度及附加质量对结构响应的影响；同时分析了非平稳非高斯激励下非线性结构参数对响应的影响。计算结果表明，指数多项式闭合法与模拟解吻合很好，平均到达速率越大，结构响应的概率密度函数曲线收缩越明显。     

随机激励下滞迟系统的稳态响应闭合解

滞迟系统属于一类典型的强非线性系统,滞迟力不仅取决于系统的瞬时变形,还与变形历程有关.虽然滞迟系统的随机振动问题已被广泛研究,但至今尚未得到滞迟系统随机响应概率密度函数的精确闭合解.本文运用迭代加权残值法获得了高斯白噪声激励下Bouc-Wen滞迟系统稳态响应概率密度函数的近似闭合解.首先,运用等效线性化法求出系统的稳态高斯概率密度函数;然后以此构造权函数,应用加权残值法求得了系统指数多项式形式的非高斯概率密度函数;最后引入迭代的过程,逐步优化权函数,提高计算所得结果的精度.以随机地震激励下钢纤维陶粒混凝土结构的稳态响应作为算例,其中Bouc-Wen模型的参数是基于拟静力学试验数据,并应用最小二乘法辨识获得.与Monte Carlo模拟结果相比,等效线性化法得到的结果精度较差;由加权残值法得到的结果能够表现出非线性特征,但其精度依然无法令人满意;采用迭代加权残值法得到的近似闭合解与Monte Carlo模拟的结果吻合非常好;对于较强随机激励情形,采用渐进迭代加权残值法具有较高的求解效率,所获得的理论解析解具有较高的精度.结果表明,所获得的近似闭合解不仅对于土木工程领域具有重要的实际应用价值,而且还可作为检验其他非线性系统随机响应预测方法的精度的标准.     

求大规模非线性随机动力系统概率密度函数解的子空间法

论文给出了一个求大规模非线性随机动力系统响应概率密度函数解的新方法,称之为子空间法.考察了此方法在求解大规模带位移项参数激励非线性随机动力系统响应概率密度函数的有效性.这里的概率密度函数解由Fokker-Planck-Kolmogorov(FPK)方程控制.该方法是基于将非线性随机动力系统状态变量空间分成两个子空间,然后在其中一个子空间上对FPK方程进行积分,采取一定措施后得到低维的FPK方程.该低维的FPK方程的维数可以人为确定,也可以取为二维,从而可以用现有的求解低维FPK方程的方法求得所需的概率密度函数解.文中给出了算例,用数值结果验证了子空间法的有效性.论文是采用作者曾提出的指数多项式闭合(EPC)法求解由子空间法降维的FPK方程.     

轴弯曲与不平衡柔性转子共振稳态响应随机分析

系统参数随机性是影响准确预测转子动力学行为的重要因素。为此,研究了具有初始轴弯曲的不平衡柔性转子在故障参数随机情况下的不确定性分析问题。首先,基于转子动力学梁元有限元理论,推导了不平衡和轴弯曲故障共同作用下的系统稳态动力学方程,并以轴心轨迹长半轴作为关键响应量建立了柔性转子共振稳态响应与输入参数间的模型函数;其次,联合广义多项式混沌展开、留一法交叉验证以及最小角回归技术实现了柔性转子共振稳态响应的非嵌入式自适应稀疏多项式混沌展开,并与基于普通最小二乘法的多项式混沌展开和基于Monte Carlo仿真的结果作了对比分析,验证了自适应稀疏展开方案的有效性、精度和效率;最后,以构建的自适应稀疏多项式混沌展开式作为近似模型,重点分析了转子圆盘处一阶共振稳态响应的随机特性,并基于Sobol指标获得了响应对各故障参数的全局灵敏度指标。     

结构可靠性分析的支持向量机响应面法

针对隐式极限状态可靠性分析问题,提出了一种支持向量机响应面法,该方法采用了与经典响应面法类似的迭代思想,由支持向量回归机替代经典响应面法中的固定多项式函数来构建响应面,并结合一次二阶矩法形成迭代过程.在此基础上,还对训练样本提出了一种改进的选取方法,从而进一步提高方法的效率.文中将所提方法与多种经典可靠性分析方法的计算结果对比分析,改进的支持向量机响应面法精度较高,调用结构分析程序的次数最少.     

非嵌入式多项式混沌法在爆轰产物JWL参数评估中的应用

介绍了非嵌入多项式混沌法的数学模型,给出了非嵌入式多项式混沌法进行不确定度量化的主要步骤。采用此方法研究了平面、散心爆轰问题数值模拟中, JWL模型参数R₁、R₂服从均匀分布的随机变量时所引起的爆轰过程计算结果的不确度性,着重分析了爆轰传播过程中压力与密度的统计特性。研究结果表明,非嵌入式多项式混沌法可以为模型输入参数不确定性的传播对输出结果响应量的影响建立一种有效不确定度评估方法,为使用JWL模型时选取参数提供参考。     

随机结构有限元分析的递推求解方法的改进

将随机结构有限元分析的递推求解方法和伽辽金投影方法相结合,提出了求解随机静力响应的改进的递推求解方法。利用随机收敛的非正交多项式展开表示由于材料、外部荷载或构件几何尺寸的随机性导致的结构随机响应。采用递推求解方法得到响应多项式展开的初始系数,并运用定义的数学算子显式地表达出来。然后,通过定义修正系数,应用伽辽金方法对随机力平衡方程在非正交多项式基上进行投影,得到了和响应展开阶次个数相同的确定的有限元方程,并进行求解得到了修正系数。数值算例表明,通过对递推求解方法中响应表达式系数的修正,以很小的计算代价较大地提高了随机响应的计算精度;与基于正交多项式展开的随机有限元方法相比,在精度相当的前提下新方法耗费的计算时间大大降低。     

随机桁架结构几何非线性问题的混合摄动-伽辽金法求解

提出应用混合摄动-伽辽金法求解随机桁架结构的几何非线性问题.将含位移项的随机割线弹性模量以及随机响应表示为幂多项式展开,利用高阶摄动方法确定随机结构几何非线性响应的幂多项式展开的各项系数.将随机响应的各阶摄动项假定为伽辽金试函数,运用伽辽金投影对试函数系数进行求解,从而得到随机桁架结构几何非线性响应的显式表达式.同已有的随机伽辽金法相比,本文所给的试函数由摄动解的线性组合而成,在求解非线性问题时,试函数的获取具有自适应性.数值算例结果表明,对于具有不同概率分布的多随机变量问题,本文方法无需对随机变量的概率分布形式进行转换,避免了转换误差,因而比同阶的广义正交多项式方法 (generalized polynomial chaos, GPC)计算精度高.同时,在结果精度相当时,和GPC方法相比,本文方法得到的试函数系数的非线性方程维度不大,方程的求解工作量小且更易求解.当随机量涨落较大时,混合摄动-伽辽金法计算所得的结构响应的各阶统计矩比高阶摄动法所得结果更逼近于蒙特卡洛模拟结果,显示了该方法对几何非线性随机问题求解的有效性.     

FPK方程的近似闭合解

讨论了FPK方程的近似闭合解问题。假设FPK方程的解具有指数多项式的形式，利用比较系数的方法确定其中待定的常数。计算表明，本方法适用于强非线性系统，在特殊情形下还能求出原方程的精确解。     

随机桁架结构几何非线性问题的混合摄动-伽辽金法求解

提出应用混合摄动$\!$-$\!$-$\!$伽辽金法求解随机桁架结构的几何非线性问题.将含位移项的随机割线弹性模量以及随机响应表示为幂多项式展开,利用高阶摄动方法确定随机结构几何非线性响应的幂多项式展开的各项系数.将随机响应的各阶摄动项假定为伽辽金试函数,运用伽辽金投影对试函数系数进行求解,从而得到随机桁架结构几何非线性响应的显式表达式.同已有的随机伽辽金法相比,本文所给的试函数由摄动解的线性组合而成,在求解非线性问题时,试函数的获取具有自适应性.数值算例结果表明,对于具有不同概率分布的多随机变量问题,本文方法无需对随机变量的概率分布形式进行转换,避免了转换误差,因而比同阶的广义正交多项式方法(generalizedpolynomial chaos, GPC)计算精度高.同时,在结果精度相当时,和GPC方法相比,本文方法得到的试函数系数的非线性方程维度不大,方程的求解工作量小且更易求解.当随机量涨落较大时,混合摄动$\!$-$\!$-$\!$伽辽金法计算所得的结构响应的各阶统计矩比高阶摄动法所得结果更逼近于蒙特卡洛模拟结果,显示了该方法对几何非线性随机问题求解的有效性.      

斜拉桥拉索在轴向窄带随机激励下的振动响应

导出了拉索在考虑垂度以及索张力沿索长变化时的参激随机微分方程,进一步给出了预测拉索在窄带随机激励下响应的近似理论解------用统计矩截断法求解矩方程,获得高斯闭合解和一阶非高斯闭合解. 以南京长江二桥约330米长的A20拉索为研究对象,对以上高斯闭合解和一阶非高斯闭合解进一步进行数值求解以获得拉索的响应,并采用Monte-Carlo数值方法对求解进行验证. 分析了拉索振动的一般特征,特别分析了激励中心频率和拉索频率比为1和2时的响应随激励带宽的变化特征,得到了一些新的结论.      

结构可靠度分析的改进BP神经网络响应面法

对功能函数不能明确表达即具有隐式功能函数的问题进行可靠度分析,常采用响应面法。其中二次多项式响应面法应用较为广泛,采用与此方法相同的思路,人们提出了BP神经网络响应面法。在此基础之上,笔者提出了改进的BP神经网络响应面法。对以上方法计算效果,通过算例进行了对比分析,其中改进BP神经网络响应面法计算精度较好,进行有限元分析次数较少。该方法用于大型复杂结构的可靠度分析,可相应提高工作效率和解题质量,具有一定实际应用价值。     

任意变厚度的旋转扁薄壳非线性稳定的幂函数解法

以幂函数为试函数,两次使用配点法成功地分离了耦合的大挠度方程,从而导出变厚度旋转扁壳非线性稳定的计算式。支座可以是弹性的。本文给出了均布或多项式分布荷载作用下,线性或多项式型变厚度的圆锥壳、球壳、余弦壳或四次多项式型旋转壳的上、下临界荷载。均布荷载作用下指数型变厚度球壳的上临界荷载同其他方法的结果作了比较。用配点法编写的程序具有收敛范围大、精度高、通用性强和计算时间少的优点。     

随机结构有限元分析的递推求解方法

提出了一种新的谱随机有限元分析方法——递推求解方法。该方法将随机结构的随机响应表示成非正交多项式展式,建立了和摄动法类似的一系列确定的递推方程,并通过确定性有限元方法对这些递推方程进行静力问题求解。算例表明,当随机量出现较大涨落时,计算结果相对于传统摄动法有不小的改进。     

多自由度强非线性耦合参激系统随机响应计算方法

本文采用非高斯矩闭合方法计算了一个含参数激振的二自由度强非线性耦合系统受随机激励的响应方差。用逐次消除法消除方程中的惯性耦合后推导出直到6阶的矩方程,采用中心累积量截断技术计算了白噪音激励下的响应方差,并与时域直接积分和数字模拟方法的结果进行比较,取得了一些有意义的结论。     

设防烈度下非比例阻尼结构地震随机响应峰值区间估计方法

建立与抗震规范设防水准相一致的地震动随机模型,给出了非比例阻尼结构体系地震反应运动方程,推导了非比例阻尼结构体系分析的复振型时域闭合解和地震随机响应的闭合解.通过随机响应的峰值因子和变异系数获得了非比例阻尼结构体系的地震随机响应峰值均值和方差,提出了设防烈度下非比例阻尼结构体系随机响应峰值区间估计方法.计算31条实际地震记录的复振型时域闭合解,并与地震随机响应的闭合解进行了比较,算例结果表明随机响应峰值区间估计方法的合理性和有效性.地震响应区间分析方法将是工程抗震设计的一个重要方法.     

简谐力激励下结构拓扑优化与频率影响分析

研究了简谐力激励下以结构指定位置稳态阶段位移响应幅值为目标函数、结构体积为约束的拓扑优化设计问题. 通过在频域上使用模态叠加法求解简谐力激励下的位移响应, 分析了激励频率和作用方向对位移响应幅值及其优化结果的影响.引入材料属性的多项式插值惩罚模型, 有效消除了动力学拓扑优化局部模态现象.分析了高频激励下位移响应幅值拓扑优化存在的稳定性差、结构不连续等问题, 并通过引入附加静位移约束, 获得了清晰合理的结构形式.理论分析和算例结果揭示了位移响应幅值优化过程中结构模态的变化规律, 验证了该拓扑优化模型的有效性.     

高斯过程响应面法研究

新型的高斯过程响应面法具有灵活性好、精度高和能量化不确定性等优点,但某些关键技术细节还存在争议。本文简要介绍了高斯过程响应面法的基本理论,结合算例详细讨论了关键技术细节的处理方法;提出了先验期望函数选用低阶多项式、相关函数选用高斯形式、粗糙度参数用边缘后验众数法固定、模型有效性必须经过验算点检验等应用建议。     

在线振动监测与故障诊断的一种新途径

文中由应力强度因子导出含裂纹单元的刚度矩阵,从而提出了一种有限元法.将该法用于单边裂纹悬臂梁,计算了应变响应.用时域法识别了模态参数,并考虑了裂纹闭合的影响.分析了应变响应对故障的敏感度,计算结果表明,应变响应相对于位移响应而言,对故障更敏感.本文还提出了一种判别函数,该函数对于损伤检测,比现有几种判别函数更有效.最后,给出了一种在线振动监测与故障诊断新方法.