一种陶瓷脆性材料动态本构参数的计算反求方法

提出了一种利用反射波和透射波响应反求陶瓷脆性材料动态材料本构参数的方法.通过数值模拟分析哑铃形试件在加载和测试中能有效地消除动态应力波的应力集中,从而得到有效的动态压缩强度,因此该文利用哑铃形得到的可靠反射波和透射波响应做反求输入信息.应用遗传算法,根据反射波和透射波响应反求材料本构参数,避免了应力波效应和应变率效应解耦的问题,不必分离结构响应和材料响应,从而快速地获取参数.     

一种脆性材料动态本构损伤参数的计算反求方法

提出一种利用平板撞击的瞬态速度剖面响应反求陶瓷脆性材料的动态本构损伤参数的方法.通过数值模拟再现平板撞击试验中的正平面波在飞片、样品和窗口材料中的传播情况,分析了"样品/窗口"界面瞬态速度剖面响应与材料参数之间的映射关系,进而建立了用响应反求本构参数的反问题,且通过一些数值试验验证该反问题是适定的.基于平板撞击速度剖面响应反求的方法,为快速地获得高应变率下一些传统方法难以测定出的本构参数提供新途径.     

确定地基固结过程中材料参数的反分析法

结合有限单元法和数学规划法提出一确定地基固结过程中材料参数的反分析法。该法可以根据固结过程中的位移值和孔隙水压力值反算地基土的弹性系数、泊松比以及渗透系数。文中通过一些算例对该法进行了验证。     

JC本构参数反求中的自适应微型遗传算法

提出自适应微型遗传算法,使种群的交叉和变异概率可以根据适应度的大小自我调节,该算法经过了测试函数的测试。本文以JC(Johnson-Cook)本构模型描述的高导无氧铜材料为例,集成运用该算法和有限元方法来反求JC动态本构参数。结果表明,该方法提高了收敛速度和搜索效率,能够快速地获取参数。     

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

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

基于反求的医疗辅助弹性护具应变能的计算

提出了基于反求工程的医疗辅助弹性护具的设计方法,推导出弹性护具应变能的计算公式. 设计中应用网格简化技术减少扫描数据以及参数化技术将三维转化为二维进行计算,减少了网格点应变能的计算量,提高了设计效率. 应用应变能密度函数建立起弹性护具几何尺寸与应变能之间的关系,为医生或患者科学地选择医疗护具提供了依据.     

基于最小二乘的工程结构边界载荷反求误差消减方法

机械装备工程结构强度计算和疲劳寿命预测等都需要准确获知结构的工作或边界载荷。然而,由于土壤和砂石等工作介质的切削及铲运等精确阻力模型难以建立,以及机械装备动力传动过程中的摩擦等损耗难以精确计算,导致准确获得装备工作结构的边界载荷非常困难。依据圣维南原理,采用有限元法建立复杂结构中局部结构应变与其边界载荷的关系,对边界载荷反求过程中误差来源和范围进行了分析和估计,提出了基于最小二乘法的复杂结构局部边界载荷误差消减方法。对推耙机横梁结构的边界载荷进行分析,验证了该误差消减方法的有效性。     

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

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

一种基于特征值反求的复合材料层合板波动特性的近似修改方法

采用一种特征值反求形式的方法用来近似层状复合板的波动特性修改.通过敏感度分析,确定层合板中能最显著地改变结构波动特性的结构参数为设计变量.采用二阶泰勒展式近似特征值和特征频率相对于设计变量的改变率.计算了不同结构变化时的色散特性,并与精确计算结果进行了比较.结果证明本方法对于层状复合板波动特性修改的近似是可行并且有效的.     

含能材料药型罩的爆炸成型及毁伤作用

为了研究由含能材料制备的药型罩爆炸成型过程及其对目标靶的终点效应,设计了球缺形药型罩的装药结构,放置与药型罩曲率相同的缓冲垫在药型罩和主装药之间,运用高速摄影系统拍摄含能材料药型罩的成型过程。实验结果表明,含能材料药型罩在爆炸作用下能够形成弹丸,弹丸速度2km/s左右。含能弹丸穿透20mm厚的装甲钢靶后反应加剧,形成大量的气体。侵彻过程含动能和化学反应的综合作用,穿孔有明显的烧蚀现象,穿孔口径0.5 D,最大穿深1.4 D。利用含能材料制备成药型罩可以实现炸药的直接驱动,这可为含能材料战斗部的工程应用提供参考。

     

Investigation of advanced strain-path dependent material models for sheet metal forming simulations

Sheet metal forming processes often involve complex loading sequences. To improve the prediction of some undesirable phenomena, such as springback, physical behavior models should be considered. This paper investigates springback behavior predicted by advanced elastoplastic hardening models which combine isotropic and kinematic hardening and take strain-path changes into account. A dislocation-based microstructural hardening model formulated from physical observations and the more classical cyclic model of Chaboche have been considered in this work. Numerical implementation was carried out in the ABAQUS software using a return mapping algorithm with a combined backward Euler and semi-analytical integration scheme of the constitutive equations. The capability of each model to reproduce transient hardening phenomena at abrupt strain-path changes has been shown via simulations of sequential rheological tests. A springback analysis of strip drawing tests was performed in order to emphasize the impact of several influential parameters, namely: process, numerical and behavior parameters. The effect of the two hardening models with respect to the process parameters has been specifically highlighted.     

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.     

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.     

覆盖件冲压仿真中的网格形态优化技术

覆盖件冲压仿真计算模型中网格密度分布的合理性与网格单元形态的优劣,对仿真结果的准确性有很大的影响.提出了一种改进的全四边形网格细分方法,使网格的密度分布适于覆盖件冲压分析计算要求,可保证细分后网格的协调性,并将算法推广以处理非结构化四边形网格和三角形四边形混合网格的细分.提出的网格细分策略,有助于提高细分后网格的质量.提出了适用于细分后四边形网格和非结构四边形网格的拓扑形态优化操作,可有效的提高网格模型的形态质量.     

Numerical study for influences of material parameters on hydrostatic bulging of metal sheet

In this paper, the influences of various material parameters, the hardening exponent (n), the rate sensitivity (m). the thickness anisotropy parameter (R) and the index M in the Hosford and Hill yield function, on the hydrostatic bulging of a circular clamped sheet of ductile metal materials are analysed by introducting a rigid-viscoplastic finite element method. By numerical studies, an empirical relationship within the average limit thickness strain – ₃ ^* and the material parameters (n andm) is obtained. Besides, it has been found that the influences of surface shapes of the yield function on the average limit thickness strain can be reflected by the Barlat'sP value which represents the effects ofR andM values.This work is supported by the National Natural Science Foundation and Natural Science Foundation of the Youth of China.     

基于有限元分析的特征值反问题求解的逆摄动方法

本文研究特征值反问题的求解方法,根据广义特征值反问题理论和有限元法的特点,以转子系统平面梁单元有限元模型结构分析的特征值反问题求解为例,给出一种新的逆摄动方法,给出了本逆摄动法较完整的理论基础,给出了其逆摄动参数的显式计算公式及相应的取值方法,本逆摄动法也可推广到其他单元类型的有限元模型特征值反问题的求解。     

A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming

In this paper an anisotropic material model based on non-associated flow rule and mixed isotropic–kinematic hardening was developed and implemented into a user-defined material (UMAT) subroutine for the commercial finite element code ABAQUS. Both yield function and plastic potential were defined in the form of Hill’s [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. R. Soc. Lond. A 193, 281–297] quadratic anisotropic function, where the coefficients for the yield function were determined from the yield stresses in different material orientations, and those of the plastic potential were determined from the r-values in different directions. Isotropic hardening follows a nonlinear behavior, generally in the power law form for most grades of steel and the exponential law form for aluminum alloys. Also, a kinematic hardening law was implemented to account for cyclic loading effects. The evolution of the backstress tensor was modeled based on the nonlinear kinematic hardening theory (Armstrong–Frederick formulation). Computational plasticity equations were then formulated by using a return-mapping algorithm to integrate the stress over each time increment. Either explicit or implicit time integration schemes can be used for this model. Finally, the implemented material model was utilized to simulate two sheet metal forming processes: the cup drawing of AA2090-T3, and the springback of the channel drawing of two sheet materials (DP600 and AA6022-T43). Experimental cyclic shear tests were carried out in order to determine the cyclic stress–strain behavior and the Bauschinger ratio. The in-plane anisotropy (r-value and yield stress directionalities) of these sheet materials was also compared with the results of numerical simulations using the non-associated model. These results showed that this non-associated, mixed hardening model significantly improves the prediction of earing in the cup drawing process and the prediction of springback in the sidewall of drawn channel sections, even when a simple quadratic constitutive model is used.     

基于泰勒级数展开的区间反演方法

实际结构或构件的几何与材料参数总包含不确定性,在对结构计算模型进行精确分析时,有时需要对参数不确定性进行量化。本文提出了一种用于区间参数识别的反演方法,即基于泰勒级数展开式分别建立参数与响应的区间中值、区间半径的对应函数关系,并通过构建两个反演问题来分步识别参数区间中值和半径,以避免区间扩张现象和简化优化反演过程。通过数值质-弹系统初步验证了方法的可行性,然后基于一组钢板的动测数据,识别了钢板的几何及材料特性参数的区间范围。研究结果表明,本文方法具有良好的区间反演精度,能有效地避免区间扩张现象,可以用于实际工程区间问题的求解。     

An inverse algorithm for the identification and the sensitivity analysis of the parameters governing micropolar elasto-plastic granular material

The article concerns the complex determination process of the material parameters governing micropolar granular material with elasto-plastic material properties. Proceeding from a gradient-based method, we split the total set of parameters and the overall identification procedure into two major categories. These are, firstly, the identification of the parameters of a standard non-polar elasto-plastic continuum, and, secondly, the determination of the remaining parameters governing the micropolar part of the constitutive model. While the first set of parameters can be obtained from homogeneous triaxial tests on, e. g., granular, cohesive-frictional materials like sand, the second set can only be determined from inhomogeneous tests, such as biaxial tests including the onset and the development of shear bands. Following this, one can obtain the first part of the identification process from a simple inverse algorithm applied to the elasto-plastic material model of non-polar solids, while the second part requires a fully inverse computation in the sense of a back analysis of the underlying boundary-value problem. In the present article, this procedure is carried out by use of the semi-discrete sensitivity analysis. Finally, the whole model is applied to the data of Hostun sand taken at the universities of Grenoble and Stuttgart. Dedicated to Professor Franz Ziegler on the occasion of his 70th birthday.     

A new method for circular grid analysis in the sheet metal forming test

Computer vision systems are employed to determine the major and minor lengths of deformed elliptic grids while determining a sheet metal's workability. The existing method identifies the ellipse using the least squares analysis. It suffers two drawbacks: assumptions in direct conflict with the observed real-world processes and an undesirable property of orientation dependence. For the remedy, this paper presents a new method that, in addition to achieving the desired property of orientation invariance, discards assumptions that conflict with real-world processes. The proposed method is implemented and tested using simulated and real-world data. Results are reported and compared with those obtained by the existing method.