弹性圆板的热过屈曲

基于精确的圆薄板轴对称大挠度变形几何方程，建立了均匀加热圆板轴对称热弹性屈曲问题位移形式的控制方程，采用打靶法和解析延拓法获得了连续依赖于温度载荷的圆板过屈曲状态解，给出了相应的数值结果。     

基于有限差分法薄板激光冲击响应的数值模拟

在激光冲击载荷作用下, 薄板变形速度快, 诱导产生的应力波的传播过程较为复杂. 传统的测量工具难以对薄板变形过程中的动态响应进行有效的测量. 本文采用理论与实验相结合的方法, 构建了薄板在激光冲击下二维轴对称平面模型, 建立其拉格朗日运动方程, 利用有限差分法求其显式解, 分析薄板在激光冲击载荷作用下薄板的变形过程和应力波的传播过程, 并研究不同工艺参数对薄板动态响应特性的影响. 结果表明, 薄板变形初期的速度为振荡式增加, 在快速的拉胀式变形过程中会出现明显的回弹现象, 在光斑边界处产生向内和向外传播的应力波, 载荷的压力-空间分布以及边界约束条件也对薄板的动态响应结果有显著的影响. 激光冲击实验得到的结果与数值结果和预测结果基本吻合. 研究方法与所得结论可为薄板激光冲击成形过程中的参数优化提供参考.      

弹性半空间地基与四边自由异性矩形板相互作用的研究

选用弹性半空间地基模型分析四边自由各向异性矩形地基板的弯曲和稳态振动解析解。将异性薄板控制微分方程与基于弹性半空间地基位移解建立的板与地基变形协调方程相结合,先按对称性分解,然后采用三角级数法得出了弹性半空间地基上四边自由各向异性矩形薄板的弯曲和稳态振动解析解,包括地基反力(幅值)、板的挠度(幅值)、板的内力(幅值)的解析表达式。克服了数值法的弊端,取消了对地基反力的假设,得到了板的内力(幅值)及地基反力(幅值)更切实际的分布规律。算例结果不但与文献结果吻合良好,而且表明对于异形板,对称载荷能引起反对称的内力和变形。该方法使得半空间地基上各向异性矩形薄板这一复杂的接触问题的求解统一化、简单化、规律化。     

弹性半空间地基与四边自由异性矩形板相互作用的研究

选用弹性半空间地基模型分析四边自由各向异性矩形地基板的弯曲和稳态振动解析解。将异性薄板控制微分方程与基于弹性半空间地基位移解建立的板与地基变形协调方程相结合,先按对称性分解,然后采用三角级数法得出了弹性半空间地基上四边自由各向异性矩形薄板的弯曲和稳态振动解析解,包括地基反力(幅值)、板的挠度(幅值)、板的内力(幅值)的解析表达式。克服了数值法的弊端,取消了对地基反力的假设,得到了板的内力(幅值)及地基反力(幅值)更切实际的分布规律。算例结果不但与文献结果吻合良好,而且表明对于异形板,对称载荷能引起反对称的内力和变形。该方法使得半空间地基上各向异性矩形薄板这一复杂的接触问题的求解统一化、简单化、规律化。     

基于变形分解的四边固支单向与双向薄板界定分析

板是工程结构中的常用构件,但需区分单向板和双向板分别进行设计,故对两者的界定研究具有重要意义。本文通过构造8节点正方体单元24种基本位移与变形向量,形成用于变形分解的正方体单元完备正交基矩阵,及均布荷载作用下板的主要基本变形分解图。基于变形分解分析,总结出单向板三大变形特征,即短边弯曲变形为主特征、浅梁变形特征和长边弯曲变形为可略特征;进而确定以长宽比1.50、2.20、2.40作为不同精度下的单、双向薄板类型划分界限。工程设计人员可依据具体工程需求,采用不同的分类方法对任意薄板进行类型判别,进而采用单向薄板或双向薄板理论进行薄板设计。     

箱型结构的解析分析

从弹性薄板理论出发,将箱型结构分为六块相互约束的简支薄板。每一块板在板内承受任意局部荷载,在四边承受待定的分布弯矩。通过板与板的边界转角协调分析,求解出每块板四边所受的分布弯矩,从而得到任意荷载作用下箱型结构的解析解。     

板材多点成形过程的有限元分析

多点成形过程采用静力隐式格式进行数值模拟是比较合适的。本文建立了用于多点成形过程分析的静力隐式弹塑性大变形有限元方法 ,给出了对稳定迭代收敛过程效果较好的板壳有限单元模型、处理多点不连续接触边界的接触单元方法以及增量变形过程中应力及塑性应变计算的多步回映计算方法。基于这些方法编制了计算软件 ,应用该软件进行了矩形板的液压胀形过程及球形模具拉伸成形过程的有限元分析 ,数值计算结果与典型的实验结果及计算结果吻合很好。最后给出了球形、圆柱形目标形状的实际多点成形过程的数值模拟结果。     

功能梯度矩形板的非线性自由振动

研究了功能梯度矩形薄板的非线性自由振动问题.采用幂律分布规律描述功能梯度材料沿厚度的梯度性质,基于von Kámán理论,建立了功能梯度薄板的非线性振动控制方程.应用Bubnov-Galerkin法得到了功能梯度矩形薄板的单模态非线性振动的时域常微分方程,借助其势能函数分析了系统的周期振动状态.采用Lindstedt-Poincaré法和Runge-Kutta法分别获得了功能梯度矩形薄板单模态非线性周期振动的摄动解和数值解.研究表明:功能梯度薄板非线性振动控制方程中包含表征拉弯耦合效应的控制项,这导致其常微分方程中出现二次项;系统振幅在板横向的正负两个方向上是不相等的,其振动存在关于板中面的不对称性.     

横观各向同性弹性半空间地基上四边自由各向异性矩形薄板弯曲解析解

选用更具广泛性的横观各向同性弹性半空间地基模型,来分析四边自由各向异性矩形地基板的弯曲解析解．将异性薄板的弯曲控制方程,与基于横观各向同性弹性半空间地基位移解建立的板与地基变形协调方程相结合,先按对称性分解,然后用三角级数法,得出横观各向同性弹性半空间地基上四边自由各向异性矩形薄板的弯曲解析解,包括地基反力、板的挠度及内力的解析表达式．该解析解克服了数值法的弊端,取消了对地基反力的假设,板的内力及地基反力求解更切实际．算例结果与文献结果吻合良好,证明本文方法的可行性．     

含冷却孔镍基合金次级取向效应的应变梯度晶体塑性有限元研究

单晶镍基合金具有优异的耐高温、高强、高韧等性能, 这些力学性能受制造过程引入的次级取向和冷却孔的影响. 已有研究大多关注单孔薄板的变形机理和力学性能, 而工程中应用的往往是多孔薄板, 当前亟需阐明多孔的塑性滑移带变形机理、次级取向效应以及冷却孔引起的应变梯度效应. 文章采用基于位错机制的非局部晶体塑性本构模型对含冷却孔镍基单晶薄板的单拉变形进行了数值模拟. 此模型基于塑性滑移梯度与几何必需位错的关系引入了位错流动项, 因此可有效刻画非均匀变形过程中的应变梯度效应. 为了全面揭示含孔镍基薄板的次级取向效应, 系统研究了[100]和[110]取向(两种次级取向)下镍基薄板的单拉变形行为, 并重点探究了在两种次级取向下冷却孔数量对薄板塑性行为的影响. 此外, 还分析了镍基合金板变形过程中各个滑移系上分切应力变化、主导滑移系开动以及几何必需位错密度的演化过程, 并讨论了塑性滑移量及其分布特征对不同次级取向镍基合金板强度的影响. 研究表明, 单孔和多孔的[110]薄板抗拉强度均低于[100]薄板, 多孔薄板的塑性变形过程比单孔薄板更为复杂且受次级取向影响更大, 并且发生滑移梯度位置主要位于冷却孔附近以及塑性滑移带区域. 研究结果可为工程中镍基合金的设计和服役提供理论指导.      

Nonproportional loading effects on elastic-plastic behavior based on stress resultants for thin plates of strain hardening materials

A stress resultant constitutive law in rate form is constructed for power-law hardening materials. The change of plate thickness is considered in the constitutive law. The elastic-plastic behavior of a plate element based on the stress resultant constitutive law under uniaxial combined tension and bending is determined under a limited number of nonproportional and unloading paths. The results based on the stress resultant constitutive law and the through-the-thickness integration method are compared within the context of both the small-strain and finite deformation approaches. The results indicate that the selection of the normalized equivalent stress resultant and the corresponding work-conjugate normalized equivalent generalized strain is appropriate for describing the hardening behavior in the stress resultant space. However, the hardening rule in a power law form must be modified for low hardening materials at large plastic deformation when finite deformation effects are considered.     

Pseudo plane stress plastic field for mode I crack growth

The pseudo plane stress field for a mode I crack growth is analyzed for both perfectly plastic and power law hardening plastic materials. When finite strain is taken into account, it is found that for perfectly plastic materials, the plastic domain is a narrow strip ahead of the crack tip. For power law hardening plastic material, the plastic domain contains a strip and a region ahead of the strip. The fracture criterion is discussed. The energy dissipated in the plastic strip is found to be proportional to the square of the thickness. Singular solutions to the field are ruled out by analysis.     

Numerical simulation of roll forming for channel section with outer edge

ＩｎｔｒｏｄｕｃｔｉｏｎＲｏｌｌｆｏｒｍｉｎｇｉｓｏｎｅｏｆｔｈｅｉｍｐｏｒｔａｎｔｍａｎｕｆａｃｔｕｒｉｎｇｍｅｔｈｏｄｓｏｆｍｅｔａｌｐｒｏｄｕｃｔｓ.Ｓｏｍｅｒｅｓｅａｒｃｈｅｒｓｓｔｕｄｉｅｄｒｏｌｌｆｏｒｍｉｎｇｐｒｏｃｅｓｓｂｙｄｉｆｆｅｒｅｎｔａｎａｌｙｓｉｓｍｅｔｈｏｄｓ .Ｔ .Ｒ .ＷａｌｋｅｒａｎｄＲ .Ｊ.Ｐｉｃｋ[1]ｃｏｎｄｕｃｔｅｄｔｈｅｇｅｏｍｅｔｒｉｃｍｏｄｅｌｉｎｇｏｆａｎＥＲＷｐｉｐｅｓｋｅｌｐ ,ａｎｄｐｒｅｓｅｎｔｅｄａｎｅｘｐｅｒｉｍｅｎｔａｌｍｅｔｈｏｄｔｏｃｈｅ…     

Boundary layers of hierarchical beam and plate models

The boundary behavior of a family of hierarchical models of linearly elastic, isotropic plates is studied. The hierarchical models are obtained by spectral semidiscretization of the displacement fields in the transverse direction and strain energy projection. The well known Reissner-Mindlin model is contained in the hierarchy as a special case. A decomposition of the boundary layers of any model in the hierarchy into a bending and a torsion layer, both of which are model dependent, is given. It is shown further that the bending and torsion layers arise as Galerkin approximations of certain (nonlinear) eigenvalue problems in the plate cross section with the subspaces used to derive the hierarchical model. It is shown that the bending layers converge, as the order of the plate model tends to infinity, to the so-called Papkovich functions on an elastic strip.The regularity of the solution on polygonal plates is investigated for the whole hierarchy of plate models and shown to equal the regularity of the plane elasticity problem.Partially supported by AFOSR Grant No. F49620-92-J0100.     

Crack tip field in a linear elastic–plastic strain-hardening material

The strip necking model for strain-hardening materials is studied in this paper, in which the stress distributed over the strip necking zone is assumed to be ultimate stress. The bi-linear stress–strain relation which can model certain features of plastic flow is adopted in this model. The stress and strain fields are calculated based on this model in this paper. The size of the strip necking region is determined by balancing the stress intensity factor due to remote loading with that due to assumed closing forces equal to the ultimate tensile strength of the material distributed over the strip necking zone. It is interesting that the strip necking region size and the crack tip opening displacement depend not only on the remote load, but also the material hardening parameters, which is different from the results of strip yield model. The results agree with experiments well, and the model has wider application.     

Material Characterization at High Strain by Adapted Tensile Tests

The strength of materials at high strain levels has been determined using the so-called Continuous-Bending-under-Tension (CBT) test. This is a modified tensile test where the specimen is subjected to repetitive bending at the same time. This test enables to create high levels of uniform strain. A wide variety of materials has been tested this way. The strength of the material after CBT testing has been measured in different ways: by secondary tensile tests, by interrupted CBT tests, and directly from the fracture in the CBT test. All methods yield similar results: the strength is largely unaffected by the cyclic pre-deformation and mainly depends on the overall increase in length. Only for multi-phase materials the strength shows a minor influence of CBT test conditions. The hardening follows the extrapolated hardening observed in a conventional tensile test, except for brass. This test method can potentially be used for measuring hardening curves at high strain levels.     

椭圆截面截卵形刚性弹体正贯穿加筋板能量耗散分析

为获得椭圆截面截卵形刚性弹体正贯穿加筋板的剩余速度,根据椭圆截面弹体贯穿靶板的破坏特征,认为贯穿过程中靶板的能量耗散方式主要为塞块剪切变形功与塞块动能、扩孔塑性变形功、花瓣动力功、花瓣弯曲变形功、靶板整体凹陷变形功、加强筋侧向凹陷变形功。推导了每种能量计算方法,计算中定量考虑了靶板扩孔、花瓣弯曲、凹陷变形的应变率效应。根据能量守恒关系,得到了椭圆截面弹体剩余速度和弹道极限速度预测公式。并通过实验结果对模型进行了验证。结果表明：考虑靶板应变硬化、应变率效应的贯穿模型可以准确预测弹体剩余速度;随着椭圆截面弹体长短轴之比的增大,靶板的弹道极限速度近似线性增大;长短轴之比小于3时,加筋板的主要耗能为花瓣弯曲变形能、整体凹陷变形能。

     

Determination of material intrinsic length and strain gradient hardening in microbending process

Microbending experiments of pure aluminum show that the springback angles increase with the decrease of foil thickness, which indicates obvious size effects and attributes to plastic strain gradient hardening. Then a constitutive model, taking into accounts both plastic strain and plastic strain gradient hardening, is proposed to analyze the microbending process of thin foil. The model is based on the relationship between shear yield stress and dislocation density, in which the material intrinsic length is related to material properties and average grain numbers along the characteristic scale direction of part. It is adopted in analytical model to calculate the non-dimensional bending moment and predict the springback angle after microbending. It is confirmed that the predictions by the proposed hardening model agree well with the experimental data, while those predicted by the classical plasticity model cannot capture such size effects. The contribution of plastic strain gradient increases with the decrease of foil thickness and is independent on the bending angle.     

Effects of specimen thickness, hardening and crack closure for the plastic strip model

The triaxial stress constraint T_z and the effective yield stress distribution in the plastic zone for strain hardening materials are considered for the strip where the thickness effect is investigated by introducing a plastic constraint factor α. This factor depends on the specimen thickness, crack length, load level and hardening exponent. A simple expression of the plastic zone length and an expression involving α are obtained. Application of the strip model to crack closure shows that the specimen thickness has an influence and the results are compared with those found finite element.     

A combined phenomenological model for the representation of anisotropic hardening behavior in high strength steel line pipes

Line pipes have anisotropic mechanical properties, such as tensile strength, ductility and toughness. These properties depend on both prestrain during the cold forming process and on the anisotropy of the mother plates. In this study, a phenomenological model combining isotropic and kinematic hardening is developed to represent anisotropic hardening behavior of high strength steel line pipes. The model is adjusted on experiments carried out on smooth and notched axisymmetric bars and plane strain specimens. The model is used to simulate bending tests carried out on large pipes containing a geometric imperfection. Numerical results suggest that prestraining in pipe forming process significantly affects the bending capacity of pipes.