利用ansysAPDL进行优化设计的例子

ANSYS优化设计步骤解析本文介绍了ANSYS优化设计的相关步骤。

共有两种方法实现ANSYS优化设计：批处理方法和通过GUI交互式地完成。

这两种方法的选择取决于用户对于ANSYS程序的熟悉程度和是否习惯于图形交互方式。

如果对于ANSYS程序的命令相当熟悉，就可以选择用命令输入整个优化文件并通过批处理方式来进行优化。

对于复杂的需用大量机时的分析任务来说(如非线性)，这种方法更有效率。

而另一方面，交互方式具有更大的灵活性，而且可以实时看到循环过程的结果。

在用GUI方式进行优化时，首要的是要建立模型的分析文件，然后优化处理器所提供的功能都可以交互式的使用，以确定设计空间，便于后续优化处理的进行。

这些初期交互式的操作可以帮助用户缩小设计空间的大小，使优化过程得到更高的效率。

优化设计通常包括以下几个步骤，这些步骤根据用户所选用优化方法的不同(批处理GUI方式)而有细微的差别。

1. 生成循环所用的分析文件。

该文件必须包括整个分析的过程，而且必须满足以下条件：a. 参数化建立模型(PREP7)。

b. 求解(SOLUTION)。

c. 提取并指定状态变量和目标函数(POST1/POST26)。

2. 在ANSYS数据库里建立与分析文件中变量相对应的参数。

这一步是标准的做法，但不是必须的(BEGIN或OPT)。

3. 进入OPT，指定分析文件(OPT)。

4. 声明优化变量。

5. 选择优化工具或优化方法。

6. 指定优化循环控制方式。

7. 进行优化分析。

8. 查看设计序列结果(OPT)和后处理(POST1/POST26)。

优化设计步骤的细节在下面列出。

批处理方式和交互方式的区别也同时指出。

第一步：生成分析文件分析文件生成是ANSYS优化设计过程中的关键部分。

ANSYS程序运用分析文件构造循环文件，进行循环分析。

分析文件中可以包括ANSYS提供的任意分析类型(结构，热，电磁等，线性或非线性)。

(注：ANSYS/LS-DYNA的显式分析不能进行优化。

用ANSYS对液压支架运动仿真及优化设计高耀东;严鹏贺【摘要】四连杆机构的尺寸对液压支架的性能存在有至关重要的影响.用ANSYS 对液压支架进行了运动仿真,计算出各构件的位移、速度和加速度等运动特性,分析结果可用于评估液压支架特性和指导液压支架设计,并在此基础上可以进行优化设计.再以改善液压支架工作性能和减轻重量为优化目标,对液压支架进行多目标优化,确定了四连杆机构的最优结构参数.从实例结果看,用ANSYS对液压支架进行运动仿真和优化设计是可行的,结果是能够满足工程需要的.【期刊名称】《机械设计与制造》【年(卷),期】2018(000)008【总页数】4页(P202-205)【关键词】ANSYS;液压支架;优化设计;运动仿真【作者】高耀东;严鹏贺【作者单位】内蒙古科技大学机械工程学院,内蒙古包头 014010;内蒙古科技大学机械工程学院,内蒙古包头 014010【正文语种】中文【中图分类】TH16;TH1221 引言液压支架是综采工作面的关键设备，用于支撑顶板和形成工作空间，并推动支架和刮板输送机向工作面移动。

四连杆机构是支撑掩护式以及掩护式液压支架中最重要的联接装置，可以使支架承受较大的横向推力，从而使立柱只受轴向力作用。

由于对液压支架的受力状态、稳定性以及工作性能都存在至关重要的影响，所以确定四连杆机构的最优尺寸是液压支架设计必须优先考虑的问题[1-3]。

液压支架四连杆机构优化设计的优化目标有：使顶梁前端水平摆幅最小、使液压支架总重量最轻、使液压支架受力最合理等，而这些特性的评估可以通过运动仿真来实现。

ANSYS是目前应用最普遍的有限元软件，在多数的液压支架设计单位都有配备，熟练掌握该软件的工程技术人员有很多。

ANSYS可以依靠本身工具创建任意复杂的几何模型，也可以从其他CAD软件导入几何模型。

可以在ANSYS中进行结构静力学分析、多体动力学分析、结构优化设计等。

所以，用ANSYS分析液压支架的动力学特性、强度和刚度特性是十分方便、可行的。

计算机技术与应用文章编号:1009 6825(2008)03 0362 02基于ANSYS 参数化语言APDL 的结构优化设计收稿日期:2007 05 08作者简介:赵长勇(1978 ),男,长江大学城建学院硕士研究生,湖北荆州 434023张系斌(1956 ),男,硕士生导师,教授,长江大学城建学院,湖北荆州 434023翟晓鹏(1980 ),男,长江大学城建学院硕士研究生,湖北荆州 434023赵长勇 张系斌 翟晓鹏摘 要:指出了传统优化设计方法的不足,概述了A NSYS 和A PDL 的基本特征,综述了基于AN SYS 参数化设计语言A PDL 的结构优化设计的步骤及其应用现状,从而给建筑结构的优化设计注入新的活力。

关键词:结构优化设计,ANSYS,A PDL 中图分类号:T P 391文献标识码:A结构优化设计理论已有近四十年的发展历史,目前在一些重要的结构(如飞机结构)上已经得到了应用,这也引起了土木和建筑工程界人士的广泛关注,寻求建筑结构优化设计的理论、方法一直在紧张有序的进行当中。

由于传统的优化方法,例如准则法、数学规划法以及两者的结合(即所谓的混合法)等静态优化方法都是基于代数方程模型的;最优控制理论中的动态规划优化方法是基于微分方程或差分方程模型的。

而这些传统数学模型的描述能力和求解方法有相当的局限性,使得最优化理论和方法在实际应用中受到了很大的限制,存在着局部最优解、维数灾难、不确定性等问题,这些困难需要寻求新的优化设计方法,才能得到最终解决。

随着有限元理论的迅猛发展和日趋成熟,特别是计算机技术的广泛应用,基于A NSYS 参数化设计语言A PDL 的结构优化设计越来越体现出它强大的生命力,这无疑给建筑结构的优化设计注入了新的活力。

1 ANSYS 和APDL 概述及其关系A NSYS 是一种运用广泛的通用有限元分析软件,其有限元分析过程主要包括:建立分析模型并施加边界条件、求解计算和结果分析3个步骤。

ANSYS高级分析-优化设计（一）本篇作为ANSYS最常用也是最有用的高级分析技术优化设计技术的开篇，主要从概念上讲述ANSYS优化设计以及在进行ANSYS优化设计分析时通用的基本步骤和考虑方法（本篇所讲述优化分析同样是基于APDL工具，其它优化设计技术具体实现过程随着对优化设计的逐渐深入会逐步描述）。

1 前言优化设计是一种寻找确定最优设计方案的技术。

所谓“最优设计”，指的是一种方案可以满足所有的设计要求，而且所需的支出（如重量，面积，体积，应力，费用等）最小。

也就是说，最优设计方案就是一个最有效率的方案。

设计方案的任何方面都是可以优化的，比如说：尺寸（如厚度），形状（如过渡圆角的大小），支撑位置，制造费用，自然频率，材料特性等。

实际上，因为ANSYS具有专业的和强大的分析能力，具有开放的使用环境，所以所有可以参数化的ANSYS选项都可以作优化设计。

ANSYS参数化设计语言APDL用建立智能分析的手段为用户提供了自动循环的功能，也就是说，程序的输入可设定为根据指定的函数、变量以及选出的分析标准决定输入的形式。

APDL允许复杂的数据输入，使用户实际上对任何设计或分析有控制权。

例如尺寸、材料、载荷、约束位置和网格密度等。

APDL扩展了传统有限元分析之外的能力，并扩展了更高级运算，包括灵敏度研究、零件库参数化建模，设计修改和设计优化。

APDL是一种为ANSYS 二次开发专门设计开发的解释性文本语言，其内容包括参数、数组参数、表达式和函数、分支和循环、重复功能和缩写以及宏和用户程序等。

2 优化设计步骤基于APDL的参数化设计的方法采用ANSYS的批处理方法进行优化的，其主要的优化设计过程通常包括以下几个步骤，这些步骤根据用户所选用优化方法的不同（批处理GUI 方式）而有细微的差别。

生成分析文件：生成分析文件与通常的ANSYS分析并没有很大的不同，唯一的区别就是：优化的分析文件必须是参数化的。

也就是说，必须定义一些参量，至少要把所有的设计变量，状态变量和目标变量定义为参量形式。

ANSYS 优化设计1.认识ANSYS优化模块1.1 什么时候我需要它的帮忙?什么是ANSYS优化？我想说明一个例子要比我在这里对你絮叨半天容易理解的多。

注意过普通的水杯吗？底面圆圆的，上面加盖的哪一种。

仔细观察一下，你会发现比较老式的此类水杯有一个共同特点：底面直径＝水杯高度。

图1 水杯的简化模型为什么是这样呢？因为只有满足这个条件，才能在原料耗费最少的情况下使杯子的容积最大。

在材料一定的情况下，如果水杯的底面积大，其高度必然就要小；如果高度变大了，底面积又大不了，如何调和这两者之间的矛盾？其实这恰恰就反应了一个完整的优化过程。

在这里，一个水杯的材料是一定的，所要优化的变量就是杯子底面的半径r和杯子的高度h，在ANSYS的优化模块里面把这些需要优化的变量叫做设计变量（DV）；优化的目标是要使整个水杯的容积最大，这个目标在ANSYS的优化过程里叫目标函数（OBJ）；再者，对设计变量的优化有一定的限制条件，比如说整个杯子的材料不变，这些限制条件在ANSYS 的优化模块中用状态变量（SV）来控制。

下面我们就来看看ANSYS中怎么通过设定DV、SV、OBJ，利用优化模块求解以上问题。

首先参数化的建立一个分析文件（假设叫volu.inp），水杯初始半径为R＝1，高度为H ＝1（DV），由于水杯材料直接喝水杯的表面积有关系，这里假设水杯表面积不能大于100，这样就有S＝2πRH＋2πR2<100（SV），水杯的容积为V＝πR2H（OBJ）。

File:volu.inp (用参数直接定义也可或者在命令栏内直接写)R=1H=1S=2*3.14*R*H+2*3.14*R*RV=10000/(3.14*R*R*H)然后再建一个优化分析文件（假设叫optvolu.inp），设定优化变量，并求解。

/clear,nostart/input,volu,inp/optopanl,volu,inpopvar,R,dv,1,10,1e-2opvar,H,dv,1,10,1e-2opvar,S,sv,,100,1e-2opvar,V,obj,,,1e-2opkeep,onoptype,subpopsave,optvolu,opt0opexec最后，在命令输入框中键入“/input,optvolu,inp”，整个优化过程就开始了。

filename,contact_thermal/title,contact_thermal example/prep7et,1,13,4,,1! 选择直接耦合单元PLANE13，单元自由度为ux,uy,temp ! 定义为轴对称et,2,48! 定义结构接触单元keyopt,2,1,1! 设定接触单元的相应选项keyopt,2,2,1keyopt,2,7,1r,2,2e11,0,0.0001,,,0.1! 定义接触单元实常数mp,ex,1,2e11! 定义管材结构及热属性mp,alpx,1,5e-5mp,kxx,1,10mp,dens,1,7880mp,c,1,500rect,0.11,0.12,0,0.02! 建模rect,0.121,0.131,0,0.02 amesh,allnsel,s,loc,x,0.11! 将内管内壁的X方向位移及温度耦合cp,1,ux,allcp,2,temp,allnsel,s,loc,x,0.12! 将内管外壁的X方向位移及温度耦合cp,3,ux,allcp,4,temp,allnsel,s.loc,x,0.121! 将外管内壁的X方向位移及温度耦合cp,5,ux,allcp,6,temp,allnsel,s,loc,x,0.131! 将外管外壁的X方向位移及温度耦合cp,7,ux,allcp,8,temp,allnsel,s,loc,y,0.02! 将内管顶部节点的Y方向位移及温度耦合nsel,r,loc,x,0,0.12cp,9,uy,allnsel,s,loc,y,0.02! 将外管顶部节点的Y方向位移及温度耦合nsel,r,loc,x,0.121,0.131cp,10,uy,allnsel,s,loc,x,0.12! 创建接触单元cm,cont,nodensel,s,loc,x,0.121cm,targ,nodetype,2real,2gcgen,cont,targ,3/soluantype,trans! 瞬态分析tunif,20! 初始平均温度tref,20! 参考温度sfl,4,conv,300,,300! 内管内壁对流边界sfl,6,conv,10,,20! 外管外壁对流边界nsel,s,loc,y,0! 约束所有底边单元的Y向位移d,all,uy,0time,180! 载荷步时间deltime,10,5,15! 定义时间步长outres,all,all kbc,1autots,on! 自动时间步长allsel solve! 求解/post1plnsol,temp! 显示温度分布plnsol,s,eqv! 显示等效应力。

大桥全长2996.8m，其中主桥采用跨度为101.5+188.5+580+217.5+159.5+116m的钢桁梁斜拉桥；非通航孔正桥采用6孔跨径64m预应力混凝土简支箱梁；东引桥采用16孔梁长32.6m预应力混凝土简支箱梁；跨大堤桥采用48.9+86+48.8m预应力混凝土连续箱梁；西引桥采用15孔梁长32.6m预应力混凝土简支梁及2孔梁长24.6m预应力混凝土简支梁，其中宁安线采用箱梁，阜景线采用T梁。

主桥采用103+188.5+580+217.5+159.5+117.5m两塔钢桁斜拉桥方案，全长1366m。

主梁为三片主桁钢桁梁，桁间距2x14m，节间长14.5m，桁高15m。

主塔为钢筋混凝土结构，塔顶高程+204.00m，塔底高程-6.00m，斜拉索为空间三索面，立面上每塔两侧共18对索，全桥216根斜拉索。

所有桥墩上均设竖向和横向约束，4#塔与主梁之间设纵向水平约束，3#塔与梁间使用带限位功能的粘滞阻尼器。

主梁为”N”字型桁式，横向采用三片桁结构，主桁的横向中心距各为14m，桁高15m，节间距14.5m[2]。

结构构造主桥采用两塔钢桁斜拉桥方案，主梁为三片主桁钢桁梁，主桁上下弦杆均为箱型截面，上弦杆内高1000mm，内宽1200mm，板厚20～48mm。

下弦杆内高1400mm，宽1200mm，板厚20～56mm。

下弦杆顶板向桁内侧加宽700mm与整体桥面板焊接。

腹杆主要采用H型截面。

H型杆件宽1200mm，高720和760mm，板厚20～48mm。

根据不同的受力区段选用不同的杆件截面，在辅助墩附近的压重区梁段，腹杆采用箱型截面杆件。

主桁采用焊接杆件，整体节点。

在节点外以高强度螺栓拼接的结构形式，上下弦杆四面等强对接拼装。

H型腹杆采用插入式连接。

箱型腹杆采用四面与主桁节点对拼的连接形式。

主桁拼接采用M30高强螺栓。

表1.2 安庆长江大桥主要工程数量表表 2.1 设计构件节点和单元数表 2.2安庆长江铁路大桥模型主要参数列表表2.3 模型实常数/nerr,0!将参数读入数据库finish/clear/PREP7!单元类型编号et,1,beam44et,2,shell63et,3,link8!材料类型编号Q420/Q370/拉索-1，混凝土C40/C50-2，刚臂-3 mp,ex,1,2.1e11mp,prxy,1,0.3mp,dens,1,7850mp,alpx,1,1e-5mp,ex,2,3.5e7mp,prxy,2,0.16mp,dens,2,2500mp,alpx,2,1e-5mp,ex,3,1e15mp,prxy,3,1mp,dens,3,0mp,alpx,3,0!读取自定义截面并编号*dim,name,char,24name(1)='Sec_1a'name(2)='Sec_2a'name(3)='Sec_3a'name(4)='Sec_4a'name(5)='Sec_5a'name(6)='Sec_6a'name(7)='Sec_7a'name(8)='Sec_1'name(9)='Sec_2'name(10)='Sec_3'name(11)='Sec_4'name(12)='Sec_5'name(13)='Sec_6'name(14)='Sec_7'name(15)='Sec_8'name(16)='Sec_9'name(17)='Sec_10'name(18)='Sec_11'name(19)='Sec_12'name(20)='Sec_13'name(21)='Sec_14'name(22)='Sec_15'name(23)='Sec_16'name(24)='Sec_x2'*do,i,1,24,1sectype,i,beam,mesh,secread,name(i),,,mesh*enddo!将ta_sec参数读入数据库!读取自定义截面并编号*dim,name2,char,15name2(1)='ta_1'name2(2)='ta_2'name2(3)='ta_3'name2(4)='ta_4'name2(5)='ta_5'name2(6)='ta_6'name2(7)='ta_7'name2(8)='ta_8'name2(9)='ta_9'name2(10)='ta_10'name2(11)='ta_11'name2(12)='ta_12'name2(13)='ta_13'name2(14)='ta_14'*do,i,1,14,1sectype,40+i,beam,mesh,secread,name2(i),,,mesh*enddoSECTYPE, 26, BEAM, HREC, , 0 SECOFFSET, USER, 0.1, 0.266SECDATA,0.2,0.266,0.006,0.006,0.006,0.006,0,0,0,0 SECTYPE, 27, BEAM, I, , 0SECOFFSET, USER, 0, 0.584SECDATA,0.240,0.240,0.6,0.016,0.016,0.012,0,0,0,0 SECTYPE, 28, BEAM, T, T_nin, 0 SECOFFSET, USER, 0, 1.4SECDATA,0.6,1.4,0.032,0.018,0,0,0,0,0,0SECTYPE, 29, BEAM, T, T_max, 0SECOFFSET, USER, 0, 2.1SECDATA,0.6,2.1,0.032,0.016,0,0,0,0,0,0 SECTYPE, 31, BEAM, HREC, , 0SECOFFSET, CENTSECDATA,0.548,0.500,0.016,0.016,0.016,0.016,0,0,0,0 SECTYPE, 32, BEAM, I, , 0SECOFFSET, CENTSECDATA,0.560,0.560,0.500,0.028,0.028,0.016,0,0,0,0 !刚臂的截面SECTYPE, 100, BEAM, CSOLID, , 0 SECOFFSET, CENTSECDATA,0.06,0,0,0,0,0,0,0,0,0r,1,r,2,0.016,!拉索截面控制!7-349r,3,134.31e-4,!7-301r,4,115.84e-4,!7-283r,5,108.91e-4,!7-241r,6,100.000e-4,finish!开始建立模型/prep7!********************************************* !!!xiaXianGan!!!********************************************* *SET,j,-14*do,i,1,(94*25+1)*2,94*25+1n,i,0,j,0n,94*25+i,1363.000,j,0FILL,i,94*25+i,94*25-1,,1,,,!中桁*if,i,eq,94*25+2,then!7截面type,1mat,1secnum,14*do,a,1,3*25,1E,i-1+a,i+a*enddo*do,b,1,11*25,1E,35*25-1+i+b,35*25+i+b *enddo*do,c,1,8*25,1E,77*25-1+i+c,77*25+i+c *enddo!6截面type,1mat,1secnum,13*do,a,1,3*25,1E,i+3*25-1+a,i+3*25+a*enddo*do,b,1,3*25,1E,32*25-1+i+b,32*25+i+b *enddo*do,c,1,3*25,1E,46*25-1+i+c,46*25+i+c *enddo*do,d,1,1*25,1E,i+76*25-1+d,i+76*25+d *enddo*do,e,1,7*25,1E,87*25-1+i+e,87*25+i+e *enddo!5截面type,1mat,1secnum,12*do,a,1,25,1E,8*25-1+i+a,8*25+i+a*enddo*do,b,1,2*25,1E,30*25-1+i+b,30*25+i+b *enddo*do,c,1,2*25,1*enddo*do,d,1,2*25,1E,72*25-1+i+d,72*25+i+d *enddo!4截面type,1mat,1secnum,11*do,a,1,3*25,1E,9*25-1+i+a,9*25+i+a*enddo*do,b,1,4*25,1E,26*25-1+i+b,26*25+i+b *enddo*do,c,1,3*25,1E,51*25-1+i+c,51*25+i+c *enddo*do,d,1,4*25,1E,68*25-1+i+d,68*25+i+d *enddo!3截面type,1mat,1secnum,10*do,a,1,1*25,1E,7*25-1+i+a,7*25+i+a*enddo*do,b,1,5*25,1E,12*25-1+i+b,12*25+i+b *enddo*do,c,1,4*25,1E,22*25-1+i+c,22*25+i+c *enddo*do,d,1,4*25,1E,54*25-1+i+d,54*25+i+d *enddo*do,e,1,5*25,1E,63*25-1+i+e,63*25+i+e *enddo!2截面type,1mat,1secnum,9*enddo*do,b,1,2*25,1E,17*25-1+i+b,17*25+i+b*enddo*do,c,1,1*25,1E,21*25-1+i+c,21*25+i+c*enddo*do,d,1,1*25,1E,58*25-1+i+d,58*25+i+d*enddo*do,e,1,2*25,1E,61*25-1+i+e,61*25+i+e*enddo*do,f,1,2*25,1E,74*25-1+i+f,74*25+i+f*enddo*do,g,1,1*25,1E,85*25-1+i+g,85*25+i+g*enddo!1截面type,1mat,1secnum,8*do,a,1,2*25,1E,19*25-1+i+a,19*25+i+a*enddo*do,b,1,2*25,1E,59*25-1+i+b,59*25+i+b*enddo*do,c,1,1*25,1E,86*25-1+i+c,86*25+i+c*enddo!左边桁*elseif,i,eq,(94*25+1)*2+1,then !7截面type,1mat,1secnum,7*do,a,1,3*25,1E,i+a,i-1+a*enddo*do,b,1,11*25,1*do,c,1,8*25,1E,77*25+i+c,77*25-1+i+c *enddo!6截面type,1mat,1secnum,6*do,a,1,3*25,1E,i+3*25+a,i+3*25-1+a*enddo*do,b,1,3*25,1E,32*25+i+b,32*25-1+i+b *enddo*do,c,1,3*25,1E,46*25+i+c,46*25-1+i+c *enddo*do,d,1,1*25,1E,i+76*25+d,i+76*25-1+d *enddo*do,e,1,7*25,1E,87*25+i+e,87*25-1+i+e *enddo!5截面type,1mat,1secnum,5*do,a,1,25,1E,8*25+i+a,8*25-1+i+a*enddo*do,b,1,2*25,1E,30*25+i+b,30*25-1+i+b *enddo*do,c,1,2*25,1E,49*25+i+c,49*25-1+i+c *enddo*do,d,1,2*25,1E,72*25+i+d,72*25-1+i+d *enddo!4截面type,1mat,1secnum,4*do,a,1,3*25,1*do,b,1,4*25,1E,26*25+i+b,26*25-1+i+b *enddo*do,c,1,3*25,1E,51*25+i+c,51*25-1+i+c *enddo*do,d,1,4*25,1E,68*25+i+d,68*25-1+i+d *enddo!3截面type,1mat,1secnum,3*do,a,1,1*25,1E,7*25+i+a,7*25-1+i+a*enddo*do,b,1,5*25,1E,12*25+i+b,12*25-1+i+b *enddo*do,c,1,4*25,1E,22*25+i+c,22*25-1+i+c *enddo*do,d,1,4*25,1E,54*25+i+d,54*25-1+i+d *enddo*do,e,1,5*25,1E,63*25+i+e,63*25-1+i+e *enddo!2截面type,1mat,1secnum,2*do,a,1,1*25,1E,6*25+i+a,6*25-1+i+a*enddo*do,b,1,2*25,1E,17*25+i+b,17*25-1+i+b *enddo*do,c,1,1*25,1E,21*25+i+c,21*25-1+i+c *enddo*do,d,1,1*25,1E,58*25+i+d,58*25-1+i+dE,61*25+i+e,61*25-1+i+e *enddo*do,f,1,2*25,1E,74*25+i+f,74*25-1+i+f *enddo*do,g,1,1*25,1E,85*25+i+g,85*25-1+i+g *enddo!1截面type,1mat,1secnum,1*do,a,1,2*25,1E,19*25+i+a,19*25-1+i+a *enddo*do,b,1,2*25,1E,59*25+i+b,59*25-1+i+b *enddo*do,c,1,1*25,1E,86*25+i+c,86*25-1+i+c *enddo!右边桁*else!7截面type,1mat,1secnum,7*do,a,1,3*25,1E,i-1+a,i+a*enddo*do,b,1,11*25,1E,35*25-1+i+b,35*25+i+b *enddo*do,c,1,8*25,1E,77*25-1+i+c,77*25+i+c *enddo!6截面type,1mat,1secnum,6*do,a,1,3*25,1E,i+3*25-1+a,i+3*25+aE,32*25-1+i+b,32*25+i+b *enddo*do,c,1,3*25,1E,46*25-1+i+c,46*25+i+c *enddo*do,d,1,1*25,1E,i+76*25-1+d,i+76*25+d *enddo*do,e,1,7*25,1E,87*25-1+i+e,87*25+i+e *enddo!5截面type,1mat,1secnum,5*do,a,1,25,1E,8*25-1+i+a,8*25+i+a*enddo*do,b,1,2*25,1E,30*25-1+i+b,30*25+i+b *enddo*do,c,1,2*25,1E,49*25-1+i+c,49*25+i+c *enddo*do,d,1,2*25,1E,72*25-1+i+d,72*25+i+d *enddo!4截面type,1mat,1secnum,4*do,a,1,3*25,1E,9*25-1+i+a,9*25+i+a*enddo*do,b,1,4*25,1E,26*25-1+i+b,26*25+i+b *enddo*do,c,1,3*25,1E,51*25-1+i+c,51*25+i+c *enddo*do,d,1,4*25,1E,68*25-1+i+d,68*25+i+d *enddomat,1secnum,3*do,a,1,1*25,1E,7*25-1+i+a,7*25+i+a*enddo*do,b,1,5*25,1E,12*25-1+i+b,12*25+i+b *enddo*do,c,1,4*25,1E,22*25-1+i+c,22*25+i+c *enddo*do,d,1,4*25,1E,54*25-1+i+d,54*25+i+d *enddo*do,e,1,5*25,1E,63*25-1+i+e,63*25+i+e *enddo!2截面type,1mat,1secnum,2*do,a,1,1*25,1E,6*25-1+i+a,6*25+i+a*enddo*do,b,1,2*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,21*25-1+i+c,21*25+i+c *enddo*do,d,1,1*25,1E,58*25-1+i+d,58*25+i+d *enddo*do,e,1,2*25,1E,61*25-1+i+e,61*25+i+e *enddo*do,f,1,2*25,1E,74*25-1+i+f,74*25+i+f *enddo*do,g,1,1*25,1E,85*25-1+i+g,85*25+i+g *enddo!1截面secnum,1*do,a,1,2*25,1E,19*25-1+i+a,19*25+i+a*enddo*do,b,1,2*25,1E,59*25-1+i+b,59*25+i+b*enddo*do,c,1,1*25,1E,86*25-1+i+c,86*25+i+c*enddo*endif*SET,j,j+14*enddo!*********************************************!!!T_beam AND shell63's nodes!!!*********************************************/prep7ngen,2,7053,1,2351,1,0,0,0.7,ngen,4,2351,7054,7054+2350,1,0,0.57,0,ngen,4,2351,7054+3*2351,7054+3*2351+2350,1,0,0.60,0, ngen,3,2351,7054+6*2351,7054+6*2351+2350,1,0,0.45,0, ngen,2,2351,7054+8*2351,7054+8*2351+2350,1,0,0.60,0, ngen,3,2351,7054+9*2351,7054+9*2351+2350,1,0,0.45,0, ngen,2,2351,7054+11*2351,7054+11*2351+2350,1,0,0.70,0,local,11,0,0,-7,0nsym,y,2351*12,7054+2351,32914+2*2351,1csys,0nsym,y,58775,7054,65828,1ngen,2,122252,2352,2352+2350,1,0,0,0.7,!生成边上4条U肋type,1*do,h,6,9,1*do,k,1,2350,1e,2351*h+k,2351*h+1+k *enddo*enddo!生成T梁1type,1mat,1secnum,27*do,k,1,2350,1e,2351*10+k,2351*10+1+k *enddo!生成U肋2条type,1mat,1secnum,26*do,h,11,12,1*do,k,1,2350,1e,2351*h+k,2351*h+1+k *enddo*enddo!生成T梁2type,1mat,1secnum,27*do,k,1,2350,1e,2351*13+k,2351*13+1+k *enddo!生成U肋另外2条type,1mat,1secnum,26*do,h,14,15,1*do,k,1,2350,1e,2351*h+k,2351*h+1+k *enddo*enddo!镜像四分之一纵向的单元!（最好用*get 最大节点号单元号）esym,,2351*12,7051,30550,1!生成横梁type,1mat,1secnum,29*do,h,7054,9404,25*do,k,h,h+11*2351,2351e,k,k+2351*enddo*enddo*do,h,37617,39967,25*do,k,h,h+10*2351,2351e,k,k+2351*enddo*enddo*do,h,35266,37616,25e,h,h+2351*12*enddo!生成横肋type,1mat,1secnum,28*do,ah,7059,9384,25*do,ac,0,15,5*do,ak,ah+ac,ah+ac+11*2351,2351 e,ak,ak+2351*enddo*enddo*enddo*do,ah,37622,39947,25*do,ac,0,15,5*do,ak,ah+ac,ah+ac+10*2351,2351 e,ak,ak+2351*enddo*enddo*enddo*do,ac,0,15,5e,ah+ac,ah+ac+2351*12 *enddo*enddo!镜像一半单元csys,0esym,,58775,7051,65354,1!补充横梁type,1mat,1secnum,29*do,h,37617,39967,25 e,h,h+86987*enddo*do,h,96392,98742,25 e,h,h+28212*enddo!补充横肋type,1mat,1secnum,28*do,ah,37622,39947,25*do,ac,0,15,5e,ah+ac,ah+ac+86987*enddo*enddo*do,ah,96397,98722,25*do,ac,0,15,5e,ah+ac,ah+ac+28212*enddo*enddo!刚臂type,1mat,3secnum,100*do,h,1,2351,1e,h,h+7053*enddo*do,h,4703,7053,1*enddo*do,h,2352,4702,1e,h,h+124604-2352*enddo!********************************************* !!!shell!!!********************************************* type,2mat,1real,2*do,i,0,11,1*do,h,7054,7054+2349,1e,h+i*2351,h+1+i*2351,h+2352+i*2351,h+2351+i*2351 *enddo*enddo*do,h,35266,35266+2349,1e,h,h+1,h+63479-35266,h+63478-35266*enddo*do,i,0,-10,-1*do,h,63478,63478+2349,1e,h+i*2351,h+1+i*2351,h+i*2351-2350,h+i*2351-2351 *enddo*enddo*do,h,37617,37617+2349,1e,h,h+1,h+1+124604-37617,h+124604-37617*enddocsys,0esym,,58775,131654,190403,1*do,h,124604,124604+2349,1e,h,h+1,h+1-(124604-96392),h-(124604-96392)*enddo!!!ShangXianGan!!!********************************************* /prep7real,1*SET,Nmax1,300000*SET,j,-14*do,i,Nmax1+1,Nmax1+(92*25+1)*2,92*25+1n,i,14.50,j,15.00n,92*25+i,1348.50,j,15.00FILL,i,92*25+i,,,,,,!中桁*if,i,eq,Nmax1+92*25+2,then!13截面type,1mat,1secnum,20*do,a,1,2*25,1E,i-1+a,i+a*enddo*do,b,1,6*25,1E,3*25-1+i+b,3*25+i+b*enddo*do,c,1,3*25,1E,70*25-1+i+c,70*25+i+c*enddo*do,d,1,10*25,1E,74*25-1+i+d,74*25+i+d*enddo*do,e,1,3*25,1E,85*25-1+i+e,85*25+i+e*enddo*do,f,1,1*25,1E,91*25-1+i+f,91*25+i+f*enddo!12截面type,1mat,1secnum,19*do,a,1,1*25,1E,i+2*25-1+a,i+2*25+a*do,b,1,1*25,1E,9*25-1+i+b,9*25+i+b *enddo*do,c,1,8*25,1E,35*25-1+i+c,35*25+i+c *enddo*do,d,1,1*25,1E,i+69*25-1+d,i+69*25+d *enddo*do,e,1,1*25,1E,73*25-1+i+e,73*25+i+e *enddo*do,f,1,3*25,1E,88*25-1+i+f,88*25+i+f *enddo!11截面type,1mat,1secnum,18*do,a,1,25,1E,10*25-1+i+a,10*25+i+a *enddo*do,b,1,1*25,1E,18*25-1+i+b,18*25+i+b *enddo*do,c,1,1*25,1E,26*25-1+i+c,26*25+i+c *enddo*do,d,1,4*25,1E,31*25-1+i+d,31*25+i+d *enddo*do,e,1,1*25,1E,43*25-1+i+e,43*25+i+e *enddo*do,f,1,5*25,1E,45*25-1+i+f,45*25+i+f *enddo*do,g,1,1*25,1E,59*25-1+i+g,59*25+i+g *enddo*do,h,1,1*25,1E,68*25-1+i+h,68*25+i+h *enddo!10截面mat,1secnum,17*do,a,1,1*25,1E,11*25-1+i+a,11*25+i+a *enddo*do,b,1,1*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,20*25-1+i+c,20*25+i+c *enddo*do,d,1,1*25,1E,25*25-1+i+d,25*25+i+d *enddo*do,e,1,4*25,1E,27*25-1+i+e,27*25+i+e *enddo*do,f,1,3*25,1E,50*25-1+i+f,50*25+i+f*enddo*do,g,1,1*25,1E,57*25-1+i+g,57*25+i+g *enddo*do,h,1,1*25,1E,60*25-1+i+h,60*25+i+h *enddo*do,hh,1,1*25,1E,67*25-1+i+hh,67*25+i+hh *enddo!9截面type,1mat,1secnum,16*do,a,1,5*25,1E,12*25-1+i+a,12*25+i+a *enddo*do,b,1,1*25,1E,19*25-1+i+b,19*25+i+b *enddo*do,c,1,4*25,1E,21*25-1+i+c,21*25+i+c *enddo*do,d,1,4*25,1E,53*25-1+i+d,53*25+i+dE,58*25-1+i+e,58*25+i+e*enddo*do,g,1,6*25,1E,61*25-1+i+g,61*25+i+g*enddo*do,h,1,1*25,1E,84*25-1+i+h,84*25+i+h*enddo!8截面type,1mat,1secnum,15*do,a,1,1*25,1E,44*25-1+i+a,44*25+i+a*enddo!左边桁*elseif,i,eq,Nmax1+(92*25+1)*2+1,then !13截面type,1mat,1secnum,20*do,a,1,2*25,1E,i-1+a,i+a*enddo*do,b,1,6*25,1E,3*25-1+i+b,3*25+i+b*enddo*do,c,1,3*25,1E,70*25-1+i+c,70*25+i+c*enddo*do,d,1,10*25,1E,74*25-1+i+d,74*25+i+d*enddo*do,e,1,3*25,1E,85*25-1+i+e,85*25+i+e*enddo*do,f,1,1*25,1E,91*25-1+i+f,91*25+i+f*enddo!12截面type,1mat,1E,i+2*25-1+a,i+2*25+a*enddo*do,b,1,1*25,1E,9*25-1+i+b,9*25+i+b *enddo*do,c,1,8*25,1E,35*25-1+i+c,35*25+i+c *enddo*do,d,1,1*25,1E,i+69*25-1+d,i+69*25+d *enddo*do,e,1,1*25,1E,73*25-1+i+e,73*25+i+e *enddo*do,f,1,3*25,1E,88*25-1+i+f,88*25+i+f *enddo!11截面type,1mat,1secnum,18*do,a,1,25,1E,10*25-1+i+a,10*25+i+a *enddo*do,b,1,1*25,1E,18*25-1+i+b,18*25+i+b *enddo*do,c,1,1*25,1E,26*25-1+i+c,26*25+i+c *enddo*do,d,1,4*25,1E,31*25-1+i+d,31*25+i+d *enddo*do,e,1,1*25,1E,43*25-1+i+e,43*25+i+e *enddo*do,f,1,5*25,1E,45*25-1+i+f,45*25+i+f *enddo*do,g,1,1*25,1E,59*25-1+i+g,59*25+i+g *enddo*do,h,1,1*25,1*enddo!10截面type,1mat,1secnum,17*do,a,1,1*25,1E,11*25-1+i+a,11*25+i+a *enddo*do,b,1,1*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,20*25-1+i+c,20*25+i+c *enddo*do,d,1,1*25,1E,25*25-1+i+d,25*25+i+d *enddo*do,e,1,4*25,1E,27*25-1+i+e,27*25+i+e *enddo*do,f,1,3*25,1E,50*25-1+i+f,50*25+i+f*enddo*do,g,1,1*25,1E,57*25-1+i+g,57*25+i+g *enddo*do,h,1,1*25,1E,60*25-1+i+h,60*25+i+h *enddo*do,hh,1,1*25,1E,67*25-1+i+hh,67*25+i+hh *enddo!9截面type,1mat,1secnum,16*do,a,1,5*25,1E,12*25-1+i+a,12*25+i+a *enddo*do,b,1,1*25,1E,19*25-1+i+b,19*25+i+b *enddo*do,c,1,4*25,1E,21*25-1+i+c,21*25+i+c*do,d,1,4*25,1E,53*25-1+i+d,53*25+i+d *enddo*do,e,1,1*25,1E,58*25-1+i+e,58*25+i+e *enddo*do,g,1,6*25,1E,61*25-1+i+g,61*25+i+g *enddo*do,h,1,1*25,1E,84*25-1+i+h,84*25+i+h *enddo!8截面type,1mat,1secnum,15*do,a,1,1*25,1E,44*25-1+i+a,44*25+i+a *enddo!右边桁*else!13截面type,1mat,1secnum,20*do,a,1,2*25,1E,i-1+a,i+a*enddo*do,b,1,6*25,1E,3*25-1+i+b,3*25+i+b *enddo*do,c,1,3*25,1E,70*25-1+i+c,70*25+i+c *enddo*do,d,1,10*25,1E,74*25-1+i+d,74*25+i+d *enddo*do,e,1,3*25,1E,85*25-1+i+e,85*25+i+e *enddo*do,f,1,1*25,1E,91*25-1+i+f,91*25+i+f *enddotype,1mat,1secnum,19*do,a,1,1*25,1E,i+2*25-1+a,i+2*25+a*enddo*do,b,1,1*25,1E,9*25-1+i+b,9*25+i+b *enddo*do,c,1,8*25,1E,35*25-1+i+c,35*25+i+c *enddo*do,d,1,1*25,1E,i+69*25-1+d,i+69*25+d *enddo*do,e,1,1*25,1E,73*25-1+i+e,73*25+i+e *enddo*do,f,1,3*25,1E,88*25-1+i+f,88*25+i+f *enddo!11截面type,1mat,1secnum,18*do,a,1,25,1E,10*25-1+i+a,10*25+i+a *enddo*do,b,1,1*25,1E,18*25-1+i+b,18*25+i+b *enddo*do,c,1,1*25,1E,26*25-1+i+c,26*25+i+c *enddo*do,d,1,4*25,1E,31*25-1+i+d,31*25+i+d *enddo*do,e,1,1*25,1E,43*25-1+i+e,43*25+i+e *enddo*do,f,1,5*25,1E,45*25-1+i+f,45*25+i+f *enddo*do,g,1,1*25,1*enddo*do,h,1,1*25,1E,68*25-1+i+h,68*25+i+h *enddo!10截面type,1mat,1secnum,17*do,a,1,1*25,1E,11*25-1+i+a,11*25+i+a *enddo*do,b,1,1*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,20*25-1+i+c,20*25+i+c *enddo*do,d,1,1*25,1E,25*25-1+i+d,25*25+i+d *enddo*do,e,1,4*25,1E,27*25-1+i+e,27*25+i+e *enddo*do,f,1,3*25,1E,50*25-1+i+f,50*25+i+f*enddo*do,g,1,1*25,1E,57*25-1+i+g,57*25+i+g *enddo*do,h,1,1*25,1E,60*25-1+i+h,60*25+i+h *enddo*do,hh,1,1*25,1E,67*25-1+i+hh,67*25+i+hh *enddo!9截面type,1mat,1secnum,16*do,a,1,5*25,1E,12*25-1+i+a,12*25+i+a *enddo*do,b,1,1*25,1E,19*25-1+i+b,19*25+i+b*do,c,1,4*25,1E,21*25-1+i+c,21*25+i+c*enddo*do,d,1,4*25,1E,53*25-1+i+d,53*25+i+d*enddo*do,e,1,1*25,1E,58*25-1+i+e,58*25+i+e*enddo*do,g,1,6*25,1E,61*25-1+i+g,61*25+i+g*enddo*do,h,1,1*25,1E,84*25-1+i+h,84*25+i+h*enddo!8截面type,1mat,1secnum,15*do,a,1,1*25,1E,44*25-1+i+a,44*25+i+a*enddo*endif*SET,j,j+14*enddo!********************************************* !!!FuGan!!!********************************************* /prep7*SET,j,-14*do,i,1,(94*25+1)*2,94*25+1k,i,0,j,0k,94*25+i,1363.000,j,0kFILL,i,94*25+i,93,,,,,j=j+14*enddo*SET,Nmax1,10000*SET,j,-14k,i,14.50,j,15.00k,92*25+i,1348.50,j,15.00 kFILL,i,92*25+i,91,,,,,j=j+14*enddo!连接右斜线1j=Nmax1+1*do,n,1,(94*25+1)*2,94*25+1 kk=0*do,i,n,n-1+20*25,25lstr,i,jkk=kk+1*if,kk,eq,1,thenlatt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,4,and,kk,le,6,then latt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,8,and,kk,le,9,then latt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,eq,7,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elselatt,1,1,1,,,,23esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(20-0)*25*enddo!连接右斜线2j=Nmax1+1+40*25*do,n,1,(94*25+1)*2,94*25+1 kk=0*do,i,n+40*25,n-1+60*25,25 lstr,i,jkk=kk+1latt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,3,and,kk,le,9,then latt,1,1,1,,,,21esize,,25lmesh,all*elselatt,1,1,1,,,,22esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(60-40)*25*enddo!连接左斜线1j=Nmax1+19*25+1*do,n,1,(94*25+1)*2,94*25+1 kk=0*do,i,n+21*25,n-1+41*25,25 lstr,i,jkk=kk+1*if,kk,eq,1,thenlatt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,13,and,kk,le,18,then latt,1,1,1,,,,21esize,,25lmesh,all*elselatt,1,1,1,,,,22esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(41-21)*25*enddo!连接左斜线2j=Nmax1+59*25+1*do,n,1,(94*25+1)*2,94*25+1 kk=0kk=kk+1*if,kk,ge,1,and,kk,le,12,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,21,and,kk,le,22,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,15,and,kk,le,16,then latt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,ge,24,and,kk,le,28,then latt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,ge,33,and,kk,le,34,then latt,1,1,1,,,,23esize,,25lmesh,all*elselatt,1,1,1,,,,21esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(95-61)*25*enddo!连接shu线j=Nmax1+1*do,n,1,(94*25+1)*2,94*25+1 kk=0*do,i,n+1*25,n-1+94*25,25 lstr,i,jkk=kk+1*if,kk,ge,1,and,kk,le,3,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,22,and,kk,le,58,then latt,1,1,1,,,,22*elseif,kk,ge,62,and,kk,le,74,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,78,and,kk,le,83,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,89,and,kk,le,93,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,eq,7,thenlatt,1,,1,,,,23esize,,25lmesh,all*elseif,kk,eq,60,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,eq,75,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,eq,86,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,eq,87,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elselatt,1,1,1,,,,21esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-93*25*enddomodmsh,detackldele,all!********************************************* !!!2X+上平联!!!********************************************* /prep7k,1,14.5,-14,15.00k,2,14.5,-7,15.00k,3,14.5,-14,10.2k,4,14.5,-7,10.2lstr,1,4lstr,2,3lstr,3,4lgen,2,all,,,,7,nummrg,kpnumcmp,kplstr,2,4lgen,93,all,,,14.5,,lgen,2,all,,,,14,,latt,1,,1,,,,24esize,,10lmesh,allmodmsh,detackldele,allkdele,all*SET,j,-14*do,i,1,(92*25+1)*2,92*25+1k,i,14.50,j,15.00k,92*25+i,1348.50,j,15.00kFILL,i,92*25+i,91,,,,,j=j+14*enddo*do,i,1,92*25+1,25lstr,i,i+92*25+1lstr,i+92*25+1,i+(92*25+1)*2*enddolatt,1,,1,,,,31esize,,20lmesh,allmodmsh,detackldele,all*do,i,1,92*25,25lstr,i+92*25+1,i+(92*25+1)*2+1*25*enddo*do,i,1+1*25,92*25+1,25lstr,i,i+92*25+1-25lstr,i+92*25+1,i+(92*25+1)*2-25*enddonummrg,kpnumcmp,kplatt,1,,1,,,,32esize,,20lmesh,allmodmsh,detackldele,allkdele,all!************************************************! !!1#ta!!!************************************************ /prep7nn=400000n,1+nn,290,0,158.9n,13+nn,290,0,132.5n,20+nn,290,0,115n,22+nn,290,0,110.5n,63+nn,290,0,110.5n,52+nn,290,-19,-6n,62+nn,290,-11.5,-34.5n,93+nn,290,19,-6n,103+nn,290,11.5,-34.5n,104+nn,290,-19,-6n,107+nn,290,19,-6fill,1+nn,13+nn,11,fill,13+nn,20+nn,6,fill,20+nn,22+nn,1,fill,22+nn,52+nn,29,fill,52+nn,62+nn,9,fill,63+nn,93+nn,29,fill,93+nn,103+nn,9,fill,104+nn,107+nn,2,real,1type,1mat,2*do,a,1+nn,12+nn,1 E,a,a+1*enddotype,1mat,2secnum,52*do,b,13+nn,19+nn,1 E,b,b+1*enddotype,1mat,2secnum,51*do,b,20+nn,21+nn,1 E,b,b+1*enddotype,1mat,2secnum,46*do,b,22+nn,51+nn,1 E,b,b+1*enddotype,1mat,2secnum,42*do,b,52+nn,61+nn,1 E,b,b+1*enddotype,1mat,2secnum,46*do,b,63+nn,92+nn,1 E,b,b+1*enddotype,1mat,2secnum,42*do,b,93+nn,102+nn,1 E,b,b+1*enddotype,1mat,2secnum,54*do,b,104+nn,106+nn,1 E,b,b+1!************************************************ !!2#ta!!************************************************* nn=450000n,1+nn,870,0,158.9n,13+nn,870,0,132.5n,20+nn,870,0,115n,22+nn,870,0,110.5n,63+nn,870,0,110.5n,52+nn,870,-19,-6n,62+nn,870,-11.5,-34.5n,93+nn,870,19,-6n,103+nn,870,11.5,-34.5n,104+nn,870,-19,-6n,107+nn,870,19,-6fill,1+nn,13+nn,11,fill,13+nn,20+nn,6,fill,20+nn,22+nn,1,fill,22+nn,52+nn,29,fill,52+nn,62+nn,9,fill,63+nn,93+nn,29,fill,93+nn,103+nn,9,fill,104+nn,107+nn,2,real,1type,1mat,2secnum,53*do,a,1+nn,12+nn,1E,a,a+1*enddotype,1mat,2secnum,52*do,b,13+nn,19+nn,1E,b,b+1*enddotype,1mat,2secnum,51*do,b,20+nn,21+nn,1E,b,b+1*enddosecnum,46*do,b,22+nn,51+nn,1E,b,b+1*enddotype,1mat,2secnum,42*do,b,52+nn,61+nn,1E,b,b+1*enddotype,1mat,2secnum,46*do,b,63+nn,92+nn,1E,b,b+1*enddotype,1mat,2secnum,42*do,b,93+nn,102+nn,1E,b,b+1*enddotype,1mat,2secnum,54*do,b,104+nn,106+nn,1E,b,b+1*enddo!************************************************ !!link!!************************************************ /prep7nn=300000mm=400000!*****************************ES*do,j,nn+1,nn+2301*2+1,2301*do,a,j,10*25+j,25type,3mat,1real,3E,a,mm+3mm=400000*enddo*do,j,nn+1+11*25,nn+2301*3,2301 *do,a,j,2*25+j,25type,3mat,1real,4E,a,mm+14mm=mm+1*enddomm=400000*enddo*do,j,nn+1+14*25,nn+2301*3,2301 *do,a,j,2*25+j,25type,3mat,1real,5E,a,mm+17mm=mm+1*enddomm=400000*enddo*do,j,nn+1+17*25,nn+2301*3,2301 type,3mat,1real,6E,j,mm+20mm=400000*enddo!****************************EM mm=400000*do,j,nn+1+21*25,nn+2301*3,2301 type,3mat,1real,6E,j,mm+20mm=400000*enddo*do,j,nn+1+22*25,nn+2301*3,2301 *do,a,j,2*25+j,25type,3mat,1real,5mm=mm-1*enddomm=400000*enddo*do,j,nn+1+25*25,nn+2301*3,2301*do,a,j,2*25+j,25type,3mat,1real,4E,a,mm+16mm=mm-1*enddomm=400000*enddo*do,j,nn+1+28*25,nn+2301*3,2301*do,a,j,10*25+j,25type,3mat,1real,3E,a,mm+13mm=mm-1*enddomm=400000*enddo!***************************************WM mm=450000*do,j,nn+1+40*25,nn+2301*3,2301*do,a,j,10*25+j,25type,3mat,1real,3E,a,mm+3mm=mm+1*enddomm=450000*enddo*do,j,nn+1+51*25,nn+2301*3,2301*do,a,j,2*25+j,25type,3mat,1real,4E,a,mm+14mm=mm+1*enddo*enddo*do,j,nn+1+54*25,nn+2301*3,2301*do,a,j,2*25+j,25type,3mat,1real,5E,a,mm+17mm=mm+1*enddomm=450000*enddo*do,j,nn+1+57*25,nn+2301*3,2301type,3mat,1real,6E,j,mm+20mm=450000*enddo!***************************************WS mm=450000*do,j,nn+1+61*25,nn+2301*3,2301type,3mat,1real,6E,j,mm+20mm=450000*enddo*do,j,nn+1+62*25,nn+2301*3,2301*do,a,j,2*25+j,25type,3mat,1real,5E,a,mm+19mm=mm-1*enddomm=450000*enddo*do,j,nn+1+65*25,nn+2301*3,2301*do,a,j,2*25+j,25type,3mat,1real,4E,a,mm+16mm=mm-1mm=450000*enddo*do,j,nn+1+68*25,nn+2301*3,2301 *do,a,j,10*25+j,25type,3mat,1real,3E,a,mm+13mm=mm-1*enddomm=450000*enddonummrg,nodenumcmp,nodeFINISH。

作业1实验设计Design ExplorationTraining Manual 1•Goal目标–演示Design ExplorationDesign Exploration中进行DOE分析的流程，并且建立响应图；–边界条件如图所示Design Modeler ds_cutout–在Design Modeler中建立模型输入参数：ds_cutout –从中得到的几何参数Bearing load（轴承载荷）输出参数：Mass（质量）Equivalent stress（等效应力）Total deformation（整体变形）Training Manual12Parameter Set 检查所有1. File>Open>Link1.dsdb1. 双击Parameter Set ，检查所有输入和输出参数。

2输入参数输出参数3. 返回到项目中31Training Manual4. 双击“Response Surface”启动DOE分析45. 双击DOEDOE 大纲给出了输入和输出DOE参数51Training Manual6. 在Outline of Design of Experiments中点击选择参数67. 在特性中定义设计变量的类型和上下限。

Ds-cutout，4.5-5-5.5上的连续变量。

78. Bearing load（负载），9-10-1111 上的连续变量8Training Manual19. 选择DOE –默认的DOE 类型是中心组合设计（Central Composite Design ）910. 查看和更新设计点1011. 点击Show Progress 展开状态栏11Training Manual1DOE 的表中的给出了9个设计点131312. 选择显示整体变形对应设计点的关系曲线12. 点击Design Points vs Parameter1Training Manual1414. 返回到项目页（ProjectPage）15. 双击ResponseSurface1516. 更新ResponseSurface17. 双击Response17 16Training Manual118. 选择二维模式，如图设置X ，Y 轴1819.选择三维模式，如图设置X ，Y 和Z 轴191Training Manual 20. 点击spider 和local sensitivities 显示图表20Training Manual121. 在响应面上点击鼠标右键选择Insert as Response Point ，将其插入到响应点22. 在需要的响应点上点击鼠标右键选择Insert as Design Point ，将其插入到设计点21其插到设计点22响应点不能和此处给出的一样！Training Manual123. 返回到Project Page2324. 双击Parameter set25. 更新所有设计点242526. 在DP1上点击鼠标右键选择Copy inputs to current 和Updated selected Design Point回到2627. 返回到Project Page271Training Manual28. 双击Solution28检查力学结果作业2What if分析多目标优化设计Training Manual•Goal （目标）–使用参数管理器探索如图所示结构的应力、质量和变形行为因为在垂直载荷的作用下几何参数是在发形行为，因为在垂直载荷的作用下，几何参数是在发生变化的。