镁合金板料热拉深成形实验与数值模拟研究--《南昌大学》2008年博士论文
镁合金板料热拉深成形实验与数值模拟研究
【摘要】：
镁合金作为目前实际应用中最轻的金属结构材料,在3C行业、汽车行业、国防、航空航天产品、人们日常的生活应用前景广泛。镁合金板料的冲压产品具有较好的力学性能和表面质量成为镁合金材料应用的一个趋势,而在冲压工艺中又以拉深工艺应用较广。由于拉深工艺过程是一个同时包含几何、材料、边界条件和力学等非线性条件于一体的极其复杂的过程,对镁合金拉深成形工艺的研究成为能否获得复杂镁合金塑性成形构件关键。
本文针对目前镁合金板料拉深成形工艺研究现状与实际需求,对镁合金板料热拉深成形过程的关键技术进行了实验研究、理论研究、数值模拟及二次开发研究。研究结果对指导镁合金产品设计、缩短产品开发周期、减少材料消耗和保证产品质量等具有科学意义和重要的实际应用价值。
主要内容如下:
设计制造了一套能自动控制拉深温度、拉深速度的镁合金板料热拉深成形的实验装置,为实现镁合金板料热拉深成形件产业化生产装置提供了参考;并用其对运用前景广阔的AZ31B、ME20M镁合金板料进行了热拉深成形物理实验,得到了两种镁合金板料热拉深成形最优工艺参数,为提出优化的镁合金板料热拉深成形工艺方案提供了依据。
利用光学金相显微镜对不同工艺参数下的AZ31B、ME20M镁合金板料热拉深成形圆筒件的金相组织结构进行了观察分析。结果表明:合理控制镁合金板料生产过程及镁合金板料热拉深成形过程工艺参数可以很好地控制镁合金板料拉深成形件金相组织结构的变化,使其具有较小晶粒组织,从而使镁合金拉深成形件具有良好的力学性能。
通过对镁合金板料进行热拉伸力学物理实验,首次提出了一个适合于镁合金板料热成形的含非常应变软化因子的高温流变应力数学模型;利用BP神经网络预测了镁合金板料热成形高温流变应力,将其预测结果与提出的数学模型预测结果、热拉伸力学物理实验结果进行了比较,进一步验证了本文提出的镁合金板料热成形高温流变应力数学模型的科学性、精确性;并用FORTRAN语言编程,在ABAQUS中进行了用户材料子程序VUMAT二次开发,将提出的高温流变应力数学模型加入到了数值模拟软件ABAQUS中,实现了在ABAQUS中加入并调用新的材料模型,使建立的材料模型更接近实际,拓展了ABAQUS软件的应用范围,提高了模拟精度。
合作设计制造了一套镁合金板料热拉深成形过程中凸凹模圆角处摩擦系数在线测量的实验装置,在综合考虑镁合金板料拉深成形过程中成形温度、压边力、润滑剂等因素的基础上,进行了镁合金板料热拉深成形凸凹模圆角处摩擦物理实验,分别提出了与成形工艺参数(成形温度、压边力)有关的凸凹模圆角处新的动态变化摩擦系数模型;并用FORTRAN语言编程,在ABAQUS中进行了用户子程序VFRIC二次开发,将提出的动态变化摩擦系数模型、修正的库仑摩擦模型相结合,加入到了数值模拟软件ABAQUS中,实现了在ABAQUS中加入真实摩擦边界条件的目的。
最后在ABAQUS软件中,分别通过不同方式获得了镁合金板料热拉深成形材料模型、摩擦边界条件,对镁合金板料热拉深成形规律进行了数值模拟研究,并与物理实验结果相比较。研究结果表明:本论文作者提出的镁合金高温流变应力模型、热拉深成形过程凸凹模圆角处的动态摩擦系数模型具有较强的科学意义与工程实用价值;利用二次开发将二者加入到ABAQUS模拟软件中去的方法是完全可行的;本文的研究成果在镁合金板料拉深成形件、乃至其它塑性成形产品生产实际中起到重要的指导意义,大大降低生产成本。
【学位授予单位】：南昌大学【学位级别】：博士【学位授予年份】：2008【分类号】：TG386.32
欢迎：、、)
支持CAJ、PDF文件格式
【相似文献】
中国期刊全文数据库
张庭芳;黄菊花;邓敏;董懿琼;游步东;兰志文;;[J];锻压技术;2006年06期
张庭芳;黄菊花;杨雪春;严萍;董懿琼;游步东;;[J];兵器材料科学与工程;2007年03期
张庭芳;谢世坤;;[J];材料导报;2009年04期
姜久光;;[J];热加工工艺;2011年15期
;[J];;年期
;[J];;年期
;[J];;年期
;[J];;年期
;[J];;年期
;[J];;年期
中国重要会议论文全文数据库
魏丽;付忠亮;李昌坤;王心朋;;[A];2009海峡两岸机械科技论坛论文集[C];2009年
魏丽;付忠亮;李昌坤;王心朋;;[A];十三省区市机械工程学会第五届科技论坛论文集[C];2009年
魏丽;付忠亮;李昌坤;王心朋;;[A];自主创新与持续增长第十一届中国科协年会论文集（2）[C];2009年
中国博士学位论文全文数据库
张庭芳;[D];南昌大学;2008年
中国硕士学位论文全文数据库
李彩霞;[D];哈尔滨理工大学;2006年
&快捷付款方式
&订购知网充值卡
400-819-9993上传列表：
- 弹塑性损伤本构用于动力学显式分析，可以计算损伤场，残余应力应变, 17:36:47,下载5次
近期下载：
- [断裂与失效] 使用VUMAT如何设置状态变量表征材料失效 - abaqus中生死单元问题，还在学习，仅供参考，有待提高 - 适用于ABAQUS的自定义材料子程序，适用于显式分析中描述材料失效的基于J2的材料模型 - ABAQUS 显式计算时 用户自定义子材料程序 VUMAT源代码， - 金属切削仿真ansys/lsdynaKwenjian - [断裂与失效] 2D 单元失效内部单元删除的小例子
- 基于VUMAT的粘弹性本构模型应力更新编程和单元失效准则以及失效单元删除 - 用于abaqus软件的单元删除，是应用在abaqus软件后处理时，删除已经失效的单元 - abaqus Johnson-cook UMAT VUMAT subroutine - ABAQUS有限元金属切削仿真模型(源自ABAQUS经典例题集三) - 采用子程序实现单元删除。子程序采用定义场变量、状态的方式，通过判断单元应力是否达到规定的限值，从而实现单元的删除。有需要的可以下载quenchplate_xpl_vumat
说明：&&最简单的一个vumat子程序，不同温度对应不用屈服，包含应力更新(The simplest VUMAT subroutine does not yield to different temperatures, including stress updates.)
文件列表：
quenchplate_xpl_vumat.f, 6075 ,
近期下载者：
相关文件：vumat_abqtest___一个模型中两种用户自定义材料
说明：&&一个模型中两种用户自定义材料，定义在平面应力状态下的材料属性(Two user custom materials in a model)
文件列表：
vumat_abqtest___一个模型中两种用户自定义材料.f, 7379 ,
近期下载者：
相关文件：您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
5.本构模型-应力更新专题-UMAT和VUMAT资料.pptx 43页
本文档一共被下载：
次 ,您可全文免费在线阅读后下载本文档。
下载提示
1.本站不保证该用户上传的文档完整性，不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值，立即自动返金币，充值渠道很便利
你可能关注的文档：
·········
··········
计算固体力学应力更新专题
庄茁
liuzhanli@tsinghua.edu.cn
;导出有限元方程（完全，更新、共旋）;应力更新;;;;;2.几种客观率;2.几种客观率;例4.1;
对于不同客观率采用了相同的材料常数，其差别是非常大的。事实上，这是误用了材料模型。材料模型必须根据不同的率转换。这是变形体，若是刚体转动，Jaumann率与Truesdell率是一致的。;Green-Naghdi率;切线模量之间的关系 ;2.几种客观率的关系;2.几种客观率的关系;2.几种客观率的关系;Comparison of different o2.几种客观率的关系;S后拉、前推;后拉、前推;Lie导数 ;后拉、前推和Lie导数 ;4.常用客观率应力更新流程;速度梯度;4.常用客观率应力更新流程;4.常用客观率应力更新流程;刚体转动求解;4.常用客观率应力更新流程;4.常用客观率应力更新流程;4.常用客观率应力更新流程;4.常用客观率应力更新流程;1
ABAQUS的UMAT;ABAQUS用户子程序;
由于主程序与UMAT之间存在数据传递，甚至共用一些变量，因此必须遵守有关UMAT的书写格式，UMAT中常用的变量在文件开头予以定义，通常格式为：
?
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,DRPLDE,
1 DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
2 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
3 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME
DIMENSION STRESS(NTENS),STATEV(NSTATV),
1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),
2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),
3 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)
?
user coding to define DDSDDE, STRESS, STATEV, SSE, SPD, SCD
and, if necessary, RPL, DDSDDT, DRPLDE, DRPLDT, PNEWDT
?
UMAT中的应力矩阵、应变矩阵以及矩阵 DDSDDE，DDSDDT，DRPLDE 等，都是直接分量存储在前，剪切分量存储在后。直接分量有 NDI个，剪切分量有 NSHR 个。各分量之间的顺序根据单元自由度的不同有一些差异，所以编写UMAT时要考虑到所使用单元的类别。;用于存储状态变量的矩阵，在增量步开始时将数值传递到UMAT中。也可在子程序USDFLD或UEXPAN中先更新数据，然后在增量步开始时将??新后的数据传递到UMAT中。在增量步结束时必须更新状态变量矩阵中的数据。
和应力张量矩阵不同的是：对于有限应变问题，除了材料本构行为引起的数据更新以外，状态变量矩阵NSTATEV中的任何矢量或者张量都必须通过旋转来考虑材料的刚体运动。;2
ABAQUS的VUMAT;2
ABAQUS的VUMAT;UMat:
For geometrically nonlinear analysis the strain increment and incremental rotation passed into the routine are based on the Hughes-Winget formulae.
– Linearized strain and rotation increments are calculated in the mid-increment con?guration.
– Approximations are made, particularly if rotation increments are large: more accurate measures can be obtained
正在加载中，请稍后...}