已阅读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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 光纤通信技术应用案例
- 企业经营管理优势分析
- 会展调查报告案例分析
- 山东省沂源县市级名校2021-2022学年中考冲刺卷化学试题含解析
- 2022年江苏省常州市金坛区白塔中学初中化学毕业考试模拟冲刺卷含解析
- 2024-2029年通信设备行业市场现状供需分析及市场深度研究发展前景及规划战略投资分析研究报告
- 公司年度总结与计划
- 2024-2029年起重机(装置)行业市场现状供需分析及市场深度研究发展前景及规划投资研究报告
- 2024-2029年贵铁合金行业市场现状供需分析及重点企业投资评估规划分析研究报告
- 2024-2029年豪华乙烯基地板(LVT)行业市场现状供需分析及重点企业投资评估规划分析研究报告
- 保安服务项目背景及需求分析
- 物业管理服务项目投标技术方案(投标方案)
- 新时代劳动教育教程(中职版劳动教育)全套教学课件
- 《婴幼儿感觉统合训练》课件-触觉
- 新生儿黄疸个案护理
- 2024中国邮政集团限公司甘肃省分公司校园招聘高频考题难、易错点模拟试题(共500题)附带答案详解
- 区域检验中心建设方案
- 加快建设教育强国科技强国人才强国 (模板)
- 电气仪表部分施工方案培训资料
- 反渗透膜行业市场分析
- 人教版英语中考模拟试题与答案解析(共三套)
评论
0/150
提交评论