比奥固结有限元中参数正交敏感性分析

用邓肯-张模型求解比奥平面固结有限元涉及到的参数较多,本文通过正交试验研究参数共同随机变化对软土地基固结沉降的影响,发现其中的k、c、G、φ和F对整个固结沉降过程的影响是高度显著的; Pc在固结的初期影响是显著的;在其他参数变异较大时,Rf对固结沉降的影响相对较小,常渗透系数Kz在固结初期影响显著。本文结果对模型参数的取值、反演、随机分析以及固结沉降的估计具有参考价值。     

基于孔压静力触探(CPTU)测试的连云港海相粘土前期固结压力确定方法研究

首先回顾了基于孔压静力触探(CPTU)测试确定前期固结压力的方法,通过连云港海相粘土场地进行的CPTU试验资料,以室内固结试验得到的前期固结压力作为参考值评估了经验方法预测前期固结压力的有效性。最简单的方法是直接建立前期固结压力和净锥尖阻力的关系,同时也是最有效的方法。     

求解饱和半空间上弹性圆板固结沉降的积分方程

本文采用解析方法分析了弹性圆板在饮和半空间上的固结沉降。考虑弹性圆板与饮和半空间的接触面上无摩擦力,且饱和半空间表面为全部透水的。运用Ｂｉｏｔ固结理论和积分方程技术,在Ｌａｐｌａｃｅ变换域上建立了弹性圆板固结沉降的对偶积分方程,并化此对偶积分方程为第二类Ｆｒｅｄｈｏｌｍ积分方程。通过对其核函数的有效数值发得到第二类Ｆｒｅｄｈｏｌｍ积分方程的解,再利用Ｌａｐａｃｅ反演技术获得弹性板在时间域中的固结沉     

某软土深基坑工程时间效应有限元分析

软土地区的深基坑因土体的固结作用和流变性而具有了时间效应。本文以Biot固结有限元法为基础,用三元件模型中的第一个线性弹簧模拟固结作用,弹性模量考虑了开挖应力路径和应力历史的影响;另外一部分(KELVIN模型)来模拟土体的流变性,以实际变形的反演来得到两个参数的大致取值,再对基坑的变形情况以及进一步开挖进行分析。假定为正常固结饱和粘土,平面应变问题。通过对某饱和软粘土地基深基坑开挖工程实例的分析,得到的挡墙水平位移曲线与实测曲线很吻合,表明程序较好地反映出基坑的时间效应。     

基于非保守体力的大变形固结有限元法

考虑土体固结过程中体力的非保守特点,完善了T.L.描述的大变形固结平衡方程,并以伴生体力代替非保守体力,提出了相应的大变形固结增量有限元列式. 当忽略土体体积应变的影响时,改进的大变形固结平衡方程及其有限元法均可退化为基于保守体力的传统形式. 采用自主开发的程序对大变形自重固结算例进行分析,结果初步表明:传统方法低估了地基固结度,但高估了地基沉降量和超静孔隙水压力. 基于非保守体力的大变形固结分析在理论上能够更合理地预测地基的固结性状.      

弹性动力学反问题的非线性及其迭代反演

分析了弹性动力学反问题的非线性及在迭代反演过程中表现出的复杂非线性现象。迭代反演的结果依赖于反演系统参数和迭代初值,而系统参数对应的Mandelbrot集和迭代初值对应的Julia集都是复杂的分形集。随反演系统状态参数的变化,完全确定性的反演系统却可能产生一系列无规则的,不可预测的迭代输出序列。反演迭代过程中出现的分形和混沌现象反映了表面简单的反演迭代后隐藏的复杂性,正是这种复杂性给迭代系统参数的合理选择带来困难,进而使反演迭代不总能给出满意的结果。     

岩石蠕变模型的比较和修正

对三种岩石在不同的应力下进行了蠕变实验,对蠕变实验曲线用粘弹模型进行了模型参数的反演,结果认为标准线性体、伯格斯和对数模型都是实验曲线的很好近似.为了进一步了解蠕变实验时间对模型参数的影响,研究了蠕变模型的预测功能.将实验蠕变曲线分成二部分,第一部分用于模型参数的反演,第二部分用于模型预测数据的检验.时间外延所得的模型的预测曲线都明显偏离了实测曲线,伯格斯和对数模型的预测曲线比实测曲线偏高,而标准线性体模型则偏低.根据不同时间的蠕变实验数据的比较和分析,认为实验时间对蠕变模型参数的反演有较为明显的影响,但只要对标准线性体模型作适当修正,引入与时间尺度有关的函数C(t),则该模型的预测能力就有较大的提高.     

基于进化多项式回归方法的土体本构关系

进化多项式回归方法是近几年来在国际上引起广泛关注的一种新颖的数据挖掘方法。本文在土样固结排水试验数据的基础上,首次使用该方法建立了土体本构关系。采用上海嘉定区某地软土的固结排水试验结果对模型进行了验证。结果表明,该本构关系能够更精确地描述土体的应力-应变特性,其预测精度比修正剑桥模型高19.6%~57.7%。另外,本文还用切线刚度法,对基于进化多项式回归方法的土体本构关系进行了ABAQUS二次开发,对软土固结排水试验进行了有限元数值模拟。     

基于五次样条函数的基坑土层等效参数反演方法

根据五次样条函数的性质,利用弹性支点法,将土压力视为不连续的分布载荷集度,对支护结构的变形进行五次样条函数拟合,并通过逐次求导反演支护结构上的受力,最后结合工程情况反演基坑土层的等效参数。本文基于上述研究思路,成功反演了支护结构的分布载荷集度,提出了一种基于五次样条函数的基坑土层等效参数反演方法,反演的土层等效参数可用于基坑工程的动态化设计、施工方案优化和变形预测,并为邻近类似土层区域的基坑设计提供土层抗剪强度指标参考。     

超固结土的一维固结理论及其试验研究

基于超固结饱和土固结理论,针对变荷载、自重应力分布等条件下提出3了相应的一维超固结土问题的控制方程,总结出六个超固结土的一维线性固结计算参数.并采用半解析法编制了求解该问题的Fortran计算程序.利用GDS高级固结仪对萧山黏土进行超固结土的一维固结渗透联合试验以及一维固结压缩试验,将测定的一维超固结土计算参数引入理论程序中,分别对比了位移、平均固结度、底部孔压这三者随时间变化的理论和试验曲线.结果表明:超固结土的一维线性固结理论比正常固结土的一维线性固结理论更接近试验结果,验证了超固结土的固结理论,表明了考虑土体应力历史对固结计算的重要性.     

A quadratic hierarchical Bayesian dynamic prediction model for infrastructure maintenance

The main objective of this study is the development of a correlation model in dynamic Bayesian belief networks (DBBNs) followed by an inverse economic analysis. This is based on a quadratic hierarchical Bayesian inference prediction method using Markov chain Monte Carlo simulations. The developed model is implemented to predict the future degradation and maintenance budget for a suspension bridge system. Bayesian inference is applied to find the posterior probability density function of the source parameters (damage indices and serviceability), given 10 years maintenance data. The simulated risk prediction under decreased serviceability conditions gives posterior distributions based on a prior distribution and likelihood data updated from annual maintenance tasks. Compared with a conventional linear prediction model, the proposed quadratic model provides highly improved convergence and closeness to the measured data. Finally, the developed inverse DBBN analysis method allows forecasts of future performance and the financial management of complex infrastructures by providing the sensitivity of serviceability and risky factors to the maintenance budgets of structural components and the overall system.     

Topology identification of complex networks from noisy time series using ROC curve analysis

As it is known, there is various unknown information in most real world networks, such as uncertain topological structure and node dynamics. Thus how to identify network topology from dynamical behaviors is an important inverse problem for physics, biology, engineering, and other science disciplines. Recently, with the help of noise, a method to predict network topology has been proposed from the dynamical correlation matrix, which is based on a general, one-to-one correspondence between the correlation matrix and the connection matrix. However, the success rate of this prediction method depends on the threshold, which is related to the coupling strength and noise intensity. Different coupling strength and noise intensity result in different success rate of prediction. To deal with this problem, we select a desirable threshold to improve the success rate of prediction by using Receiver Operating Characteristic (ROC) curve analysis. By the technique of ROC curve analysis, we find that the accuracy and efficiency of topology identification is mainly determined by coupling strengths. The success rate of estimation will be reduced if the coupling is too weak or too strong. The presence of noise facilitates topology identification, but the noise intensity is not always crucial to the effectiveness of topology identification.     

Parameter identification for improved viscoplastic model considering dynamic recrystallization

The successful application of viscoplastic model considering dynamic recrystallization depends on how well the parameters are identified. However, it is difficult to obtain satisfactory parameters using conventional parameter identification methods. The reasons are due to difficulties in obtaining homogeneous deformation, high complexity of physical process described by the model and large number of parameters. In this paper, the material parameters are identified by inverse analysis. Global information on objective function is firstly studied by an improved uniform random sampling method; secondly, a hybrid global optimization method, which combines the genetic algorithm, the Levenberg–Marquardt algorithm, the augmented Gauss–Newton algorithm and the flexible tolerance method, is constructed and an inverse analysis numerical procedure, which combines the proposed optimization method with the finite element analysis, is proposed; at last, a set of satisfactory material parameters for 26Cr2Ni4MoV is obtained by the proposed inverse analysis numerical procedure.     

用于板料成形反向模拟的特殊应力更新算法

提出了一种改进的反向模拟法,以最终构型为研究对象,采用Euler坐标系,基于虚功原理获得有限元列式. 改进的反向模拟法采用了一种基于塑性流动理论的本构方程,可以充分考虑应变历史对塑性变形的影响. 为了避免流动理论应力更新算法过程中关于未知量\Delta\lambda 的非线性方程的求解,引入等效应力思想,无需Newton-Raphson迭代直接计算未知量\Delta \lambda . 盒形件的拉深实例中,传统的基于塑性形变本构方程的反向模拟法和改进的基于塑性流动本构方程的反向模拟法计算结果,分别与基于增量有限元法的正向数值模拟求解器LS-DYNA计算结果进行对比. 通过获得的坯料轮廓、成形极限图、等效应变分布、计算效率等的比较,验证了所提出的基于塑性流动理论本构模型的应力更新算法的有效性.      

基于混合网格的低声爆反设计方法研究

低声爆设计方法已成为新一代军民用超声速飞机研制过程中必须解决的关键难题之一。针对传统SGD低声爆外形反设计方法无法对声爆近场非线性效应进行描述和分析的缺点,提出了利用CFD方法求解得到的声爆近场压力分布代替F函数进行低声爆反设计的方法。声爆近场预测采用点-点对接的结构/非结构混合网格,充分利用非结构网格对复杂外形适应性强和结构化网格计算效率高的优点。结果分析表明,基于改进后的低声爆反设计方法得到的方案在声爆超压以及感觉噪声级等方面都比基于原始SGD方法得到的方案有较大改善。     

平面结构可靠性预测的随机边界元逆分析法

在结构可靠性计算的随机边界元理论和逆分析法基础上,用随机边界元方程结合改进的卡尔曼滤波算法来预测结构的可靠性,提出了用于平面机械结构可靠性预测的边界元逆分析法理论;并将方法应用于同时存在未知缺陷和随机边界载荷模型,进行了结构可靠性预测,在识别缺陷的同时预测出随机边界载荷的分布参数及结构体上各点的结构可靠度指标,给出了平面结构模型的数值示例。     

An interval perturbation method for exterior acoustic field prediction with uncertain-but-bounded parameters

This paper proposes two interval analysis methods, called the first-order interval parameter perturbation method (FIPPM) and the modified interval parameter perturbation method (MIPPM), for use in exterior acoustic field prediction when there are uncertainties in both the material properties and the external load. Interval variables are used to quantitatively describe the uncertain parameters in the face of limited information. The conventional first-order Taylor expansion and perturbation terms are employed in the FIPPM, while the MIPPM introduces modified Taylor series to approximate the non-linear interval matrix and vector. The high-order terms of the Neumann expansion are retained to calculate the interval matrix inverse. A numerical example is given by comparing the results with a Monte Carlo simulation to demonstrate the feasibility and effectiveness of the proposed methods at evaluating the sound pressure ranges in an exterior acoustic field.     

Interaction solutions for a reduced extended $$\mathbf{(3}\varvec{+}{} \mathbf{1)}$$ ( 3 + 1 ) -dimensional Jimbo–Miwa equation

Heavy-duty industrial robots have great advantages in the manufacturing industry. Considering the heavy process load and low stiffness of the robot, an accurate and efficient dynamic model plays an important role in the behavior analysis and performance improvement in the robot. This paper presents a novel methodology for the inverse dynamic analysis of the heavy-duty industrial robot with elastic joints. In particular, high-order kinematics and dynamics are concisely deduced using Lie group to deal with elastic joints for the robot inverse dynamic analysis. Meanwhile, position errors of the end-effector due to elastic joints are evaluated through the inverse dynamic analysis when the robot is in heavy-duty applications. Compared with previous approaches, the advantage of proposed method is that new formulas for inverse dynamic analysis are shown to be more concise and computationally efficient using Lie group. Moreover, the position error evaluation method considering dynamic forces is proved to be more accurate than the traditional method when the robot is in the high-speed application. Because of the high computational efficiency and accurate evaluation results, the proposed approach is applicable to trajectory optimization and position error compensation, especially for the robot in heavy-load and high-speed applications.     

基于Kriging代理模型的迭代更新高效反演方法

为提高混凝土坝等大体积结构参数反演效率和精度,减少由于应用有限元进行大量正分析而产生的计算机时,建立了一种结合Kriging代理模型和粒子群优化(PSO)算法的迭代更新反演方法。通过拉丁超立方抽样(LHS)方法确定初始样本点的空间分布,并使用有限元正分析获取对应的响应值,构建粗糙的初始代理模型,结合具有全局寻优能力的PSO算法,反演大体积结构的分区弹性模量,随之再代入有限元模型中,计算获取新的位移响应,并将其作为新样本加入到样本集中,通过迭代更新获得局部更高精度的代理模型。工程实际算例表明,该方法对混凝土坝等大体积结构参数反演精度较高和适用性好,且能大幅减少传统有限元模型反演方法所需消耗的正分析机时,提高反演效率。     

共轭梯度法求解非线性多宗量稳态传热反问题

应用共轭梯度法求解非线性多宗量稳态热传导反问题。采用八节点的等参单元在空间上进行离散，建立了便于敏度分析的非线性正演和反演的有限元模型，可直接求导进行敏度分析。给出了相关的数值验证，对测量误差及测点数目的影响作了初步探讨，结果表明，采用的算法能够对非线性稳态热传导中导热系数和边界条件联合反问题进行有效的求解，并具有较高精度。