HyperStudy软件-结构优化.ppt

A,1,HyperStudy优化,A,2,HyperStudy分析时所支持的求解器:,AbaqusANASYSLS-DYNANastranOptiStructPAM-CRASH2GRADIOSSHyperFormMADYMOADAMS,A,3,一通过hypermesh可以定义为设计变量的参数:,Shellthickness（壳单元料厚）Springstiffness（弹簧刚度）Concentratedmass（集中质量）Compositeplythickness（复合结构的组成结构板厚）Compositeplyangle（复合结构中各组成结构间的角度）Shapes(ascreatedbytheHyperMorphmodule,orbytheperturbationspanel)（通过HyperMorph或摄动面板下产生的结构的形状）,可为设计变量的参数名称,形状优化,应用范围,结构的灵敏度分析、优化,结构的灵敏度分析、优化,结构的灵敏度分析、优化,结构的灵敏度分析、优化,简化车身模型中，接头对车身的灵敏度及优化,设计变量的定义,A,4,1）Hyperform可用的设计变量：,二用于制造工艺优化的设计变量,Sheetthickness（钣金厚度）Friction（摩擦系数）Forces（载荷）Shapes（ascreatebyHyperMorph）（通过hypermorph定义的形状）,2）HyperXtrude可用的设计变量:,Controlpointsofthebearingprofilecurve（轴承轮廓曲线的控制点）Shapes(ascreatedbyHyperMorph)（通过hypermorph定义的形状）,设计变量的定义,A,5,三各求解器所支持的设计变量及对应的参数类型,设计变量的定义,A,6,三各求解器所支持的设计变量及对应的参数类型,设计变量的定义,A,7,HyperStudy的输入文件,HyperStudy认可的输入文件为:,Studyfiles(.xml)Modelfiles(.tpl，.hm，.hf，.mdl，.xls）Preferencefiles（preference.mvw）求解器脚本文件,A,8,针对不同求解器输出结果所支持的文件类型,Solver,Result,File,Remarks,Abaqus,Modelmass（质量）,centerofgravity（重心）momentsofinertia（转动惯量）Frequency（模态频率）,.dat,如果要在该文件中输出mass信息，则需要在计算输入文件中定义一个动态载荷步。,Nodalnumber（节点数）,Displacement（位移）,reactionforce（反作用力）,.fil,Elementnumber（单元数）,Stress（应力）,Strain（应变）,.fil,Nodalposition（节点位置）,.inp,Nodaldata（节点数据），elementdata（单元数据）,.odb,如果要得到向量中的具体元素的数值，需要使用resvector()函数。,节点向量包含x、y、z三个方向的数据。,HyperStudy的输出文件,内容：模型定义、历史定义，警告与出错信息提示,为计算输入文件，主要内容为模型定义,二进制文件，内容为模型信息与计算结果,计算结果文件，存放计算所得的单元、节点等结果。,*STEP*DYNAMIC1.0,0.01*ELFILE*ENDSTEP,2.如果要得到具体编号的单元或节点的数据，需要使用getval2函数。,1.所有单元与节点的数据以向量的形式存放。,A,9,Solver,Result,File,Remarks,OptiStructRadioss（Bulkformat）,Modelvolume（模型体积），Mass（质量），Frequency（模态频率），Bucklingfactor（屈曲因子）,.out,Nodaldisplacement（节点位移），reactionforce（反作用力），elementstress,strain（单元应力、应变），strainenergydensity（应变能密度）,.res,.h3d,针对不同求解器输出结果所支持的文件类型,HyperStudy的输出文件,Nodaldisplacement（节点位移），reactionforce（反作用力），elementstress,strain（单元应力、应变），strainenergydensity（应变能密度）,压缩的二进制结果文件，包含模型与结果信息,ASCII格式的结果文件，包含分析的注释,如果要得到向量中的具体元素的数值，需要使用resvector()函数,如果要得到向量中的具体元素的数值，需要使用resvector()函数,二进制结果文件,A,10,Solver,Result,File,Remarks,RADIOSS（blockformat）,Nodeldisplacement（节点位移），Velocity（速度），Acceleration（加速度）,.T01,elementstress、strain（单元应力、应变）,.T01,.A00,Nodaldata（节点数据），elementdata（单元数据）,Nastran,Model,mass（模型质量）,.f06,Nodeldata（节点数据），Elementdata（单元数据）Frequency（模态频率）,.op2,Towritethemasstothe.fo6file,usePARAM,GRDPNT,0intheinput.,Vectorcontainstimehistory,Vectorcontainstimehistory,针对不同求解器输出结果所支持的文件类型,HyperStudy的输出文件,如果要得到向量中的具体元素的数值，需要使用resvector()函数,二进制结果文件,如果要得到向量中的具体元素的数值，需要使用resvector()函数,输出数据的文件,结果存放文件,A,11,Solver,Result,File,Remarks,LS-DYNA,Modelandpartmass（模型、部件质量），Momentsofinertia（转动惯量），Centerofgravity（重心）,d3hsp,Modelenergies（模型能量）,glstat,nodout,Nodaldisplacement（节点位移），Velocity（速度），Acceleration（加速度）,Elementstress、strain（单元应力、应变）,elout,Nodaldata（节点数据），Elementdata（单元数据）,d3plot,Nodalposition（节点位置）,Keywordfile,文件存放历史数据。,Vectorcontainsx,yorzcomponentofthenode.,reactionforce（反作用力）,spcforc,针对不同求解器输出结果所支持的文件类型,HyperStudy的输出文件,如果要得到向量中的具体元素的数值，需要使用resvector()函数,该文件存放历史数据。,ASCII文件，存放全局数据，模型大部分类型的能量通过该文件输出。,ASCII文件，存放节点数据。,ASCII文件，存放单元（梁、壳）数据。,文件存放历史数据。,SPC反作用力。该文件存放力、力矩。,A,12,RADIOSS(BlockFormat),$file-both,RADIOSS(BulkDataFormat),Optistruct,$file/$filescrC:temp,Solverinputarguments（运行参数）,Solver（求解器）,Ls-DYNA,i=$file(MEMORY=5000000),Abaqus,job=.inp（memory=200Mb）interactive,MotionSolve,Excel,在Excel表格中输入设计变量并计算响应。,Nastran,$file,Excel,无,Inputfile（输入文件类型）,.fem文件,.fem文件,.bdf文件,/.dat文件,.key文件,.inp文件,.csv文件,运行控制,A,13,SizeOptimization尺寸优化,尺寸优化时，和尺寸相关的属性都可以作为尺寸优化的变量。常用的尺寸优化设计变量有：1）结构料厚thickness2）集中质量mass,尺寸优化的变量,A,14,一在excel表格基础上进行优化,1.目的是创建一个近似值，然后对这个近似值进行优化。2.该表格包含五列数据，如右图所示。第一列表示设计的次数，第二和第三列表示每个设计的两个设计变量，第四和第五列表示先前DOE分析所得到的结果。,1.打开excel表格，并将其另存为.txt格式的文本文件（制表符分隔）。,2.创建model时，modeltype选择HyperStudy。,3.创建response响应时，建立名称为index的响应用来表示设计的次数，其Responseexperssion可表示为“convert(getenv(“STUDY_RUN_NUMBER”)-1”,4.创建结果response响应时，响应矢量Vector的Vectorresourecefile选择referencefile类型的文件，并且选择前面得到的.txt格式的文本文件。Type和Request采用默认的Unknown和Block1,Component则根据该矢量的作用选择表格中相应数据表示的列(columnx),5.创建DOE分析时，Controlledfactors中的DOEClass选择RunMatrix，并相应选择先前得到的.txt文件。,A设计变量与结果在excel表格中已经给出，以此作为优化的依据。,由于Excel中已经列出了设计变量与相应的设计结果，因此前面建模时不需要进行初始计算，而仅在response中设置相应的响应结果即可。,SizeOptimization尺寸优化,A,15,B：在excel表格中输入设计变量，同时有响应的计算方法，以此作为优化的依据。,一在excel表格基础上进行优化,Exemple：HS-1070（Theobjectiveistofindthecross-sectionaldimensions,width,andheightofabeamthatminimizesthebeamvolumewhilekeepingthetipdeflectionbelow0.35mm.）如右图示,1.创建HyperStudy可以评估的矩阵输入。打开excel文件，选择tools菜单下的加载宏命令，将宏命令hw_hst_genpdd.xla加载入，如果没有，则在Altair的安装目录下templateshst找到hw_hst_genpdd.xla，并加载即可。,2.创建Model时，选择spreadsheet作为创建的Model的类型，并且前一步中加载入宏命令的excel表格作为Model的输入。,3.创建设计变量l时，通过AddModelParameter命令，在excel表格中选择设计变量(DesignVariable)和响应(Response)。,SizeOptimization尺寸优化,A,16,ShapeOptimization形状优化,形状优化的设计变量是结构的形状，因此在创建设计变量时需要先定义形状。通常使用HyperMorph定义用于形状优化的设计变量形状。,形状优化的变量,A,17,建立分析模型(Studysetup),步骤：,1定义分析,2建立模型参考输入文件类型,3定义设计变量参考设计变量定义,4初次分析运行参考求解器定义、运行控制,5建立响应参考HyperStudy输出文件,6设计变量间的相互关系定义,7灵敏度定义,A,18,设计变量的交互作用(ControlledDesignVariableInteraction),DOE分析(Postprocessing),左图中为设计变量C0与C1的交互作用，其中x轴表示其中一个设计变量C0，y轴表示分析的响应。图中的两条曲线则分别表示设计变量C1在其极限值时（最大与最小设计值）响应随设计变量C1变化而产生的变化。如果两条曲线相互平行，说明两个设计变量之间没有交互作用，也即两变量改变对响应的影响是相互独立的。如果两曲线不平行，则说明该二变量之间有交互作用。,A,19,响应拟合(Approximation),误差方差(Residuals),表示通过求解器计算得到的响应值与使用回归方程预测得到的响应值之间的差异，该差异值越小，所得到的回归方程越接近实际。,根据DOE分析结果建立响应关于设计变量的关系曲线，在设计变量发生更改时，根据拟合得到的关系曲线即可得到相应的响应结果，而不需要对整个模型进行分析来获得结果，这样可以节约计算时间，从而提高成本。,A,20,方程拟合(Approximation),1定义一个拟合方程,步骤：,曲线方程拟合的类型：1）LeastSquaresRegression,如图中所示：定义曲线拟合方程时，程序会在方程内自动为每个设计变量分配一个变量名，用于拟合曲线。,2）MovingLeastSquares,3）HyperKriging,A,21,方程拟合(Approximation),2输入矩阵,通过DOE分析计算，得到一系列不同设计变量时，相应的响应结果，以此作为输入，建立拟合方程。,A,22,方程拟合(Approximation),3创建拟合方程,4残差,拟合模型,拟合模型的阶次,拟合模型的精度,拟合方程及系数,DOE分析得到的响应结果与通过拟合方程得到的相应响应结果之间的差别对比。,A,23,方程拟合(Approximation),5trade-off,6Anova,通过拟合方程得到的响应结果与设计变量之间的关系图表。,评估误差及各因子对结果的影响作用。,A,24,优化分析(Optimizationstudy),典型的目标函数类型(TypicalObjectiveFunctions),1结构质量与体积(StructuralMassandVolume)2位移(Displacements)3频率(Prequency)4能量(Energies)forDYNA/ABAQUS5伤害指标(HeadInjuryCriteriaHIC,FemurLoads,etc),典型的约束函数类型(TypicalConstraints),1位移(Displacements)2应力(stress)3塑性应变(PlasticStrain)4频率(Prequency),典型的设计变量类型(TypicalConstraints),1尺寸优化(结构厚度，梁截面参数)2形状优化(修改结构的外部边界)3材料响应(不同密度的泡沫材料、钢材),A,25,优化分析(Optimizationstudy),1定义优化,优化算法类型：1）AdaptiveResponseSurfaceMethod自适应响应面法2）SequentialQuadraticProgramming序列二次规划法3）MethodofFeasibleDirections模式搜索法/可行方向法4）GeneticAlgorithm遗传算法5）UserDefined-Xopt用户自定义优化方法,4）遗传算法基于生物进化上“适者生存”的原则。每次迭代优化时，保留适应度高的因子，淘汰适应度低的因子，进入下次优化，经过数次迭代得到最好的结果，即为问题的最优解或次优解。,步骤：,A,26,优化分析(Optimizationstudy),2选择用于优化的设计变量,3定义约束,A,27,优化分析(Optimizationstudy),程序重启动后，是否接着指定的优化或DOE分析步继续运行。,优化类型选择,优化目标：,1）Minimize；结果最小化2）Maximize；结果最大化3）MinMax；最大值最小化4）MaxMin；最小值最大化5）SystemIdentification；系统识别,最大迭代次数,结果绝对收敛条件,结果相对收敛条件,约束容差,设计变量收敛条件,优化初次移动步长,步长界限,初始变量波动,约束的显示,出现失败分析的处理,创建优化目标,选择既非约束又非目标的响应的评估方式。,一般定义优化的目标时，采用默认的优化参数即可。如果有特殊需求，则可根据需要修改相应的优化参数。,4定义目标,A,28,优化分析(Optimizationstudy),5结果处理,a优化曲线,b优化数据表,所有响应（包括用于优化的约束和目标）、设计变量在优化过程中的数值变化曲线。,所有响应（包括用于优化的约束和目标）、设计变量在优化过程中的数值列表。,失效的优化方案,可接受的优化方案,接近成功的优化方案：约束超差在容差范围内。,A,29,随机性研究(Stochasticstudy),A,30,随机性研究(Stochasticstudy),1.定义一个随机分析,1)SimpleRandomSampling简单的随机抽样；2)LatinHypercubeSampling；3)HammersleySampling。,数据采样方法：,步骤,A,31,随机性研究(Stochasticstudy),2.定义随机变量,随机变量确定时需注意：,1）随机变量的离散与连续性（与设计变量的离散与连续性相同，在studysetup时确定）；2）随机变量的分布：aNormalDistribution正态分布bUniformDistribution平均分布cTriangularDistribution三角分布dExponentialDistribution指数分布eWeibullDistribution威布尔分布,随机变量分布类型选择,根据统计规律，正态分布是自然界多数现象的存在状态，因此确定随机变量时一般采用正态分布形式。,A,32,随机性研究(Stochasticstudy),3.定义随机变量的相关性,4.选择响应,5.计算抽样数据,使用二维、三维以及统计方法列出随机变量的抽样数据分布。其中Histogram柱状图，表示随机变量取值频率的分布；pdf概率分布图，表示随机变取值与均值接近程度的概论；cdf累积分布图，表示随机变量取值小于或等于均值的概论。,A,33,随机性研究(Stochasticstudy),6.后处理,对抽样的随机变量进行计算得到的响应的分布,使用二维、三维以及统计方法列出响应与随机变量或响应之间的数据分布。左图为响应结果的统计学分布其中Histogram柱状图，表示响应结果值出现频率的分布；pdf概率分布图，表示响应结果值与目标值接近程度的概率；cdf累积分布图，表示响应结果值小于或等于目标值的概率。,A,34,随机性研究(Stochasticstudy),6.后处理,响应结果统计,A,35,随机性研究(Stochasticstudy),6.后处理,分析结果可靠性,Reliability：可以满足该响应约束条件的设计的概率；Probabilityof