基于等几何分析方法的带孔结构形状优化设计

基于等几何分析方法的带孔结构形状优化设计Author:SnailAnalysisGeometry目录

目录引言在计算机辅助几何设计〔CAGD〕中，常用Non-UniformRationalB-Splines(NURBS)实现任意曲面造型。在国际标准化组织〔ISO〕1991年公布的工业产品几何定义STEP标准中，NURBS被定义为唯一的自由型曲线曲面的表示方法。计算机辅助几何设计1引言有限元分析1使用插值多项式分段逼近原模型模型转换上耗费大量时间和工作量单元细分时需要不断与原模型交互C0连续性，难以处理高阶连续问题引言等几何分析〔IGA-IsogeometricAnalysis〕方法用NURBS基函数作为形函数，代替有限元法的插值多项式，实现对精确几何模型的分析。NURBS基函数可以构造任意高阶连续的近似函数，克服了有限元分析方法通常仅有C0连续性的弊端。等几何分析1直接使用原模型的NURBS信息无需模型转换，不损失几何精度单元细分简便且能保持模型不变可方便用于求解高阶连续问题引言等几何分析方法在实现计算机辅助几何设计〔CAGD〕和有限元分析〔FEA〕结合的同时，自然在结构优化中拥有了独到的优势：—直接将模型的NURBS控制点及其权值作为优化对象—根据优化后的控制点信息简便精确地得到优化后的形状等几何分析1NURBSNURBSNURBS目录

目录简单结构的等几何分析方法B样条2B样条是NURBS的根底，p次B样条曲线的定义为：式中,是控制点，是B样条基函数，定义为：其定义在节点矢量上。简单结构的等几何分析方法B样条基函数的性质：—局部支撑性—非负性—标准性—可微性给定节点矢量和控制点信息〔坐标和权值〕即可确定一条B样条曲线以下图是节点矢量为的二次B样条曲线B样条2非均匀有理B样条〔NURBS〕在B样条方法的根底上引入了权因子与分母，弥补了B样条难以描述标准解析形状的缺陷，NURBS曲线定义为：NURBS曲面定义为：式中，双变量有理基函数为：

NURBS2简单结构的等几何分析方法在有限元分析方法中，常用等参单元来离散具有复杂形体的结构，防止在整体坐标系中直接构造插值函数，也即形函数。依据等参单元的思想，结构离散后的单元分为物理空间的单元和参数空间的单元，它们之间可通过雅克比矩阵进行坐标转换而联系起来。参数空间形状规那么且形函数统一。等参单元分析方法2八节点等参曲边单元的形函数节点坐标简单结构的等几何分析方法等几何分析方法在等参单元思想的根底上，用NURBS基函数代替等参单元的形函数，使模型建立和模型分析拥有了统一的“语言〞NURBS。基于NURBS的等几何分析方法2双二次NURBS曲面的有理基函数控制点坐标基函数随参数空间的改变而改变，也即不同单元的基函数不同。简单结构的等几何分析方法为了提高分析精度，等几何分析方法常在参数空间中插入节点以进行单元细分〔h改进方法〕，或者提高NURBS基函数的阶数〔p改进方法〕。几何模型在CAGD软件中建立，然后通过程序读取模型信息并直接进行单元细分和分析，中间无需对模型进行处理。基于NURBS的等几何分析方法2目录

目录带孔结构的等几何分析方法在分析边界复杂或者内部带孔的模型，很难采用单片NURBS来实现建模，即使建模成功，在分析过程中产生的物理空间单元常是扭曲不均匀的。多片NURBS模型的等几何分析方法3带孔结构的等几何分析方法采用多片NURBS来拼装成复杂的模型，多片NURBS的控制点、基函数和参数空间节点等由连通数组(ConnectivityArrays)联系起来。采用四片NURBS进行拼装，并进行单元细分。多片NURBS模型的等几何分析方法3控制点网格物理空间单元带孔结构的等几何分析方法在边界施加均布压力时，需要依据虚功等效原那么将受力单元边界上的压力等效移置到相关控制点上。均布压力的移植需要对此边界进行积分，压力常数可提取到积分符号外。NURBS基函数的阶数越高，每个控制点影响的单元边界越多，同理，与某单元边界相关的控制点也越多，载荷移置的工作量也就越大。多片NURBS模型的等几何分析方法3控制点上等效载荷受载边方向上单变量NURBS基函数受力点上载荷相关控制点个数带孔结构的等几何分析方法施加均布压力后，经过等几何分析得出的变形示意图〔放大50倍〕和VonMises应力云图。多片NURBS模型的等几何分析方法需要对分析模型人为地进行划分，这样CAGD文件中的NURBS信息就会被舍弃，违背了等几何分析方法将CAGD与FEA结合起来的初衷。而且还要消耗时间来构造多片NURBS模型，不利于程序的集成化和自动化多片NURBS模型的等几何分析方法3带孔结构的等几何分析方法3物理空间参数空间带孔结构的等几何分析方法由于其参数空间经过了一组曲线的修剪，有些单元〔被修剪单元〕已不再是矩形单元，所以关键问题在于区分被修剪单元并对其进行精确的数值积分。3带孔结构的等几何分析方法先进行全局单元细分，再区分被修剪单元并将之分解成三角形单元，然后筛选出有作用控制点〔依据基函数的局部支撑性〕。3带孔结构的等几何分析方法3目录

目录基于等几何方法的优化等几何分析方法的一个重要优点就是把CAGD、FEA和结构优化统一起来。在结构优化中，可以直接将控制点的坐标和权系数作为优化变量，根据优化后的控制点坐标和权系数，可以简便精确地得到优化后的几何形状，省去了一系列模型转换过程。由于在参数空间不改变的情况下〔即单元划分不变〕，物理空间的几何模型会随着控制点和权系数的变化而发生变化，这使得基于等几何分析的优化方法在优化过程中不用重新划分单元就可以得到优化后的形状，极大地节约了机时。优化手段4基于等几何方法的优化孔形优化设计4圆圈标记的控制点记为变化控制点〔variablecontrolpoint〕，其坐标为优化变量。矩形框住的控制点记为联合控制点〔conjunctcontrolpoint〕，其坐标改变规律与变化控制点对称。模型最外边界上的控制点记为固定控制点〔fixedcontrolpoint〕,其坐标不变，以保证模型边界固定。其余控制点记为连接控制点〔linkedcontrolpoint〕,按照与变化控制点和联合控制点的远近进行等比例调节，以保证物理空间单元均匀。基于等几何方法的优化以单元最大VonMises应力最小为优化目标，保证圆孔面积不变小，优化数学模型用下式表达为孔形优化设计4变化控制点的坐标向量带孔平板的面积带孔平板初始面积圆心到边界的距离目录

目录结论与展望结论5尚不能简单灵活地对复杂模型进行分析缺乏结论与展望展望5等几何分析还可将平面问题方便地转化为三维壳问题，区别只在于控制点是否位于同一平面内，而用于分析计算的参数空间依旧为二维平面。由于NURBS的光滑特性，等几何分析无需像有限元分析方法为保证单元外表的协调一致性而引入旋转自由度，极大简化了薄壳的分析过程。下一阶段可采用等几何分析方法对壳体中的孔进行形状优化。

Thankyouforyourattention!区分被修剪单元为了区分被修剪单元，需要先依据修剪曲线的方向法那么来判断参数空间中的任意一点是否在分析区域。规定将内部区域修剪掉的曲线为顺时针方向，将外部区域修剪掉的曲线为逆时针方向，曲线上任意一点的切向量方向与曲线方向一致。通过判断参数空间中的点所对应的叉积方向，就可以得到点与分析区域的关系。区分被修剪单元通过判断单元的中心点和四个顶点与分析区域的关系，就可以区分该单元是否为被修剪单元。区分过程分为两步：—判断单元中心点到修剪曲线的最近距离—找出修剪曲线上与单元顶点距离最近的四个点

被修剪单元上的数值积分需要将其分解为三角形单元，被修剪单元按照其顶点在分析区域的个数分为三种类型，每种类型对应着不同的分解形式。如被修剪单元同时被多条修剪曲线修剪，或者修剪曲线只对被修剪单元的某一边进行了修剪等，这时就要采用四叉树法先把单元细化成上述三种类型。

被修剪单元上的数值积分先计算出三角形在修剪曲线上的顶点坐标再按照三角形的高斯积分方法和修正的积分方法NEFEM分别对直边三角形单元和曲面三角形单元进行数值积分直边三角形单元的数值积分曲边三角形单元的数值积分

其余需要注意的问题对于几何外形比较复杂的分析模型，其物理空间的边界由参数空间的修剪曲线映射而成。NURBS曲面与修剪曲线的使用参数没有直接的解析关系。如果边界条件直接施加在修剪曲线上，需要做一定处理使之等效转换在被修剪NURBS曲面上。如果边界条件只是施加在修剪曲线上的一局部时那么要进行节点反求和插值细化。附录B：等几何分析方法一般过程

与等参元分析方法类比等几何分析方法的控制点对应于有限元中的节点；等几何分析的NURBS基函数对应于有限元中的形函数；等几何分析方法在物理空间中的每个区域对应于有限元的单元，它由参数空间中的区域映射而成；等几何分析过程可参照等参元的求解步骤；由于NURBS曲面中每个区域的基函数是不同的，所以等几何分析的求解过程要比等参元繁琐。附录B：等几何分析方法一般过程

等几何分析过程以第k个单元的计算为例参考文献[1]T.J.R.Hughes,J.A.Cottrell,Y.Bazilevs,Isogeometricanalysis:Cad,finiteelements,NURBS,exactgeometryandmeshrefinement,Comput.Methods.Appl.Mech.Engrg.194(39–41)(2005)413–4195.[2]J.A.Cottrell,T.J.R.Hughes,Y.Bazilevs,IsogeometricAnalysis:TowardIntegrationofCADandFEA,Wiley,ChiChester(2021).[3]S.Lipton,J.A.Evans,Y.Bazilevs,T.Elguedj,T.J.R.Hughes,Robustnessofisogeometricstructuraldiscretizationsunderseveremeshdistortion,Comput.MethodsAppl.Mech.Engrg.199(2021)357–373.[4]H.J.Kim,Y.D.Seo,S.K.Youn,IsogeometricanalysisfortrimmedCADsurfaces,Comput.MethodsAppl.Mech.Engrg.198(2021)2982–2995.[5]K.Reed,Jr.,D.Harrod,W.Conroy,TheInitialGraphicsExchangeSpecification(IGES)(Version5.0),USDepartmentofCommerce,1990.[6]W.A.Wall,M.A.Frenzel,C.Cyron,Isogeometricstructuralshapeoptimization,Comput.MethodsAppl.Mech.Engrg.197(2021)2976–2988.[7]A.P.Nagy,M.M.Abdalla,Z.Gurdal,Isogeometricsizingandshapeoptimizationofbeamstructures,Comput.MethodsAppl.Mech.Engrg.199(2021)1216–1230.