薄板剪应力起皱数值模拟研究.pdf

第 12 卷 4 期 应用基础与工程科学学报 Vol. 12 ,No. 4 　　 　 　　 　 2004 年 12 月 JOURNAL OF BASIC SCIENCE AND ENGINEERING December 2004 ( ) 文章编号 2004　　中图分类号 :TG386 　　文献标识码 :A 薄板剪应力起皱数值模拟研究 孙振忠 , 　杨玉英 , 　郑明复 ( 哈尔滨工业大学材料科学与工程学院 ,黑龙江 哈尔滨 150001) 摘要 :利用有限元法模拟了剪皱试样分叉后继变形的畸变与起皱过程 ,皱纹的位 置、形状、起皱方向与实验结果一致. 分析了剪应力失稳起皱机理 ,并给出了临界 失稳起皱的判据和剪应力起皱的评价指标. 最后还讨论了屈服强度 σ、硬化指 s 数 n 值、板厚方向性指数 r 值对抗剪皱性能的影响 ,为板材抗剪皱性研究提供 了一种有效的数值模拟方法. 关键词 :薄板 ;剪应力 ;起皱 ;数值模拟 在大型薄板冲压成形中 ,起皱现象是经常出现的主要质量问题之一. 而且与破裂现象 相比 ,起皱现象更具多态性. 皱纹的形成和增长受应力状态、板材的力学性能、工具的几何 形状、接触条件等诸多因素的影响. 根据引起起皱的外力条件不同 ,可以分为压应力起皱、 不均匀拉应力起皱、剪应力起皱及弯曲力起皱等[1 ] . 对压应力作用下的毛坯起皱问题 ,其 典型代表是圆筒件拉深时法兰起皱 ,许多学者进行了大量的研究 ,给出了多种判断起皱的 公式和临界压边力公式[25 ] . 对于不均匀拉应力下的起皱 , 日本学者吉田清太对大量的汽 ( 车覆盖件冲压成形中的起皱现象进行了研究 ,总结提出的方板对角拉伸实验 即 YBT 实 ) 验 ,对不均匀拉应力下的起皱机理、起皱规律、影响起皱高度的因素、消皱措施等进行了 大量研究[67 ] . 最近又有很多学者利用有限元法对不均匀起皱机理、影响因素及相关数值 模拟进行了研究[811 ] . 而对剪应力下的起皱 ,尤其对板材的抗剪应力起皱评价方法及相关 的数值模拟研究较少[1213 ] . 1 　有限元模型基础理论 1. 1 　单元模型 ( ) x i i = 1, 2 , 3 为空间固定直角坐标系, ei 为x i 的单位向量. 考虑三维的板壳单元, 由 ( ) 中面上各点沿其法向 厚度 伸长而张成的空间定义. 将板壳离散成二维单元, 不妨以四 ξ 节点单元为例. 为描述方便定义位于中面各离散点处局部坐标系 i , 板壳单元内任一点 p 的空间位置可表示为 4 ξ ( ξ ξ) 3 ( ) x i = N m 1 , 2 x + δV 1 ∑ i , m 2 m i , 3 m m = 1 式中 , ξ为厚向坐标, N