已阅读5页,还剩1页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
JournalofMaterialsProcessingTechnology140(2003)616621EfficiencyenhancementinsheetmetalforminganalysiswithameshregularizationmethodJ.H.Yoon,H.HuhDepartmentofMechanicalEngineering,KoreaAdvancedInstituteofScienceandTechnologyScienceTown,Daejeon305-701,SouthKoreaAbstractThispapernewlyproposesameshregularizationmethodfortheenhancementoftheefficiencyinsheetmetalforminganalysis.Theregularizationmethodsearchesfordistortedelementswithappropriatesearchingcriteriaandconstructspatchesincludingtheelementstobemodified.Eachpatchisthenextendedtoathree-dimensionalsurfaceinordertoobtaintheinformationofthecontinuouscoordinates.Inconstructingthesurfaceenclosingeachpatch,NURBS(non-uniformrationalB-spline)surfaceisemployedtodescribeathree-dimensionalfreesurface.Onthebasisoftheconstructedsurface,eachnodeisproperlyarrangedtoformunitelementsascloseastoasquare.Thestatevariablescalculatedfromitsoriginalmeshgeometryaremappedintothenewmeshgeometryforthenextstageorincrementalstepofaforminganalysis.Theanalysisresultswiththeproposedmethodarecomparedtotheresultsfromthedirectforminganalysiswithoutmeshregularizationinordertoconfirmthevalidityofthemethod.2003ElsevierB.V.Allrightsreserved.Keywords:Meshregularization;Distortedelement;NURBS;Patch;Finiteelementanalysis1.IntroductionNumericalsimulationofsheetmetalformingprocessesenjoysitsprosperitywithaburstofdevelopmentofthecom-putersandtherelatednumericaltechniques.Thenumericalanalysishasextendeditscapabilitiesforsheetmetalformingofcomplicatedgeometrymodelsandmulti-stageforming.Inthecaseofacomplicatedgeometrymodel,however,severelocaldeformationoccurstoinducetheincreaseofthecom-putingtimeanddeterioratetheconvergenceoftheanalysis.Distortionandseveredeformationofthemeshgeometryhasaneffectonthequalityofforminganalysisresultsespeciallyinthecaseofmulti-stageforminganalysiswhenthemeshgeometryformedbytheforminganalysisatthefirststageisusedfortheforminganalysisatthenextstage.Thisillbehaviorofthedistortedmeshcanbeavoidedbytherecon-structionofthemeshsystemsuchasthetotalortheadaptiveremeshingtechniques.Theadaptiveremeshingtechniqueisknowntobeanefficientmethodtoreducedistortionofelementduringthesimulation,butitstillneedstremen-douscomputingandputsrestrictionsamongsubdividedelements.Correspondingauthor.Tel.:+82-42-869-3222;fax:+82-42-869-3210.E-mailaddress:hhuhkaist.ac.kr(H.Huh).Effectivemethodstoconstructameshsystemhavebeenproposedbymanyresearchers.Typicalmethodscouldber-method1inwhichnodalpointsareproperlyrearrangedwithoutthechangeofthetotaldegreesoffreedomofthemeshsystem,h-method2inwhichthenumberofmeshesisincreasedwithelementsofthesamedegreesoffreedom,andp-method3inwhichthetotaldegreesoffreedomofthemeshsystemisincreasedtoenhancetheaccuracyofso-lutions.SluiterandHansen4andTalbertandParkinson5constructedtheanalysisdomainasacontinuousloopandcreatedelementsinsub-loopsdividedfromthemainloop.Lo6constructedtriangularelementsinthewholedomainandthenconstructedrectangularelementsbycom-biningadjacenttriangularelements.Inthispaper,ameshregularizationmethodisnewlyproposedinordertoenhancetheefficiencyoffiniteele-mentanalysesofsheetmetalforming.Themeshregular-izationmethodautomaticallyfindsoutdistortedelementswithsearchingcriteriaproposedandcomposespatchestobemodified.Eachpatchisthenextendedtothree-dimensionalsurfacesinordertoobtaintheinformationofthecontinuouscoordinatesonthethree-dimensionalsurface.Thesurfaceenclosingeachpatchisdescribedasathree-dimensionalfreesurfacewiththeuseofNURBS(non-uniformrationalB-spline).Onthebasisoftheconstructedsurface,eachnodeisproperlyarrangedtocomposeregularelementsclosetoasquare.Thestatevariablescalculatedfromitsoriginalmesh0924-0136/$seefrontmatter2003ElsevierB.V.Allrightsreserved.doi:10.1016/S0924-0136(03)00801-XJ.H.Yoon,H.Huh/JournalofMaterialsProcessingTechnology140(2003)616621617geometryaremappedintothenewmeshgeometryfortheforminganalysisatthenextstage.Numericalresultscon-firmtheefficiencyoftheproposedmethodandtheaccuracyoftheresult.Itisalsonotedthatthepresentmethodiseffec-tiveinthecrashanalysesofsheetmetalmembersobtainedfromtheformingsimulation.2.RegularizationofthedistortedelementTheregularizationproceduretomodifydistortedele-mentsisintroducedinordertoenhancetheefficiencyofanalysisforthenextfiniteelementcalculation.Thedistortedelementsareselectedwithappropriatesearchingcriteriaandallocatedtoseveralpatchesforregularization.Thepatchesareextendedtothree-dimensionalsurfaceswiththeuseofNURBSforfullinformationofthecontinuouscoor-dinatesonthethree-dimensionalsurface.Onobtainingthenewcoordinatesofeachnode,thedistortedelementsareregularizedtoaregularelementthatisclosetoasquare.2.1.ThecriterionofmeshdistortionDistortedmeshesareselectedwiththetwogeometricalcriteria:oneistheinnerangle;andtheotheristheaspectratiooftheelementside.2.1.1.InnerangleTheinnerangleofaquadrilateralelementshouldbeclosetotherightangleforgoodresultsfromfiniteelementcalcu-lation.Zhuetal.7definedthereasonableelementwhenthefourinneranglesareformedwiththeangleof9045whileLoandLee8proposedtheinnerangleof9052.5asthesamecriterion.ThecriterionofmeshdistortionfortheinnerangleisdeterminedbyconstitutingEq.(1).AmeshisregardedasdistortedwhenEq.(1)islessthan/3or(i)maxinEq.(3)9isgreaterthan/6.Thecriterionisratherstrictinordertoavoidthegeometricallimitationincaseofapplyingtheregularizationmethodinconfinedregions:vectorfQ=1e1+2e2+3e3+4e4(1)|vectorfQ|=radicaltpradicalvertexradicalvertexradicalbt4summationdisplayi=1(i)2(2)i=vextendsinglevextendsinglevextendsingle12ivextendsinglevextendsinglevextendsingle(3)2.1.2.AspectratiooftheelementTheidealaspectratiooftheelementsideshouldbeunitywhenthefoursidesofanelementhavethesamelength.TheaspectratioisdefinedasEq.(4)andthenthedistortionisdefinedwhenitislessthan5thatcouldbemuchlessforastrictcriterion:maxr12,r23,r34,r41minr12,r23,r34,r41(4)whererijisthelengthofeachelementside.Fig.1.Processforconstructionofapatch.2.2.Domainconstruction2.2.1.ConstructionofthepatchDistortedelementsselectedbythecriteriaofmeshdistor-tionaredistributedinvariousregionsaccordingtothecom-plexityoftheshapeofformedgeometry.Theseelementsareallocatedtopatchesconstructedfortheefficiencyofthealgorithm.Theshapeofpatchesismadeupforrectangu-larshapesincludingalldistortedelementsforexpandingtheregionofregularizationandapplyingtoNURBSsurfaceex-plainedinnextsection.ThisprocedureisshowninFig.1.Whenholesandedgesarelocatedbetweendistortedele-ments,theregionsarefilleduptomakepatchesarectangularshape.Thepatchisthenmappedtoathree-dimensionalfreesur-facebyusingNURBSsurface.Theprocedureisimportanttoobtainentireinformationofthecontinuouscoordinatesonthethree-dimensionalsurface.NURBSsurfacecande-scribethecomplexshapequicklybyusinglessdatapointsanddoesnotchangetheentiredomaindataduetothelocalchange.2.2.2.NURBSsurfaceNURBSsurfaceisgenerallyexpressedbyEq.(5)asthep-orderintheu-directionandtheq-orderinthev-direction10:S(u,v)=summationtextni=0summationtextmj=0Ni,p(u)Nj,q(v)wi,jPi,jsummationtextni=0summationtextmj=0Ni,p(u)Nj,q(v)wi,j(5)wherePi,jisthecontrolpointsastheu-,v-direction,wi,jtheweightfactorandNu,p(u),Nj,q(u)thebasisfunctionthatareexpressedbyEq.(6):Ni,0=braceleftBigg1ifuiuui+1,0otherwise,Ni,p(u)=uuiui+puiNi,p(u)+ui+p+1uui+p+1ui+1Ni+1,p1(u)(6)Inordertomapthenodesfromthepatchesontothecon-structedsurface,anumberofpointsarecreatedfortheircoordinatesontheNURBSsurface.Thelocationofeach618J.H.Yoon,H.Huh/JournalofMaterialsProcessingTechnology140(2003)616621Fig.2.Selectingdirectionofdistortedelements.movingnodebyapplyingaregularizationmethodisdeter-minedsuchthatthelocationofapointhastheminimumdistancebetweennodesonNURBSsurface.Theinforma-tiononthecoordinatesofthenodalpointstobemovedisstoredtoconstructanewmeshsystem.2.3.RegularizationprocedureTheregularizationmethodiscarriedoutwiththeunitofapatchthatformsarectangularshape.FiniteelementstoberegularizedisselectedbytheorderofFig.2.Eachselectedelementisdividedbytwotriangularelementsandthenthedividedelementismadeofarighttriangularelementbyrelocatingthevertexonthecirclehavingthediameterfromvectorx1tovectorx2asshowninEq.(7)andFig.3.Whentheprocedureterminates,thesameprocedureisrepeatedintheoppositedirection:vectorx1+vectorx22=vectorxcen,|vectorx1vectorx2|2=r,vectorxcurvectorxcen=vectorxdir,vectorxnew=vectorxdir|vectorxdir|rfactor+vectorxcen(7)Thefinallocationofanoderelocatedbyusingthereg-ularizationmethodissubstitutedforthelocationofapointonNURBSsurface.Aftertheregularizationproce-dureisfinished,asimplesoothingprocedureiscarriedoutbyEq.(8)fortheroughregiongeneratedduringtheFig.3.Rcedure:PN=summationtextNi=1AiCisummationtextNi=1Ai(8)wherePNisthecoordinateofanewnode,AitheareasofadjacentelementsandCithecentroidoftheadjacentelements.2.4.LevelofdistortionAsadistortionfactor,levelofdistortion(LD)isnewlyproposed.LDcanbeusedtoevaluatethedegreeofimprove-mentintheelementquality:LD=AB(9)whereA=summationtext4i=1|sini|4,B=tanhparenleftbigkBprimeparenrightbig(10)Bprime=minr12,r23,r34,r41maxr12,r23,r34,r41,k=tanh1()(11)LDhasthevaluebetween0and1;whenLD=1,theele-mentisanidealelementofasquareandwhenLD=0,thequadrilateralelementbecomesatriangularelement.iarethefourinneranglesofanelement,soAisthefactorfortheinnerangle.BisthefactorfortheaspectratioofelementsidesandisdefinedbythehyperbolictangentfunctioninordertomakeLDlesssensitivetothechangeofB.Forex-ample,whenthereasonableaspectratiooftheelementsideis1:4,thevalueofBcanbeadjustedbyapplying=0.25and=0.6suchthattheslopeofthefunctionBischangedabruptlyaroundthevalueofBprime=0.25.Consequently,thevalueofLDdecreasesrapidlywhentheaspectratioBprimeislessthan0.25whilethevalueofLDincreasesslowlywhentheBprimeisgreaterthan0.25.Thisschemecanregulatethein-nerangleandtheaspectratiotohavetheequaleffectontheLD.2.5.MappingofthestatevariablesWhentheregularizedmeshsystemisusedforthenextcalculationoftheforminganalysisorthestructuralanaly-sis,mappingofthestatevariablesisneededformoreaccu-rateanalysisconsideringthepreviousforminghistory.Themappingprocedureistomapthecalculatedstatevariablesintheoriginalmeshsystemontotheregularizedmeshsys-tem.AsshowninFig.4,asphereisconstructedsurround-inganewnodesuchthatthestatevariablesofnodesinthespherehaveaneffectonthestatevariablesofthenewnode.ThestatevariablesofthenewnodearedeterminedfromthestatevariablesoftheneighboringnodesinthespherebyimposingtheweightingfactorinverselyproportionaltothedistancebetweenthetwonodesasshowninEq.(12):Vc=summationtextmj=1Vj/rjsummationtextmi=11/ri(12)J.H.Yoon,H.Huh/JournalofMaterialsProcessingTechnology140(2003)616621619Fig.4.Controlsphereformappingofthestatevariables.whereVjisthestatevariablecalculatedontheoriginalmeshsystem,andrjthedistancebetweenthenewnodeandtheneighboringnodes.3.Numericalexamples3.1.ForminganalysisofanoilpanWhileoilpansareusuallyfabricatedwithatwo-stageprocessinthepressshop,thepresentanalysisiscarriedoutwithasingle-stageprocessasshowninFig.5thatdescribesthepunchanddieset.Theregularizationmethodcanbeappliedtothefiniteele-mentmeshsystemwheneverneededforenhancementofthecomputationefficiency.Inthisexamplefordemonstration,themethodisappliedtotheanalysisofoilpanformingattwoformingintervalsforregularizationofdistortedmeshesasdirectedinFig.6.Fig.7explainstheprocedureoftheregularizationmethod.Fig.7(a)showsthedeformedshapeatthepunchstrokeofFig.5.Punchanddiesetforoilpanforming.60%forming6080%forming100%formingRegularizationRegularizationFig.6.Applyingtheregularizationmethodtotheforminganalysis.Fig.7.Procedureofregularization:(a)searchingdistortedelements;(b)constructingpatchesfordistortedelements;(c)regularizationofdistortedelements.60%andthreepartsofmeshdistortionbytheformingpro-cedure.Itindicatesthatthenumberofpatchestobecon-structedis3.Distortedmeshesareselectedaccordingtothetwogeometricalcriteriaformeshdistortio
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 图书馆改造装修工程施工组织设计
- 四川省阿坝藏族羌族自治州松潘县2023-2024学年小升初易错点语文检测卷含答案
- 冀教版五年级数学下册课件《异分母分数加减 第2课时》
- 《分布式生活垃圾中转站臭气处理技术规程》征求意见稿
- 2024年中考数学总复习:几何初步及三角形-知识讲解(提高)
- 林业与沿海防护林建设
- 黑色金属矿行业的技术水平与生产效率
- 2024年为战士喝彩作文
- 2024年中职招生工作总结(经典12篇)
- 2024年《金刚川》观后感12篇
- 2024年04月上海市药品和医疗器械不良反应监测中心2024年公开招考工作人员笔试参考题库附带答案详解
- 积极心理学:塑造刚健自信的中国青年智慧树知到课后章节答案2023年下上海思博职业技术学院
- 愚公移山英文 -中国故事英文版课件
- 2023储能柜通用技术规范
- 红色喜庆国风小孩周岁宴教育课件ppt模板
- 压铸件的缺陷分析及检验
- UMAT符号说明
- 幼儿园教师如何备课 (新).ppt
- 个人竞聘演讲稿ppt课件.ppt
- 血型课件.ppt
- 毒、麻、精神药品,高危药品管理.ppt
评论
0/150
提交评论