《应用LSDYNA进行薄板成形仿真.ppt》由会员分享，可在线阅读，更多相关《应用LSDYNA进行薄板成形仿真.ppt（65页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、应用应用LS-DYNA进行薄板进行薄板成形仿真成形仿真高级培训高级培训主要内容主要内容薄板成形工艺和术语薄板成形仿真的准备工作前处理阶段的工作仿真分析结果及其用途输入文件的常用模板其它常见的问题及处理方法建议的参考资料薄板成形工艺及术语薄板成形工艺及术语说明：请参看DYNAFORMAPPLICATIONMANUAL成形工艺包括拉延（DRAW）切边（TRIM）翻边（FLANGE）拉延正拉延(toggle)反拉延(inverted)压边(binderwrap)拉延筋(drawbead)回弹（springback）程序中使用的重要术语减薄（thinning）拉裂(fracture)起皱(wrinkl
2、e)成形极限图(FLD)薄板成形工艺及术语薄板成形工艺及术语正拉延：第一个动作第一个动作-上压边圈闭合，压紧板料。上压边圈固定于外上方的支架，由固定压力控制，通过压力调整可以减小或增大固定压力。压边圈闭合时，不对板料做下垂控制或反吹。第二个动作第二个动作-上凸模将板料拉延成形。上凸模固定于内压头，将板料从压边圈下拖进凹模型腔。下凹模是实心的并固定不动，反拉延第一个动作第一个动作上压边圈闭合以紧固板料并置拉延筋。板料由下凸模支撑以得到一个好的外形。当压边圈闭合时，与正拉延相比，对板料做下垂控制和反吹。第二个动作第二个动作-上压边圈压紧下压边圈，从而使板料绕凸模拉伸。该动作的目的在于用拉延筋将板料
3、固定于压边圈之间，防止板料移动。所有外凸形状均在该动作内实现。第三个动作第三个动作-上凹模闭合，对板料做二次拉伸，并拉延出最终的制件形状。与常规拉延或拉伸成形相比，反拉延多出一个动作。并非所有的零件都适于反拉延成形极限图成形极限图成形极限图通常的获得的方法：大量实验DYNAFORM的FLD图取自经验公式：如何制作如何制作FLD以上公式确定极限图的最低点。板料紧贴模具形状成形。然而，成形(切边)结束后，将工件从模具中取出时，其形状会发生变化，这种形状的变化称之为“回弹”。回弹现象随材料屈服应力、模具间隙、弯曲半径的增大而显著，随材料厚度的增加而减弱。在模具设计过程中必须考虑回弹因素，这样工件的最
4、终形状才能满足外观的设计要求及装配要求。不同材料的回弹程度也不同。回回弹弹影响工件总体回弹的两个组成部分是：弹弹性性回回弹弹冲压成形中的弹性变形，主要是由于凹凸模之间的间隙比板料厚度大的缘故。回弹角可能会很大，这将导致严重的形状偏差。屈服应力、模具间隙、板料厚度都是影响弹性回弹的重要因素。塑塑性性回回弹弹塑性变形区残余应力的释放导致应征偏差与回弹现象。影响塑性回弹的主要因素是应力和板料厚度。“二次弯曲”或所谓的“整形”成形，是消除简单工件回弹的常用方法。回弹示意图回弹的计算LS-NIKE求解：输入文件nikin;LS940以前LS-DYNA隐式求解器：与拉延同在一个文件中（seamlesssw
5、itchsolving）回弹量：任意节点在拉延构形和回弹构形中的位移。Resultantdisplacementplot.回弹方向：由原构形指向新构形.displacementvectorplot.薄板成形仿真的准备工作薄板成形仿真的准备工作-板料的材料应变的选择？（工程、真实）各向异性？（正交、纵向）屈服准则？（与各向异性性质相关）硬化准则？（曲线、幂指数、多项式）应变率效应？（是否明显？）选择LS-DYNA何种材料本构？*mat_3-parameter_barlat*mat_transversely_anisotropic_elastic_plastic*mat_fld_anisotrop
6、ic薄板成形仿真的准备工作薄板成形仿真的准备工作-板料的材料1.试验单向拉伸2.获得应力应变曲线3.拟和选定的LS-DYNA本构参数？由产品商提供？完全由试验测得常用的方法：薄板成形仿真的准备工作薄板成形仿真的准备工作Springback分析中的材料?-板料的材料薄板成形仿真的准备工作薄板成形仿真的准备工作-摩擦系数摩擦性能试验装置（？国内用户）能够充分反映操作过程中的涂油条件结构简单精确控制主要参数与模拟结果拟和薄板成形仿真的准备工作薄板成形仿真的准备工作-拉延筋单位长度的阻力*contact_drawbeadLCIDRF:确定板料在拉延筋弯曲展直变形的阻力,Fbending=B()LCID
7、NF:确定板料在下压过程中的法向力,Fnormal=N(),据此确定摩擦力,Ffriction=Fnormal拉延筋阻力F=Ffriction+Fbending薄板成形仿真的准备工作薄板成形仿真的准备工作-拉延筋单位长度的阻力如何确定两条曲线？1.DBFP？2.建立自己的数据库试验装置不同形式的拉延筋如何拟和试验与仿真的结果薄板成形仿真的准备工作薄板成形仿真的准备工作-压机参数压机的类型，模具结构压边力or压边间隙行程曲线（速度-时间）前处理阶段前处理阶段-几何模型的读入几何模型最好为IGES带工艺补充部分标明切边线对局部的修改分析模型较CAD/CAM精度高对丢失、缝隙、形态差的曲面进行修复减
8、少划网格的工作量前处理阶段前处理阶段应遵守基本的有限元建模准则，下面是一些具体要求：1.纵横比小于10；2.内角不超过30；3.翘曲角不超过20；4.尽量用四边形单元而避免使用三角形单元；5.部件单元的法线方向一致；6.部件边界检查；7.重单元检查；8.避免划分太小的单元。尽量使用DYNAFORM的网格EDIT功能获得理想的网格-单元的划分规则前处理阶段前处理阶段-自动定位，确定各部件之间的距离部件的行程测量时应该考虑到壳单元的厚度自动生成速度-时间运动曲线压边力曲线前处理阶段前处理阶段-其它问题是否有后续工序？是否考虑重力影响？是否使用自适应网格？何时使用？仿真结果及其利用仿真结果及其利用仿
9、真结果：给定板料在给定的模具和工艺条件下的形状及成形状态。厚度分布（减薄率）应变分布（FLD）起皱缺料Skidmarks(划伤)回弹量残余应力成形力，压边力涉及到的人员包括：车身零件设计工程师工艺设计工程师模具设计工程师分析工程师仿真结果及其利用仿真结果及其利用在设计工作的早期阶段评价覆盖件及其模具设计、工艺设计的可行性(assessfeasibility)；在试冲试模阶段进行故障分析，解决问题(troubleshooting)；在批量生产阶段用于缺陷分析，改善覆盖件生产质量(improvequality)，同时可用来调整材料等级。例例1：设计拉延筋位置：设计拉延筋位置设计1设计2单位长度拉压
10、筋阻力相同，拉延筋的布置不同单位长度拉压筋阻力相同，拉延筋的布置不同板料厚度减薄率对比板料厚度减薄率对比设计1：减薄率50%设计2：减薄率40%成形极限图对比成形极限图对比设计1：未通过设计2：通过厚度分布成形极限图例2微型车顶盖的可成形性-其它应用仿真结果及其利用仿真结果及其利用成形力-选择压机压边力-控制阻力是否存在局部开裂？板料的坯料形状？板料在后续工序中的成形剩余空间？如果达到极限，在后续工序中容易开裂。输入文件的常用模板输入文件的常用模板*KEYWORD*TITLEtemplate for metal forming$use this input deck as a guide fo
11、r setting up metal forming$simulations with LS-DYNA$-$-$set the termination time,mass scaling time step size$*CONTROL_TERMINATION$ENDTIM$*CONTROL_TIMESTEP$DTINIT TSSFAC ISDO TSLIMT DT2MS$-$-$The adaptive frequency should be small enough to adapt every 2-4 mm of punch$travel.*CONTROL_ADAPTIVE$ADPFREQ
12、 ADPTOL ADPOPT MAXLVL TBIRTH TDEATH LCADP IOFLAG x.x 4.000 2 3 0.0 0.0 0 0$ADPSIZE ADPASS IREFLG ADPENE*CONTROL_HOURGLASS$IHQ QH*CONTROL_SHELL$WRPANG ITRIST IRNXX ISTUPD THEORY BWC MITER 20.000 2 -1 1 2 2 1*CONTROL_CONTACT$SLSFAC RWPNAL ISLCHK SHLTHK PENOPT THKCHG ORIEN 0.010 2 1 4 2 3$USRSTR USRFAC
13、 NSBCS INTERM XPENEN$select an output frequency,additional databases.Add strain to plot file.$*DATABASE_RCFORC$dt *DATABASE_BNDOUT$dt *DATABASE_BINARY_D3PLOT$dt *DATABASE_EXTENT_BINARY$neiph neips maxint strflg sigflg epsflg rltflg engflg 0 0 0 1 0 0 0 0$cmpflg ieverp beamip dcomp shge stssz n3thdt
14、0 0 0 0 0 0 0$-$-$-$Create the output file dynain,including part#1 only.Add threeCreate the output file dynain,including part#1 only.Add three$nodes for constraining rigid body motion during springback.Code 7=xyz,$nodes for constraining rigid body motion during springback.Code 7=xyz,$code 5=yz,and c
15、ode 3=z-direction constrained.$code 5=yz,and code 3=z-direction constrained.$*INTERFACE_SPRINGBACK_DYNA3D*INTERFACE_SPRINGBACK_DYNA3D$psid$psid 6 6$nid code$nid code xxx 7 xxx 7 xxx 5 xxx 5 xxx 3 xxx 3$*SET_PART_LIST*SET_PART_LIST 6 6 1 1$-$-*CONTROL_IMPLICIT_GENERAL$imflag dt0 iefs 0 0.00 1*CONTROL
16、_IMPLICIT_NONLINEAR$nlsolvr ilimit maxref dctol ectol rctol lstol 0 0 5 0.002 0.00 0 0$dnorm divflag inistif 0 0 0*CONTROL_IMPLICIT_LINEAR$lsolvr prntflg negeig 0 0 0*CONTROL_IMPLICIT_AUTO$iauto iteopt itewin dtmin dtmax 1 0 0 1.0e-4 0*CONTROL_IMPLICIT_STABILIZATION$ias scale tstart tend 1 0.0 0 0$-
17、$enter material properties in parenthesis(xx)$sheet*MAT_TRANSVERSELY_ANISOTROPIC_ELASTIC_PLASTIC 1 (rho)(youngs)(pr)(sigy)0.0000000 (r)1$*DEFINE_CURVE 1$enter material stress vs PLASTIC strain$die*MAT_RIGID 2 (rho)(youngs)(pr)$enter initial sheet thickness$*SECTION_SHELL$SECID ELFORM SHRF NIP PROPT
18、QR/IRID ICOMP 1 16 0.0000000 7.0000000 0.0000000 0.0000000 0 (thick)(thick)(thick)(thick)$*SECTION_SHELL$SECID ELFORM SHRF NIP PROPT QR/IRID ICOMP 2 2 0.0000000 1.0000000 0.0000000 0.0000000 0$-$*PARTsheet$PID SID MID EOSID HGID GRAV ADPOPT TMID 1 1 1 0 0 0 1 0*PARTdie$PID SID MID EOSID HGID GRAV AD
19、POPT TMID 2 2 2 0 0 0 0 0$punch displacement*BOUNDARY_PRESCRIBED_MOTION_RIGID$pid dir vad lcid scale*DEFINE_CURVE 2$time punch disp$-$binder load*LOAD_RIGID_BODY$pid dir lcid scale*DEFINE_CURVE 3$time binder load$enter friction coefficient fs for each contact interface:$sheet vs die$*CONTACT_FORMING
20、_ONE_WAY_SURFACE_TO_SURFACE$ssid msid sstyp mstyp sboxid mboxid spr mpr 1 2 3 3 0 0 0 0$fs fd dc vc vdc penchk bt dt$sfs sfm sst mst sfst sfmt fsf vsf$sheet vs punch$sheet vs binder$-DRAW BEAD LINE 1-*CONTACT_DRAWBEAD$PROPNAME BEADPRO 15$SSID MSID SSTYP MSTYP SBOXID MBOXID SPR MPR 1 1 4 2 1 0 0$FS F
21、D DC VC VDC PENCHK BT DT 0.100E+00 0.000E+00 0.000E+00 0.000E+00 0.200E+02 0.000E+00 0.100E+21$SFS SFM SST MST SFST SFMT FSF VSF 0.100E+01 0.100E+01$LCIDRF LCIDNF DBDTH DFSCL NUMINT 12 0 0.600E+01 0.250E+01*DEFINE_BOX 1 0.390E+03 0.791E+03-.725E+03 0.200E+03 0.989E+03 0.142E+04*CONSTRAINED_EXTRA_NOD
22、ES_SET$PID NSID 6 1*SET_NODE_LIST*NODE$enter or include node data$*ELEMENT_SHELL$enter or include element data$use part numbering:1=sheet,2=die,3=punch,4=binder$*ENDDYNAFORM常见的问题常见的问题及处理办法及处理办法1.怎样实现多工序分析怎样实现多工序分析?两种方法:a)一个文件中定义所有PART，定义边界条件和工艺条件。同时定义不同接触界面的生、死时间。(卡片*CONTACT_xxxxx,)。这种方法使得输入文件相当复杂，容
23、易出错。建议有相当经验的用户可以使用。b)通过LSDYNA生成“dynain”文件的功能实现。这个文件包括板料前一个工序结束时所有的重要信息，包括厚度变化、应力分布、网格形态等。此时，输入文件中可将前一个工序的PART直接INCLUDE到文件中，容易操作。易于使用，避免出错。2.做回弹分析时，求解不收敛，怎么办做回弹分析时，求解不收敛，怎么办?以下几个方法用于提高收敛性a)采用更多个载荷部。由*CONTROL_IMPLICIT_GENERAL关键字的NSTEPSB参数给定。b)调整*CONTROL_IMPLICIT_STABILIZATION的参数SCALE。如果第一个载荷步即不收敛，选择较大
24、的数值;如果接近结束时不收敛，选择较小的数值。如果还不收敛，增加载荷步。c)调整*CONTROL_IMPLICIT_NONLINEAR关键字的LSTOL参数。多数情况，选择较小的数。d)改变*CONTROL_IMPLICIT_LINEAR关键字的LSOLVER。不同的求解器具有不同的收敛行为。e)改变约束条件。事实证明不同的约束影响收敛速度。可能需要尽量减少总位移。f)如果求解因为出现负的特征值，CONTROL_IMPLICIT_LINEAR中的NEGEIG给为缺省值.g)检查网格形态，是否存在切边过程产生过度差的单元。如果存在，修改单元。h)检查是否存在过分细长单元，应用表明这样的单元往往导
25、致收敛出现问题。3.当用当用seamlessmethod做做forming和和springback分析时，得不到收敛的分析时，得不到收敛的解解.而且不希望再做一次成形，怎样利用前面的结果只做回弹分析？而且不希望再做一次成形，怎样利用前面的结果只做回弹分析？两种方法:a)d3dumpXX生成dynain文件,然后写一个小的输入文件，包括此dynain文件.，同时文件中应包含适合的约束条件，然后进行回弹的求解。b)下面介绍由d3dumpXX生成dynain：c)写一个输入文件restart.k，包括*INTERFACE_SPRINGBACK_DYNA3D关键字d)执行求解器，ls950I=rest
26、art.kr=d3dumpxxmemory=?e)b)写一个输入文件，包含要修改的控制卡片，用上面的格式求解：ls950I=restart.kr=d3dumpxxmemory=?4.当按问题当按问题3的第二种方法做回弹分析时的第二种方法做回弹分析时,但求解时程序退出，但求解时程序退出，message文件说计算正常结束。这是什么原因？文件说计算正常结束。这是什么原因？当从dump重启动进行回弹分析时,如果在拉延分析过程中使用自适应网格，程序启动后首先寻找一在拉延分析中自适应文件，如果adaptive设置时其DEATH时间小于termination时间,将不会生成此文件,回弹分析也不能进行。可采用
27、如下两种方法，避免这个问题:a)在adaptive设置时使用缺省的时间设置。缺省的death是个非常大的数值，大于拉延时间。b)由d3dumpxx生成dynain文件,由这个文件用问题3第一个方法进行回弹分析。5.在在LSDYNA中有许多接触类型供用户选择。在板成形分析中常中有许多接触类型供用户选择。在板成形分析中常使用那种类型使用那种类型?在冲压仿真中，常用的接触是专用于板成形的FORMING类接触,。在这种接触里，模具的厚度不考虑，同时，降低对模具网格的连续性要求。通常用如下两种类型：a)*CONTACT_FORMING_ONE_WAY_SURFACE_TO_SURFACE,通常使用这种类
28、型。b)*CONTACT_FORMING_SURFACE_TO_SURFACE.要用更多的时间。当使用前种方式时，如果模具的网格比板料的网格还要密，则往往出现较大的接触穿透，此时，可使用这种接触类型。6.当做拉延分析时，出现很明显的接触穿透，如何控制？当做拉延分析时，出现很明显的接触穿透，如何控制？可尝试以下几种方式：a)增大*CONTROL_CONTACT关键字中的SLSFAC参数，此参数通常范围在0.1之间。数值由于稳定性原因，一般不超过0.1。b)选择其它的接触类型。如果模具的网格比板料网格密，使用CONTACT_FORMING_SURFACE_TO_SURFACE。c)使用constr
29、aint接触类型。FORMING_ONE_WAY_SURFACE_TO_SURFACE,增加OptionialCardA，其参数SOFT=4即可。7.当使用当使用constrain类型的接触时，网格产生畸变或程序崩溃，如何处理?当使用这种接触类型时，在板料与模具的接触界面没有阻尼，所以应采用一些措施如*DAMPING_GLOBAL消除高频震荡，应用表明，效果较好。这类接触有一个限制，即板料两面中只能有一面定义。另一面要采用罚函数接触算法。8.在在LSDYNA中有许多壳单元的算法，在板成形分析中那些适中有许多壳单元的算法，在板成形分析中那些适合合?如果只是进行拉延分析，选择BT(element#
30、2)很快，结果可以接受。如果拉延之后要进行回弹分析，选择全积分单元(element#16)。好处是可获得较好的回弹结果及收敛效果；缺点是时间长。其它相关问题其它相关问题1.在FEMB27（DYNAFORM）中，所有要求输入文件名的操作，都可以在提示区使用Unix的ls,pwd,cd等功能列表当前路径文件、查看当前路径、改变当前路径。当改变当前路径后，以后的操作都在此路径中进行，即已不是启动FEMB的工作目录。2.FEMB（DYNAFORM）中Part的概念来源于求解器ls-dyna。其定义为材料类型和单元类型一致的单元组合。由此可知，实体单元与壳单元不能在一个同Part中共存。3.从空数据库开
31、始作业，必须首先创建Part，FEMB（DYNAFORM）将新创建的Part置为当前Part；若无特殊提示，新生的几何模型（点、线、面）和有限元模型（节点、单元）都归属为当前Part。因此在创建各种模型前要检查当前Part是否为所需要Part，否则，可由3DMOUDULE的Part菜单中的ChangeCurrentPart改变当前Part。4.在FEMB（DYNAFORM）中，所有的曲线包括应力应变曲线都称为LoadCurve，应力-应变曲线在定义Part的材料性质时输入。除材料应力-应变曲线，其他曲线都定义参量依时间变化，即曲线的横坐标为时间。曲线可以在提示区输入，也可以生成固定格式的文件，
32、然后读入到数据库中。当然，也可以将输入的曲线写出文件到当前路径中。5.一般地，工作站不连接打印机，因此要生成图形文件然后传到微机上打印。FEMB中的打印功能的使用如下：在SETUP中选择HARDCOPYDEFAULT/PRINTERNAME菜单，选择printtofile/GIF格式，然后由3D模块的PRINT打印，给定文件名后，即可生成图形区的GIF文件，在微机WINDOWS98的画笔中可以编辑及打印。6.建议在写出ls-dyna求解的关键字文件前由MODEL CHECK和NODEOPTION菜单检查模型，避免重叠单元、重节点、自由节点、不连续网格的存在。7.Ls-dyna中 的 壳 单 元
33、 缺 省 为 Belyschok_tsay,此 单 元 计 算 速 度 快，在DYNAFORM中缺省为沿厚度3个积分点。当进行板成形或汽车碰撞分析时，板壳的翘曲过大后应增加积分点到5、7。8.在DYNAFORM后处理菜单中，拾取3D模块的save可保存任一时步的变形网格，以Nastran文件格式保存。9.在进行板成形分析的过程中，拉延筋的模拟不要用几何曲面构造，应在DRAWBEAD菜单中生成拉延筋。用几何曲面构造拉延筋往往导致计算不稳定。10.在应用DYNAFORM进行多工序或回弹分析时，若计算由于negativearea而终止，可将读入的dynain和nikin的网格可进行局部的轻微修改。1
34、1.在DYNAFORM中，由凸模（或凹模）单元法向偏移生成凹模（或凸模）时两者的thickness应是实际模具间隙，一般为板料的1.2-1.4倍。12.在DYNAFORM中，凸模的行程应是由DEFINETOOL/AUTOPOSITION将所有模具部件自动定位后，由DEFINETOOL/MINIMUMDIST.测量得到的凸、凹模之间的距离减去模具间隙。14.Ls-dyna中时间积分步长与材料密度有关。在板成形分析中，可以适当增大板料材料的密度。同时，增加模具的运动速度也是一个办法。建议初学者不要使用。此方法在碰撞分析中不能使用。15.Ls-dyna中的37号材料是一种常用的弹塑性材料，适于模拟钣
35、金成形。在输入材料参数时，除真实应力-应变曲线外，一定要正确给定各向异性指数r值；硬化指数n可以不输入，但在后处理进行FLD分析时，材料模型简化为指数硬化方式，需要输入n。16.使用自适应网格功能时，计算时间将增加。为减少计算时间，可增大*control_adaptive关键字第一个参数（网格重划分的时间间隔值），并将第二排第二个参数置为1。建议的参考文献建议的参考文献1.MetalformingSimulationwithExplicitTimeIntegration-ApproachingSolutionsforIndustrywithLS-DYNAByKarlSchweizerhof.C
36、AD-FEM,Germany.1stInternationalLS-DYNAConference,1993.2.AnElaborateStudyoftheStampingofaCarBodyComponentUsingLS-DYNA3DByASrange,PThilderkvist,KMattiassonstInternationalLS-DYNAConference,1993.3.AsensitivetyStudyofsimulationParametersinSheetMetalFormingandSpringbackSimulationsusingLS-DYNAByMingF.Shi,D
37、ouglasG.PrinceandWenhuithInternationalLS-DYNAUserConference.4.ImplicitSpringbackCaculationusingLS-DYNAByBradleyMaker,LSTC,USA.5thInternationalLS-DYNAUserConference.1998.5.InputParametersforMetalformingsimulationsusingLS-DYNAByBradleyMakerandHaixinZhu,LSTC,USA.6thInternationalLS-DYNAUserConference.20
38、00.6.AnisotropicBehaviorandItsApplicationonSheetMetalStampingByWeilongHuandJeanneHe.ETA.USA.6thInternationalLS-DYNAUsersConference.2000.7.EquivalentDrawbeadanditsApplicationinOptimizationofAutobodyFormingProcess.ByXionghuiZhou,et.al.ShanghaiJiaotongUniv.China.6thInternationalLS-DYNAUsersConference.2000.Thankyou!