一种预测圆形金属薄板低速冲击凹坑尺寸的方法
【技术领域】
[0001] 本发明提供一种预测圆形金属薄板低速冲击凹坑尺寸的方法，属于力学设计领 域。
【背景技术】
[0002] 金属薄板经常容易受到沙石、维修工具等意外低速冲击事件，导致不可恢复的塑 性变形，形成冲击凹坑，直接影响其静强度和疲劳等性能，因此，需要对这类低速冲击碰撞 现象进行评估金属薄板受低速冲击下的冲击凹坑尺寸对于设计人员是非常关键的输入参 数之一，尤其对材料选取和尺寸设计至关重要，因此，如何获取金属薄板在低速冲击下的冲 击凹坑尺寸一直以来都是工程人员需要解决的课题。通过金属薄板低速冲击试验可以直接 测定其冲击凹坑尺寸，但是试验成本比较高，并且费时费力周期长，尤其在初始设计阶段， 把任意材料、尺寸和冲击速度情况下的冲击凹坑尺寸都通过试验确定是不切实际的；有限 元数值模拟方法需要建立复杂的有限元模型，计算复杂，计算效率低；现有的一些解析方法 比较复杂，理论求解也不方便。在评估金属薄板受低速冲击下的冲击凹坑尺寸时，常采用经 典的圆形金属薄板作为受冲击的靶板。因此，本发明提出了一种预测圆形金属薄板低速冲 击凹坑尺寸的解析方法，该方法非常简单实用，仅仅需要少量金属薄板材料参数和几何参 数，就可以很容易得到与不同冲击速度下的冲击凹坑尺寸，可见本发明具有重要学术意义 和工程应用价值。
【发明内容】
[0003] 本发明提供了一种预测圆形金属薄板低速冲击凹坑尺寸的新方法，该方法具有计 算简便，精度高等优点，其技术方案如下：
[0004] 步骤一、提出预测圆形金属薄板低速冲击凹坑尺寸的新方法的假设条件。
[0005] 假设条件包括：
[0006] (1)金属薄板的冲击凹坑形状为旋转对称的，并且不考虑回弹变形对凹坑尺寸的 影响；
[0007] (2)满足基尔霍夫-乐甫薄板假设，因此可忽略面外法向应力和横向剪切应力的 影响，主要为径向拉伸应力、径向弯曲应力和周向弯曲应力；
[0008] (3)金属薄板的材料为弹塑性线性强化材料，其应力-应变关系如图3所示，弹性 阶段，线和为塑性阶段，其对应的斜率分别为E和E% %为弹性极限应变；
[0009] (4)忽略空气阻力、冲击过程摩擦力等能量耗散，认为冲击能量全部转化为应变 能。
[0010] 步骤二、根据圆形金属薄板低速冲击时的受载形式和凹坑特征，确定冲击凹坑的 变形函数。
[0011] 根据前面基本假设（1)，由于冲击凹坑形状具有旋转对称性，因此，可在圆柱坐标 系r0z下来描述冲击凹坑形状，其中，r为径向距离坐标，0方位角坐标。冲击凹坑区域中 面的Z向变形用W表示，可知W仅为关于r的函数，而与0无关。圆形金属薄板边缘固支， 中心作用有集中载荷P，则距中心为r的同心圆处有总剪力需与集中载荷p平衡，具体可表 示为
[0012] 2JrrQr=p⑴
[0013] 式中，Qr为距板中心为r的同心圆处的剪力。
[0014] 根据剪力的极坐标公式：
[0015]
(2)
[0016] 将式⑵代入式（1)中，可得
[0019] 对式（4)积分三次，可得
[0020]
(5)
[0021] 式中，Q、(：2和C3为待定常数；D为弯曲刚度，具体可表示为
[0022]
供）
[0023] 式中，E#为屈服段模量，v为泊松比，t。为板的厚度。
[0024]圆形金属薄板固支约束时的边界条件需满足：
[0025]当r=0变形函数需满足如下边界条件：
[0026]
(7)
[0027] 当r=r。变形函数需满足如下边界条件：
[0028]
(S)
[0029] 将式（5)分别代入式（7)和式⑶中，可得
[0030]
(9)
[0031] 将式（9)代入式（5)中，可得
[0032]
(10)
[0033] 步骤三、利用步骤二提出的圆形金属薄板冲击凹坑变形函数，进一步建立对应的 应变关系。
[0034]根据前面基本假设（2)，圆形金属薄板在受低速冲击的变形过程中，主要靠径向拉 伸变形、径向弯曲变形和周向弯曲变形吸收冲击能量。
[0035] 圆形金属薄板径向拉伸应变￡1^可表示为
[0041] 式（13)的变量替换相当于对变量r进行了归一化处理，因此0彡x彡1。
[0042] 将式（13)代入式（12)中，可得
[0043]
(14)
[0044] 式中，H为常数，可表示为
[0045]
(1.5):
[0046] 薄板的径向弯曲曲率K,周向弯曲曲率K^可表示为
[0047]
〇6):
[0048] 将式（13)代入式（16)中，可得
[0049]
(17)
[0050] 薄板的径向弯曲应变\,和周向弯曲应变eeb可表示为
[0051]
<18)
[0052] 步骤四、根据功能原理，建立圆形金属薄板在低速冲击条件下控制方程，再利用数 值求解方法求解控制方程的待定参数，最终确定冲击凹坑尺寸。
[0053] 根据前面基本假设（3)，圆形金属薄板径向拉伸变形的应变能可表示为
[0057] 为使式（20)表述更为清楚简洁，进行变量替换
[0058] Ut=A1p4+A2p2+A3 (21)
[0059] 式中，A。AjPA3为中间变量，可表示为
[0063] 同理可得，圆形金属薄板弯曲变形的应变能
[0064]
[0065] 式中，A4、A#PA6为中间变量，可表示为
[0069] 圆形金属薄板总应变能可表示为
[0070] U=Ut+Ub=Ajp4+(A2+A4)p2+A5p+(A3+A6) (29)
[0071] 冲击物的冲击能量为
[0072] Q=mgh(30)
[0073] 式中，m为冲击物质量，g为重力加速度，h为冲击高度。习惯上通常用Q表征冲击 能量。
[0074] 总的冲击能量还需要考虑凹坑深度的影响，即
[0075] Q*=mg(h+ 8 ) (31)
[0076] 式中，S为凹坑深度。
[0077] 冲击凹坑深度可表示为
[0078]
(32)
[0079] 根据前面基本假设（4)，由于冲击能量全部转化为应变能，即
[0083] Ajp4+ (A2+A4)p2+ (A5-Bj)p+A3+A6-Q= 0 (37)
[0084] 其中，为中间变量，可表示为
[0085]
(38)
[0086] 通过数值方法很容易求出方程式（37)中的未知量p，再将p的解代入式（10)和式 (32)中，即可确定与冲击能量Q对应的冲击凹坑变形和冲击凹坑深度5。此外，当冲击物 为水平冲击时，即凹坑深度不会导致额外冲击能能量时，只需令式（37)中的队为0即可。
【附图说明】
[0087] 图1为固支约束的圆形金属薄板受低速冲击的示意图。
[0088] 图2为固支约束的圆形金属薄板受低速冲击时简化的受力形式。
[0089] 图3为刚塑性线性强化材料应力-应变本构关系。
[0090] 图4为是本发明所述方法的流程框图。
[0091] 图中符号说明如下：
[0092] 图1中的x为直角坐标系下的坐标，y为直角坐标系下的坐标，z为直角坐标系下 的距离坐标，〇为直角坐标系下的原点，r为圆柱坐标系下径向距离坐标，0为圆柱坐标系 下方位角坐标。
[0093] 图2中的p为圆形金属薄板中心作用的集中载荷，Q1^为距圆形金属薄板中心为r 的同心圆处的面外剪力。
[0094]图3中的E为线性强化材料模型中应力-应变曲线在弹性阶段的曲线斜率，f为 线性强化材料模型中应力-应变曲线在塑性阶段的曲线斜率，〇为应力， 〇s为线形强化材 料模型的屈服应力，e为应变，％为最大弹性应变。
【具体实施方式】
[0095] 图4为本发明所述方法的流程框图，本发明份4步实现，具体为：
[0096] 步骤一、提出预测圆形金属薄板低速冲击凹坑尺寸的新方法的假设条件。
[0097] 假设条件包括：
[0098] (1)金属薄板的冲击凹坑形状为旋转对称的，并且不考虑回弹变形对凹坑尺寸的 影响；
[0099] (2)满足基尔霍夫-乐甫薄板假设，因此可忽略面外法向应力和横向剪切应力的 影响，主要为径向拉伸应力、径向弯曲应力和周向弯曲应力；
[0100] (3)金属薄板的材料为弹塑性线性强化材料，其应力-应变关系如图3所示，其中 线