资源描述:
《Matlab Implementation of the Finite Element Method in Elasticity英文学习材料》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、Computing69,239–263(2002)DigitalObjectIdentifier(DOI)10.1007/s00607-002-1459-8MatlabImplementationoftheFiniteElementMethodinElasticityJ.Alberty,Kiel,C.Carstensen,Vienna,S.A.Funken,KielandR.Klose,KielReceivedJune18,2001;revisedFebruary25,2002Publishedonline:December5,2002ÓSpringer-Ve
2、rlag2002AbstractAshortMatlabimplementationforP1andQ1finiteelements(FE)isprovidedforthenumericalsolutionof2dand3dproblemsinlinearelasticitywithmixedboundaryconditions.Anyadaptationfromthesimplemodelexamplesprovidedtomorecomplexproblemscaneasilybeperformedwiththegivendocumentation.Num
3、ericalexampleswithpostprocessinganderrorestimationviaanaveragedstressfieldillustratethenewMatlabtoolanditsflexibility.AMSSubjectClassification:65N30,65R20,73C50.Keywords:Finiteelementmethod,elasticity,Matlab.1.IntroductionUnlikecomplexblack-boxcommercialcomputercodes,thispaperprovides
4、asimpleandshortopen-boxMatlabimplementationofP1andQ1finiteelements(FE)forthenumericalsolutionsoflinearelasticityproblems.Insteadofcoveringallkindsofpossibleproblemsinonecode,theproposedtoolaimstobeplain,easytounderstand,andtomodify.Therefore,onlysimplemodelexamplesareincludedtobeada
5、ptedtowhateverisneededinthespiritof[1].Wepresentanerrorestimatorwhichillustratestheaccuracyofthenumericalcalculation.AdifferentapproachisrealizedbythecompanyComsolwiththefiniteelementpackageFemlab[7].Itisatoolforthesolutionofpartialdifferentialequationsofmanyphysicalphenomena,includin
6、gstructuralmechanics.Aspecialfeatureisthecomfortablemeshgeneratorandthepossibilitytocombinedifferentphysicalproblems.ThecodeiswrittenforthemostpartinMatlabandcanbemodified,exceptforthemeshgeneratorandthesolverfornonlinearproblems.Thebasisfunctionsinthefiniteelementmethodcanbechosentob
7、elinear,quadratic,cubicoroffourthorder.Thesolvercanuseadaptivityaswell.Theplanofourpaperisasfollows.Themodelproblem,theNavier-Lame´equations,withgeneralboundaryconditionsisdescribedinSect.2;theweakformulationanddiscretisationinSect.3.Theheartofthiscontributionisthedatarepresentatio
8、nofthetriangulation,theDir