MSC.NASTRAN--MSC.NASTRAN--结构构优化化 第一部分结构优化简介 结构优化分类结构优化分类 •设计优化•形状优化•动响应优化•超单元优化•空气弹性变形优化 设计优化设计优化 概念概念•什么是设计优化 在计算机上实现自动修改分析模型参数以达到预期目标并满足设计要求•基本优化问题的数学描述*寻找一组设计变量 x1,x2,...*使得函数F(X)最小*并且满足 不等式约束、等式约束、副边界条件 设计模型设计模型•SOL 200——支持设计灵敏度和优化•分析类型：静力分析、模态分析、bucking、直接频率*、模态频率*、模态瞬态*、静线弹性、线弹性颤动（*包含声振）•设计变量：尺寸特性（包括超单元）、形状（超单元只有节点可变） Bulk Data: DESVAR•设计变量与特性的关系 Bulk Data: DVPREL1, DVPREL2, DEQATN•定义目标函数 执行控制：DESOBJ Bulk Data: DRESP1, DRESP2, DEQATN•定义约束条件 Bulk Data: DCONSTR, DCONADD Case Control: DESSUB, DESGLB Version69新功能概括新功能概括•Beam截面库•单元及特性输入检查•格式化的灵敏度值输出•Mode tracking•多边界条件：bucking, normal modes, flutter•新的自然频率响应类型 NASTRAN结构优化的优点结构优化的优点•有效解决小到大规模问题•可靠的收敛特性•用户界面、用户定义方程的多样性•近似概念的完善•逐步加强•NASTRAN分析的可靠性•NASTRAN的一部分——经济水平的支持•方便利用NASTRAN中熟悉的分析工具 结构优化的难题结构优化的难题•大量的设计变量•大规模的约束条件•指导设计的结构响应值与设计变量是隐式函数。
函数式难以确定，需要梯度信息以确定寻查方向•对设计工程师的要求： *设计目标函数的表达，尽可能用最少的设计变量； *对约束条件的适当描述； *分析模型的完整性 结构优化的解决方法近似技术近似技术•执行细节有限元分析•计算所有的约束条件，删除不很关键的条件(DESREEN)•计算剩余约束条件的梯度•产生与设计变量有关的高可靠性近似响应值•解决近似问题•修正分析变量•执行设计的具体分析•计算所有的约束条件•检查优化的收敛性•必要时重复以上过程 第二部分结构优化设计模型 设计模型与分析模型关系 设计模型与分析模型比较•分析模型中的单元特性是设计模型中设计变量的函数例如：在工字梁截面参数优化中， 分析模型中的截面面积、惯性积等是设计变量宽、高和厚度的函数 设计模型如何定义？•设计变量定义•目标函数定义•特性或几何形状与设计变量的关系•判别设计响应•定义设计约束条件•优化过程控制参数设置 设计变量定义DESVAR 说明：说明： 字段 内容•ID 设计变量识别号•LABEL 用户定义输出名（字符）•XINIT 变量初始值•XLB 下限•XUB 上限•DELXV 目标函数定义DESOBJ格式格式•执行控制卡: DESOBJ(NIN)=N 定义目标函数，必须结合以下命令的使用：•Bulk Data: DRESP1 N ... or DRESP2 N ... •定义目标函数的例子 使重量最小： DESOBJ(MIN) = 10 DPESP1, 10, w, weight 设计变量与特性关系的定义DVPREL1 说明：说明： 字段 内容•ID 识别号•TYPE 特性卡类型•PID 特性卡识别号•FID 分析模型中单元特性在特性卡中的字段位置•PMIN 特性的最小值•PMAX 特性的最大值•C0 关系式的常数(缺省值为0)•DVIDi 设计变量卡DESVAR的识别号•COEFi 线性关系式的系数 判别设计响应DRESP1 说明：说明： 字段 内容•ID 识别号•LABEL 用户定义输出名•RTYPE 响应类型(WEIGHT,EIGN,STRESS等)•PTYPE 单元特性名(PBAR,PSHELL等)•REGION 用以筛选约束条件•ATTA ATTB 响应属性 ATTi 定义设计约束条件 DCONSTR结合以下命令的使用以定义约束条件：•执行控制卡 DESSUB 或 DESGLB•Bulk Data DRESP1 或 DRESP2 DCONSTR •定义设计约束条件的例子 定义单元特性组1和2类型的单元应力值范围 SUBCASE 20 ANALYSIS = STATICS DESSUB = 100 … DRESP1,11,SAMAX,STRESS,PBAR, ,7, ,1, + +, 2 DRESP1,12,SBMAX,STRESS,PBAR, ,14, ,1, + +, 2 DCONSTR, 100, 11, -85., 85. DCONSTR, 100, 12, -85., 85. 优化过程控制参数设置 DOPTPRM •DOPTPRM卡中常用参数 DESMAX——设计容许的最大迭代次数 P1——输出控制 P2——输出控制（缺省为1，输出目标和设计变量值） DELP——两次迭代单元特性容许差值 （缺省0.02） DELX——两次迭代设计变量容许差值 （缺省1） CONV2——收敛准则（缺省值为0.01，两次迭代差值小于CONV2时，优化迭代结束） IPRINT——结果输出控制（缺省值为0，不输出；2输出迭代过程结果及寻查方向；5输出迭代过程目标函数及设计变量的值） METHOD——数字优化方法（1改进的可行方向法；2序列线性规划；3序列二次规划 ） 定义不独立设计变量与独立设计变量的关系DLINK 说明说明 DLINK定义了以下关系式:•字段 内容•ID 识别号•DDVID 特性卡类型•C0 关系式的常数(缺省值为0)•CMULT 乘子•IDVi 独立设计变量DESVAR的识别号•Ci 线性关系式中IDVi的系数 •用DLINK定义设计变量关系的例子: 设计变量t2,t3可写作独立变量t1,t4的函数,DLINK可实现此关系式 … DESVAR, 1, t1, 1., 0.01, 5. DESVAR, 2, t2, 1., 0.01, 5. DESVAR, 3, t3, 1., 0.01, 5. DESVAR, 4, t4, 1., 0.01, 5. DLINK, 11 , 2, 0., 0.333, 1, 2., 4, 1. DLINK, 12 , 3, 0., 0.333, 1, 1., 4, 2. DVPREL1, 21, PSHELL, 101, 4, 0.01, 5., , , + +, 1, 1.0 DVPREL1, 22, PSHELL, 102, 4, 0.01, 5., , , + +, 2, 1.0 DVPREL1, 23, PSHELL, 103, 4, 0.01, 5., , , + +, 3, 1.0 DVPREL1, 24, PSHELL, 104, 4, 0.01, 5., , , + +, 4, 1.0 ... 定义设计方程式DEQATN•定义一个或多个方程式,用以设计灵敏度分析或P-单元分析•格式: DEQATN EQID EQUATION 定义方程中的常数项 DTABLE•定义DEQATN中的常数项•格式 DTABLE LABL1 VALU1 LABL2 VALU2 -etc.- 设计变量与特性关系的定义DVPREL2 •用DVPREL2定义变量与特性关系的例子: 形心C1,C2可写作独立变量B,H的线性关系式,可用DVPREL1实现惯性积I1,I2与独立变量B,H的函数关系,结合DEQATN,用DVPREL2实现… DESVAR, 10, B, 0.3, 0.1, 1. DESVAR, 11, H, 0.4, 0.1, 1. DVPREL2, 250, PBAR, 120, 4, , , 501, , + +, DEVAR, 10, 11 DVPREL2, 251, PBAR, 120, 5, , , 502, , + +, DEVAR, 10, 11 DVPREL2, 252, PBAR, 120, 6, , , 503, , + +, DEVAR, 10, 11 DEQATN 501 AERA(B,H) = B*H DEQATN 502 I1(B,H) = B*H**3/12. DEQATN 503 I1(B,H) = H*B**3/12. DVPREL1, 260, PBAR, 120, 12, -0.5, , , , + +, 10, -0.5 DVPREL1, 261, PBAR, 120, 13, 0.05, , , , + +, 11, 0.5... 判别设计响应DRESP2 说明：说明： 字段 内容•ID 识别号•LABEL 用户定义输出名•EQID DEQATN识别号•LABLi DTABLE卡中定义的常数符号•NRi DRESP1的识别号•Gm 节点号•Cm 节点自由度号 用用DRESP2定义设计响应的例子定义设计响应的例子 节点1,2,3,4的位移u1,u2,u3,u4可用DRESP1卡定义, 结合DEQATN和DRESP1, 自定义响应DEF可用DRESP2实现:… DESOBJ(MIN) = 21... DRESP1, 11, U1,DISP, , , 2, , 1 DRESP1, 12, U2,DISP, , , 2, , 2 DRESP1, 13, U3,DISP, , , 2, , 3 DRESP1, 14, U4,DISP, , , 2, , 4 DRESP2, 21, DEF, 100, , , , , , + +, DRESP1, 11, 12, 13, 14 DEQATN 100 F1(A,B,C,D) = (2.*A+D)/3. ; + F2 = (A+2.*D)/3. ; + F = SQRT( (B-F1) * * 2 + (C-F2)* *2 )... 结构优化的Bulk Data 命令DESVARDLINKDVPREL1或DVPREL2DRESP1或DRESP2DCONSTRDCONADDDOPTPRM 结构优化中情况控制命令•DESOBJ•DESSUB•DESGLB 第三部分灵敏度分析 设计灵敏度分析概念•设计灵敏度分析 是目标函数和约束条件对设计变量的偏导数以至二阶偏导数的计算,即结构响应量对设计变量的变化率. •设计灵敏度分析为近似技术提供梯度信息•位移灵敏度•特征值灵敏度•为节省内存,NASTRAN中不形成总刚度矩阵偏导数或总质量矩阵偏导数,而在总刚度矩阵偏导数迭加过程中与相应的位移项相乘,只形成虚载荷的列向量,该乘积在设计变量变化单元的相应自由度有变化. 设计灵敏度分析的实现计算代价昂贵•减少设计变量•减少约束条件数•减少方程右式的虚载荷向量数目 设计灵敏度系数矩阵 设计灵敏度输出参数DSAPRT说明:START=i ——从第i次迭代开始输出BY=j ——定义输出迭代的间隔为jEND=k ——输出到第k次迭代为止ALL——输出所有的DRESP1和DRESP2定义的响应n——只有在SET命令定义的响应 设计灵敏度分析和优化的参数•有限插分法选择控制CDIF•将设计变量和节点修正过程写入PUNCH文件的控制DESPCH（0输出最后一次迭代结果；小于0，不输出；大于0，输出每一次迭代结果）•DSNOKD（1Buckling灵敏度分析包括刚度微分；0不包括）•NASPRT（输出迭代过程控制）•OPTEXIT（用户预先定义退出点）•SOFTEXIT（Yes,收敛则停止迭代；缺省时No，不停止迭代）•UPDTBSH（形状优化参数：Yes,修正边界形状；缺省值No，不修正边界形状） 静态响应灵敏度•定义设计模型•定义静态分析类型 SOLUTION 200 ANALYSIS = STATICS•用OPTEXIT定义灵敏度系数输出 PARAM,OPTEXIT,4•调用设计优化分析,静态响应灵敏度在设计优化过程中自动执行. 特征值灵敏度•定义设计模型•定义分析类型 ANALYSIS = MODES•用OPTEXIT定义灵敏度系数输出 PARAM,OPTEXIT,4•调用设计优化分析,特征值灵敏度在设计优化过程中自动执行 Buckling Load Factor灵敏度•定义设计模型•定义屈曲分析类型 ANALYSIS = BULK•用OPTEXIT定义灵敏度系数输出 PARAM,OPTEXIT,4•调用设计优化分析, Buckling Load Factor灵敏度在设计优化过程中自动执行 第四部分•形状优化 形状优化方程描述 •SOL 200•分析类型（执行控制命令 ANALYSIS）•定义设计变量（ Bulk Data: DESVAR）•设计变量与分析模型中特性的关系 Bulk Data: DVPREL1, DVPREL2, DEQATN （形状）（形状）DVBSHAP,DVSHAP,BNDGRID•定义设计响应 Bulk Data: DRESP1, DRESP2, DEQATN•定义目标和约束条件 Bulk Data: DCONSTR, DCONADD 执行控制：DESOBJ，DESSUB, DESGLB•提供必要的控制参数 Bulk Data: DOPTPRM, DSCREEN形状优化的设计模型输入准备 形状优化的解决方法•建立辅助模型•收集位移解用作基本向量•定义设计任务 形状优化的辅助模型•目的:在以下方法中结合辅助模型的使用帮助产生形状基本向量 *直接输入形状法 *几何边界形状法 *分析边界形状法•什么是辅助模型? 辅助模型的变形用以产生形状基本向量,虽然和原模型几何条件一样,但边界条件、载荷条件、材料类型可能不同。
辅助模型•辅助边界模型CASE CONTROL命令 AUXCASE： 分解符,表示辅助边界模型在Case Control Section中的开始 AUXMODEL=n•BEGIN BULK AUXMODEL = n ：开始建立辅助边界模型 n 辅助模型举例AUXCASEAUXMODEL=1SUBCASE=10 SPC=100 LOAD=160 TEMP(LOAD)=162 DISPLACEMENT=ALLSUBCASE=20 SPC=100 LOAD=161 TEMP(LOAD)=162 DISPLACEMENT=ALLBEGIN BULK…….BEGIN BULK AUXMODEL=1……. 定义节点位置与设计变量的关系DVGRID 定义形状基本向量DVSHAP•将位移矩阵的列与设计变量识别号联系起来，定义形状基本向量•格式 DVSHAP,DVID,COL1,SF1,COL2,SF2,COL3,SF3•说明： DVID——DESVAR的识别号 COLi——位移矩阵的列数 Sfi——对应位移矩阵的列数的系数•注意 位移矩阵必须已经由NASTRAN产生，可由数据库中得到。
ASSIGN DISPMAT=‘file.master’ DBLOCATE DATABLK=(UG/UGD,GEOM1/GEOM1D, GEOM2/GEOM2D), LOGICAL=DISPMAT 直接输入形状法•预先进行辅助模型分析，DBLOCATE形成位移矩阵文件，然后用DVSHAP产生形状基本向量数据•举例：圆孔管道开口优化（设计变量为椭圆a,b） 辅助模型的几何、连接等与原型一样，只是加载情况不一样 几何边界形状法•不需预先的辅助模型，依据边界几何数据，自动产生形状基本向量•每次迭代结束，随形状变化，形状基本向量也在更新•用BNDGRID定义边界，DVGRID定义节点相对于变量的变化量 分析边界形状法•建立许多针对边界的边界辅助模型•分析边界形状法使用的卡AUXCASE——辅助模型执行控制命令分隔符AUXMODEL——形状优化产生边界形状的参考辅助模型（定义识别号）BNDGRID——形状优化中设计边界或面中节点的识别号DVBSHAP——设计变量与边界形状的关系，从辅助模型得到设计变量与边界形状基本向量的线性关系 指定设计边界节点BNDGRID•格式： BNDGRID,C,GP1,GP2,-etc.-•说明： C指自由度号，GPi 指形状边界节点的识别号•举例 BNDGRID,123,41,42,43,44,45,46 •借助辅助模型，将边界形状与设计变量联系起来•格式 DVBSHAP,DVID,AUXMOD,COL1,SF1,COL2,SF2,COL3,SF3•说明 DVID指设计变量识别号 AUXMOD指辅助模型识别号 COLi指辅助模型执行控制卡中载荷序列识别号 SFi比例定义边界形状与变量的关系DVBSHAP 结束结束 。