拓扑优化简介及在ansys软件中的实现

拓扑优化简介及在ANSYS软件中的实现,拓扑优化概述,ANSYS中拓扑优化过程,实例讲解,主要内容,拓扑优化概述,拓扑优化概述,工程结构优化,尺寸优化：以几何尺寸为设计变量,而材料的性质,结构的拓扑和几何形状保持不变,形状优化：以连续体几何区域的边界线或边界面为设计变量，拓扑保持不变。,拓扑优化：寻求结构的最优拓扑在计区域内寻求材料的最优分布问题。,拓扑描述方式,拓扑优化求解算法,均匀化方法密度法变厚度法拓扑函数描述法,有限元法优化准则法序列凸规划法无网格法,均匀化方法,优化准则法序列凸规划法,拓扑优化概述,可用来解决以下问题：（1）体积约束下的最大刚度设计：以柔顺度为目标函数，体积为约束函数；（2）刚度约束下的最小体积优化：以体积为目标函数，刚度为约束函数；（3）体积约束下的最大动刚度设计：以n阶自振频率为目标函数,体积为约束函数；（4）以上多种工况的组合优化问题；,依赖于单元的伪密度来决定材料去留，0(去掉)、1(保留),拓扑优化概述,ANSYS中拓扑优化过程,批处理的方式-APDL语言一,求解步骤,操作方式,GUI交互方式一般用户,建立几何模型直接在ANSYS中建立导入CAD三维建模软件中中性格式PRO/E软件与ANSYS软件无缝连接定义单元类型拓扑优化的有效单元类型为:PLANE2*或PLANE82*.SOLID92或SOLID95SHELL93定义实常数(壳厚度等),ANSYS中拓扑优化过程,定义材料属性要求输入杨氏模量及泊松比；注意泊松比不是缺省的0.3，需要定义；对重力、旋转或惯性载荷需要定义密度；记住使用统一的单位制；最方便的是使用材料库(MPREADwiththeLIBoption,orPreprocessorMaterialPropsMaterialLibrary).,ANSYS中拓扑优化过程,划分网格对清楚的拓扑结果建议采用细而均匀的网格。然而优化执行多次,因此网格过细运行时间也会增加；将不优化的单元类型设置为2或更大。,不参与优化的,优化的,TYPE2elements,TYPE1elements,只对类型号为1的单元进行优化，若想排除优化的区域，将该区域单元类型设为2或更大。,ANSYS中拓扑优化过程,施加载荷约束固定点，对称边界条件等等外载荷力、压力、温度及惯性载荷诸如重力及角速度,约束对结果有直接的影响,底边铰支UX=UY=0,底边简支UY=0,ANSYS中拓扑优化过程,施加尽可能少的约束；多种载荷，选择性施加:在一个载荷步中施加所有载荷结果形状对所有载荷共同作用提供最大刚度；产生非保守逼近，因为它假定所有载荷同时作用并为整个幅值。在独立的载荷步中分别施加各自的载荷(LSWRITEorSolutionWriteLSFile.)并指定多载荷步优化求解结果形状为对每一载荷的等刚度,ANSYS中拓扑优化过程,求解拓扑优化问题:指定优化控制(TOVAR及TODEF命令)体积减少量(作为百分比)载荷步数目收敛容差开始优化(TOLOOP命令)指定循环次数,体积减少量对结果影响显著,ANSYS中拓扑优化过程,从主菜单选择TopologicalOpt-SetUp-BasicOpt然后指定体积减少量；接下来，选择TopologicalOptRun输入迭代数，并开始优化。,ANSYS中拓扑优化过程,通用后处理器中查看伪密度等值图“topoplot”PLNSOL,TOPOorGeneralPostprocPlotResultsNodalSolution红色表示要保留的材料(pseudo-density1.0)；蓝色表示可以去掉的材料(pseudo-density0.0)。,ETAB,EDENS,TOPOPLETAB,TOPOPLETAB,EDENSPRETAB,EDENSESEL,S,ETAB,EDENS,0.9,1EPLOT,提取保留单元,ANSYS中拓扑优化过程,三个有用的图形操作（对所有结果图均适用）：设等值图线数为2:/CONTOUR,2或UtilityMenuPlotCtrlsStyleContoursUniformContours.关闭位移缩放:/DSCALE,OFF或UtilityMenuPlotCtrlsStyleDisplacementScaling.关闭图例栏:/PLOPTS,INFO,OFF或UtilityMenuPlotCtrlsWindowControlsWindowOptions.GUI操作对应命令流的输出单步查看最终整体输出,Helpisveryhelpful！,ANSYS中拓扑优化过程,对一长正方形平板零件，底边中部受到均匀的压力6.5MPa，顶部两侧受到集中载荷3.3KN。本问题的目标是在体积减少70%的条件下，结构的柔顺度最小。,实例一力载荷下的拓扑优化,3.3KN,6.5MPa,长宽=160120,弹性模量70MPa，泊松比为0.33,实例讲解,1定义分析类型PreferenceStructure2定义单元类型、平板厚度以及材料属性（1）定义单元类型与参数选项PreprocessorElementTypeAdd/Edit/DeleteAdd（2）定义实常数PreprocessorRealConstantsAdd/Edit/DeleteAdd（3）定义材料属性PreprocessorMaterialPropsMaterialModelsStructuralLinearElasticIsotropic3.利用尺寸变量建立模型PreprocessorModelingCreateAreasRectangleByDimensions4划分网格PreprocessorMeshingMeshTool,实例讲解,5施加载荷和边界条件SolutionDefineLoadsApplyStructuralDisplacementSymmetryB.C.OnNodesSolutionDefineLoadsApplyStructuralPressureOnNodes,节点集合的建立,实例讲解,6求解（1）指定体积减少量TopologicalOpt-SetUp-BasicOpt（2）开始优化TopologicalOptRun（输入迭代数）7查看求解结果GeneralPostProcPlotResultsContourPlotNodalSoluNodalSolutionTopologicalOptimizationDensitiesTopologicalOptPlotDensities,实例讲解,在实例一中的模型上施加热边界条件如下：,实例二力热载荷综合作用下的拓扑优化,3.3KN,6.5MPa,长宽=160120,实例讲解,热-结构耦合分析,采用顺序耦合分析的方法，即首先进行整机温度场分析，然后利用热分析结果即节点温度作为“体载荷”施加到随后的结构分析中。,耦合方法,温度场边界条件,转换单元,转换材料属性,转换接触算法,设置参考温度,清除物理环境,设置边界条件,分析流程,保存温度场物理环境,温度场计算,包括节点温度！,实例讲解,二维实体：PLANE55四节点四边形单元PLANE77八节点四边形单元PLANE35三节点三角形单元三维实体：Solid70一阶六面体单元（可退化为四面体）Solid87二阶四面体单元Solid90二阶六面体单元（可退化为四面体）,热分析单元,热分析单元转换为结构单元PreprocessorElementTypeSwitchElemType,二维实体：PLANE77PLANE183PLANE82三维实体：Solid87Solid187Solid92Solid90Solid186Solid95,推荐采用Solid87、Solid90，避免使用六面体单元,实例讲解,实例三ANSYSWORKBENCH中的拓扑优化,实例讲解,一个两端固