ANSYSLS-DYNA数值模拟霍普金森压杆试验资料
- 格式:docx
- 大小:173.51 KB
- 文档页数:11
简述“霍普金斯”杆测量材料动态应力应变曲线的原理;选用一种大型软件对其进行计算模拟,并对模拟结果进行分析。
答:选用ABAQUS大型有限元软件一、“霍普金斯”压杆理论:Hopkinson压杆技术源于1914年B.Hopkinson测试压力脉冲的试验工作,后来R.M.Davies对它进行了改进。
1949年,H.Kolsky在这些基础上建立了进行材料单轴动态压缩性能试验的试验方法,测试了高应变率下金属材料的力学性能,这个方法称为分离式Hopkinson压杆(或Kolsky杆)技术。
其原理是将试样夹持于两个细长弹性杆(入射杆与透射杆)之间,由圆柱形子弹以一定的速度撞击入射弹性杆的另一端,产生压应力脉冲并沿着入射弹性杆向试样方向传播。
当应力波传到入射杆与试样的界面时,一部分反射回入射杆,另一部分对试样加载并传向透射杆,通过贴在入射杆与透射杆上的应变片可记录入射脉冲,反射脉冲及透射脉冲,由一维应力波理论可以确定试样上的应力、应变率、应变随时间的变化,以及应力、应变曲线。
5O多年来,此技术广泛用在高变形速率下材料力学性能的测试。
研究人员也对Hopkinson压杆试验方法进行了系统深入的研究,使该技术不断地改善和发展。
J.Harding 等在1960年将用于单轴压缩试验的Hopkinson压杆推广到了单轴拉伸试验,在此基础上,1983年又提出至今被广泛使用的Hopkinson拉杆试验方法。
W.E.Backer等、J.D.Campbell 等、J.Dully等又提出了Hopkinson扭杆技术,可对于试样施加高应变速率的纯扭转载荷。
为提高试验精度,前人在应力波的弥散效应、三维效应、应力波分离、试样中的瞬态平衡对试验结果的影响等方面做了大量工作。
分离式Hopkinson压杆实验的示意图如下:图8-1 分离式Hopkinson压杆示意图上图表示了压杆、试件和测试仪器等的位置安排。
压杆由高强度合金钢制成。
压杆与试件的接触面需要加工得很平并且保持平行。
霍普金森压杆实验中的脉冲整形技术——笔记1、霍普金森压杆实验中的一些问题的现状:关于霍普金森压杆技术有效性的讨论过去主要集中试件的尺寸效应,波在杆中的二维弥散修正等。
实验过程中试件是否处于应力均匀状态以及试件是否以恒应变率变形这两个问题上所给予的关注并不多,或者说还没有找到一个非常可行的方法来解决这两个问题。
2、常规霍普金森压杆技术所遇到的问题:要得到有效并精确的数据,下列霍普金森压杆的假设必须得到满足:1)压杆中的波传播必须是平面一维的,因为应变片所测得的杆的表面应变通常代表压杆整个横截面上的轴向应变。
2)试件中的应力和应变均处于均匀状态3)此外,为保证得到有效的应力—应变数据,还应该使试件中的应变随时间变化的历史也是均匀的,即试件的变形是在恒应变率的条件下进行的。
所对应的问题:1)二维波动效应(或称为波的弥散效应)2)在高应变率霍普金森压杆实验中,加载的上升时间在10μs左右,高速撞击将导致明显的应力波传播效应(纵向惯性效应),低应变率下的试件中应力均匀性的结论不再成立,因而这时的试件也不可能处于实际的应力均匀状态3)在常规霍普金森压杆实验中,子弹的撞击在入射杆中产生一个梯形的入射脉冲。
由于试件横截面的增加和试件材料的硬化,应变率则会随时间减小以致于不能在整个实验中保持为一恒定值3、常规霍普金森压杆对不同材料测试时存在的主要问题:3.1金属类材料:因为金属的弹性行为发生在非常小的变形下,在这样的小变形下,要得到精确的实验数据,因弥散效应引起的波的振荡问题和试件中应力均匀性是必须要考虑的敏感问题。
3.2脆性材料:首先,在霍普金森压杆实验中必须保持脆性材料试件两个端面严格平行以增加实验数据的精度,因为试件端面的不平行或不平整都可能导致局部失效和应变测量的不精确。
其次,常规霍普金森压杆实验中陡峭的梯形脉冲也导致脆性试件在小变形下的严重应力不均匀。
此外,经典的梯形入射脉冲还会导致脆性材料试件非恒应变率变形。
霍普金森压杆实验流程下载温馨提示:该文档是我店铺精心编制而成,希望大家下载以后,能够帮助大家解决实际的问题。
文档下载后可定制随意修改,请根据实际需要进行相应的调整和使用,谢谢!并且,本店铺为大家提供各种各样类型的实用资料,如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等,如想了解不同资料格式和写法,敬请关注!Download tips: This document is carefully compiled by theeditor. I hope that after you download them,they can help yousolve practical problems. The document can be customized andmodified after downloading,please adjust and use it according toactual needs, thank you!In addition, our shop provides you with various types ofpractical materials,such as educational essays, diaryappreciation,sentence excerpts,ancient poems,classicarticles,topic composition,work summary,word parsing,copy excerpts,other materials and so on,want to know different data formats andwriting methods,please pay attention!霍普金森压杆实验是一种用来研究材料的压缩性质的实验方法,下面是霍普金森压杆实验的流程:1. 准备实验所需材料和设备:实验所需的材料包括被测样品和压杆,设备包括载荷传感器、位移传感器、数据采集系统等。
基于霍普金森压杆系统的动态压痕实验张新;候兵;李玉龙【摘要】设计了基于分离式Hopkinson压杆系统(SHPB)的动态压痕实验装置,使其能够获得动态压痕实验中试样的压痕位移和所受冲击载荷的实验数据,从而得到一种新的材料动态性能测试技术.对3种不同的材料进行了动态压痕实验,得到了材料的动态硬度和率敏感性等动态性能,与采用其他实验技术所测得的性能数据具有很好的一致性,并实现了利用动态压痕的方法获取材料的动态性能.%A new technology for testing dynamic properties of materials was proposed by developing an experimental setup based on a split-Hopkinson pressure bar system. The experimental setup can be used in dynamic indentation experiments to obtain the indentation displacements of materials and the impact load suffered by them. To validate the experimental setup, the ABAQUS/Explicit FE code was used to numerically simulate dynamic indentation experiments. Dynamic indentation experiments were conducted on aluminum alloy LY12, pure titanium and oxygen-free copper. The dynamic properties such as dynamic hardness and rate sensitivity were obtained for these three kinds of materials. The dynamic properties obtained by dynamic indentation experiments are consistent with those by other experimental methods.【期刊名称】《爆炸与冲击》【年(卷),期】2011(031)003【总页数】7页(P256-262)【关键词】固体力学;动态硬度;压痕;SHPB;率敏感性【作者】张新;候兵;李玉龙【作者单位】西北工业大学航空学院,陕西西安710072;西北工业大学航空学院,陕西西安710072;西北工业大学航空学院,陕西西安710072【正文语种】中文【中图分类】O347压痕硬度实验是一种应用广泛的材料性能测试技术,主要思想是强制压迫“较硬”的压头进入“较软”材料的表面,通过记录压入过程中载荷和压痕深度的关系曲线,研究材料抵抗外物侵入的能力,即材料硬度[1]。
A N S Y S L S-D Y N A数值模拟霍普金森压杆试验ANSYS/LS-DYNA数值模拟霍普金森压杆试验1 功能概述大多数材料在强度等力学性质方面都表现出某种程度的加载率或应变率敏感性,高幅值短持续时间脉冲和荷载所引起材料力学性质的应变率效应,对于抗动载的结构设计和分析是非常重要的。
这些动载来至常规武器侵彻与爆炸、偶然爆炸和高速撞击等许多军事和民用事件,对于这些事件的理论分析和数值模拟必须知道材料的高应变率强度、断裂特性和应力-应变关系等本构性质。
要研究材料在脉冲动载作用下的力学性质的实验设备和实验必须模拟类似现场的应变率条件,分离式霍普金森杆被公认为是最常用最有效的研究脉冲动载作用下材料力学性质的实验设备。
数值模拟是一种依靠电子计算机对工程问题和物理问题乃至自然界各类问题进行研究的技术。
它利用材料的本构函数,结合有限元或有限容积的概念,采用数值计算和图像显示的方法,因此具有如下优势:(1)检验理论结果是否正确;(2)弥补实验与观测得不足;(3)利用模拟结果,了解非线性过程中的因果关系与主要物理机制;(4)预测在不同初始条件与边界条件下非线性过程的发展情形;(5)数值模拟成本低,可以带来巨大社会经济效益。
由于很多材料的本构性质已经知道,因此在设计产品时,可以利用材料的本构性质通过仿真来模拟复杂的系统。
ANSYS/LS-DYNA数值模拟霍普金森压杆试验,就是通过ANSYS/LS-DYNA软件来模拟霍普金森压杆实验,通过设置弹丸不同速度,对试件进行研究。
霍普金森压杆实验分为自由式和分离式两种,本仿真采用分离式的办法。
2 原理简介2.1 霍普金森压杆实验简介霍普金森杆实验装置的基本原型最早是由Hopkinson提出的,它可用于测量冲击载荷的脉冲波形。
1949年Kolsky将压杆分成两段,试件置于输入杆和输出杆中间,通过加速的质量块、短杆撞击或炸药爆轰产生加速脉冲,利用这一装置可测量材料在冲击载荷作用下的应力-应变关系。
Kolsky的工作是一项革命性改进,现代的分离式霍普金森杆都是在其基础上发展而来,所以分离式霍普金森杆也称之为Kolsky杆。
在这半个多世纪的时间里,分离式霍普金森杆实验技术得到了大力的发展,由最初的压缩实验系统发展到拉伸和扭转实验系统甚至是压扭、拉扭复合系统;其测试材料的种类己由金属发展到非金属,由韧性材料到脆性材料;从常温实验发展到高、低温实验;从较低应变率实验发展到较高应变率实验。
另外,由于自由式霍普金森杆技术由于能够实验高过载,已经成为高g值加速度传感器的标准标定实验技术。
2.2 ANSYS/LS-DYNA简介在求解冲击、爆炸问题时,只有采用非线性数值算法才能解决。
广泛使用的有限元程序如LS-DYNA、DYTRAN、ABAQUS和AUTODYN等专长于求解非线性问题。
DYTRAN可用于爆炸、高速侵彻、船体撞击毁损等分析领域,但是在处理冲击问题的接触算法上远不如LS-DYNA全面。
ABAQUS可以分析复杂的固体力学和结构力学系统,特别是能够驾驭非常庞大复杂的问题和模拟高度非线性问题,但是对爆炸和冲击过程的模拟相对不足。
AUTODYN可用于处理几何和材料大变形的非线性瞬态动力分析数值模拟,尤其在弹药工程领域应用广泛,可对聚能射流现象、破甲穿甲弹侵彻靶板的作用机理进行仿真。
LS-DYNA作为世界上最著名的通用显示动力分析程序,能够模拟真实世界的各种复杂问题,特别适合于求解各种二维、三维非线性结构的高速碰撞、爆炸和金属成型等非线性动力冲击问题,同时可以求解传热、流体及流固藕合问题。
1996LSTC公司与ANSYS公司合作推出ANSYS/LS-DYNA,大大增强了LS-DYNA的分析能力。
ANSYS/LS-DYNA 程序是功能齐全的几何非线性(大位移、大转动和大应变)、材料非线性(140多种材料动态模型)和接触非线性(50多种)程序。
它以Lagrange算法为主,兼有ALE和Euler算法;以显式求解为主,兼有隐式求解功能;以结构分析为主,兼有热分析、流体-结构耦合功能;以非线性动力分析为主,兼有静力分析功能(如动力分析前的预应力计算和薄板冲压成型后的回弹计算);军用和民用相结合的通用结构分析非线性有限元程序。
运用ANSYS/LS-DYNA进行数值模拟的一般步骤包括:前处理、求解及后处理。
前处理工作主要包括:算法及有限元单元的选择、材料模型的确定、模型的建立、网格的划分、接触的定义、约束及载荷的施加等步骤,其中前处理工作由ANSYS/PREP7完成。
前处理操作完成后,系统将生成相应的K文件。
此时可根据数值模拟的实际情况来决定是否修改K文件内容。
由于K文件为二进制文件,需采用文本编辑器将K文件打一开、编辑。
K文件修改完毕后,调用LS-DYNA970版求解器进行求解,生成图形文件和时间历程文件。
求解过程结束后,启动后处理程序LS-PREPOST。
LS-PREPOST可读取LS-DYNA的计算结果d3Plot文件,进行计算数据的整理以及二次运算。
3 仿真材料及模型尺寸的选择分离式霍普金森压杆实验包括子弹、输入杆、输出杆及试件。
本仿真采用钨钢作为子弹、输入杆以及输出杆的材料,试件选用5A06-C铝。
数值模拟的模型尺寸如表3.1,数值模拟材料参数如表3.2,表3.3。
表3.1 数值模拟的模型尺寸模型长度l/m 外径∅/m 钨钢子弹0.2 1.5×10−2钨钢输入杆 1 1.5×10−2钨钢输出杆0.8 1.5×10−25A06-C铝试件8×10−3 1.3×10−2表3.2 子弹、输入杆、输出杆的材料参数和模型材料密度ρ/kg/m3杨氏模量E/Pa泊松比μ切线模量E T/Pa失效模量E S/Pa模型钨钢子弹7678 2.1E11 0.29 - - Isotropic 钨钢输入杆7678 2.1E11 0.29 - - Isotropic 钨钢输出杆7678 2.1E11 0.29 - - Isotropic表3.3 试件的材料参数和模型材料密度ρ/kg/m3杨氏模量E/Pa泊松比μA/Mpa B/Mpa n c m5A06-C铝试件2680 7.9E10 0.33 235.4 622.3 0.58 0.174 1.05 接上表:材料Ef P1SRSpecificHeatFailStressRoomTemp(K)MeltTemp(K)模型5A06-C 铝试件1e-6 4.77e-6 -9 294 1050Johnson-Cook上述3个表中采用的是kg-m-s单位制,在建模过程中采用的是cm-g-us单位制,因此需进行单位换算。
4 霍普金森压杆实验LS-DYNA有限元模型建立4.1 创建单元及材料类型首先选择单元类型,由于霍普金森压杆实验中子弹、输入杆、输出杆及试件均为圆柱体且共轴,因此,有限元模型单元类型选用三维实体SOLID164单元,并采用拉格朗日算法。
由于该实验过程时间短,为了方便输入参数以及仿真,本仿真采用cm-g-us单位制进行建模。
因此,在建模的过程中,对单位进行了换算,将通用的kg-m-s单位制参数换算成cm-g-us单位制。
另外,为了防止模型尺寸对数值模拟结果的影响,数值模拟过程中实体建模采用全比例尺寸。
为了节省数值计算工作量,采用给弹丸赋值初速度的办法忽略了炮膛和支架,并通过约束输出杆末端的办法忽略吸收杆,同时圆柱体的对称结构,建立了1/4三维实体模型,在后处理的时候还原成圆柱体。
4.2 划分网格在划分网格前,需要给实体模型赋予网格属性,即给实体模型选择之前创建的材料类型和单元类型。
在数值模拟过程中,网格密度太小会产生虚波现象,网格密度太大,对改善波形起的作用不是很大,相反会耗费大量机时。
因此,应该给实体模型选择合适的网格密度。
本仿真选择映射式网格划分有限元模型,并对前人的结果进行分析,不同网格密度的计算结果进行比较,设置了合适的网格数量。
4.3 定义接触LS-DYNA程序中处理接触-碰撞界面主要采用三种不同的算法,即:节点约束法、对称罚函数法和分配参数法。
在此,我采用了最常用的对称罚函数法。
对称罚函数法中接触刚度值K是个重要参数,根据应力波相关知识,弹性波在两个截面相等、波阻抗相同以及互相接触的弹性杆中传播时,应力波在杆的接触处无反射。
这里选择接触刚度值为1.2,此时入射杆上的反射波很微弱。
满足要求。
在数值模拟中,子弹与输入杆之间的接触类型为表面-表面自动接触类型,子弹为接触面,输入杆为目标面。
输入杆与试件之间的接触类型为节点-表面自动接触类型,试件与输出杆之间的接触类型也为节点-表面自动接触类型,由于试件的网格密度较大,所以试件为接触面,输入杆和输出杆为目标面。
4.4 定义对称和子弹速度由于采用了1/4三维实体模型,所以需要对对称面才用对称约束。
另外,由于省略了吸收杆,因此对输出杆末端进行了约束,建立了其z轴无位移约束,x轴和y轴方向的平动约束和转动约束。
另外,前处理中定义了子弹的速度,并可以通过修改K文件中的*INITIAL-VELOCITY-GENRATION关键字来修改子弹的速度,从而进行不同速度的比较。
4.5 定义求解条件为了求得实验的完整结果,所定义的求解时间应略大于应力波在压杆中完整传播一个来回的时间。
通过查询资料,求解时间设为600us,最终结果达到了预想的要求。
5 实验结果分析图5.1 分离式霍普金森压杆LS-DYNA模型图如图5.1所示,改图为分离式霍普金森压杆LS-DYNA模型,从左到右依次是输出杆、试件、输入杆以及子弹。
通过LS-DYNA建立模型后,输入相关求解条件,然后生成K文件,用LS-DYNA求解器求解,并获得名为3dplot的文件,然后通过LS-DYNA后处理软件LS-Prepost软件分析求解结果。
图5.2 子弹速度为10m/s时试件靠近输入杆一侧的应变弹性波速是材料的基本参数,它可由材料的弹性模量和密度计算得出。
如图5.2所示,当t=184.94us时,试件开始有应变,说明弹性波已经通过输入杆进入试件了,可以计算出数值模拟中弹性波在压杆中的传播速度C0C0=(1/184.94)×106=5407.16m/s图5.3 不同子弹速度下输入杆中点处应变图5.4 不同子弹速度下输入杆中点处应力图5.5 不同子弹速度下输入杆中点出加速度在分离式霍普金森杆数值模拟中,不同子弹速度下,输入杆上应变、应力以及加速度与时间关系如图5.3、图5.4、图5.5所示。
如图5.3,图5.4所示,可以发现:(1)随着子弹速度的增大,入射波应力幅值随之增大;(2)随着子弹速度的增大,入射波的第一峰值点在时间域内的位置不发生变化;(3)入射波近似为矩形波,入射波在刚开始稍有振荡现象,表明波形弥散现象存在,但不是很严重。