基于静力响应面的结构有限元模型修正方法

提出了基于静力响应面的结构有限元模型修正方法.运用响应面方法,将结构静力响应和结构参数之间复杂的隐式关系用显式函数近似表达出来;在此响应面模型(函数)基础上,通过优化计算对结构有限元模型参数进行修正.阐述了基于静力响应面的结构有限元模型修正方法的基本理论和一般实现过程.对两跨连续梁结构的静力模型修正数值算例分析结果表明:基于静力响应面的有限元模型修正方法可以减少结构有限元计算的次数、提高模型修正的优化效率,结构有限元模型修正结果具有可接受的精度.     

基于贝叶斯信息融合的岸桥有限元模型修正研究

为建立精确的岸桥有限元模型,研究了基于贝叶斯信息融合的模型修正方法．通过方差分析,确定待修正参数,利用中心复合试验设计获取样本点,根据有限元计算结果与实测的结果残差为目标函数获得响应样本．拟合样本点和响应样本值构建二阶多项式响应面模型,并检验响应面模型的精度．基于贝叶斯理论更新融合系数来优化响应面参数,从而获得修正模型．以宁波大榭3号岸桥为工程背景,对比修正后的模态频率和实测频率,最大频率相对误差不超过5%,进而验证了基于贝叶斯信息融合的动力学有限元模型修正方法的有效性．修正后的有限元模型可进一步应用于岸桥的健康监测和安全评估．     

结构动力学有限元模型修正的目标函数及算法

结构动力学有限元模型修正是结构动力学领域的一个热点问题,回顾了结构动力学有限元模型修正研究的发展历史和现状,简要评述了结构动力学有限元模型修正所使用的设计变量,着重阐述了各种结构动力学有限元模型修正方法中所使用的目标函数及修正算法,讨论了工程结构动力学有限元模型修正的一些策略,最后对结构动力学有限元模型修正技术的发展进行了总结和展望.     

基于频响函数提升小波总能量的桁架结构随机模型修正

为了解决模型修正问题中的随机性,构建了一种基于提升小波总能量的随机模型修正方法.首先,将结 构的加速度频响函数进行提升小波变换,并提取提升小波总能量来代替加速度频响函数;然后,以待修正参数作为输入,提升小波总能量为输出构建响应面代理模型代替原来的有限元模型;接着,运用蒙特卡洛抽样抽取响应样本,并设定阈值筛选响应样本;最后,以代理模型预测得到的响应和抽样所得真实响应之间的差值最小为 目标函数,通过布谷鸟优化算法寻优求解待修正参数的均值.算例表明,所提方法修正后参数的最大误差小于3.3％,相应的频响函数曲线重合度高.     

A COMPLEX PARAMETER MODEL UPDATING METHOD FOR DAMPING CHARACTERISTICS OF FINITE ELEMENT MODELS

阻尼对于结构动力学响应具有重要的影响,但有限元模型一般很难对阻尼特性进行精确建模.基于实测频响函数,研究了一种有限元模型阻尼特性的复参数修正方法.以待修正区域各单元质量、刚度矩阵的比例修正系数为复修正参数,建立了单元矩阵比例修正的灵敏度方程直接算法,并对比分析了复修正参数与不同阻尼特性之间的数学关系.以六自由度集中参数模型和25杆平面桁架模型为例,验证了复参数修正方法在阻尼特性修正中的有效性.     

有限元模型阻尼特性的复参数修正方法研究

阻尼对于结构动力学响应具有重要的影响,但有限元模型一般很难对阻尼特性进行精确建模.基于实测频响函数,研究了一种有限元模型阻尼特性的复参数修正方法.以待修正区域各单元质量、刚度矩阵的比例修正系数为复修正参数,建立了单元矩阵比例修正的灵敏度方程直接算法,并对比分析了复修正参数与不同阻尼特性之间的数学关系.以六自由度集中参数模型和25杆平面桁架模型为例,验证了复参数修正方法在阻尼特性修正中的有效性.     

有限元模型修正研究进展:从线性到非线性

有限元分析在实际工程中得到了广泛应用.然而有限元模型由于受到网格划分、边界条件和材料物理参数不确定性等的影响,与真实结构有差异. 因此须通过试验数据加以修正,使其尽可能接近实际结构,以保证之后的结构动力模拟分析和监测等具有实际意义. 经过多年发展,有限元模型修正技术已经能够成功应用于一些实际工程,但现代工程技术的进步对有限元模型修正提出了更高要求,修正后的有限元模型不仅要有较高的精确度,还需要为后续应用给出具有指导意义的置信度.而现有的有限元模型修正、确认方法多基于结构线性的假设,而未能考虑实际结构中广泛存在的非线性.因此本文以土木工程结构模型修正的一些研究成果为例,通过对传统有限元模型修正的发展历程进行全面回顾;总结评述传统有限元修正技术的主要方法,以及包括有限元模型确认在内的最新研究进展;重点探讨有限元模型修正技术向非线性发展的技术路线和目前主要研究成果,展望其未来发展方向, 并提出值得研究的问题.      

基于径向基函数的频响函数模型修正

针对使用频响函数进行有限元模型修正存在的时间成本和精度问题,结合模态参与变异系数法和模态动能法分别优化激励点和测点,使获得的模态信息更完整;然后,引入径向基函数(RBF)模型减少原有限元模型计算次数,并根据均方根误差准则对所构建代理模型的参数(spread)进行优选,提高模型预测精度;最后,选用一个36自由度的二维桁架模型进行可行性验证,对比有限元法、Kriging模型及二阶响应面模型的修正精度和迭代时间,结果表明,本文方法具有较好的优势。     

有限元模型修正研究进展:从线性到非线性

有限元分析在实际工程中得到了广泛应用.然而有限元模型由于受到网格划分、边界条件和材料物理参数不确定性等的影响,与真实结构有差异.因此须通过试验数据加以修正,使其尽可能接近实际结构,以保证之后的结构动力模拟分析和监测等具有实际意义.经过多年发展,有限元模型修正技术已经能够成功应用于一些实际工程,但现代工程技术的进步对有限元模型修正提出了更高要求,修正后的有限元模型不仅要有较高的精确度,还需要为后续应用给出具有指导意义的置信度.而现有的有限元模型修正、确认方法多基于结构线性的假设,而未能考虑实际结构中广泛存在的非线性.因此本文以土木工程结构模型修正的一些研究成果为例,通过对传统有限元模型修正的发展历程进行全面回顾;总结评述传统有限元修正技术的主要方法,以及包括有限元模型确认在内的最新研究进展;重点探讨有限元模型修正技术向非线性发展的技术路线和目前主要研究成果,展望其未来发展方向,并提出值得研究的问题.     

Recent progress on finite element model updating: From linearity to nonlinearity

有限元分析在实际工程中得到了广泛应用.然而有限元模型由于受到网格划分、边界条件和材料物理参数不确定性等的影响,与真实结构有差异. 因此须通过试验数据加以修正,使其尽可能接近实际结构,以保证之后的结构动力模拟分析和监测等具有实际意义. 经过多年发展,有限元模型修正技术已经能够成功应用于一些实际工程,但现代工程技术的进步对有限元模型修正提出了更高要求,修正后的有限元模型不仅要有较高的精确度,还需要为后续应用给出具有指导意义的置信度.而现有的有限元模型修正、确认方法多基于结构线性的假设,而未能考虑实际结构中广泛存在的非线性.因此本文以土木工程结构模型修正的一些研究成果为例,通过对传统有限元模型修正的发展历程进行全面回顾;总结评述传统有限元修正技术的主要方法,以及包括有限元模型确认在内的最新研究进展;重点探讨有限元模型修正技术向非线性发展的技术路线和目前主要研究成果,展望其未来发展方向, 并提出值得研究的问题.     

基于WENO-THINC/WLIC模型的水气二相流数值模拟

水气二相流与诸多领域的实际工程问题密切相关. 对二相流运动进行高精度的数值模拟是计算流体力学研究的难点和热点. 针对开敞水域的自由表面流运动问题, 将水和空气均视为不可压缩流体, 采用五阶加权基本无震荡(weighted essentially non-oscillatory, WENO)格式求解描述流体运动的纳维斯托克斯(Navier-Stokes, NS)方程, 利用以加权线性界面算法改进的多维双曲正切函数界面捕捉法(tangent of hyperbola for interface capturing with weighed line interface calculation, THINC/WLIC)追踪水气界面, 建立WENO-THINC/WLIC水气二相流运动数值模型. 模型采用分步计算法离散求解控制方程, 通过压力投影法求解压强场, 并利用三阶总变差递减(total variation diminishing, TVD)龙格库塔(Runge-Kutta, RK)法对时间项进行离散求解. 通过对环境速度场下Zalesak's disk和shearing vortex界面运动问题, 线性液舱晃荡问题以及溃坝问题的模拟结果与理论分析或试验结果的比较, 对所建立的水气二相流数值模型的适用性及模拟精度进行了验证. 结果表明, 本模型的模拟结果与物理模型或理论分析结果吻合良好, 能较为准确地再现不可压缩水气二相流运动现象. 鉴于WENO格式和THINC法本身在算法及应用等方面仍在不断改进, 本研究提出的WENO-THINC耦合模型为后续更高精度的二相流计算模型开发提供了一种研究思路.      

基于神经网络对扁壳结构载荷位置识别问题的研究

在建立扁壳结构识别问题的力学模型基础上,运用神经网络方法构造了三层结构的神经网络,并以此对扁壳在集中力作用下的载荷位置的识别问题进行了研究,结果表明,采用有限元方法获取基本样本,运用神经网络方法对扁壳结构载荷位置的识别具有较好的效果。     

固有频率法评估损伤的阈值研究

通过分析损伤可测性的影响因素,对固有频率法判断损伤阈值的估算方法进行了研究,并以 悬臂梁模型为例,应用该方法研究了悬臂梁的损伤阈值. 数值及理论计算的结果表明：在已 知频率测量精度及方法偏差时,该方法可以判别固有频率法损伤识别的阈值,在早期损伤阶 段,该方法判别的阈值与数值模拟结果吻合良好.     

基于一类非局部宏-微观损伤模型的裂纹模拟

结合近场动力学和统一相场理论的基本思想, 最近提出了一类非局部宏-微观损伤模型, 为固体裂纹扩展模拟提供了新途径. 本文在此基础上改进了微观损伤准则, 并给出损伤的$\bar{\lambda}-\ell$语言以刻画固体破坏过程中位移场的不连续程度. 在改进模型中, 首先根据两物质点(即物质点对)之间的变形量, 基于相对临界伸长量的历史最大超越程度, 给出表征物质键性能退化的微细观损伤. 进而, 对影响域内的物质键损伤进行空间局部加权平均, 获得宏观拓扑损伤. 通过引入能量退化函数, 建立基于能量的损伤与宏观拓扑损伤之间的关系, 由此将其嵌入连续损伤力学基本框架, 形成了问题求解的基本方程. 该模型是一类非局部化模型, 可采用有限单元法进行离散求解, 避免了经典局部损伤力学所面临的网格敏感性问题. 文中, 进一步将其应用于具有强非线性回弹特性的裂纹扩展模拟问题. 实例分析表明, 本文方法不仅可以把握裂纹扩展模式, 而且能够定量刻画裂纹扩展过程中的载荷-变形关系. 最后指出了需要进一步研究的问题.      

????????????????????????????

基于群体与适应度的概念，应用改进的PSO算法，从随机解出发，提出了基于PSO 算法的本构模型参数识别方法. 该方法解决了橡胶类材料大应变时硬化现象的本构模型参数 的确定这一难题. 首先通过单轴拉伸本构模型实验，在针对硬化实验曲线存在拐点的情 况下，应用PSO算法进行拟合，最后利用简单剪切实验进行验证. 结果表明该方法科学可行， 且具有速度快、精度高、易于收敛等优点. 并且有效地解决了本构模型参数识别的困难，可 广泛应用于各种复杂材料.     

利用三维随机模型分析沥青混合料的蠕变行为

将沥青混合料看作由粗骨料和沥青砂组成的两相复合材料,根据给定的级配生成凸多面体骨料,然后利用随机投放算法建立沥青混合料试样的三维随机模型.采用广义Maxwell模型刻画沥青砂的本构行为,其参数通过单轴蠕变实验获得.在对三维随机模型的有效性进行验证之后,采用参数化建模方法建立包含不同骨料分布、含量和级配的沥青混合料有限元模型,通过数值模拟研究骨料分布、含量和级配对沥青混合料蠕变行为的影响.结果表明:骨料分布对沥青混合料力学性质的影响较小;沥青混合料瞬时弹性模量随骨料含量的增加呈近似线性增大,而且弹性模量与骨料含量之间的关系曲线处于Paul上下边界内;在计算考虑的粒径范围内,骨料平均粒径越小,沥青混合料抵抗变形的能力越强.     

一种基于湍动能方程的转捩判定方法简

转捩现象是阻碍阻力高精度求解的主要问题之一. Menter 和Langtry 所提出的γ-θ转捩模型通过引入涡量雷诺数和间歇因子输运方程来驱动转捩,但是其中很多经验公式的理论立足点有待商榷. 驱使层流转变到湍流依赖的仍然是平均速度的一阶和二阶相关量,它们组合构成了湍动能方程的耗散尺度. 在湍动能方程中做合适的耗散平衡后,仅仅依靠湍动能方程可以有效地捕捉转捩现象. 采用自然转捩和旁路转捩测试算例进行了验证,结果证明该方法与试验值匹配较好,具有一定的工程实用价值.     

基于径向基点插值法的旋转Mindlin板高次刚柔耦合动力学模型

将无网格径向基点插值法(radial point interpolation method,RPIM)用于中心刚体?旋转柔性板的动力学分析.基于浮动坐标系方法和一阶剪切变形理论即Mindlin板理论,考虑剪切变形的影响,并计入板面内变形的非线性耦合变形项,采用径向基点插值法描述板的变形场,保留动能中有关非线性耦合变形项...     

基于单元破裂的岩石裂纹扩展模拟方法

传统离散元方法在处理破裂问题时, 采用界面上的准则进行判断, 裂纹只能沿着单元边界扩展. 当物理问题存在宏观或微观裂隙时, 在界面上应用准则具有其合理性; 而裂纹沿着单元边界扩展, 使得裂纹路径受网格影响较大, 扩展方向受到限制. 针对上述情况, 可以基于单元破裂的方式, 构建连续- 非连续单元法, 并应用于岩石裂纹扩展问题的模拟. 该方法在连续计算时, 将单元离散为具有物理意义的弹簧系统, 在局部坐标系下由弹簧特征长度、面积求解单元变形和应力, 通过更新局部坐标系和弹簧特征量, 可进一步计算块体大位移、大转动, 连续问题计算结果与有限元一致, 同时提高了计算效率. 在此基础上, 引入最大拉应力与莫尔—库伦的复合准则, 判断单元破裂状态和破裂方向, 并采用局部块体切割的方式, 在单元内形成初始裂纹. 裂纹两侧相应增加新的计算节点, 同时引入内聚力模型描述裂纹两侧的法向、切向作用与张开度及滑移变形之间的关系. 按此方式, 裂纹尖端处的扩展路径可穿过单元内部和单元边界, 在扩展方向的选取上更为准确. 最后, 通过三点弯曲梁、单切口平板拉伸、双切口试样等典型数值试验, 模拟裂纹在拉伸、压剪等各种应力状态下的扩展问题, 并对岩石单轴压缩试验的破坏过程进行模拟, 分析裂纹形成与应力—应变曲线各阶段之间的对应关系. 结果表明: 连续—非连续单元法通过单元内部破裂的方式, 可以显示模拟裂纹萌生、扩展、贯通直至形成宏观裂缝的过程.      

一种基于湍动能方程的转捩判定方法

转捩现象是阻碍阻力高精度求解的主要问题之一. Menter 和Langtry 所提出的γ-θ转捩模型通过引入涡量雷诺数和间歇因子输运方程来驱动转捩,但是其中很多经验公式的理论立足点有待商榷. 驱使层流转变到湍流依赖的仍然是平均速度的一阶和二阶相关量,它们组合构成了湍动能方程的耗散尺度. 在湍动能方程中做合适的耗散平衡后,仅仅依靠湍动能方程可以有效地捕捉转捩现象. 采用自然转捩和旁路转捩测试算例进行了验证,结果证明该方法与试验值匹配较好,具有一定的工程实用价值.