多种群遗传算法在PBX本构模型参数识别中的应用

利用多种群并行结构对标准遗传算法SGA进行并行化处理,引入移民算子和精华种群形成多种群遗传算法MPGA,并设计了自适应交叉和变异概率对算法的收敛速度进行改进。结合ABAQUS软件和改进的多种群遗传算法,建立了材料本构模型参数识别方法。采用该方法对PBX炸药黏弹性损伤本构模型参数进行了模拟识别,并同基于标准遗传算法的参数识别方法进行了比较。结果证明,基于改进多种群遗传算法IMPGA的方法对克服算法未成熟收敛有显著的效果,识别结果更稳定。同时该方法的收敛速度更快,寻优能力更强,适合复杂非线性问题的优化,此方法可以被应用到其他材料本构模型的参数识别中。     

基于爆炸切割试验的有机玻璃本构模型参数反演

为了获取爆炸切割数值模拟中有机玻璃（PMMA）的材料本构模型参数,建立了一种基于神经网络的有机玻璃Johnson Holmquist ceramics (JH-2)本构模型参数反演方法：基于从爆炸切割试验和现有研究得到的JH-2本构模型经验参数,确定本构模型参数的调整区间;使用LS-DYNA数值模拟软件对2.5 mm宽爆炸切割索切割14 mm PMMA平板过程进行数值模拟并收集平板损伤数据集;建立PMMA平板本构模型参数与损伤数据之间的神经网络模型;通过训练完成的神经网络模型对PMMA平板的JH-2本构模型参数进行反演。为验证通过反演参数的可靠性,进行了4.2 mm宽爆炸切割索切割19 mm PMMA平板试验和有限元数值模拟,计算结果中的平板损伤情况与实验结果相差较小,表明通过反演获得的JH-2本构模型参数能较好地应用于PMMA平板爆炸切割数值模拟。传统材料参数获取方法,该参数反演方法相较于可以通过较少的试验及测试,获得比较准确的材料本构模型参数。     

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

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

基于机器学习软骨细胞的时间依赖性力学行为及本构参数反演

探究软骨细胞机械负载下的力学特性对于理解软骨细胞的正常和病理状态以及骨性关节炎的病因至关重要. 基于软骨细胞有限元计算模型的力学响应与其本构参数之间的高度复杂非线性, 本文提出了分别利用双向深度神经网络TW-Deepnets模型和随机森林RF模型并结合有限元方法来识别软骨细胞本构参数的两种反演方法. 首先, 建立了软骨细胞的无侧限压缩实验有限元模型, 收集MSnHS本构参数空间点与对应的有限元计算模型的压缩反作用力响应数据集. 其次, 结合贝叶斯超参数优化算法搭建了用于软骨细胞本构参数反求的TW-Deepnets模型和RF模型, 对有限元收集的数据进行训练, 并利用单个软骨细胞受到50%压缩程度下的实验数据对软骨细胞的MSnHS本构参数进行了反求. 最后, 通过与实验曲线的对比验证了所提出的反演方法的有效性, 并引入决定系数R2对两种模型的预测准确性进行了对比评估, 检验了模型对各本构参数的预测性能, 分析了MSnHS本构模型中各参数影响软骨细胞力学响应的重要性占比. 结果表明, 本研究提出的本构参数反演方法能够有效获取软骨细胞的本构参数值, 从而准确描述软骨细胞的时间依赖性力学特性, 该方法也可进一步推广到生物细胞在静态或动态负载条件下的复杂参数反演问题.      

弹性或弹塑性土体中桩基的大变形分析

采用弧坐标,首先建立了位于弹性地基或弹塑性地基上并具有初始位移的桩基大变形行为的非线性微分方程组,并采用Winkeler模型来模拟地基对桩基的抗力;其次,应用微分求积方法离散非线性微分方程组,得到一组离散化的非线性代数方程,并给出了利用Newn-Raphson方法求解非线性代数方程的步骤;作为应用给出了数值算例,得到了桩顶受组合载荷作用时,变形后桩基的构形、弯矩和剪力,考察了土的弹性和弹塑性性质、桩基初始位移、荷载等参数对桩基力学行为的影响.最后将模型进行简化,得到了小变形理论的解析解,并比较了由大变形理论与小变形理论所得结果的差别.     

主应力轴旋转下黏土累积变形的分数阶元件模型

在长期交通载荷作用下土体塑性累积变形本构模型对路基沉降计算至为关键.元件组合模型可以计算岩土体循环累积应变,但现有的各类元件模型未能反映饱和软黏土的主应力轴循环旋转现象.在对饱和软黏土进行等向固结条件下的主应力轴循环旋转加载试验及非等向固结下的循环扭剪试验基础上,将Abel黏壶代替Burgers模型中的Newton黏壶,得到分数阶Burgers模型;利用遗传算法优化循环塑性累积应变的Burgers模型和分数阶Burgers模型的参数,通过对比两组模型的计算值与试验值,发现分数阶模型更适合模拟计算循环载荷下饱和软黏土的累积变形.     

基于遗传算法的弹性地基加肋板肋梁无网格优化分析

基于遗传算法及一阶剪切理论,提出一种弹性地基上加肋板肋条位置优化的无网格方法.首先,通过一系列点来离散平板及肋条,并用弹簧模拟弹性地基,从而得到加肋板的无网格模型;其次,基于一阶剪切理论及移动最小二乘近似原理导出位移场,求出弹性地基加肋板总势能;再次,根据哈密顿原理导出结构的弯曲控制方程,并通过完全转换法处理边界条件;最后,引入遗传算法和改进遗传算法,以肋条的位置为设计变量、弹性地基板的中心点挠度最小值为目标函数,对肋条位置进行优化达到地基板控制点挠度最小的目的.以不同参数、载荷布置形式的弹性地基加肋板为例,与ABAQUS有限元结果及文献解进行比较.研究表明,采用所提出的无网格模型可有效求解弹性地基上加肋板弯曲问题,结果易收敛,同时基于遗传算法与改进混合遗传算法所提出的无网格优化方法均可有效优化弹性地基加肋板肋条位置,后者计算效率相对较高,只进行了三次迭代便可获得稳定的最优解,此外在优化过程中肋条位置改变时只需要重新计算位移转换矩阵,又避免了网格重构.     

基于能量密度等效的超弹性压入模型与双压试验方法

针对超弹性材料压入问题, 本文基于能量密度中值等效原理, 提出了描述球、平面、锥3类压头独立压入下载荷、深度、压头几何尺寸和Mooney-Rivlin本构关系参数之间关系的半解析超弹性压入模型(semi-theoretical hyperelastic-material indentation model, SHIM), 进而提出了球、平面、锥压入组合的双压试验方法(indentation method due to dual indenters, IMDI). 正向验证表明, 基于系列超弹性材料的本构关系参数, 由SHIM分别预测的球、平面、锥3类压入下的载荷-位移曲线与有限元分析(finite element analysis, FEA)结果之间密切吻合; 反向验证表明, 基于系列超弹性材料的FEA条件本构关系下3类压入的载荷-位移曲线, 由双压试验方法预测的Mooney-Rivlin本构关系与FEA条件本构关系密切吻合. 针对3种超弹性橡胶, 完成了球、平面、锥压入试验, 应用双压试验方法获得的3组Mooney-Rivlin本构关系均与单轴拉伸试验结果吻合良好.      

基于能量密度等效的超弹性压入模型与双压试验方法

针对超弹性材料压入问题,本文基于能量密度中值等效原理,提出了描述球、平面、锥3类压头独立压入下载荷、深度、压头几何尺寸和Mooney-Rivlin本构关系参数之间关系的半解析超弹性压入模型(semi-theoretical hyperelastic-material indentation model, SHIM),进而提出了球、平面、锥压入组合的双压试验方法 (indentation method due to dual indenters, IMDI).正向验证表明,基于系列超弹性材料的本构关系参数,由SHIM分别预测的球、平面、锥3类压入下的载荷-位移曲线与有限元分析(finite element analysis, FEA)结果之间密切吻合;反向验证表明,基于系列超弹性材料的FEA条件本构关系下3类压入的载荷-位移曲线,由双压试验方法预测的Mooney-Rivlin本构关系与FEA条件本构关系密切吻合.针对3种超弹性橡胶,完成了球、平面、锥压入试验,应用双压试验方法获得的3组Mooney-Rivlin本构关系均与单轴拉伸试验结果吻合良好.     

基于遗传算法的弹性地基加肋板肋梁无网格优化分析

基于遗传算法及一阶剪切理论, 提出一种弹性地基上加肋板肋条位置优化的无网格方法. 首先, 通过一系列点来离散平板及肋条, 并用弹簧模拟弹性地基, 从而得到加肋板的无网格模型; 其次, 基于一阶剪切理论及移动最小二乘近似原理导出位移场, 求出弹性地基加肋板总势能; 再次, 根据哈密顿原理导出结构的弯曲控制方程, 并通过完全转换法处理边界条件; 最后, 引入遗传算法和改进遗传算法, 以肋条的位置为设计变量、弹性地基板的中心点挠度最小值为目标函数, 对肋条位置进行优化达到地基板控制点挠度最小的目的. 以不同参数、载荷布置形式的弹性地基加肋板为例, 与ABAQUS有限元结果及文献解进行比较. 研究表明, 采用所提出的无网格模型可有效求解弹性地基上加肋板弯曲问题, 结果易收敛, 同时基于遗传算法与改进混合遗传算法所提出的无网格优化方法均可有效优化弹性地基加肋板肋条位置, 后者计算效率相对较高, 只进行了三次迭代便可获得稳定的最优解, 此外在优化过程中肋条位置改变时只需要重新计算位移转换矩阵, 又避免了网格重构.      

粘弹塑性统一本构模型的参数评估

针对粘弹塑性统一本构模型参效众多的特点和传统评估参效方法的缺点，应用非线性最小二秉法与敏感系数分析方法相结合，提出了一个比较系统的对粘弹塑性统一本构模型的参数进行评估的方法。对铜合金在室温下不同加载速率的试验曲线进行参数适配，表明本方法的合理性与正确性。为粘弹塑性统一本构模型的工程应用奠定了基础。     

基于遗传算法的设计地震反应谱标定方法

基于自然进化原理的遗传算法是一种解决复杂非线性规划问题十分有前景的方法。对遗传算法中一些有代表性的算法作了简要分析和评价。介绍了用地面峰值加速度(PGA)、峰值速度(PGV)和峰值位移(PGD)三个参数标定设计地震反应谱的模型。针对此分段非线性函数规划问题,运用MATLAB编写遗传算法程序进行了计算。计算表明遗传算法在全局优化、非线性优化、多参数优化等方面表现出了传统算法无法比拟的优势。     

黄土损伤与流变耦合模型及参数研究

以试验为基础,对不同含水量Q3黄土的蠕变特征进行了分析研究,结合黄土结构损伤特性及其对黄土力学特征影响,通过考虑结构损伤对变形模量的劣化来反映黄土结构性对蠕变规律的影响,提出了黄土损伤与流变耦合作用的本构模型。该模型不仅具有参数少的优点,而且能够准确描述黄土的突发性破坏特征。利用结构性软土流变的研究成果,给出了黄土流变模型参数的确定方法,并对模型参数随非饱和黄土含水量变化的规律进行了研究。数值模拟结果表明,该模型能够很好地描述黄土的加速蠕变阶段的变形特征。     

Parameter determination and response analysis of viscoelastic material

The parameter determination of viscoelastic material is a multi-variable, multi-aim nonlinear optimization problem, which made the optimization process very complicated. In this paper a hybrid optimal algorithm was proposed to determine the viscoelastic parameters in the constitutive relation according to the experimentally obtained mechanical properties. This algorithm merges the Broydon–Fletcher–Goldfarb–Shanno search into a genetic algorithm framework as a basic operator in order to enhance the local search capability. The proposed hybrid algorithm not only can reduce the iterative times greatly but can abolish the limitation of initial parameter values. Nonlinear material characteristic curve-fitting was carried out using the proposed algorithm and other existing approaches. And the comparison results show this algorithm is accurate and effective. The numerical simulation and experimental study of viscoelastic cantilever beam also indicates that the finite element formulation and the calculative viscoelastic model parameters are reliable. The proposed optimization method can be extended to further complex parameter estimation researches.     

岩石三轴蠕变试验黏弹性解析及参数识别

基于岩石三轴蠕变试验,推导了相应的黏弹性解析解,并对三参量H-K流变模型给出 了黏弹性解析解的显式以及流变参数的识别方法. 在此基础上,针对花岗岩及片岩两种岩石 的室内蠕变分级加载试验曲线进行了流变参数的识别与反演. 分析表明,同一岩样在不同加 载级别下的弹性及黏弹性参数在量值上均有所不同,其中黏弹性模量和黏滞系数的量值随 载荷的增加出现损伤衰减,基于试验反演所获参数的解析蠕变曲线与试验蠕变曲线吻合较好.     

土结构性本构模型研究现状综述

土本构模型的建立是一个重要而又复杂的问题,到目前为止,国内外学者们已提出数以百计的土本构模型,诸多文献也对这些模型进行了评述和归纳。然而这些土本构模型多是在扰动土或砂土的基础上发展和建立起来的,它们难以描述由于土结构性引起的各种非线性行为,其计算结果与实际情况相差甚远。天然土体一般都具有一定的结构性,所以有必要建立考虑土结构性影响的土本构模型。针对这个现实,目前有些学者已基于各种理论和方法,提出了一些可以考虑土结构性影响的土本构模型,并得了较好的应用。但在目前的文献中还很少有对土的结构性本构模型研究进行归纳,基于此,本文简要介绍了一下目前土的结构性本构模型研究现状,并提出了这些本构模型在应用中所存在的问题。     

CONSTITUTIVE MODELS AND HYPERGRAVITY PHYSICAL SIMULATION OF SOILS$^{\bf 1)}$

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.     

土体的本构模型和超重力物理模拟

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.