空爆载荷作用下固支弹塑性圆板的动力学模型

吴迪 米国 郭香华 张庆明

彭克锋, 潘昊, 赵凯, 郑志军, 虞吉林. 铜粉末动态压缩行为的多颗粒有限元分析[J]. 高压物理学报, 2019, 33(4): 044102. doi: 10.11858/gywlxb.20180665
引用本文: 吴迪, 米国, 郭香华, 张庆明. 空爆载荷作用下固支弹塑性圆板的动力学模型[J]. 高压物理学报, 2022, 36(5): 054202. doi: 10.11858/gywlxb.20220525
PENG Kefeng, PAN Hao, ZHAO Kai, ZHENG Zhijun, YU Jilin. Dynamic Compaction Behaviors of Copper Powders Using Multi-Particle Finite Element Method[J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 044102. doi: 10.11858/gywlxb.20180665
Citation: WU Di, MI Guo, GUO Xianghua, ZHANG Qingming. Dynamic Model of Clamped Elastoplastic Circular Plate under Air Blast Loading[J]. Chinese Journal of High Pressure Physics, 2022, 36(5): 054202. doi: 10.11858/gywlxb.20220525

空爆载荷作用下固支弹塑性圆板的动力学模型

doi: 10.11858/gywlxb.20220525
详细信息
    作者简介:

    吴 迪(1997-),女,硕士研究生,主要从事结构冲击动力学研究. E-mail:15010059608@163.com

    通讯作者:

    郭香华(1974-),男,博士,副教授,主要从事爆炸与冲击仿真、材料与结构动力学研究.E-mail:guoxh@bit.edu.cn

  • 中图分类号: O383

Dynamic Model of Clamped Elastoplastic Circular Plate under Air Blast Loading

  • 摘要: 基于等效单自由度模型理论,考虑变形过程中弯矩和膜力的联合承载,建立了空中近爆作用下固支弹塑性圆板的动力学模型,实现了圆板加载和卸载的全过程描述。基于文献中的试验工况建立有限元模型,利用有限元数值模拟对空爆作用下固支圆板的动态响应进行分析,并将动力学模型计算结果与试验结果及有限元数值模拟结果进行对比,验证了动力学模型计算结果的准确性。结果表明,理论计算结果与试验结果及数值模拟结果吻合较好,建立的动力学模型可以应用于固支圆板在空中近爆载荷作用下的大变形预测,为结构抗爆提供技术支持。

     

  • 颗粒金属材料具有多尺度、孔隙度、流动性及一定强度等特性,已被广泛应用于许多重要工业领域,例如航空航天、自动化、电子[12]。颗粒金属经过冲击压缩可以制备出不同性能的材料,因此金属粉末压实越来越受到工业制造领域的关注[3]。相较于传统的压实方法,冲击压缩可以将粉末快速熔化、快速固化,压缩得更加密实[45]。国内外学者从20世纪80年代开始对颗粒金属材料的冲击压缩行为进行研究,但其冲击压缩过程的力学响应尚未得到很好的理解,而这恰恰是颗粒金属材料设计和应用的关键问题,尤为必要进行深入研究。

    颗粒金属材料的冲击压缩过程历时短、反应迅速,且粒子直径很小,准确捕捉其响应特征较为困难[68]。一些学者以单质金属粉末、金属粉末混合物作为研究对象,借助温度压力传感器等瞬态测试工具,通过高速飞片冲击试验对材料冲击前、后的初态和终态进行测定[910]。然而,实验测试难以获得压缩过程的响应信息。近来,大量数值模拟方法从细观角度出发研究颗粒材料的冲击压缩行为,如欧拉流体动力学方法、有限元法、离散元法,主要研究的不均匀颗粒材料包括铝[11]、铜[12]、多晶铁[13]、铍[14]等。连续介质有限元方法可以表征局部相对密度分布、应力分布和粉末流动,但由于粉末是不均匀的,不能从粒子尺度表征颗粒接触和相互作用。离散元方法能弥补上述缺点,可从颗粒水平观察力学响应,但由于简化了粒子接触算法,所以不适用于压实过程。多颗粒有限元方法[1517]集成以上两种方法的优点,适用于完全压实过程且能很好地阐明颗粒局部变形及力链传递。

    已有的一些理论模型可以用来描述颗粒金属材料的动态压缩行为。压力-孔隙比(P-α)模型[18]仅需要定义参数α,即孔隙材料与基体材料体积之比,就可以捕捉到压力随孔隙比的变化关系。P-α模型通过引入α的方式成功分离孔穴压溃和基体压缩的影响,但是α的演化关系需要额外补充。经验性的球壳压溃模型[6]随后被发展,该模型提出包含惯性项的空穴压溃方程,成功应用于静水压加载,压力为准静态和惯性相关项之和。基于球壳压溃模型,许多学者对其进行修正,分别引入偏应力[7]、温度[19]、黏性[20]、微惯性[21]的影响。然而,大部分模型都是基于经验提出的,直接得到动态的应力-应变关系在实验上是困难的,而细观有限元模型可以得到局部应力-应变信息。因此,本研究采用细观有限元模拟冲击压缩过程,研究颗粒金属材料的冲击压缩模型。

    本研究采用多颗粒有限元方法,从粒子层次建立颗粒金属材料的二维细观有限元模型,研究在恒速冲击载荷作用下铜粉末的力学行为,分析冲击波的传播特性,观察变形及速度分布特征,表征冲击速度与冲击波速度的关系及波后动态应力-应变状态。

    模型简化为二维等径圆形颗粒,颗粒直径为30 μm,初始颗粒堆积方式大致分为4类:随机堆积、简单立方堆积、密排六方堆积和蜂窝堆积[2224]。其中密排六方堆积及随机堆积的准静态压缩测试结果表明,压力-平均相对密度关系的变化趋势与实验结果吻合较好,随机堆积可以较好地描述颗粒间的摩擦作用机制。为此,采用随机堆积的建模方式,颗粒在空间随机分布,颗粒之间不重叠。

    采用Matlab软件编写颗粒堆积模型,建立细观分布模型的步骤为:(1)计算出在一定空间内设定孔隙率ϕ(孔隙部分所占比重)下等径颗粒的个数,如表1所示;(2)随机撒点的坐标作为颗粒核点的初始位置,通过调整核点位置确保任意两点之间的距离大于或等于2倍颗粒半径;(3)实时更新颗粒位置坐标,直到所有颗粒都不重叠为止。例如图1颗粒材料生成试件的孔隙率分别为0.50、0.35、0.25。

    图  1  颗粒金属材料的几何模型
    Figure  1.  Geometric models of granular materials
    表  1  不同孔隙率下细观模型所含颗粒数
    Table  1.  Particle numbers of mesoscale model at different porosities
    Porosity (ϕ) Number of particles (N)
    0.60 142
    0.50 177
    0.45 195
    0.40 213
    0.35 230
    0.30 248
    0.28 255
    0.25 265
    下载: 导出CSV 
    | 显示表格

    数值模拟采用ABAQUS/Explicit 6.11,基体材料选择实验数据丰富的铜颗粒。对每个颗粒单独划分网格后导入有限元软件,颗粒网格划分如图2(b)所示,边界处网格进行加密处理,保证内部颗粒相互接触作用的精度,内部采用粗网格能有效提高计算效率,以节约计算时间。经过网格收敛性分析,每个颗粒划分成100个六面体单元,包含226个节点,采用体单元可以增加计算精度,但需约束节点在厚度方向的位移,限制粒子在二维平面内的变形。

    图  2  铜粉末恒速压缩的细观有限元模型(a)以及颗粒初始网格划分(b)
    Figure  2.  A finite element model of copper powders under constant-velocity compression (a) and the initial mesh of a particle (b)

    铜粉末在高速冲击加载时表现出应变率硬化、应变硬化、温度软化等特性,在此采用典型的应变率、温度相关的Johnson-Cook模型表征[25],写作

    σ=(σ0+Bεn)(1+Cln˙ε˙ε0)[1(TTrTmTr)m]
    (1)

    式中:屈服应力σ0为90 MPa,应变硬化系数n为0.31,材料参数B为292 MPa,C为0.025,参考应变率˙ε0取为1 s–1,熔点温度Tm为1360 K,参考温度Tr为300 K,温度软化系数m为1.09。

    本研究采用的基体材料模型可以较好地表征真实铜粉末在冲击加载时的温度响应特征,粉末在压缩变形后,部分塑性功转化为热能,转化系数为0.9,从而带来温度升高的现象[26]

    颗粒试件位于长和宽均为0.5 mm的正方形容器内,上端、下端及支撑端刚性板均固定不动,冲击端刚性板以速度V向右恒速压缩,如图2(a)所示。颗粒之间的摩擦系数设为0.2,颗粒与刚性板无摩擦[15]。颗粒之间的所有接触面定义为通用接触。法向接触属性采用硬接触,切向选用罚函数接触。

    通过对初始孔隙率为0.28的样本进行冲击压缩模拟,观察发现在动态冲击下,颗粒金属材料呈现出高度局部化的变形特征。随着冲击速度的提高,变形局部化的特征更加明显,波阵面变得更加平整,波阵面后方的颗粒也随着冲击速度的增大而变形更加急剧,由均匀变形渐渐转变为不均匀变形。应力波沿加载方向以近似平面的形式逐层传播,波阵面宽度大约为2~4个颗粒直径。当冲击速度为95 m/s时,波阵面后方还未完全压实,仍存在一些孔隙,孔隙率的减小依赖于颗粒的平动、转动和塑性变形。当速度增加到一定值时,波后几乎完全压实,没有孔隙,此时波后应变为最大压实应变,如图3所示。

    图  3  不同时刻试件变形图
    Figure  3.  Deformation patterns at different times

    单轴冲击压缩下,应力波以近似一维应变平面波的形式从冲击端向支撑端传播,因此颗粒材料的冲击压缩问题可以简化为一维问题进行分析。引入一维速度场分布定量化地表征冲击波的传播特征。从有限元结果可以直接得到在Lagrange坐标系下每个节点的速度信息,将各截面上的节点在加载方向的速度分量进行平均,得到铜粉末在不同时刻t的一维速度场分布。速度场结果很好地捕捉到高速冲击下铜粉末呈冲击波模式的变形特征,表明速度急剧变化集中在一条狭窄的区域内。冲击速度越高,此区域越窄;而冲击速度较低时,波阵面越宽,如图4所示。

    图  4  不同冲击速度下试件内的一维速度(v)分布
    Figure  4.  One-dimensional velocity distribution in specimens under different impact velocities

    塑性冲击波波阵面的位置可由局部应变率最大值的位置确定。基于一维应力波理论,连续性方程为

    v/X=ε/t
    (2)

    可知速度梯度与应变率相等。一维速度场分布特征表明,高速冲击时,铜粉末的速度分布可以分为3个区域:波前的未变形区域、波后的压实区和状态发生突变的区域。依据速度梯度的变化趋势可以看出速度的变化程度,速度梯度最大处,速度变化程度最剧烈,将该位置作为此时刻冲击波波阵面的位置,如图5所示。

    图  5  一维速度分布及其速度梯度 (V=225 m/s, t=0.4 μs, ϕ=0.28)
    Figure  5.  One-dimensional velocity distribution and the corresponding velocity gradient (V=225 m/s, t=0.4 μs, ϕ=0.28)

    不同时刻的波阵面位置与时间近似呈线性关系,此直线的斜率即为冲击波波速。冲击速度越高,斜率越大,冲击波波速越大,如图6所示。

    图  6  不同冲击速度下冲击波位置随时间变化的关系
    Figure  6.  Variation of the shock front position with impact time at different impact velocities

    将多颗粒细观有限元模拟所得结果与文献中气枪冲击实验结果[27]及欧拉有限元数值模拟结果[28]进行比较,如图7所示。Benson[28]所采用的欧拉有限元方法适用于速度更高的情形,而在较低的冲击速度下,该方法对波阵面位置的捕捉不准确,导致冲击波波速的结果偏高。而采用多颗粒有限元方法所得结果与实验结果吻合较好,二维简化模型所得冲击波波速略低于实验结果。Borg等[29]通过对陶瓷颗粒进行二维简化,数值模拟分析其冲击Hugoniot关系,所得结果类似。原因是相同孔隙率下,二维模型的堆积配位数比三维少,较真实粉末有更多的自由空间,粒子连接链较短,在一定程度上妨碍冲击波的传播,从而导致冲击波波速较实验结果偏低。此外,在较低的冲击速度下,颗粒材料存在孔隙和分布不均匀性,波阵面不平整,选取波阵面位置和计算波速时会带来一定的误差。随着冲击速度的增大,波阵面越来越狭窄,位置较易确定,波速计算精度提高,误差逐渐减小。

    图  7  冲击波速度-粒子速度关系与已有文献结果对比
    Figure  7.  Comparison of the shock wave speed-particle velocity relations with the data in the literature

    由第2节可知,在较高速度冲击下,应力波在颗粒试件中近似以塑性冲击波的形式沿加载方向逐层传播,基于应力波理论[30]分析,冲击波波阵面上的质量守恒和动量守恒方程分别为

    vBvA=Vs(εBεA)
    (3)
    σBσA=ρ0Vs(vBvA)
    (4)

    式中:Vs为冲击波速度,{vA, εA, σA}、{vB, εB, σB}分别为波阵面前后的速度、应变、应力状态。

    波前颗粒材料由于处于较疏松的状态,颗粒间存在间隙,弹性波不能向前传播,因此应力、应变、速度都假设为零,即波前状态为vA=0、εA=0、σA=0,与冲击速度无关。

    在动态冲击载荷作用下,波后应变与冲击速度是有关的,会随着冲击速度的增大而逐渐增大,直到达到理论上的最大压实应变。当冲击速度增大到一定程度时,波后应变逐渐逼近理论上能达到的最大压实应变,写作

    εB=εm
    (5)

    式中:εm=ϕ,为最大压实应变。

    恒速加载时波后速度为vB=V,因此可求得冲击波波速

    Vs=V/ϕ
    (6)

    将(6)式代入(4)式可得波后应力

    σB=ρ0V2/ϕ
    (7)

    观察发现,在速度较高时波前状态为vA=0、εA=0、σA=0,波后状态为vB=VεB=ϕσB=–ρ0V2/ϕ,仅含有一个材料参数初始孔隙率ϕ。基于此特征,发展以动态锁定应变为唯一参数的冲击波模型,如图8所示。图8中箭头表示从波前状态沿Rayleigh线突跃到波后状态。该模型未考虑材料的率敏感性和温度的影响,仅以初始孔隙率ϕ表征较高冲击速度下铜粉末的动态应力-应变关系。

    图  8  动态应变锁定冲击波模型
    Figure  8.  Dynamic shock model with strain locked

    通过已有文献实验结果验证了多颗粒有限元方法数值模拟的有效性,继续通过此方法计算不同孔隙率的铜粉末在不同冲击速度下粒子速度与冲击波速的关系。孔隙率分别为0.25、0.28、0.30、0.35、0.40、0.45、0.50、0.60的Vs-v关系,如图9所示。而通过动态应变锁定的冲击波模型得出的Vs-v关系与细观有限元模拟结果吻合极好,由此证明动态应变锁定的冲击波模型可以很好地表征不同孔隙率铜粉末的Hugoniot关系。

    图  9  不同孔隙率的铜粉末粒子速度与冲击波速度的关系
    Figure  9.  Variation of the shock wave speed with the particle velocity at different porosity of copper powders

    不同冲击速度下,冲击载荷不同;冲击速度越大,冲击端应力越大。高速恒速加载情形下,冲击端应力近似等于波后应力,由细观有限元模拟提取出的冲击端应力如图10所示。颗粒材料的不均匀性以及孔隙的存在导致冲击力振荡较大,但基本维持在某一特定值上下。由动态应变锁定的冲击波模型计算得出的波后应力可以很好地捕捉到这一水平。

    图  10  不同冲击速度下冲击端应力与冲击波模型所得波后应力比较
    Figure  10.  Comparison of the stress at the impact end obtained from finite element models and the stress behind the shock wave front predicted by the shock model at different impact velocities

    对于不同孔隙率的颗粒材料,波后应力随着孔隙率的增大逐渐减小,在某一恒定速度加载情形下,对塑性冲击波传到支撑端前这一时间范围内的冲击端应力进行平均,将所得应力与通过理论模型计算得出的波后应力进行比较。结果表明,冲击端平均应力与冲击波模型所得波后应力相对误差较小(小于6%),如图11所示。相对误差由下式求得

    图  11  不同孔隙率下冲击端应力与冲击波模型所得波后应力比较
    Figure  11.  Comparison of the stress at the impact end and the stress behind the shock wave front obtained by the shock model under different porosities
    δ=|σavgσs|σavg×100%(8)
    (8)

    式中:σavg为冲击端平均应力,σs为理论模型计算所得波后应力。

    采用多颗粒有限元方法建立二维颗粒金属材料的冲击压缩模型。基于不同孔隙率下等径圆形颗粒的随机堆积模型,采用恒速压缩的加载方式,得到铜粉末的一维速度场分布,确定了塑性冲击波波阵面位置,进而得到了粒子速度与冲击波波速的关系。波后应变随着冲击速度增大,最终保持为最大压实应变不变,而波前应力、应变、速度均为零。基于颗粒材料的这一特征,发展了以动态锁定应变为唯一参数的冲击波模型,该模型可以很好地表征Vs-v关系并确定波后应力振荡的平均值。同时需要指出的是,该冲击波模型仅在速度较大、波后应变为最大压实应变时有较好的适用性,而在速度较小时会有较大的误差,在中低速下还需要发展更加精确的模型。颗粒金属材料的细观特性对其宏观力学性能有重要影响,本研究虽然从细观尺度建立模型,可以分析其宏观和微观力学响应,但是并未研究细观参数颗粒形状、尺寸、分布等的影响。基体材料选用的是Johnson-Cook模型,考虑了基体材料的应变率效应,但是在高速冲击加载时,应变率效应与惯性效应难以解耦,理论模型并未考虑应变率效应,即使这样仍可以较好地捕捉其主要特征。

  • 图  载荷等效示意图

    Figure  1.  Equivalent load diagram

    图 (4)和式(5)的计算结果对比

    Figure  2.  Comparison of calculation results between Eq. (4) and Eq. (5)

    图  板元受力分析

    Figure  3.  Force analysis of plate microelement

    图  弯矩和膜力联合承载的屈服条件

    Figure  4.  Yield condition of combined bending moment and membrane force

    图  固支圆板发生弹塑性大变形后的应力场分布

    Figure  5.  Stress field of clamped circular plate after large elastic-plastic deformation

    图  加载阶段的等效质量-弹簧-阻尼模型

    Figure  6.  Equivalent mass-spring-damper model in loading phase

    图  卸载阶段的等效质量-弹簧-阻尼模型

    Figure  7.  Equivalent mass-spring-damper model in unloading phase

    图  试验布置

    Figure  8.  Test setup

    图  有限元模型

    Figure  9.  Finite element model

    图  10  各工况下有限元数值模拟得到的中心点位移-时间曲线

    Figure  10.  Displacement-time curves of center point by finite element simulation under various calculation cases

    图  11  峰值压力点及其拟合函数

    Figure  11.  Peak pressure points and fitting functions

    图  12  各工况下动力学模型计算结果与试验及有限元数值模拟结果的对比

    Figure  12.  Comparison of dynamic model results with test and finite element simulation results under various cases

    图  13  各工况下动力学模型与试验得到的靶板轮廓线的对比

    Figure  13.  Comparison of circular plate deformations between dynamic model and test under various cases

    表  1  3种高强度钢的试验参数[12]

    Table  1.   Test parameters of three kinds of high strength steel[12]

    CaseMaterialThickness/mmYield strength/MPaMass of TNT/kgDetonation distance/mm
    1Weldox700E steel88006250
    2B1 steel810706250
    3B2 steel816006250
    4Weldox700E steel1280010250
    5B1 steel12107010250
    6B2 steel12160010250
    下载: 导出CSV

    表  2  各工况下试验与有限元模拟得到的中心点位移对比

    Table  2.   Comparison of center point displacements between test and finite element simulation under various calculation cases

    CaseMaximum displacement of center point Final displacement of center point
    Test/mmSim./mmError/%Test/mmSim./mmError/%
    1145148.842.65 129137.426.53
    2128134.675.21102111.889.69
    310399.52−3.388579.07−6.98
    4166152.44−8.17153141.95−7.22
    5148133.87−9.55118120.121.80
    6115106.74−7.188379.28−4.48
    下载: 导出CSV

    表  3  各工况下动力学模型计算结果与试验结果的对比

    Table  3.   Comparison between dynamic model results and test results under various cases

    CaseMaximum displacement Final displacement
    Test/mmSDOF model/mmError/%TestSDOF model/mmError/%
    1145141.48−2.43 129129.220.17
    2128126.03−1.54102110.067.90
    3103101.49−1.478586.211.42
    4166160.24−3.47153148.57−2.90
    5148141.27−4.55118127.668.19
    6115114.19−0.708395.0314.49
    下载: 导出CSV
  • [1] 王芳, 冯顺山, 俞为民. 爆炸冲击波作用下靶板的塑性大变形响应研究 [J]. 中国安全科学学报, 2003, 13(3): 58–61. doi: 10.16265/j.cnki.issn1003-3033.2003.03.016

    WANG F, FENG S S, YU W M. Study on large plastic deformation response of target plate under explosive blast wave [J]. China Safety Science Journal, 2003, 13(3): 58–61. doi: 10.16265/j.cnki.issn1003-3033.2003.03.016
    [2] 吴成, 金俨, 李华新. 固支方板对水中爆炸作用的动态响应研究 [J]. 高压物理学报, 2003, 17(4): 275–282. doi: 10.3969/j.issn.1000-5773.2003.04.006

    WU C, JIN Y, LI H X. A study on square plate dynamic response under underwater explosion [J]. Chinese Journal of High Pressure Physics, 2003, 17(4): 275–282. doi: 10.3969/j.issn.1000-5773.2003.04.006
    [3] 何建, 肖玉凤, 陈振勇, 等. 空爆载荷作用下固支矩形钢板的塑性极限变形 [J]. 哈尔滨工业大学学报, 2007, 39(2): 310–313, 329. doi: 10.3321/j.issn:0367-6234.2007.02.036

    HE J, XIAO Y F, CHEN Z Y, et al. Plastic limited deformation analysis of the clamped rectangular steel plate subjected to air non-contact explosions [J]. Journal of Harbin Institute of Technology, 2007, 39(2): 310–313, 329. doi: 10.3321/j.issn:0367-6234.2007.02.036
    [4] DUFFEY T A. Large deflection dynamic response of clamped circular plates subjected to explosive loading [R]. Albuquerque: Sandia National Laboratories, 1967.
    [5] NURICK G N, PEARCE H T, MARTIN J B. Predictions of transverse deflections and in-plane strains in impulsively loaded thin plates [J]. International Journal of Mechanical Sciences, 1987, 29(6): 435–442. doi: 10.1016/0020-7403(87)90004-X
    [6] JONES N. Impulsive loading of a simply supported circular rigid plastic plate [J]. Journal of Applied Mechanics, 1968, 35(1): 59–65. doi: 10.1115/1.3601174
    [7] FLORENCE A L. Circular plate under a uniformly distributed impulse [J]. International Journal of Solids and Structures, 1966, 2(1): 37–47. doi: 10.1016/0020-7683(66)90005-9
    [8] CLOETE T J, NURICK G N. On the influence of radial displacements and bending strains on the large deflections of impulsively loaded circular plates [J]. International Journal of Mechanical Sciences, 2014, 82: 140–148. doi: 10.1016/j.ijmecsci.2014.02.026
    [9] CAMPBELL T I, CHARLTON T M. Finite deformation of a fully fixed beam comprised of a non-linear material [J]. International Journal of Mechanical Sciences, 1973, 15(5): 415–422. doi: 10.1016/0020-7403(73)90040-4
    [10] KARAGIOZOVA D, YU T X, SHI S Y, et al. On the influence of elasticity on the large deflections response of circular plates to uniform quasi-static pressure [J]. International Journal of Mechanical Sciences, 2017, 131/132: 894–907. doi: 10.1016/j.ijmecsci.2017.07.032
    [11] BIGGS J M. Introduction to structural dynamics [M]. New York: McGraw-Hill, 1964.
    [12] 钟贤哲. 高强度钢Weldox700E的动态本构关系及其抗爆特性 [D]. 北京: 北京理工大学, 2019.

    ZHONG X Z. Dynamic constitutive relations and explosion resistance properties of high strength steel Weldox700E [D]. Beijing: Beijing Institute of Technology, 2019.
    [13] RIGBY S E, TYAS A, BENNETT T. Single-degree-of-freedom response of finite targets subjected to blast loading: the influence of clearing [J]. Engineering Structures, 2012, 45: 396–404. doi: 10.1016/j.engstruct.2012.06.034
    [14] TAYLOR G I. The distortion under pressure of diaphragm which is clamped along its edge and stressed beyond the elastic limit [J]. Underwater Explosion Research, 1950.
    [15] KYOHEI K, PIAN T H H. Large deformations of rigid-plastic circular plates [J]. International Journal of Solids and Structures, 1981, 17(11): 1043–1055. doi: 10.1016/0020-7683(81)90012-3
    [16] WANG X D, YU T X. Parkes revisited: effect of elastic deformation at the root of a cantilever beam [J]. International Journal of Impact Engineering, 1991, 11(2): 197–209. doi: 10.1016/0734-743X(91)90006-2
    [17] 汤安民, 许明达, 赵蕾. 自由变形结构弹塑性分析的新方法 [J]. 广西大学学报(自然科学版), 2009, 34(1): 36–39. doi: 10.13624/j.cnki.issn.1001-7445.2009.01.008

    TANG A M, XU M D, ZHAO L. New method for elasto-plastic analysis of free-form deformation structure [J]. Journal of Guangxi University (Natural Science Edition), 2009, 34(1): 36–39. doi: 10.13624/j.cnki.issn.1001-7445.2009.01.008
    [18] 许明达. 自由变形结构弹塑性分析的新方法及应用 [D]. 西安: 西安理工大学, 2009.

    XU M D. The structure of free-form deformation elasto-plastic analysis and application of new method [D]. Xi’an: Xi’an University of Technology, 2009.
  • 加载中
图(13) / 表(3)
计量
  • 文章访问数:  259
  • HTML全文浏览量:  177
  • PDF下载量:  44
出版历程
  • 收稿日期:  2022-03-07
  • 修回日期:  2022-03-21
  • 刊出日期:  2022-10-11

目录

/

返回文章
返回