机械可靠性英文文献及翻译.docx

机械可靠性英文文献及翻译.docx

《机械可靠性英文文献及翻译.docx》由会员分享，可在线阅读，更多相关《机械可靠性英文文献及翻译.docx（10页珍藏版）》请在冰豆网上搜索。

机械可靠性英文文献及翻译

MechanicalReliabilitySimulation

JianfengXu,YalianXie,DanXu

ShanghaiInstituteofProcessAutomationInstrumentationShanghai,PRC

Abstract—Thisarticleprovidesanumericalmethodtosimulatethereliabilityofacomplicatedmechanicalsystemwithmultipleinputrandomparameters.ThemethodbasesonthetheoryofmechanicalreliabilityandusesacommercialFEMsoftware,ANSYS.ThequantitativereliabilityofthemechanicalsystemisobtainedbyMonteCarlosimulation.Themethodisverifiedbytheoreticalanalysisandisappliedtoarealstructure.Thismethodhaspreventedthetediousmathematicalcalculationandismoreeconomicalthanexperimentalmethods.Itcanbeusedasacomplementoftraditionalexperimentalmethods.

Keywords—MechanicalReliability;FiniteElement

Method;MonteCarloSimulation

I.INTRODUCTION

BasedonthenationalstandardGB/T3187-1994[1],reliabilityhasbeendefinedas“theprobabilitythataproductwillsatisfactorilyperformitsintendedfunctionundergivencircumstancesforaspecifiedperiodoftime”.Reliabilityhasbeenviewedasoneofthemostcriticalparametersfortheproducts.Forexample,itissaidinUSthatreliabilitywillbethefocusofcompetitionforthefutureproducts.MostJapanesecompaniestakemorereliableproductastheirmainobjectiveforcompanydevelopment.Therefore,toaccuratelypredictthereliabilityoftheproduct,andtooptimizethedesignbasedonthereliabilityrequirementsisverymeaningful.

Generallyamechanicalsystemiscomposedbymanyparts.Asystem’sreliabilityisdeterminedbymanyparameters,suchasthenumberofcomponents,thequalityofthecomponentsandalsohowthosecomponentsareassembled.Itisalsorelatedtothematerialproperties,manufacturingtechnologies,dimensiontolerancesandsoon.Mostofthereferences

currentlyavailablearefocusedonthetheoreticalderivationforreliabilitycalculation,however,duethecomplexity,itrequirestediousmathematicalcalculation.Anothermethodtoobtainthereliabilityofasystemisbyexperiments,whichrequiresenoughsamplesandtestingtime,thustheexperimentalmethodcanonlybeapplicabletothosemostcriticalsystems.

Thisarticleusesmechanicalstrengthreliabilityasanexample,illustrateshowFEMcanbeappliedtothereliabilitysimulation.ThemethodemployscommercialFEMsoftware,ANSYS,andtakesadvantageoftheintegratedMonteCarlosimulationfunction.Asimpleexampleisanalyzedtoverifythecorrectnessofthemethod.Theverifiedmethodisthenappliedtoarealstructuretoobtainthequantitativereliabilityofacomplexsystemwithmultipleinputrandomparameters.

II.STRENGTH-STRESSINTERFERENCETHEORY

Normallythestructuralmaximumstressandthematerialstrengtharedistributedwithcertainprobabilitydensity.Figure1showsanexampleofthedistribution,in

whichfσandgsrepresentthemaximumstress

andmaterialstrength,respectively.InFigure1（a）,since

gsisalwayslargerthanfσ,thereliabilityofthestructureis1.However,inFigure1（b）,itisnot

guaranteedthatgsbelargerthanfσ.Thereisan

interferencezoneinthedistribution,asshownintheshadedareaandzoomedinFigure2.Thusthereliabilityofthestructureislessthan1.Basedonthemathematicalderivation,thereliabilitycanbecalculatedbyintegrationwiththedistribution.

978-1-61284-666-8/11$26.00

2011IEEE1147

RPσs−∞fσσgsdsdσ

fσgs

III.METHODDESCRIPTION

ThecommercialFEMpackageANSYShasaPDS（ProbabilisticDesignSimulation）module[2].Themoduledefinesdimension,materialpropertiesandotherparametersasrandomvariableswithcertaindistributions.Theoutputstressanddisplacementcanbedefinedasresultsvariables.SamplesarecollectedbasedonMonte

（a）

fσ

（b）

gs

Carlomethodtoobtainthedistributionofoutputvariables.

WiththePDSmodule,asuitableanalysisfileshouldbefirstlygeneratedtodefinethesimulationprocessandtodefinetheinputandoutputvariables.Ifnecessary,correlationbetweeninputparameterscanalsobedefined.MonteCarlomethodisthenemployedtocollectenoughsamplesforsimulation.Sinceeachrunisindependent,the

Figure1Stressandstrengthdistribution

gs

fσ

Figure2Stress-strengthinterferencezone

Tocreateanalysisfileandtodefineparameters

Todefinetheinputrandomparametersanddistributions

Todefinecorrelationsbetweeninput

Todefineoutputrandomvariables

Todefinemethodforprobabilisticsimulation

Responsesurfacesimulation

Resultsreview

Figure3.Flowchartforstrength

methodissuitableforparallelruntosavetime.Generally50~200samplesarecollectedinMonteCarlomethodtoinsurethecorrectnessoftheresults.Insomecases,theoutputvariableisasmoothfunctionoftherandominputvariables.Insuchcases,theresponsecanbeapproximatelyevaluatedbyamathematicalfunction.Byassessingthecoefficientsinthefunctionwiththesamplingpointsandruns,responsesurfacecanbedeterminedandnumerousrunscanbeconductedtoobtainbetterdistributionswithlimitedtime.BasicallytheflowcharttosimulatethestrengthreliabilitywithANSYScanbeshowninFigure3,inwhichthedashedstepsindicatetheyareoptional.

IV.VERIFICATIONSANDAPPLICATIONS

Thecorrectnessofthemethodisverifiedbyatensionrodwhichisanexampleinthereference[3].The

tensionP~N29400,441N,materialstrength

S~N1054.48,41.36N

mm2,radiusoftherod

r~N3.195,0.016mm,theoreticalcalculation

showsthereliabilityoftherodis0.999.

1148

reliabilityofapressurevessel.TheshapeandparameterdistributionsareshowninFigure6.

Barelement

P

Figure4.FEMmodel

（a）MonteCarloSimulation

（b）ResponseSurfaceSimulation

Figure5.Histogramofstrength-stressdifference

WithANSYS,therodismodeledasabarelementwhichisfixedinonend.Thetensionisappliedattheotherend.TheaxialstressiscalculatedinANSYS.Thetensionforce,radiusandmaterialstrengtharedefinedasrandominputvariablesinPDS.Theoutputvariableisdefinedasthedifferencebetweenthestrengthandmaximumstress.Avaluelessthan0indicateastrengthfailure.OnehundredsamplesarecollectedwithMonteCarlomethod.Resultsshowthatthedifferenceapproximatelyfollowsanormaldistributionwithnormal137.73Nandvariance45.42N.Lateron,responsesurfaceareobtainedbasedontheprevioussamplesand10000runsareconducted.Thenormalandvariancecanbebetterobtainedas137.65Nand44.52N,respectively.Fromwhichitisconcludedthatthereliabilityoftherodis0.9988and0.9990withthetwomethods.Thisresultmatcheswellwiththetheoreticalcalculationwhichdemonstratesthecorrectnessofthemethod.Figure4showstheFEMmodel,andFigure5isthehistogramofthestrength-stressdifference,whereKPaisunit.

Themethodisthenappliedtoanalyzethe

r

L

Fixed

r~U1990,2010mmL~U9990,10010mm

Thickness~U7.5,8.5mm

Pressure~U2800,3200KPaYoung'sModule~N200,10GPaYieldStres~N1000,50MPa

Figure6PressureVessel

InPDSmodule,therandominputparametersaredefinedinaccordancewithFigure6.With50MonteCarlosamplingpoints,thesimulationshowsthatthedifferencebetweenstrengthandstresscanbeapproximatelydescribedasanormaldistribution

N229.22,67.251MPa.Afurtherresponsesurface

simulationdemonstratesthatthedifferencematchthedistributionasN228.86,65.823MPa,asshownin

Figure7.Fromthesedistributions,thereliabilityofthe

1149

pressurevesselcanobtainedas0.99967and0.99974,

respectively.

（a）MonteCarloSimulation

（b）ResponseSurfaceSimulation

Figure7.Strength-stressdeviationhistogramofpressurevessel

V.CONCLUSIONS

FEMandMonteCarlosimulationmethodcanbeintegratedtoobtainthemechanicalreliabilityofacomplexstructurewithmultipleinputrandomparameters.Althoughthisarticleconsidersonlythestaticcases,themethodcanalsobeappliedtodynamiccaseswithsuitableanalysisfile.

Themethoddescribedinthisarticlepreventstediousmathematicalcalculation.Thestructurereliabilitycanbeobtainedwithlesstimeandmoney.Thismethodcanbeacomplementofthetraditionalexperimentalmethods.

REFERENCES

[1]ReliabilityandMaintainabilityterms,GB/T3187-94

[2]ANSYSHelp,ANSYSCorp

[3]LiuP.,BasicofReliabilityEngineering,MetrologyPress,2002

中文翻译：

机械可靠性仿真

JianfengXu,YalianXie,DanXu

上海工业自动化仪表研究所

上海，中华人民共和国

摘要：

本文提供了一种多输入随机参数模拟复杂机械系统的可靠性计算方法。

该方法基于可靠性理论，采用商业有限元软件ANSYS。

机械系统的可靠性定量是通过蒙特卡罗模拟。

该方法是通过理论分析验证和实际应用的结构。

该方法避免了繁琐的数学计算和比实验方法更经济。

它可以作为传统实验方法的一种补充。

关键词：

机械可靠性；有限元；蒙特卡罗仿真方法。

一、引言

根据国家标准GB／T3187-1994[1]，可靠性被定义为“在一定的条件下在指定的一段时间内，一个产品很好地履行其预期的功能的概率”。

可靠性一直被看作是一个产品的最重要的参数。

例如，据说在美国，可靠性将是今后产品竞争的焦点。

大多数日本企业以更可靠的产品作为公司发展的主要目标。

因此，为了准确地预测产品的可靠性，并优化基于可靠性要求的设计是非常有意义的。

一般的机械系统是由许多部分。

一个系统的可靠性取决于许多参数，如组件的数量，质量的部件和这些部件组装。

这也与材料的性质，相关的制造技术，尺寸公差等有关。

大多数现有的文献主要集中在可靠性理论推导计算，然而，由于它的复杂性，需要繁琐的数学计算。

还有一种实验法来获得系统的可靠性试验，需要足够的样本和测试时间，因此，实验方法只能适用于那些最重要的系统。

本文采用机械强度可靠性为例，说明了有限元法可以应用于可靠性仿真。

该方法采用商业有限元软件，ANSYS，并利用MonteCarlo模拟功能集成。

一个简单的实例分析验证了该方法的正确性。

这种验证的方法随后便应用到实际结构，来获得多输入随机参数的一个复杂系统的可靠性定量。

二、拉压干涉理论。

通常情况下，结构的最大应力和材料强度的分布有一定的概率密度。

图1显示的分布的一个例子，在fσ和gs分别代表最大应力和材料强度的情况下。

图1（a）中因为gs总是大于fσ，结构的可靠度是1。

然而，在图1（b），它是不能保证gs大于fσ。

图中分布着一个干扰区，如图所示，在阴影区域和放大的图2。

因此，结构的可靠度小于1。

基于数学推导的可靠性，可以按照分布做一体化计算。

RPσs−∞fσσgsdsdσ

fσgs

（a）

fσgs

（b）

图1拉压分布

gs

fσ

图2拉压干扰区

创建分析文件并定义参数

定义输入随机变量以及分配

定义输入数据之间的相关性

定义输出随机变量

定义概率模拟方法

响应面模拟

结论综述

图3拉力流程图

三、描述方法

商业有限元软件ANSYS有PDS（概率设计仿真）模块[2]。

该模块定义尺寸，材料特性和其他参数作为具有