基于遗传算法的混凝土三维非稳态温度场反分析

大体积混凝土结构施工期间的合理温度控制问题非常重要，而精确进行温度计算所需的一些材料参数往往不易直接测得，需要根据一些易测得的量进行反求。本文基于三维瞬态温度场有限元求解理论与反问题理论，建立了混凝土三维瞬态温度场反问题求解数值模型。运用遗传算法寻求非线性反演问题全局最优解，只需要若干点温度实测值便可实现混凝土多个热学参数如绝热温升、导温、导热系数及热交换系鼓等的同时反演，算例对本文反演方法的反演精度及数值稳定性给出了满意的证明。     

自适应无网格热弹塑性接触模型研究

提出一种自适应无网格热弹塑性接触求解模型,求解接触问题的线性规划-增量初应力法与基于应变能梯度的自适应无网格法相结合,给出了模型计算理论和算法实现.通过圆柱体与弹塑性平面热弹塑性接触算例对模型进行验证.对是否考虑材料应变硬化,是否考虑摩擦力和热输入,是否考虑材料屈服强度温度相关等情况的两种算例进行了讨论.结果表明,该模型能有效地求解考虑不同情况下的热弹塑性接触问题,在较真实地模拟接触状况的同时,具有较高的计算精度和计算效率.     

弹塑性结构优化的并行算法

本文给出了弹塑性结构优化设计的序列二次规划算法，分别考虑了弹塑性结构分析，灵敏度分析和二次规划的Ｌｅｍｋｅ算法的并行计算。针对弹塑性结构灵敏度的不连续性，比较了连续模型和间断模型的计算结果，结果是接近的。算例表明此算法有很好的并行性。     

基于变形场测量数据主元压缩的模型参量反求方法

利用主元分析法对数字图像相关技术所测结构变形信息进行数据压缩,并进一步用于力学模型未知参量的反求计算。首先,为降低数字图像相关技术所测庞大数据的应用成本,提出利用主元分析法对结构表面变形场数据进行压缩,实现在保留结构表面变形信息主要特征的前提下显著降低数据量的目的;其次,针对压缩后的数据建立了基于最小二乘法的力学模型参量反求模型,并利用高斯牛顿法进行求解;最后,以具体算例从计算精度、收敛速度和抗噪性等方面验证了数据压缩对模型参量反求的效果。研究结果表明,所提方法在显著降低使用数据量的前提下,能够有效提高力学模型参量反求计算的收敛速度,特别是对于包含多个模型参数的反求问题,具有较高的精度和较好的稳定性。     

基于同伦技术的偶应力反问题求解

引入Bregman距离构造同伦函数,建立偶应力反问题的一种求解模式,可以对偶应力理论相关参数进行识别.利用有限元技术,建立了偶应力正/反问题数值求解模型,该模型不仅考虑了非均质的影响,而且也便于反问题的敏度分析.文中给出了相关的数值算例,并对信息误差和不同的函数形式计算结果做了初步探讨,得到满意的结果.数值算例表明该模...     

三维板料成形过程的显式有限元分析

本文采用基于随动坐标系的假定应变域壳单元及显式有限元格式求解三维板料成形问题。板材料服从Ｈｉｌ各向异性弹塑性准则，板料与模具之间的接触界面由主仆面接触搜寻法处理，接触力由罚参数法计算。文中给出了几个三维成形过程的计算实例。数值算例表明，本文方法具有较高的计算精度和计算效率，可在微机上分析中等复杂程度的成形过程     

各向异性弹塑性中厚度板壳的有限元分析

本文在文献［２，３］的基础上，提出了一个解各向异性弹塑性中厚度板壳问题的有限元方法。考虑材料各向异性的特点，采用了Ｈｉｌｌ推广的Ｈｕｂｅｒ－Ｍｉｓｅｓ屈服准则；借用Ｏｗｅｎ的剪切修正系数，正确计及了叠层复合材料壳体的横向剪切效应；为了避免“自锁”现象，文中采用了９节点的Ｈｅｔｅｒｏｓｉｓ二次壳单元；特别是本文利用插值外推的思想，提出了一个带预测的弧长增量控制法，显著提高了确定变形路径的计算效率。几个数值算例表明本文给出的有限元方法对于各向异性中厚度板壳的弹塑性分析有较好的精度，尤其是对具有复杂变形路径的结构计算，收敛速度提高更快。     

弹塑性杂交/混合有限元法与实施

本文提出两种结构弹塑性分析的有限元方法——杂交/混合非协调元法。该法采用场变量分解,文[1]对增量形式的杂交应力元能量泛函进行修正,简化了高阶矩阵的求逆运算,从而提高了多变量非线性有限元分析的计算效率和收敛精度。文中按文[2]给出的第二种能量泛函的修正形式建立了弹塑性平面问题中的杂交/混合非协调四边形单元。最后给出算例,并与实验解及文[3]解进行比较。表明该方法分解弹塑性问题是十分有效的。     

混凝土材料动态本构参数的分阶段计算反求技术

基于不同应变率下的分离式Hopkinson压杆实验,提出了一种混凝土材料动态本构参数的分阶段反求方法。该分阶段反求方法中,结合不同的实验模型对未知参数进行敏感性分析,并依据敏感程度对参数进行分类,每类参数对应于1组实验,再利用反求方法分阶段确定各类参数,在每一类参数反求时均将上一次反求结果作为已知条件。结果表明,该分阶段反求方法充分利用所有实验的测量响应信息,能快速获取混凝土本构中较难确定的关键参数,为参数的确定提供了一种有效的途径。     

基于有限元法-正则化的弹性模量反求算法研究

生物组织的病变往往导致其物理特性的变化,研究生物组织的物理特性变化对诊断生物组织病变与否能够提供参考.本文提出了根据已知生物组织部分边界位移信息、外部载荷、边界条件和有限元模型,利用正则化方法反求弹性模量的方法.首先将求解弹性模量问题描述成一个力学反问题,建立逆系统辨识的数学模型,并应用正则化原理建立生物组织弹性模量反求算法.最后利用仿真算例验证了算法的有效性,并结合图像识别技术试验获取对象部分边界位移信息,反求其弹性模量,证明了本文反求方法的可行性.     

蚁群算法求解二维拉压不同模量反问题

利用光滑函数技术对二维拉压不同模量本构关系进行光滑化处理,采用初应力方法求解二维拉压不同模量正问题的有限元方程。在此基础上,建立了基于连续域蚁群算法的二维拉压不同模量反问题的数值求解模型,考虑了区域非均质的影响,实现了对拉压弹性模量和泊松比的单一/组合识别。通过两个数值算例,对所提算法进行了数值验证,分别探讨了蚁群算法相关参数、测点分布和数据噪音等对识别结果的影响。数值验证表明,所提算法可有效地求解二维拉压不同模量反问题,并具有较好的计算精度。     

Non-probabilistic information fusion technique for structural damage identification based on measured dynamic data with uncertainty

Based on measured natural frequencies and acceleration responses,a non-probabilistic information fusion technique is proposed for the structural damage detection by adopting the set-membership identification(SMI) and twostep model updating procedure.Due to the insufficiency and uncertainty of information obtained from measurements,the uncertain problem of damage identification is addressed with interval variables in this paper.Based on the first-order Taylor series expansion,the interval bounds of the elemental stiffness parameters in undamaged and damaged models are estimated,respectively.The possibility of damage existence(PoDE) in elements is proposed as the quantitative measure of structural damage probability,which is more reasonable in the condition of insufficient measurement data.In comparison with the identification method based on a single kind of information,the SMI method will improve the accuracy in damage identification,which reflects the information fusion concept based on the non-probabilistic set.A numerical example is performed to demonstrate the feasibility and effectiveness of the proposed technique.     

高维非线性振动系统参数识别

将增量谐波平衡非线性识别推广到高维振动系统, 推导了基于增量谐波平衡的多自由度非线性系统的识别方程. 针对一个两自由度系统进行了数值模拟计算, 讨论了系统在单周期、倍周期和混沌运动状态下的参数识别, 以及噪声对识别结果的影响, 验证了增量谐波平衡非线性识别在多自由度系统中的有效性. 结果表明, 该方法具有较高的计算效率和识别精度, 以及良好的抗噪能力.      

悬臂梁中开闭裂纹的损伤识别

利用裂纹诱导弦挠度函数,建立了悬臂Euler-Bernoulli中开闭裂纹位置、深度、初始张开角等损伤参数的识别方法.为此,首先将梁中开闭裂纹等效为单向扭转弹簧,给出了考虑裂纹缝隙效应的裂纹梁等效抗弯刚度,并得到悬臂Euler-Bernoulli开闭裂纹梁弯曲挠度的显式闭合解及裂纹诱导弦挠度函数,证明了裂纹诱导弦挠度的...     

Identification of the Ten Inertia Parameters of a Rigid Body

Standard experiments for identifying inertia parameters of a rigid body only provide a subset of the inertia parameters of the body [1–10]. In addition, they do not use in the estimation process the complete information included in the equations of motion of the rigid test body. The objective of the work described in this paper is the simultaneous, automatic experimental identification of the ten inertia parameters of a rigid body using the complete information hidden in the nonlinear model equations of the test body. This task has been solved in several steps:– mathematical modelling of the special motions of a rigid body in space. These model equations have been mapped into a form suitable for identification purposes (identification hypothesis)– design of a special measurement robot for performing the identification experiments– laboratory experiments providing test data used for the identification experiments– identification of the inertia parameters and accuracy tests.The accuracy of the identified parameters is satisfactory.     

基于模态数据的长细结构小损伤精确识别研究

本文针对现有的损伤识别方法不能满足部分结构损伤识别精度要求的现状,对结构的小损伤精确识别方法开展研究．以长细结构为研究对象,对具有不同损伤位置和损伤程度的圆柱形的轻阻尼梁结构进行了数值分析和实验研究,应用数值计算方法和实验确定的特征向量和特征频率对长细结构裂缝参数进行识别计算．本文在研究过程中编制了一个创新性的预测程序,通过其一次性生成目标函数图来选择合适的初始参数,从而对识别结果进行分析．研究结果表明,应用本文提出的识别方法,裂缝位置的识别误差可以控制在0.05 %~0.28 %范围内,裂缝深度识别误差低于7 %．     

基于子结构解耦的连接特性识别新方法

连接部件动态特性的准确辨识对预测装配式机械结构的动力学行为有重要意义. 针对传统基于子结构解耦的连接结构动力学特性识别方法难以直接利用可测量数据进行辨识及辨识结果受噪声影响显著等问题, 本文提出了一种新方法. 首先, 提取子结构解耦基本方程在测试自由度上的分量, 并经矩阵变换得到显含连接动刚度矩阵的形式, 而后由真实连接动刚度矩阵分解为已知的初始矩阵与待求的增量矩阵, 推导了具有收敛性质的增量型方程以增强界面自由度较多时辨识的数值稳定性, 并采用多项式拟合动刚度将其转化为了拟合系数的频域估计方程, 按给定准则选取合适的频率点联立方程后, 得到了只需装配体测试自由度上的频响函数来辨识连接特性的迭代公式. 最后, 以若干算例说明了算法的具体流程. 对10自由度弹簧?质量块系统进行了数值仿真, 验证了所提方法的正确性及抗噪性; 对包含一处胶接连接的T形梁结构和包含两处螺栓连接的L形梁结构进行了试验, 所辨识连接结构与残余结构重组的装配体有限元模型计算的频响函数与测量值在较宽频带内吻合较好, 表明了该方法能有效识别实际装配体结构中的连接特性.      

A Newton type iterative method for heat-conduction inverse problems

An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.     

A Newton type iterative inverse method for heat-conduction problems

An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.     

In-situ mechanical properties identification of 3D particulate composites using the Virtual Fields Method

This paper presents an identification procedure based on the Virtual Fields Method (VFM) for identifying in situ mechanical properties of composite materials constitutive phases from 3D full-field measurements. The new procedure, called the Regularized Virtual Fields Method (RVFM), improves the accuracy of the VFM thanks to the imposition of mechanical constraints derived from an appropriate homogenization model. The developed algorithms were validated through virtual experiments on particulate composites. The robustness of both the VFM and the RVFM was assessed in the presence of noisy strain data for various microstructures. A study was also carried out to investigate the influence of the size of region of interests on the reliability of the identified parameters. Accordingly, the optimum size of region of interest was determined based on full-field measurement requirements and accuracy of the identified parameters. This study enables determining, a priori, the required magnification level of 3D images for composites of any mechanical and morphological characteristics.