patran错误日志及解决方法

1Introduction这篇文章我本人在学习patran/natran过程中遇到的问题，及后来找到的解决方法，这篇文章也在逐步更新中，希望这篇文档能给那些学习用patran/nastran的一点帮助。YuanchongxinDelft2011/10/3121.THESOLUTIONFORTHERESIDUALSTRUCTUREANDTHEAPPLIEDLOADSFORTHECURRENTSUBCASEAREZERO.后来将边界条件由１２３４５６改为１２３，即将位移约束变成simplifiedsupported，就没有此问题了。2.USERWARNINGMESSAGE4124(IFS3P)THESPCADDORMPCADDUNIONCONSISTSOFASINGLESET在图中用了RB3的MPC，其中dependentnode(ux,uy,uz),independent(ux,uy,uz,rx,ry,rz),有可能是这里的问题。不过这个倒不影响计算结果。3.建立夹层结构的有限元网络对于meshing,sweepcanproducethesolidelementonthebasisoftheshellelement.另外sweep下的loft，可以在两个shellmesh之间，创建solidelement.对于夹层结构的modeling,shellandsolidelementshouldsharecommonnode.否则算出来的结果，solidstress为零。至于如何共用节点，则需用到element下的sweep命令。4.NoPARAMvaluesweresetintheControlFile不管失败还是成功的f06文件中，都会出现这句话。这句话是正常的。如果出现结果为０，并且还没有fatal，另一种可能是因为mpc没有包括全部边界节点，在对夹层结构圆柱进行分析时，发现由于没有完全把边界节点进行displacement界定，以及没有完全进行mpc限制，f06文件中为０，后来改正一下，就好了。5.RBE3加载与使用可以用RBE３进行axialload的加载，但不能用于bendingmoment的加载，曾尝试过，发现f06文件中全部为零。后改成直接对边界曲线进行加载，问题解决。对比下将bendingmoment加载到曲线和加载到节点上的区别：曲线：bendingmoment为40,EIGENVALUE=2.187715E-02节点：bendingmoment为40,EIGENVALUE=2.187715E-02结果完全一样，这说明将载荷加到曲线或曲线上的节点上，得到的结果完全相同。Axialload加载到RBE３点上和加载到节点上的区别：RBE3：axialload=30,EIGENVALUE=-2.379083E+00节点：USERFATALMESSAGE4683(LNNRIGL)KDIFMATRIXNEEDEDFOREIGENVALUEANALYSIS由此可以看出，对于圆柱体来说，轴向载荷只能依靠rbe３来加载，而弯曲载荷则靠RBE2加载。出现此问题，还有可能没有载荷加载。另外，奇怪的事，当弯曲载荷施加于平板时，用linearstatic分析可以成功，但是用105计算时却出现4683的错误。后来将载荷改为轴向压力时，可使上述错误消失。这说明载荷方向不对，也会出现类似错误。那至于圆柱应该如何呢？3RBE3的加载，dependent的dof应该是123,indepdendent应该是123456.不然造成无解。不过有个参考书上正好相反，但用nastran计算有误。不知为何，这里需要提到的是计算结果与参考书相差为2，现在尚未找出原因。6.USERFATALMESSAGE9050(SEKRRS)^^^RUNTERMINATEDDUETOEXCESSIVEPIVOTRATIOSINMATRIXKLL.^^^USERACTION:CONSTRAINMECHANISMSWITHSPCIORSUPORTIENTRIESORSPECIFYPARAM,BAILOUT,-1TOCONTINUETHERUNWITHMECHANISMS.以前也遇到这种情况，这次遇到后，又在一节点加载了位移约束，就解决了。看来这种错误主要是由于约束不够，线性方程组无解造成的。还有可能是没有equivalence，Thisshouldsolveyourproblemorreducethenumberoffailedratios.有一次就遇到此情况。7.常用材料定义对比表MAT1MAT2MAT3MAT8MAT9isotropicanisotropic(2)orthotropic(3)orthotropic(2)anisotropic(3)•UseaFORCEentryifyouwanttodefineastatic,concentratedforceatagridpointby一个点的力specifyingavector.•UseaFORCE1entryifthedirectionisdeterminedbyavectorconnectingtwogridpoints.两个点的力•UseaFORCE2entryifthedirectionisspecifiedbythecrossproductoftwosuchvectors.以上的乘积SPC1,42,123456,1,2,3,4,6,7,8,942为SID。SID:loadsetidentificationnumber载荷序号G:Gridpoint加载点CID：Coordinatesystemidentificationnumber坐标序号F：Scalefactor放大因子Ni：三维方向的微量，由CID决定例子：FORCE3802000.1.0.0.加载在节点8上的力，大小为2000N.向量为100，分析坐标为coord0.SID=3。Force1和force2的例子在《nastranusergood》P142.48.MPCThedisplacementconstraincanaffecttheoutputresultsharply.Whenanalysisthebucklingofthecylinder,Isetthetopendas12456,whichmeansthatitcanmoveatzdirection.Whilethebottomissetas123456,whichmeansfixed.Theresultisequaltothereference,butthedeformationdoesn’tresemblemuch.Idon’tknowwhy.1.ThempcissetasRBE3Dependent123456Independent123IftheRBE3ischangedasDependent123Independent123456TheobtainedEigenvaluedoesn’taltermuch.Theformeris0.337,thelatteris0.332.2.Theloadcanaddononeside,andtheothersideisfixed.Itisnotgoodtosetcompressiveloadatbothend,inthiscase,thepatranispronetofail.3.Undercompressiveloadforcylinder,Thedeformationisnotverysymmetry,andsomeparthaveabnormallargedeformation,Idon’tknowwhynow.9.关于几个应力的计算方法vonMisesStress：Goodquestion.VonMisesStressisactuallyamisnomer.ItreferstoatheorycalledtheVonMises-Henckycriterionforductilefailure.Inanelasticbodythatissubjecttoasystemofloadsin3dimensions,acomplex3dimensionalsystemofstressesisdeveloped(asyoumightimagine).Thatis,atanypointwithinthebodytherearestressesactingindifferentdirections,andthedirectionandmagnitudeofstresseschangesfrompointtopoint.TheVonMisescriterionisaformulaforcalculatingwhetherthestresscombinationatagivenpointwillcausefailure.TherearethreePrincipalStressesthatcanbecalculatedatanypoint,actinginthex,y,andzdirections.(Thex,y,andzdirectionsaretheprincipalaxesforthepointandtheirorientationchangesfrompointtopoint,butthatisatechnicalissue.)5VonMisesfoundthat,eventhoughnoneoftheprincipalstressesexceedstheyieldstressofthematerial,itispossibleforyieldingtoresultfromthecombinationofstresses.TheVonMisescriteriaisaformulaforcombiningthese3stressesintoanequivalentstress,whichisthencomparedtotheyieldstressofthematerial.(Theyieldstressisaknownpropertyofthematerial,andisusuallyconsideredtobethefailurestress.)TheequivalentstressisoftencalledtheVonMisesStressasashorthanddescription.Itisnotreallyastress,butanumberthatisusedasanindex.IftheVonMisesStressexceedstheyieldstress,thenthematerialisconsideredtobeatthefailurecondition.Theformulaisactuallyprettysimple,ifyouwanttoknowit:(S1-S2)^2+(S2-S3)^2+(S3-S1)^2=2Se^2WhereS1,S2andS3aretheprincipalstressesandSeistheequivalentstress,orVonMisesStress.Findingtheprincipalstressesatanypointinthebodyisthetrickypart.OctahedralShearStress：HydrostaticStress：1st,2nd,and3rdinvariantstresses：PrincipalStresses：Themagnitudesofthetwoprincipalstressesfromthe2DMohr’scirclemethodarecalculatedaccordingthefollowingequations:TrescaShearStress：MaximumShearStress：6Magnitude：10.加载位移和载荷要加载在节点