低碳钢材料扭转和平弯问题的有限元分析.pdfVIP

Home 低碳钢材料扭转和平弯问题的有限元分析 赵慧娟、庄茁、郑泉水 北京清华大学工程力学系100084 一、前言 在现有的国家标准 （如8653-88）中，通过拉伸实验测得常温下金属的杨氏 模量、泊松比等材料常数。在许多工程应用中，也采用拉伸实验来测定诸如低碳 钢这类常用材料的幂硬化指数。在大变形的情况下，低碳钢等金属材料在拉伸实 验中会出现颈缩等局部化塑性失稳现象，此时材料均匀变形的假设不再成立。传 统的拉伸实验已经无法测得材料真正的本构关系，特别是在大变形下的塑性本构 关系。针对低碳钢这种材料，通过圆柱拉伸和扭转实验表明，在圆柱大扭转实验 过程中，试件始终保持均匀变形直到断裂破坏为止，没有任何局部化现象发生。 另外，从实验数据中可以明显地看出，扭转实验提供的真实有效的塑性区的数据 范围大概是拉伸实验的 5 倍（图 1）。因此，我们提出，能否用实心圆柱的扭转 实验作为研究低碳钢这类弹塑性材料在大变形特征下更为有效的基本实验，而以 传统的拉伸实验作为补充。 拉伸实验名义应力－名义应变关系曲线 扭转实验母线伸长率与扭矩关系曲线 图1 低碳钢材料圆柱试件拉伸与扭转试验结果的比较。 - 1 - Home 在本文中，主要讲述如果针对材料本构关系编制 ABAQUS 材料用户子程序 UMAT 的方法以及在模拟低碳钢材料圆柱扭转和平弯问题过程中，有关ABAQUS 高 级有限元软件的应用。 二、低碳钢材料本构与UMAT 通过观察实验曲线发现，低碳钢在进入塑性时的应力—应变关系表现为典型 的幂硬化曲线的特征。因此，采用对于低碳钢材料通用的幂硬化本构模型对已有 的扭转实验数据进行拟合。对于幂硬化材料 / f ( ) σ σ ε e Y e  ≤ ε/ εY ε εY (1) f ( )ε  N ( )( ) sign ε ε / εY ε εY ′ ′ ′ ′ 这里σ 3σ σ 2 为等效应力，ε 2ε ε 3 为等效应变，σ 为等效屈服应 e ij ij e ij ij Y ′ ′ 力，ε 为等效屈服应变，N 为幂硬化指数，其中σ 和ε 分别是σ 和ε 的偏分量。 Y ij ij ij ij ABAQUS 高级有限元软件提供了一些用户接口，对于特殊材料或单元，用户 可以根据需要编写自定义程序对软件进行补充。利用ABAQUS 软件的这一功能， 针对幂硬化材料的本构关系，我们编写了用户程序 UMAT。利用每一增量步的应 力、应变和应变增量，通过用户给定的本构关系求出下一步的应力和相应的刚度 矩阵∂σ/ ∂ε （在UMAT 中为DDSDDE 矩阵）。在输入文件中，需要给出UMAT 所需用 σ ε