粘弹性随机分析的虚功原理及随机有限元

利用增量法处理粘弹性本构关系中的遗传积分，将粘弹性材料的随机性、结构几何形状的随机性、外载荷的随机性引入虚功方程，应用摄动方法，研究了粘弹性随机分析的虚功原理和粘弹性随机有限元。研究发现，尽管粘弹性本构关系具有时间相依性，其随机摄动格式并不存在“长期项”的影响，算例表明，应用该方法进行粘弹性结构的随机模拟，计算效率较高、精度较高。     

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

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

基于随机有限元的非线性结构稳健性优化设计

结合结构优化技术和摄动随机有限元方法研究了非线性结构稳健设计问题。将结构稳健性优化设计问题构造为双目标优化问题。优化目标包含结构性能函数的期望值和标准差。约束函数的变异也给予考虑，并采用基于函数梯度的算法进行求解。为对具有路径相关特征的非线性结构性能及结构响应的平均值及标准差进行分析。本文采用缩减的随机变量，提出了基于增量法的摄动随机有限元计算格式。在此框架下，进一步提出以一般泛函形式表达的结构性能的平均值和方差及其灵敏度的计算格式。为显示方法的有效性。文中给出几个数值算例。     

随机结构在随机激励下的二阶摄动随机势能泛函及其应用

利用小参数摄动法,建立了随机结构在随机激励下的二阶振动随机势能泛函。并由此推导了二阶摄动随机变原理,作为应用,建立了随机有限元的计算列式。     

线性随机结构的非平稳随机响应变异性分析

对于具有随机参数的结构受到非平稳随机激励的问题，给出了结构随机响应变异系数的虚拟激励摄动算法。它应用虚拟激励法先将随机荷载转化为确定性荷载，以使随机问题精确地转化为仅结构参数具有随机性的问题，从而将问题归结为应用随机摄动法求解单随机问题。求解过程简单高效，且有较高的精度。     

结构动力分析的随机变分原理及随机有限元法

将结构动力系统的参数及激励的随机性直接引入结构的动力泛函变分表达式中，基于瞬时最小势能原理，应用小参数摄动法，建立了随机结构动力分析的随机变分列式及相应的确机有限元法。算例表明，应用此法分析随机结构动力响应，具有程序实施简便，计算效率高的优点。     

一种新型随机有限元法

将Legendre积分法应用于随机结构的有限元分析，针对多随机变量非线性问题，建立基于Legendre积分法的随机有限元算法及列式。选择不同的Legendre积分点数目进行算例分析，并用Monte-Carlo法的计算进行对比，考察该方法的有效性。计算结果显示，单随机变量问题在很少样本点的情况下，一阶矩、二阶矩既有较高的精度，在选点数较多时，多随机变量问题的一阶矩、二阶矩也有足够的精度。考虑到计算上有很高的效率，该方法在随机有限元的计算上具有一定的价值。     

Legendre积分法在随机有限元法中的应用

将Legendre积分法应用于随机结构的有限元分析，针对非线性问题，建立基于Legendre积分法的随机有限元列式。选择不同的Legendre积分点数目进行算例分析，并用Monte—Carlo法的计算进行对比研究，考察该方法的有效性。计算结果表明本文提出的Legendre积分随机有限元有很高的计算效率，在精度上，较少的积分点在一阶矩、二阶矩计算上即有较高的精度，在积分点数较多时，三阶矩、四阶矩也有较高的精度。     

架设阶段悬索桥静力问题的随机有限元分析

架设阶段的悬索桥不仅具有显著的几何非线性行为，而且存在材料特性的变异。本文建立了以增量形式表达的非线性摄动法随机有限元列式，并编制了悬索桥静力分析程序ＳＮＡＰ。针对一座国内在建的悬索桥，计算了该桥在某一施工阶段考虑缆索弹性模量随机性时主梁跨中点位移的变化。结果表明，材料随机性的影响十分明显。     

非确定性结构静动态特性稳健优化设计

本文研究了考虑参数随机性的结构静动态特性稳健性优化设计问题的数学模型和数值求解。在考虑结构设计变量和其研究了考虑参数随机性的结构静动态特性稳健性优化设计问题的数学模型和数值求 解. 在考虑结构设计变量和其他参数随机分布的前二阶矩的条件下，采用基于二阶摄动法的 随机有限元方法对结构响应的平均值和方差进行近似求解. 在摄动法有限元分析的框架下， 提出以一般函数形式表达的结构性能的平均值和标准差及其灵敏度的计算格式. 将结构 稳健性优化设计问题构造为双目标优化问题，优化目标包含结构性能函数的期望值和标准 差，约束函数的变异也给予考虑. 优化问题采用基于函数梯度的算法进行求解. 文中给出的数值算例显示了方法的有效性.     

含随机参数的非线性黏弹性有限元计算

进行了粗粒土与结构接触面单调和循环加载试验，基于宏细观测量结果, 扩展了 损伤概念以 描述该类接触面在受载过程中的物态演化, 及由于物态演化导致的力学特性从初始状态到最终 稳定状态的连续变化过程. 揭示了接触面损伤的细观物理基础主要是接触面内土的颗粒破碎 和剪切压密这两种物态演化；指出接触面的剪胀体应变可以划分为可逆性和不可逆性剪胀体 应变两部分，其中不可逆性剪胀体应变可作为接触面损伤发展的宏观量度，因此其归一化 形式可作为一种损伤因子的定义；提出了建立粗粒土与结构接触面一种损伤本构关系的基本思路.     

平纹机织碳纤维复合材料的多尺度随机力学性能预测研究

平纹机织碳纤维复合材料在结构上具有多尺度特性和空间随机性. 同时, 组分材料会因存储条件和组成相成分、批次的不同导致力学性能有所差异. 当考虑各尺度结构和组分性能参数不确定性进行随机力学性能预测时, 存在以下两个难点: 一是随机变量众多, 使得对不确定性传递方法的精度和效率提出了要求; 二是由于随机参数之间存在高维相关性, 需要建立高精度的相关性模型. 针对以上问题, 本文提出了基于混沌多项式展开和Vine Copula的平纹机织复合材料多尺度随机力学性能预测方法, 综合考虑了平纹机织碳纤维复合材料微观及介观尺度的材料、结构随机参数, 基于自下而上层级传递的策略逐尺度地研究力学性能不确定性. 该方法采用Vine Copula理论构造相关随机变量的高维联合概率分布, 并运用非嵌入式混沌多项式展开法实现不确定性传递. 结果显示, 本方法构造的相关性模型几乎与原模型一致, 且能够高效准确地实现各尺度力学性能的随机预测.      

随机结构系统的一般实矩阵特征值问题的概率分析

由于工程实际结构的复杂性和所用材料在统计上的离散性以及测量、加工、制造误差的存在，必然导致具有随机参数的随机结构振动系统，按结构参数的性质来划分，随机振动问题包括两方面内容：(1)确定结构问题；(2)随机结构问题。本文以现代数学理论为依托，研究了随机结构系统的一般实矩阵的特征值问题。根据Kronecker代数、向量值和矩阵值函数的灵敏度分析、一般二阶矩法和概率摄动技术给出了计算随机结构系统的一般实矩阵的特征值和特征向量的数值方法，可以有效地得出随机结构系统的一般实矩阵的特征向量的统计量，发展了2D矩阵值函数的随机结构系统的特征值问题概率分析理论。     

埋地管道随机振动的摄动分析

埋地管道的材料特性和沿线的土壤性质存在差异，对这种差异性采用随机参数描述有一定的合理性，因此在管道的随机振动分析中考虑结构参数的随机性是必要的，对于管道的抗震设计具有现实意义。对于在空间相关的地震地面随机激励下的埋地管道，将结构参数看作是随机变量，采用摄动分析法，推导了随机响应的相关函数和功率谱密度函数的解析表达式，大大方便了工程应用。针对某输油管道，给出了计算结果。     

随机杆系结构几何非线性分析的递推求解方法

建立了随机静力作用下考虑几何非线性的随机杆系结构的随机非线性平衡方程. 将和 位移耦合的随机割线弹性模量以及随机响应量表示为非正交多项式展开式，运用传统的摄动方法获 得了关于非正交多项式展式的待定系数的确定性的递推方程. 在求解了待定系数后，利用非 正交多项式展开式和正交多项式展开式的关系矩阵，可以很方便地得到未知响应量的二阶统计矩. 两杆结构和平面桁架拱的算例结果表明，当随机量涨落较大时，递推随机有限元方法比基于 二阶泰勒展开的摄动随机有限元方法更逼近蒙特卡洛模拟结果，显示了该方法对几何非线性 随机问题求解的有效性.     

复合材料层合板的谱随机无网格伽辽金法

针对基于摄动理论的计算方法不适用于分析随机结构参数大变异问题,提出谱随机无网格伽辽金法,该方法基于随机场正交分解理论,将随机场采用Karhunen-Loève级数展开为一系列不相关随机变量,再引入结构位移随机响应的混沌多项式分解,结合无网格伽辽金法,从而导出含随机变量的复合材料层合板的谱随机无网格伽辽金法,给出结构响应的统计特征值的计算公式,该方法既适用于随机参数大变异情况,又具有无网格法的优势;数值算例结果表明该方法是正确有效的。     

Change-of-variable interval stochastic perturbation method for hybrid uncertain structural-acoustic systems with random and interval variables

The response analysis from hybrid uncertain structural-acoustic systems with random and interval variables (HUSAS) plays an important role in the optimal design of structural-acoustic systems. In this work, a hybrid uncertain numerical method known as the change-of-variable interval stochastic perturbation method (CVISPM) is proposed to predict the interval of the response probability density function and the response confidence interval of a HUSAS. This method is based on perturbation analysis and the change-of-variable technique. In the proposed method, the response of a HUSAS is approximated as a linear function of random variables using the stochastic perturbation analysis. According to the approximated linear relationships between the response and the random variables, the change-of-variable technique is introduced to calculate the response probability density function. Based on the response probability density function, the interval perturbation approach is used to predict the interval of the response probability density function and the response confidence interval. A numerical example of a shell structural-acoustic system with random and interval variables was employed to verify the effectiveness and precision of the proposed method.     

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

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