等几何分析

等几何分析等几何分析等几何分析等几何分析研究进展摘要等几何分析就是一种刚刚兴起的数值分析方法,对现有的CAE产生了很大的影响。等几何分析法的出现于发展,缓解与消除了困扰CAE多年的难题,开启了一条结合 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 、分析与优化等三方面的途径。本文阐述了等几何分析产生的背景、意义与相关的定义,还介绍了等几何分析从首先提出到现如今的10年发展历程,包括基础理论体系的发展与完善,新型样条的构建,网格细分方法的研究,计算效率的提高,以及其她方面(如边界条件的施加、接触分析、结构优化等)的进展,展示了等几何分析相对于基于拉格朗日插值的有限元法的优势。关键字等几何分析有限元NURBS发展现状1前言有限元分析就是目前应用最广泛的一种数值分析方法,且由于结合了能够高速运算的计算机,有限元法得到了大多数人的支持。有限元法就是将连续的物体离散成有限个单元,单元之间通过节点连接在一起,并将节点处的未知量作为基本未知量,使得无限自由度问题转换成了有限自由度的问题,在利用力学原理近似的求解出未知量。这一突出优点使得有限元法得到广泛应用,各类有限元软件也层出不穷,如ABAQUS、ANSYS、LS-DYNA、HyperMesh等。不过这一突出的优点也大大的限制了有限元的进一步发展。首先,有限元法求得的结果的精确度与网格的细化程度有关,网格越细,则计算结果的精度越高,而计算时间与计算所需的内存也将随之增加,而以目前的水平来瞧,还无法做到超高精度的细化网格。Sandia国家实验室曾做过一项统计,在汽车、航空航天与造船行业,大约全部分析时间的80%用于网格划分及划分前的几何模型准备[1]。其次,网格划分使得应力不连续,且在处理大变形问题中,单元的过度扭曲导致精度严重损失。第三,网格划分工具对几何形状的识别精度较低,特别就是划分复杂高级曲面时无法精确划分,容易划分出大变形网格。再者,网格划分就是建立在几何模型的基础上,若几何模型发生改变,那么须得重新划分网格,花费大量时间。最后,在处理网格畸变、网格移动如动态裂纹扩展、冲压成型等问题时需要进行网格重构,不仅浪费计算时间,还会损害计算精度[2]。网格就是有限元分析的基础,而以上缺陷都就是网格划分造成的,就是有限元法无法避免的。基于以上原因,在2005年,Hughes等[3]提出等几何分析的思想。该方法直接结合了CAD中的几何模型,将其中的几何信息作为有限元分析的输入信息,大大地节省划分网格的时间。等几何分析与有限元法有许多相同之处,可以说就是有限元法的发展,但其具有一套独立的理论体系。该方法采用描述几何形状的NURBS函数作为基函数,具有几何精确特性,且离散的几何形状不随单元的稀疏而改变,这意味着即使就是比较稀疏的网格划分,也能精确描述研究对象的几何形状,具有很高的数值精度[4]。NURBS本身就具有网格,一个NURBS实体包含若干个NURBS单元,分析时,这些单元成为精确描述几何形状的实体单元。另外,类似于有限元的网格,NURBS单元也可以细分,基函数的次数也可提高,计算结果更加精确,但几何形状不改变。于就是,Hughes将其命名为等几何分析。2等几何分析简介2、1B样条基函数由于NURBS基函数就是B样条基函数的线性组合[5],这里先讨论B样条基函数的构造。B样条基函数由节点矢量构建,如,式中ui为节点,n与p分别就是B样条基函数的个数与阶数。基函数由Cox-deBoor递推公式定义为[6]:当p=0时,当p>0时,由基函数Ni,p(u)与控制点Pi可表示出B样条曲线:由于B样条曲线具有局部性质,因而,可将上面的p次B样条曲线方程改写为分段表示形式:若给定(m+1)×(n+1)个控制点Pi,j(i=1,2,…,m,j=1,2,…,n)的阵列,构成一张控制网格。又分别给定参数u与v的阶数p与q,以及两个节点矢量与V=,这样,就定义了一张p×q次张量积B样条曲面,方程为:由于B样条曲面也具有局部性质,可将其分段表示为:2、2NURBS曲线与曲面2、2、1NURBS曲线p次NURBS曲线可以表示为一分段有理多项式矢函数:式中,被称作权或全因子,分别与控制点相联系。2、2、2NURBS曲面类似于B样条曲面,NURBS曲面也可分段表示为:同样的,一个三变量NURBS实体也可表示为:2、3计算流程基于NURBS的等几何分析法的分析思路如下[3]:由节点向量积确定NURBS片;通过节点插值将计算域细分为单元;每个基函数的支撑域包含少量单元;由基函数的控制点定义几何模型;采用等参概念,即场变量与几何模型采用相同的基函数表示,而基函数的系数即为自由度或控制变量;通过节点插值或基函数阶数可进一步细化单元,有h型细化、p型细化与k型细化;采用类似于有限元的方法,可将等参NURBS片构建的数组组装成全局数组;施加Dirichlet边界条件有几种方法。最粗糙的方法就是加在控制变量上,这种近似法会导致比较大的误差。然而,对于一些特殊情况,如齐次边界条件,该方法能满足精确度要求。此外,Dirichlet边界条件常常通过变分近似法或几何近似法施加。3等几何分析的发展3、1基础理论体系Hughes等提出了等几何分析的概念后,Bazilevs等[7]用数学的知识对其进行分析与误差估计, 证明 住所证明下载场所使用证明下载诊断证明下载住所证明下载爱问住所证明下载爱问 了等几何分析的收敛性与稳定性等特征,这个结论为之后等几何分析的发展奠定了扎实的理论基础。Cottrell与Hughes[8]研究了等几何分析中网格的细化与近似连续性。Gomez等[9]通过等几何分析研究了Cahn-Hilliard相域建模问题。Lipton等[10]研究了等几何离散化的鲁棒性。Sevilla等[11]提出了增强的NURBS有限元方法。Shaw与Roy[12]创立了基于NURBS的参数无网格法。Hughes等[13]人发现对B样条使用高斯积分计算效率并不理想,因而基于Half-Point-rule提出了一种宏单元积分法,该方法比高斯积分成本减少一半、3、2网格细化等几何分析中主要有3种细化方法[14]:基于节点插入的h型细化方法,基于基函数升阶的p型细化方法以及升阶与节点插入相结合的k型细化方法。之后,为了提高计算精度,Xu等[15]提出了可优化计算域内部控制点的位置的r型细化方法。徐岗等[16]在r型细化方法的基础上提出了r-p型细化方法,并应用于二维热传导问题上。没多久,徐岗、朱亚光等[17]提出了基于局部误差估计的自适应r细化方法,大大提高计算域参数化的优化效率。Pilgerstorfer等[18]进一步从理论层面挖掘了计算域参数化以及节点分布对等几何分析求解的影响,发现对一般问题而言,求解误差与计算域的等参数线网格(或等参数面网格)的均匀性及正交性有关,这也为构造适合分析的计算域参数化提供了重要的理论依据。3、3新型样条最初的等几何分析采用的基函数就是CAD中最广泛应用的B样条与NURBS基函数,后来许多学者陆续采用T样条等新兴样条[19]。T样条就是由Sederberg等[20]在2003年提出的,可以实现局部加密,使单个样条表示复杂模型的能力大大提升。Bazilevs等[7]率先将T样条用于等几何分析,并用于简单地二维与三维结构问题。Wang等[21]提出一种将任意非结构化二次网格转化为标准T样条表面的方法。不过,T样条的局部细分结果依赖于控制网的拓扑结构,其复杂性也就是不确定的,另外T样条的混合函数不一定总就是线性无关的[22]。于就是,Nguyen-Thanh等[23]提出了PHT样条(多项式分层T样条),PHT样条就是定义在层次T网格上双三次C1连续的样条空间,其继承了NURBS与T样条的所有优点,解决了NURBS与T样条相互转化的困难的问题。PHT样条具有很好的局部加细性质,这使其在拟合、缝合、简化、自适应曲面重构与等几何分析中得到了广泛的应用[22]。BernardMourrai等[24]研究了一般T网格上的维数计算公式,为将PHT样条向一般情形推广奠定了基础。3、4计算效率由于等几何分析法采用了高阶基函数,使得刚度矩阵相比于有限元法要稠密的多,大大地增加了计算时间,为此,许多学者展开了相关的研究。郭利财等[25]提出了一种基于矩阵分解的并行算法,从刚度矩阵的装配上节省计算时间。而后她们[26]还提出了直接使用粗网格的精确解作为细网格的初始解,大大提高了迭代速度。Veiga等[27]提出一种基于重叠加法型施瓦茨预条件算子的重叠计算域分解方法,证明等几OAS算子关于子区域数就是可扩展的。Gahalant等[28]分析了网间迁移算子的逼近性质与松弛算子的光顺性质,以此为基础证明了等几何分析多重网格方法的单元尺寸无关的最优收敛速度。Collier等[29]研究了基函数个数n,基函数次数p与光滑阶k对串行多波直接求解器性能的影响。3、5边界条件的施加由于基函数不同于有限元中的多项式基函数,它一般不具有插值的特征,若直接在控制点处施加边界条件而不釆取任何处理方式就会由于收敛速率的恶化导致结果出现明显误差。于就是,王东东等[30]通过一种变换将控制变量与边界上的配置节点联系起来,并提出一种对边界条件增强处理的方法。之后,她们提出了基于罚函数位移边界条件施加方式的几何精确NURBS有限元方法[31]。张勇等[32]提出了比例边界等几何分析方法,并将之应用于波导本征问题的求解,求解结果比传统等几何分析精度更高。陈涛等[33]提出一种基于样条拟合的Dirichlet边界条件施加方法,可通过最小二乘拟合引入更多的边界配点以获得精度更高的近似值。后来,她们又提出采用Nitshe法施加位移边界条件的新方法[34]。3、6结构优化用传统有限元法进行优化时,网格内部几何近似会在响应分析时导致精度问题,在设计灵敏度分析时会更不利。为此,Hughes等[3]提出将等几何分析用于结构优化设计。之后,Wall等[35]给出了一种基于等几何分析的结构形状优化设计框架,并将其应用到二维线弹性结构的最优化问题。Ha等[36]在等几何结构优化中采用了T样条,在产生相同优化设计形状的情况下,可获得更高效的求解过程。Nagy等[37]将等几何分析用于曲梁的尺寸与形状优化,提出了多级设计方法。Seo等[38]采用样条表面与修剪曲线表示设计模型的外部与内部边界,提出了新的控制点灵敏度分析公式与内前面合并算法。Qian与Sigmund提[39]出了一种将等几何形状优化的类矩形NURBS面片扩展到布局复杂的几何的优化设计方法。Hassani等[40]使用了基于控制点的SIMP法与类似于SIMP法的惩罚技术,推导应用了优化准则。3、7接触分析有限元法采用非光滑边界逼近光滑边界,人为导致接触边界的不光滑性,在某些情况下导致较大的误差。另外,由于单元边界处法线方向不唯一而导致搜索困难,搜索成本高。引入等几何分析将大大减小以上原因带来的计算误差。Lu[41]首先提出了接触问题的等几何分析框架,在几何描述与分析中统一采用NURBS基函数。Temizer等[42]提出了等几何分析的点-面接触算法,定性研究了各种有限变形无摩擦热弹接触问题,得到了令人满意的结果。DeLorenzis[43]将等几何分析用于二维大变形Coulomb摩擦接触问题。Stadler等[44]在有限元框架下采用NURBS接触面,得到Cn连续的有限变形弹性接触模型,并用于分析二维问题。DeLorenzis等[45]将等几何分析用于三维无摩擦大变形光滑接触问题与有摩擦接触问题。Kim[46]提出一种新的接触匹配算法用于精确确定等几何分析框架下的接触面。4展望等几何分析目前应用在简单外形的问题较多,而对于外形复杂的几何模型,构造适合分析的的体参数化以及可能出现的样条不连续问题目前仍就是一个不易解决的难题。目前的研究主要集中在基于节点插人的局部加细样条上,如T样条、PHT样条等。若能找到基于升阶的局部加细样条,对以后的研究工作有很大的帮助。等几何分析的高阶性使得其计算效率低,如何提高计算效率,使得等几何分析能真正应用于实际中就是今后的研究重点之一。等几何分析采用样条基函数作为形函数,由于基函数不满足正交条件,其网格存在重叠或错误,引起求解错误,该方面有待研究。5结论10年的发展,使等几何分析得等到逐步完善,诸多问题(如网格细化、边界条件的施加、计算效率等)得到了解决与提高。越来越多的学者将其运用到各个领域中(如板壳问题、结构优化、流固耦合等)。等几何分析方法的提出不仅为CAD/CAE的无缝融合提供新思路,还为成熟的几何设计与计算领域提出了新的课题。虽仅有10年的发展,但等几何分析仍就是一种创新优秀的数值计算方法。参考文献[01]BazilevsY,CaloVM,CottrellJA,etal、IsogeometricanalysisusingT-splines[J]、ComputerMethodsinAppliedMechanicsandEngineering,2010,199(5-8):229-263[02]葛建立,杨国来、同几何分析研究进展[J]、力学进展,2012,46(2):771-784[03]HughesTJR,CottrellJA,BazilevsY、Isogeometricanalysis:CAD,finiteelem-ents,NURBS,exactgeometryandmeshrefinement[J]、ComputerMethodsinAppliedMechanicsandEngineering,2005,194(39-41):4135-4195[04]张汉杰、薄梁板结构NURBS几何精确有限元分析[J]、力学季刊,2010,31(4):469-477[05]RogersDF、AnIntroductiontoNURBSwithHistoricalPerspective[M]、Acade-micPress,2001[06]DeBoorC、OnCalculationwithB-spine[J]、J、appox、theory,1972(6):50~62[07]BazilevsY,BeiraodaVeigaL,CottrellJA,etal、Isogeometricanalysis:approxi-mation,stabilityanderrorestimatesforh-refinedmeshes[J]、MathematicalMod-elsandMethodsinAppliedScience,2006,16(7):1031-1090[08]CottrellJA,HughesTJR,RealiA、Studiesofrefinementandcontinuityinisog-eometricstructuralanalysis[J]、ComputerMethodsinAppliedMechanicsandEn-gineering,2007,196:4160-4183[09]GomezH,CaloVM,BazilevsY,HughesTJR、Isogeometricshellanalysis:theReissner-Mindlinshell[J]、ComputerMethodsinAppliedMechanicsandEnginee-ring,2010,199:276-289[10]LiptonS,EvansJA,BazilevsY,ElguedjT,HughesTJR、Robustnessofisogeo-metricstructuraldiscretizationsundersevermeshdistortion[J]、ComputerMeth-odsinAppliedMechanicsandEngineering,2010,199:357-373[11]SevillaR,Fernandez-MendezS,HuertaA、NURBS-enhancedfiniteelementmet-hod[J]、Int、J、Numer、Eng、,200876:56-83[12]ShawA,RoyD、NURBS-basedparametricmesh-freemethod[J]、ComputerMet-hodsinAppliedMechanicsandEngineering,2008,197:1541-1567[13]HughesTJR,RealiA,SangalliG、EfficientquadratureforNURBS-basedisoge-ometricanalysis[J]、ComputerMethodsinAppliedMechanicsandEngineering,2010,199:301-313[14]BeiraodaVeigaL,BuffaA,RivasJ,etal、Someestimatesforh-p-krefinementinisogeometricanalysis[J]、NumerischeMathatik,2011,118(2):271-305[15]XuG,MourrainB,DuvigneauR,etal、Parameterizationofcomputationaldoma-ininisogeometricanalysis:methodsandcomparision[J],ComputerMethodsinAppliedMechanicsandEngineering,2011,200(23/24):2021-2031[16]徐岗,王毅刚,胡维华、等几何分析中的r_p型细化方法[J]、计算机辅助设计与图形学学报、2011,23(12):2019-2024[17]徐岗,朱亚光,邓立山等、局部误差驱动的等几何分析计算域自适应优化方法[J]、计算机辅助设计与图形学学报、201426(10):1622-163[18]PilgerstorferE,JuttlerB、Boundingtheinfluenceofdomainparameterizationandknotspacingonnumericalstabilityinisogeometricanalysis[J]、ComputerMeth-odsinAppliedMechanicsandEngineering,2014,268:589-613[19]李伟伟、基于B样条与NURBS的等几何分析研究[D]、吉林:吉林大学、2013[20]TSederberg,JZheng,ABakenov,etal、T-splinesandT-NURCCs[C],ACMTransGraph、2003[21]WangW,ZhangY,ScottMA,etal、ConvertinganunstructuredquadrilateralmeshtoastandardT-splinesurface[J]、ComputationalMechanics、2011,48(4):477-498[22]康红梅、适合于分析与建模的若干样条的研究[D]、合肥:中国科技大学,2014[23]Nguyen-ThanhN,Nguyen-XuanH,BordasSPA,etal、Isogeometricanalysisus-ingpolynomialsplinesoverhierarchicalT-meshesfortwo-dimensionalelasticsolids[J]、ComputerMethodsinAppliedMechanicsandEngineering,2011,200(21-22):1892-1908[24]BernardMourrain、OnthedimensionofsplinespacesonplanarT-subdivisions[M]、INRIATechnicalreports,2010[25]GuoLC,HuangZJ,GuNJ、AparallelNURBS-basedisogeometricanalysisfra-meworkonmulti-coreCPUs[J]、JournalofComputationInformationSystem,2012,8(16):6549-6600[26]CollierN,PardoD,DalcinL,etal、Thecostofcontinuity:astudyoftheperform-anceofisogeometricfiniteelementsusingdirectsolvers[J]、ComputerMethodsinAppliedMechanicsandEngineering,2012,213:353-361[27]daVeigaLB,ChoD,PavarinoL,etal、OverlappingSchwarzmethodsforisogeo-metricanalysis[J],SIAMJournalonNumericalAnalysis,2012,50(3):1394-1416[28]郭利财、通过并行计算与多重网格提升等几何分析计算效率[D]、合肥:中国科技大学,2013[29]GahalautK,KrausJ,TomarS、Multigridmethodsforisogeometricdiscretization[J]、ComputerMethodsinAppliedMechanicsandEngineering,2013,253:413-425[30]WangDD,XuanJC、AnimprovedNURBS-basedisogeometricanalysiswithen-hancedtreatmentofessentialboundaryconditions[J]、ComputerMethodsinAp-pliedMechanicsandEngineering,2010,199(37-40):2425-2436[31]王东东,轩军厂,张灿辉、几何精确NURBS有限元中边界条件施加方式对精度影响的三维计算分析[J]、计算力学学报,2012,29(1):31-37[32]张勇、比例边界等几何分析方法I:波导本征问题[J]、力学学报,2012,44(2):382-392[33]陈涛,莫蓉,万能、等几何分析中Dirichlet边界条件的配点施加方法[J]、机械工程学报,2012,48(5):157-164[34]陈涛,莫蓉,万能、等几何分析中采用Nitsche法施加位移边界条件[J]、力学学报,2012,44(2):369-381[35]WallWA,FrenzelMA,CyronC、Isogeometricstructuralshapeoptimization[J]、ComputerMethodsinAppliedMechanicsandEngineering,2008,197(33):2976-2988[36]HaSH,ChoiKK,ChoS、NumericalmethodforshapeoptimizationusingT-spl-inebasedisogeometricmethod[J]、StructuralandMultidisciplinaryOptimization,2010,42(3):417-428[37]NagyAP,AbdallaMM,GadurdalZ、Isogeometricsizingandshapeoptimisationofbeamstructures[J]、ComputerMethodsinAppliedMechanicsandEngineering,2010,199(17-20):1216-1230[38]SeoYD,KimHJ,YounSK、Isogeometrictopologyoptimizationusingtrimmedsplinesurfaces[J]、ComputerMethodsinAppliedMechanicsandEngineering,2010,199(49-52):3270-3296[39]QianX,SigmundO、IsogeometricshapeoptimizationofphotoniccrystalsviaCo-onspatches[J]、ComputerMethodsinAppliedMechanicsandEngineering,2011,200(25-28):2237-2255[40]HassaniB,KhanzadiM,TavakkoliSM、Anisogeometricalapproachtostructuraltopologyoptimizationbyoptimalitycriteria[J]、StructuralandMultidisciplinaryOptimization,2012,45(2):223-233[41]LuJ、Isogeometriccontactanalysis:Geometricbasisandformulationforfriction-lesscontact[J]、ComputerMethodsinAppliedMechanicsandEngineering,2011,200(5-8):726-741[42]TemizerI,WriggersP,HughesTJR、ContacttreatmentinisogeometricanalysiswithNURBS[J]、ComputerMethodsinAppliedMechanicsandEngineering,2011200(9-12):1100-1112[43]DeLorenzisL,TemizerI,WriggersP,etal、AlargedeformationfrictionalcontactformulationusingNURBS-basedisogeometricanalysis[J]、InternationalJournalforNumericalMethodsinEngineering,2011,87(13):1278-1300[44]StadlerM,HolzapfelGA,KorelcJ、CncontinuousmodellingofsmoothcontactsurfacesusingNURBSandapplicationto2Dproblems[J]、InternationalJournalforNumericalMethodsinEngineering,2003,57(15):2177-2203[45]DeLorenzisL,WriggersP,ZavariseG、Amortarformulationfor3Dlargedefor-mationcontactusingNURBS-basedisogeometricanalysisandtheaugmentedLagrangianmethod[J]、ComputationalMechanics,2012,49(1):1-20[46]KimJY,YounSK、Isogeometriccontactanalysisusingmortarmethod[J]、Inter-nationalJournalforNumericalMethodsinEngineering,2012,89(12):1559-1581