厚壁球壳冲压成形过程的热力耦合分析

在高温下成形的金属，工件与模具、外界环境之间存在剧烈的热交换，同时塑性变形功及摩擦消耗功又转化成热能，这些都引起变形体内温度的急剧变化。在金属塑性成形过程中，当温度和应变速率对材料的屈服流动应力产生明显影响时，或高温下成形的金属，应力与应变关系必须采用粘塑性本构模型。采用热力耦合粘弹塑性模型对厚壁球壳的热拉伸成形过程进行了分析，模型如图1所示。毛坯材料为LD10锻件，凸、凹模材料用5CrNnMo，材料加热温度为470℃，并经过充分保温，凸、凹模初始温度均为120℃。润滑采用油基石墨，     

喷射成形FGH4095的热变形特征

研究了喷射成形高温合金FGH4095在温度为1010—1140℃、应变速率为0.001—10.0s^-1下的热变形行为,并使用双曲正弦模型建立了喷射成形FGH4095的本构方程.结果表明:本构方程预测结果与实验值符合较好,为喷射成形FGH4095热加工过程的数值模拟提供了较准确的模型.     

冲击载荷下混凝土本构模型构建研究

 在对混凝土动态力学性能和现有本构模型综合分析的基础上,构建了一个新的适用于冲击响应问题数值分析的混凝土本构模型。该本构模型全面考虑了压力、应力第三不变量、变形的硬化和软化、应变率强化以及拉伸损伤等各个影响因素。将其加入LTZ-2D程序,确定了本构模型参数,对混凝土靶板的穿透问题进行了数值验证分析。计算得到的弹体剩余速度同实验结果基本一致,同时得到了混凝土靶板破裂的计算图像。计算结果及其分析表明,所构建的本构模型能够较好地反映冲击载荷作用下混凝土动态响应的主要特性。     

粘弹性减阻流动的雷诺应力模型研究

针对由Giesekus本构方程描述的粘弹性槽道减阻流动,建立了低Re数雷诺应力模型.通过与DNS数据的比较,该模型对平均速度、平均变形、减阻率和应力的预测较准确,而对时均速度剖面和脉动速度强度的预测精度有待提高.另外,对雷诺应力再分配项(压强变形率相关项)的模化有待改进.     

材料参数拟合方法对弹靶侵彻仿真的影响

弹靶侵彻仿真中材料参数对计算结果有着至关重要的影响。为寻求一套适用于弹靶侵彻仿真计算的材料参数拟合方法,借助前期开展的靶板材料动态力学性能试验、靶板材料断裂试验,通过不同拟合方法依次得到不同的JC本构模型及失效模型参数,依据试验建立有限元计算模型,将数值计算结果与试验结果进行对比。结果表明:(1)对于同一材料的力学性能试验,采用不同的拟合方法可得到不同的JC本构、JC失效参数,二者会对弹靶仿真结果造成一定影响;(2)在不考虑温度软化项的前提下,采用高应变率作为参考应变率进行拟合能更加准确地表征材料在高应变率下的应力-应变关系,更加适用于弹靶侵彻强瞬态、高应变率作用过程仿真;(3)对于同一JC本构模型,采用平均应力三轴度拟合的JC失效模型较采用初始应力三轴度拟合的JC失效模型所得战斗部剩余速度计算结果偏小,仅采用拉伸试件结果拟合的JC失效模型较采用扭转、拉伸试件结果拟合的JC失效模型所得战斗部剩余速度计算结果偏小。     

单晶钨纳米线拉伸变形机理的分子动力学研究

利用分子动力学方法, 对本课题组率先采用金属催化的气相合成法制备出的高纯度单晶钨纳米线进行拉伸变形数值模拟, 通过分析拉伸应力-应变全曲线及其微观变形结构, 揭示出单晶钨纳米线的拉伸变形特征及微观破坏机理. 结果表明: 单晶钨纳米线的应力-应变全曲线可分为弹性阶段、损伤阶段、相变阶段、强化阶段、 破坏阶段等五个阶段, 其中相变是单晶钨纳米线材料强化的重要原因; 首次应力突降是由于局部原子产生了位错、孪生等不可逆变化所致; 第二次应力突降是发生相变的材料得到强化后, 当局部原子再次产生位错导致原子晶格结构彻底破坏而形成裂口、且裂口不断发展成颈缩区时, 材料最终失去承载能力而断裂. 计算模拟得到的单晶钨纳米线弹性模量值与实测值符合较好. 关键词： 分子动力学 应力应变曲线 微观机理 单晶钨纳米线     

简论金属材料JC本构模型的精确性

通过比较JC模型预测结果与6种金属(2024-T351铝合金、 6061-T6铝合金、OFHC无氧铜、4340高强钢、Ti-6Al-4V钛合金和Q235软钢)在不同应变率及温度下的实验数据,对JC本构模型的精确性进行了关键评估。为了进一步评估其精确性,采用JC本构模型和失效准则对平头弹正撞2024-T351铝合金靶板进行数值模拟,并与实验结果比较。结果表明:JC本构模型只适用于中、低应变率和温度下的Mises材料,对非Mises材料该模型预测的剪切应力-应变曲线和失效与实验结果吻合较差;同时,JC本构模型的精度随应变率和温度的提高而降低,特别是在高应变率条件下利用实验得到的动态增强因子进行相应数值模拟时,所得计算结果与弹道穿透实验结果不一致,说明其表达式(即准静态应力-应变关系×动态增强因子)是不恰当的。     

退火铝合金中Portevin-Le Chatelier效应的数值模拟研究

通过综合考虑沉淀动力学和动态应变时效机理，建立了一套具有明确物理内涵的唯象本构模型.并在对不同加载应变率下的单轴拉伸实验的数值模拟中，得到了3种类型的Portevin-Le Chatelier效应应力曲线.计算结果与实验数据的较好吻合，验证了该模型的有效性. 关键词： Portevin-Le Chatelier效应 动态应变时效 数值模拟     

高温高应变率下G550冷弯钢的动态力学行为

为研究G550冷弯钢在高温和高应变率下的动态力学性能,采用高温同步控制霍普金森拉杆装置,开展了不同温度下的高应变率拉伸试验,并在高速液压拉伸试验机上进行了室温下的中应变率拉伸试验。通过获得的应力-应变曲线,得到了材料的本构模型,结合微观形貌分析,探究了温度和应变率对流变应力的影响。结果表明：G550冷弯钢具有明显的应变率强化和温度软化效应。在特定的高应变率范围内(1 000～1 500 s^-1),温度对流变应力的影响大于应变率。基于温度软化系数随温度变化的特征,提出了G550冷弯钢的修正Johnson-Cook本构模型。该模型可以较好地描述G550冷弯钢在高温和高应变率下的动态力学行为,从而为G550冷弯钢在高温、爆炸冲击相关的有限元仿真提供参考。     

含扭转晶界位错Al金属拉伸强度第一性原理预测

本文基于密度泛函理论第一原理方法,从影响力学性能本质的电子结构计算上,对含Σ 5{001}扭转晶界位错Al金属拉伸强度进行了预测,发现其理论拉伸强度达到8.73 GPa,临界应变为 24%.拉伸强度低于文献报道(Phys. Rev. B 75, 174101 (2007))的倾斜晶界位错Al金属的理论拉伸强度9.5 GPa,但其临界应变却远大于倾斜晶界的16%.本研究结果表明,通过工艺参数控制,改变缺陷形态,可极大地改变其力学性能.进一步地,从电子结构层次上, 分析了含晶界位错Al金属拉伸断裂行为的实质,通过分析电荷密度分布、键长变化等,发现其断裂处发生在晶界处;理论计算结果将对Al金属结构设计及力学性能改善具有重要的指导作用.     

基于无网格点插值法的旋转悬臂梁的动力学分析

将基于多项式点插值的无网格方法用于旋转悬臂梁的动力学分析. 利用无网格点插值方法对柔性梁的变形场进行离散, 考虑梁的纵向拉伸变形和横向弯曲变形, 并计入横向弯曲变形引起的纵向缩短, 即非线性耦合项, 运用第二类Lagrange方程推导得到系统刚柔耦合动力学方程. 与有限元法相比, 该方法只需节点信息, 无需定义单元, 具有前处理简单的优势; 构造的形函数采用更多的节点插值, 具有高阶连续性. 将无网格点插值方法的仿真结果与有限元和假设模态法进行比较分析, 验证了该方法的正确性, 并表明其作为一种柔性体离散方法在刚柔耦合多体系统动力学的研究中具有可推广性.     

Accurate free vibration analysis of thick laminated circular plates with attached rigid core

This paper deals with the free vibration behavior of laminated transversely isotropic circular plates with axisymmetric rigid core attached at the center. The governing equations of motion are obtained based on Mindlin's first-order shear deformation plate theory. Two possible categories of vibration modes related to up-down translation of the core and wobbly rotation of the core about a diameter are studied. Accurate natural frequencies hitherto not reported in the literature are presented for a wide range of thickness-to-radius ratio, inner-to-outer radius ratio, mass and moment of inertia ratios of the core and various boundary conditions at the outer edge of the plate. Numerical results are compared with those of a three-dimensional finite element method (3-D FEM) to demonstrate the high accuracy and reliability of the current analysis.     

Thermal effects in high-power Nd:YAG disk-type solid state laser

An investigation of thermal effects in a high-power Nd:YAG disk-type solid state laser pumped with different pump beam transverse profiles is carried out by numerical simulation based on the finite element method (FEM). Impact of the heat sink on the thermal effects is included in the simulation. The distribution of first principle stress, thermally induced birefringence, including the distribution and variation of the birefringence loss, are studied. The characteristics of the phase variation are analyzed with consideration of the temperature gradient, deformation, strain and thermal stress. Thermal lensing is explored as a function of pump power and of the radius pumped with different pump beam transverse profiles. The non-parabolic part of optical phase distortion is simulated. Furthermore, the characteristics of the bi-focus of the disk laser are also studied. Experiments on the maximum tensile stress distribution and depolarization loss are carried out. The presented calculations are in qualitative agreement with the experimental observations.     

激光辐照下充压圆筒变形的相似律问题

 以充压圆筒在强激光辐照下的变形与破坏响应为研究目标，首先利用相似分析方法，主要采用了方程分析法和量纲分析法，确定了相似准则。然后利用ANSYS有限元程序，模拟了三组几何和物理相似的大小模型。通过对计算结果的分析表明，满足相似准则的大小模型，在相应的时间和空间，相对的物理参量是一致的。     

Thermal postbuckling and vibration analyses of functionally graded plates

Thermal postbuckling and vibration behaviors of the functionally graded (FG) plate are investigated. The material properties of the FG plate are assumed to vary continuously through the thickness of the plate and as temperature with the nonlinearity. The nonlinear finite element equations based on the first-order shear deformation plate theory are formulated for the FG plate. The von Karman nonlinear strain–displacement relationship is used to account for the large deflection of the plate. The incremental form considering the initial displacement and initial stress is adopted for the nonlinear temperature-dependent material properties of the functionally graded material. The numerical result shows the characteristics of the thermal postbuckling and vibration of the FG plate in the pre- and post-buckled regions.     

Analytic description of the fusion and fission processes through compact quasi-molecular shapes

Recent studies have shown that the characteristics of the entrance and exit channels through compact quasi-molecular shapes are compatible with the experimental data on fusion, fission and cluster radioactivity when the deformation energy is determined within a generalized liquid drop model. Analytic expressions allowing to calculate rapidly the main characteristics of this deformation path through necked shapes with quasi-spherical ends are presented now; namely formulas for the fusion and fission barrier heights, the fusion barrier radius, the symmetric fission barriers and the proximity energy.     

大展弦比气弹机翼柔性非线性变形的气动力效应分析与风洞试验

从飞行器刚弹耦合动力学模型出发,引入柔性机翼准定常假设,建立大柔性飞行器非线性静气动弹性气动力方程,利用非线性迭代求解思路模拟了柔性飞行器的静气动弹性响应行为,开展了大展弦比飞机静气动弹性风洞试验验证,采用气动力有限基本解与机翼的耦合计算,发现了大柔性飞机大变形状态下载荷及结构变形形式随风速的变化规律.传统基于小变形假设的线性分析方法和刚体分析由于无法考虑气动面随结构变形的曲面气动力因素和结构变形后的非线性刚度特性,均与风洞试验存在一定的误差.对于大展弦比柔性飞机的非线性静气动弹性分析十分必要.      

Single-step scanner-based digital image correlation (SB-DIC) method for large deformation mapping in rubber

The incremental digital image correlation (DIC) method has been applied in the past to determine strain in large deformation materials like rubber. This method is, however, prone to cumulative errors since the total displacement is determined by combining the displacements in numerous stages of the deformation. In this work, a method of mapping large strains in rubber using DIC in a single-step without the need for a series of deformation images is proposed. The reference subsets were deformed using deformation factors obtained from the fitted mean stress-axial stretch ratio curve obtained experimentally and the theoretical Poisson function. The deformed reference subsets were then correlated with the deformed image after loading. The recently developed scanner-based digital image correlation (SB-DIC) method was applied on dumbbell rubber specimens to obtain the in-plane displacement fields up to 350% axial strain. Comparison of the mean axial strains determined from the single-step SB-DIC method with those from the incremental SB-DIC method showed an average difference of 4.7%. Two rectangular rubber specimens containing circular and square holes were deformed and analysed using the proposed method. The resultant strain maps from the single-step SB-DIC method were compared with the results of finite element modeling (FEM). The comparison shows that the proposed single-step SB-DIC method can be used to map the strain distribution accurately in large deformation materials like rubber at much shorter time compared to the incremental DIC method.     

剪切流动条件下液滴变形和断裂的数值模拟

本文采用扩散界面法研究了剪切流动条件下悬浮液滴变形和断裂的动力学机制。控制方程采用考虑表面张力影响的Navier-Stokes-Cahn-Hilliard方程描述。计算网格采用均匀矩形交错网格。采用基于压力增量的近似投影法计算 Navier-Stokes方程,采用完全近似多重网格法计算Cahn-Hilliard方程。稳态液滴变形规律及液滴拉伸断裂过程的计算结果与试验结果符合较好,表明本文模型能够很好的研究液滴变形及断裂机理。     

Analysis of ultrasonic attenuation in particle-reinforced plastics by a differential scheme

Attenuation of ultrasonic longitudinal waves in some particle-reinforced polymer composites is studied theoretically by a micromechanical model based on a differential (incremental) scheme. A set of differential equations is established by which the attenuation spectrum of the composite can be computed from the known properties of viscoelastic matrix and elastic particles. For a composite reinforced with glass particles with radius 0.15 mm, the proposed scheme is shown to predict the attenuation in better agreement with the foregoing experimental results than the previous simplistic independent scattering model. Based on this scheme, the dependence of the longitudinal attenuation spectrum of a particulate polymer composite on the wavelength-to-particle radius ratio and the particle volume fraction is examined in detail. It is then shown theoretically that the attenuation of the composite decreases monotonically with the particle volume fraction when the particle radius is sufficiently small compared to the incident wavelength, while it shows non-monotonic particle-fraction dependence when the ratio of the particle radius to the wavelength is larger. To examine this theoretical finding from an experimental point of view, the longitudinal attenuation in a glass-particle-reinforced polyester composite with particle radius 0.0225 mm is measured for different particle volume fractions. The measured attenuation characteristics are shown to support the qualitative features of the theoretical prediction.