利用ANSYS生成变截面箱梁
用ANSYS 软件进行变截面翼板主梁优化设计
用ANSYS 软件进行变截面翼板主梁优化设计3.4.1 参数化建模在双梁箱形龙门起重机结构设计时应选择最佳结构参数，使产品达到最佳经济和技术性能指标，降低自重，提高产品的性能价格比，使梁的强度和刚度都达到许用值。
由第4 章ANSYS 分析梁的截面尺寸过程可知在梁全长范围内有时刚度有富余，有时强度有富余，因此有必要进一步优化。
3.4.2. 设定模型参数上下翼板厚度：A=30主腹板厚度：C =24主梁宽度：B = 3500主梁高度：H = 8500主梁长度：L = 120000翼板厚度变化幅值：M = 2主梁长度变化次数：N =203.4.3. 编制翼板厚度离散化的语句因为工艺制造的原因，主梁翼板变化是离散的，因此生成的关键点这一步骤必须实现参数化。
ANSYS 的参数化设计语言（APDL ）中有分支和循环的功能，利用这一功能，将关键点的标号与其坐标关联起来，实现参数化。
这一过程如下：I =20*DO,l,1,4* N + l,1* IF,MOD ( 1,2 ),NE,0,THEN //分支语言设定关键点的Z坐标。
Z =0ELSEIF,MOD ( l,4 ),EQ ,0Z =MELSE,Z ＝-MEND IFK,I ( I-l ) *L ( 4*10 ),O,ZEND DOFLST,3,4*N + 1,3DO,I,l,4*N + 1,lFITEM,3,IEND DOBSPLIN, PSlX在语句K,I,( I-1 ) * L / ( 4*N ),O,Z 中，K 表示关键点，I 表示关键点的标号，（I-1 ) * L / ( 4*N ）表示关键点的X 坐标，Y 坐标为零，Z 坐标为Z 。
用上面的这一程序就能自动画出一系列的关键点，如果需要改变翼板的形状的时候，只需要改变参数N 和M 就行了。
这一程序虽然小，却是非常关键的一步，它直接关系到优化设计能否进行。
因为模型的所有尺寸都己经在建模开始前设计了ANSYS 的参数，因此在建模中涉及模型的尺寸必须用这些参数而不能数字，在设定的梁单元常数时候就必须用参数表示，建模的其它方面与前面分析主梁时完全一致。
ANSYS中自定义梁截面分析的一个例子
ANSYS中自定义梁截面分析的一个例子finish/clear,nostart/prep7et,1,82 !注意划分截面单元需要用PLANE82单元。
此处千万别用梁单元mp,prxy,1,1/3 !定义泊桑比mp,ex,1,2.07e11 !定义弹性模量k,1,0,0k,2,400,0k,3,400,40k,4,0,40a,1,2,3,4!lsel,all!lesize,all,0,,4!smrtsize,5!其实上面这三条命令这样划分网格也可以，但是没有下面的方法好。
下面更好控制网格质量。
asel,all !选择所有的面esize,,5 !确定划分网格的大小，每边划成5份。
如果你愿意，也可以改变网格数量。
amesh,all !给截面划分网格secwrite,jm2,sect,,1 !将截面命名为jm2.sect写入当前工作目录中去。
!下面开始建立梁中轴线。
注意截面保存后，生成的关键点、面、线等都已经不存在，所以即使不重新开始，!编号也是重新开始。
但是我还是重新开始，清空了内存数据。
finish/clear,nostart/prep7et,1,188 !定义梁单元，该单元必须是支持采用自定义截面的单元类型。
mp,prxy,1,1/3 !定义泊桑比mp,ex,1,2.07e11 !定义弹性模量sectype,1,beam,mesh !指定读入的截面类型在后面使用中编号secoffset,cent !指定截面在梁纵轴上的偏移量secread,'jm2','sect',,mesh !读入截面。
如果截面保存在其他路径，可以采用绝对路径的方法确定SECPLOT,1,1 !画出截面，并显示截面的网格划分。
k,5,1,10000k,6,1,0k,7,1,0,5000k,8,5000,0,0 !前两个关键点是为了建立梁，后两个作为方向关键点使用l,5,6lsel,s,line, ,1,5,1 !选择梁单元的轴线latt,1,,1,,7,8,1 !将材料号、截面参考号、实常数（如果有的话）、方向关键点等信息分配给!上面已经选择好的还没有划分单元的梁轴线/lesize,all,,,10 !指定梁纵向划分网格的尺寸。
ANSYS变截面梁
注：此为网上下载文档，其中有许多看不太明白的地方，故按个人想法进行了更改，更改后的文档也许看的亦不明白，大概了解方法即可。
文档所在网址：/view/04e43ebbfd0a79563c1e722d.html一、前言变截面梁，即两端截面不同的梁。
在工程中通常将梁受力吃紧的一端采用较大的截面，另一端采用较小的截面，以实现等强度设计，节省原材料。
这种变截面梁生活中普遍存在，远至古代娶亲用的花轿轿杆、农村大车车辕，直至当今建筑钢结构中各种采光大蓬横梁（包括弯曲横梁）、各种连杆摇臂结构等。
特别是悬臂梁几乎处处可见，就连输电用的水泥电杆也做成根粗尖细的。
由于这种变截面梁在工程中普遍存在，在ANSYS程序结构分析中，专门设置了变截面梁单元的功能，使用起来特别方便。
具体操作过程在下面的实例分析中详细介绍。
自由度耦合即构件连接处两个节点的自由度（包括移动自由度和转动自由度）变化是一致的，主节点如何变化，从节点随着同样变化。
自由度耦合在静力分析时常用在连接件上，特别用在具有转动的连接件上。
例如：汽车挂挡手柄连接端的球铰，各种销钉与耳环的连接，各种转盘与转轴的连接等。
自由度耦合的概念与自由度释放的概念正好相反。
在ANSYS程序中没有自由度释放功能，只有自由度耦合功能，但用自由度耦合功能完全可以达到自由度释放的目的。
这种功能具体操作过程，也在下面的实例分析中详细介绍。
二、雨蓬计算分析该雨蓬结构是由11根变截面工字钢梁、1根等截面工字钢梁、3根圆管钢梁、1根槽钢梁和6根拉杆钢组成。
其结构如下：图1雨蓬结构1．原始数据（1）坐标数据关键点号XYZ1：0，0，02：0，0，63003：0，4500，0 （第一拉杆上端点）4：0，1000，0 （变截面梁Z轴方向）5：0，1000，0 （第一排工字钢Z轴方向）6：0，1000，6300（第六排槽钢Z轴方向）7: 0,0,4200 （拉杆与雨棚连接点，亦自由度耦合处）（2）材料数据主钢梁变截面工字钢：大端300×150×8×6，小端150×150×8×6第一横梁工字钢：150×100×8×6第二、第三、第四横梁圆管钢：Φ102×4 Ri=47 R0=51第五横梁槽钢：160×80×8×6拉杆：Φ102×4A=1231.5 =0.0001(3) 载荷由玻璃均布载荷计算而得:中间节点: F=5560N边节点: F=2780N角节点: F=1690N(4) 单元梁单元: BEAM188杆单元: LINK8(5) 边界条件主梁固定于墙上:位移和转角全约束；拉杆与主梁连接: 位移耦合,转角自由；拉杆与墙连接: 位移约束, 转角自由。
ansys建模实例-混凝土箱梁(命令流)
/TITLE, BRIDGE XX/PREP7! 定义单元类型ET,1,SHELL63ET,2,SOLID45ET,3,LINK8!定义1,钢，2，混凝土，3，无限刚度MP,EX,1,2.1E11MP,DENS,1,7.8E3MP,PRXY,1,0.3mp,alpx,1,1.0e-5MP,EX,2,3.5E10MP,DENS,2,2.6E3MP,PRXY,2,0.1667MP,EX,3,3.5E15MP,DENS,3,7.8E3MP,PRXY,3,0.3*set,nz,154!定义箱梁N1R,1,7,7,7,7!定义箱梁N2R,2,8,8,8,8!定义箱梁N3R,3,4,4,4,4!定义混凝土solid65的实常数R, 4!定义钢绞线的实常数R, 5, 140!降温=预加力/线膨胀系数*钢筋截面积*弹性模量tem1=-1190.4e6/(1.0e-5*2.11e11)tem2=-1339.2e6/(1.0e-5*2.11e11)!定义纵向梯度参数表*dim, z, array, nz*SET,z(1),130,260,371,482,593,702,811,920,1029,1138,1247,1356,1465,1574,1683,17 87.19,1891.38,1995.57*SET,z(19),2099.76,2203.95,2308.14,2412.33,2516.52,2620.71,2724.9,2829.09,2933. 28,3037.47,3141.66,3245.85,*SET,z(31),3350.04,3454.23,3558.42,3662.61,3766.8,3870.99*SET,z(37),3975.18,4079.37,4183.56,4287.75,*SET,z(41),4391.94,4496.13,4600.32,4704.51,4808.7,4912.89,5017.08,5121.27,5225. 46,5329.65,5433.84,5538.03,5642.22,5746.41*SET,z(55),5850.6,5954.79,6058.98,6163.17,6267.36,6371.55,6475.74,6579.93,6682. 47,6785.01,*SET,z(65),6887.55,6990.09,7092.63,7195.17,7297.71,7400.25,7502.79,7605.33 *SET,z(73),7707.87,7810.41,7912.95,8024.28,8135.61,8246.94,8358.27,*SET,z(80),8469.6,8580.93,8680.93,8780.93,8880.93,8980.93,9080.93,9180.93,9280. 93,9380.93,9480.93*SET,z(91),9580.93,9691.93,9802.93,9913.93,10024.93,10135.93,10246.93,10346.93, 10446.93,*SET,z(100),10546.93,10646.93,10746.93,10846.93,10946.93,11046.93,11146.93,1124 6.93,11358.26*SET,z(109),11469.59,11580.92,11692.25,11803.58,11914.91,12005.91,*SET,z(115),12096.91,12187.91,12278.91,12369.91,12460.91,12551.91,12642.91,1273 3.91,12824.91,12915.91,13006.91,13125.91*SET,z(127),13244.91,13363.91,13482.91,13601.91,13720.91,13839.91,13939.91,1403 9.91,*SET,z(135),14139.91,14239.91,14339.91,14439.91,14539.91,14639.91,14739.91,1483 9.91,14942.41,15044.91*SET,z(145),15147.41,15249.91,15352.41,15454.91,15557.41,15659.91,15770.91,1588 1.91,15992.91,16109.91!!定义横隔梁的梯度参数表*set,nd,18*dim, del, array, nd*SET,del(1), 2,5,15,31,47,63,76,82,92,95,98,108,114,125,132,142,*SET,del(17),150,153!!定义初始面座标参数表*set,ni,36*dim, initial, array,ni, 3*SET,initial(1,1),-633,-533,-399.6667,-266.333,-133,0,133.333,266.6666667,400,5 33.333,633,-633,-533,-399.667,-266.3333333,-133,*SET,initial(17,1),0,133.3333333,266.6666667,400,533.3333333,633,-633,-633,633, 633,-766,-699.5,-633,-566.5,-500,500,566.5,633,*SET,initial(35,1),699.5,766*SET,initial(1,2),0,0,0,0,0,0,0,0,0,0,0,0*SET,initial(12,2),60,60,60,60,60,60,60,60,60,60,60,*SET,Initial(23,2),184,276,184,276,368,368,368,368,368,368,368,368,368,368, *SET,initial(1,3),115.82 ,97.53 ,73.13 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-73 .19 ,-97.59 ,-115.82 ,*SET,initial(12,3),115.82 ,97.53 ,73.13 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-7 3.19 ,-97.59 ,-115.82 ,115.82 ,*SET,initial(24,3),115.82 ,-115.82 ,-115.82 ,140.16 ,127.99 ,115.82 ,103.66 ,91.49 ,-91.49 ,-103.66 ,-115.82 ,-127.99 ,-140.16!定义横隔梁参数表*set,nv,25*dim,vertical,array,nv,3*SET,vertical(1 ,3),97.53 ,73.13 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-73.19 ,-97.59 ,97.53 ,73.13 ,48.73 ,24.34*SET,vertical(14,3),0.00 ,-24.40 ,-48.79 ,-73.19 ,-97.59 ,73.19 ,48.73 ,24.34 , 0.00 ,-24.40 ,-48.79 ,-73.19*SET,vertical(1,2),184,184,184,184,184,184,184,184,184,*SET,vertical(10,2),276,276,276,276,276,276,276,276,276,*SET,vertical(19,2),368,368,368,368,368,368,368*SET,vertical(1,1),-533,-399.667 ,-266.333 ,-133 ,0 ,133.333 ,266.667 ,400.000 ,533.333 ,-533.000 ,-399.667 ,-266.333 ,-133*SET,vertical(14,1),0,133.333 ,266.667 ,400.000 ,533.333 ,-400.000 ,-266.333 ,-133.000 ,0.000 ,133.333 ,266.667 ,400.000!!定义混凝土的节点参数表*set,nc,135*dim,concrete,array,nc,3*SET,concrete(1,1),-1301.625 ,-1301.625 ,-1301.625 ,-1248.225 ,-1248.225 ,-1248 .225 ,-1194.825 ,-1194.825 ,-1194.825 ,-1141.425*SET,concrete(11,1),-1141.425 ,-1141.425 ,-1088.025 ,-1088.025 ,-1088.025 ,-103 4.625 ,-1034.625 ,-1034.625 ,-981.225 ,-981.225*SET,concrete(21,1),-981.225 ,-927.825 ,-927.825 ,-927.825 ,-874.425 ,-874.425 ,-874.425 ,-821.025 ,-821.025 ,-821.025*SET,concrete(31,1),-766.000 ,-766.000 ,-766.000 ,-699.500 ,-699.500 ,-699.500 ,-633.000 ,-633.000 ,-633.000 ,-566.500*SET,concrete(41,1),-566.500 ,-566.500 ,-500.000 ,-500.000 ,-500.000 ,-450.000 ,-450.000 ,-450.000 ,-400.000 ,-400.000*SET,concrete(51,1),-400.000 ,-333.167 ,-350.000 ,-333.167 ,-266.333 ,-250.000 ,-266.333 ,-199.667 ,-199.667 ,-199.667*SET,concrete(61,1),-133.000 ,-133.000 ,-133.000 ,-66.500 ,-66.500 ,-66.500 ,0. 000 ,0.000 ,0.000 ,66.500*SET,concrete(71,1),66.500 ,66.500 ,133.000 ,133.000 ,133.000 ,199.667 ,199.667 ,199.667 ,266.333 ,250.000*SET,concrete(81,1),266.333 ,333.167 ,350.000 ,333.167 ,400.000 ,400.000 ,400.0 00 ,450.000 ,450.000 ,450.000*SET,concrete(91,1),500.000 ,500.000 ,500.000 ,566.500 ,566.500 ,566.500 ,633.0 00 ,633.000 ,633.000 ,699.500*SET,concrete(101,1),699.500 ,699.500 ,766.000 ,766.000 ,766.000 ,821.025 ,821. 025 ,821.025 ,874.425 ,874.425*SET,concrete(111,1),874.425 ,927.825 ,927.825 ,927.825 ,981.225 ,981.225 ,981. 225 ,1034.625 ,1034.625 ,1034.625*SET,concrete(121,1),1088.025 ,1088.025 ,1088.025 ,1141.425 ,1141.425 ,1141.425 ,1194.825 ,1194.825 ,1194.825 ,1248.225*SET,concrete(131,1),1248.225 ,1248.225 ,1301.625 ,1301.625 ,1301.625*SET,concrete(1, 2), 468.25,443.219 ,418.188 ,468.250 ,440.716 ,413.181 ,468.250 ,443.219 ,408.175 , 468.250*SET,concrete(11, 2), 443.219 ,403.169 ,468.250 ,430.703 ,398.163 ,468.250 ,430.703 ,393.156 ,468.250 ,418.188*SET,concrete(21, 2), 388.150 ,468.250 ,418.188 ,383.144 ,468.250 ,418.188 ,378.138 ,468.250 ,418.188 ,373.131*SET,concrete(31, 2), 468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250*SET,concrete(41, 2), 418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188*SET,concrete(51, 2), 368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000*SET,concrete(61, 2), 468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250*SET,concrete(71, 2), 418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188*SET,concrete(81, 2), 368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000*SET,concrete(91, 2), 468.25,418.188 ,368.000 ,468.25,418.188,368.000 ,468.25,418.188,368,468.25*SET,concrete(101, 2), 418.188,368,468.25,418.188,368,468.250 ,418.188,373.131,468.25,418.188*SET,concrete(111, 2), 378.138,468.25,418.188,383.144,468.25,418.188 ,388.15,468.25,430.703,393.156 *SET,concrete(121, 2), 468.25,430.703,398.163,468.25,443.219,403.169,468.25,443.219,408.175,468.25 *SET,concrete(131, 2), 440.716,413.181,468.25,443.219,418.188,*SET,concrete(1, 3), 238.1649941,238.165 ,238.165 ,228.394 ,228.394 ,228.394 ,218.623 ,218.623 ,218. 623 ,208.852*SET,concrete(11, 3), 208.852 ,208.852 ,199.082 ,199.082 ,199.082 ,189.311 ,189.311 ,189.311 ,179.540,179.540*SET,concrete(21, 3), 179.540 ,169.769 ,169.769 ,169.769 ,159.998 ,159.998 ,159.998 ,150.227 ,150.227 ,150.227*SET,concrete(31, 3), 140.159 ,140.159 ,140.159 ,127.991 ,127.991 ,127.991 ,115.823 ,115.823 ,115.823 ,103.655*SET,concrete(41, 3), 103.655 ,103.655 ,91.488 ,91.488 ,91.488 ,82.339 ,82.339 ,82.339 ,73.190 ,73.190 *SET,concrete(51, 3), 73.190 ,60.961 ,64.041 ,60.961 ,48.732 ,45.744 ,48.732 ,36.534 ,36.534 ,36.534*SET,concrete(61, 3), 24.336 ,24.336 ,24.336 ,12.168 ,12.168 ,12.168 ,0.000 ,0.000 ,0.000 ,-12.168*SET,concrete(71, 3), -12.168 ,-12.168 ,-24.336 ,-24.336 ,-24.336 ,-36.534 ,-36.534 ,-36.534 ,-48.732 ,-45.744*SET,concrete(81, 3), -45.744 ,-60.961 ,-64.041 ,-60.961 ,-73.190 ,-73.190 ,-73.190 ,-82.339 ,-82.339 ,-82.339*SET,concrete(91, 3), -91.48756137,-91.488 ,-91.488 ,-103.65,-103.65,-103.655 ,-115.8232527,-115.8232 527,-115.82,-127.99*SET,concrete(101, 3), -127.991,-127.9910,-140.159,-140.159,-140.159,-150.227 ,-150.23,-150.23,-159.99 80217,-159.9980217*SET,concrete(111, 3), -159.998,-169.768,-169.768,-169.768,-179.539,-179.540 ,-179.5397,-189.3106,-189.3106364,-189.3106364*SET,concrete(121, 3), -199.0815,-199.0815,-199.0815,-208.852,-208.8523,-208.8523,-218.62,-218.623,-21 8.623251,-228.3941226*SET,concrete(131, 3), -228.3941226,-228.3941226,-238.1649941,-238.1649941,-238.1649941,!节点编号，生成箱梁的节点*do, i, 1,1,1*do, j, 1, nin, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)*enddo*enddo*do, i, nz,nz,1*do, j, 1, nin, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)*enddo*enddo*do, i, 2,nz-1,1*do, j, 1, nin, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)*enddo*enddo!生成横隔板节点*set, infl, ni*(nz+1)*set, nfl, 39*do, i, 1, 18, 1*SET,k,del(i)*do, j, 1, 25, 1n, infl+i*nfl+j, vertical(j,1), vertical(j,2), vertical(j,3)+z(k)*enddo*do, j, 26, 32n, infl+i*nfl+j, vertical(j-7,1), vertical(j-7,2), vertical(j-7,3)+z(k-1)*enddo*do, j, 33, nfln, infl+i*nfl+j, vertical(j-14,1), vertical(j-14,2), vertical(j-14,3)+z(k+1) *enddo*enddo!节点编号，生成混凝土的节点*set,incr,infl+18*nfl+1*do, i, 1, 1,*do, j, 1, ncn, incr+j+nc*i, concrete(j,1),concrete(j,2),concrete(j,3)+z(i)*enddo*enddo*do, i, 2, nz*do, j, 1, ncn, incr+j+nc*i, concrete(j,1),concrete(j,2),concrete(j,3)+z(i)*enddo*enddo!生成钢绞线的节点*set,insty1,incr+nc*(nz+1)!Y1*do, i, 37, 127, 1n, insty1+(i-36), concrete(53,1),concrete(53,2),concrete(53,3)+z(i)n, insty1+(i-36)+127-37+1, concrete(68,1),concrete(68,2),concrete(68,3)+z(i) n, insty1+(i-36)+2*(127-37)+2, concrete(83,1),concrete(83,2),concrete(83,3)+z(i) *enddo!Y2*set,insty2,insty1+3*(127-37+1)*do, i, 66,114,1n, insty2+(i-65), concrete(56,1),concrete(56,2),concrete(56,3)+z(i)n, insty2+(i-65)+141-66+1, concrete(80,1),concrete(80,2),concrete(80,3)+z(i) *enddo!生成钢绞线锚固区节点*set,instb1,insty2+2*(141-66)*set,instb2,instb1+40!Y1*do, i, 1, 4*SET,k,27+in, instb1+i, -350.000 , 298.188+12*i, 64.0413+z(k)n, instb2+i, -350.000 , 298.188+12*i, 64.0413+z(k)*enddo*do, i, 1, 4*SET,k,27+in, instb1+i+4, 0 , 298.188+12*i, z(k)n, instb2+i+4, 0 , 298.188+12*i, z(k)*enddo*do, i, 1, 4*SET,k,27+in, instb1+i+8, 350.000 , 298.188+12*i, -64.0413+z(k)n, instb2+i+8, 350.000 , 298.188+12*i, -64.0413+z(k)*enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+12, -350.000 , 418.188-12*(i+6), 64.0413+z(k)n, instb2+i+12, -350.000 , 418.188-12*(i+6), 64.0413+z(k)*enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+16, 0 , 418.188-12*(i+6), z(k)n, instb2+i+16, 0 , 418.188-12*(i+6), z(k)*enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+20, 350.000 , 418.188-12*(i+6), -64.0413+z(k)n, instb2+i+20, 350.000 , 418.188-12*(i+6), -64.0413+z(k)*enddo!Y2*do, i, 1, 4*SET,k,56+in, instb1+i+24, -250.000 , 298.188+12*i, 45.743+z(k)n, instb2+i+24, -250.000 , 298.188+12*i, 45.743+z(k)*enddo*do, i, 1, 4*SET,k,56+in, instb1+i+28, 250.000 ,298.188+12*i, -45.743+z(k)n, instb2+i+28, 250.000 ,298.188+12*i, -45.743+z(k)*enddo*do, i, 1, 4*SET,k,115+i+6n, instb1+i+32, -250.000 ,418.188-12*(i+6), 45.743+z(k) n, instb2+i+32, -250.000 ,418.188-12*(i+6), 45.743+z(k) *enddo*do, i, 1, 4*SET,k,115+i+6n, instb1+i+ni, 250.000 ,418.188-12*(i+6), -45.743+z(k) n, instb2+i+ni, 250.000 ,418.188-12*(i+6), -45.743+z(k) *enddo!生成箱梁的单元*do, i, 1,nz-1,1!底板，shell63，钢，厚度8type,1mat,1real,2*do, j, 1, 10, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j,*enddo!顶板，shell63，钢，厚度7type,1mat,1real,1*do, j, 27, 30, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j,*enddo*do, j, 32, 35, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j,*enddo!腹板、底板加劲肋，shell63，钢，厚度4type,1mat,1real,3*do, j, 1, 11, 1e, i*ni+j, (i+1)*ni+j, (i+1)*ni+j+11, i*ni+j+11*enddoe, i*ni+12, (i+1)*ni+12, (i+1)*ni+23, i*ni+23e, i*ni+23, (i+1)*ni+23, (i+1)*ni+24, i*ni+24e, i*ni+24, (i+1)*ni+24, (i+1)*ni+29, i*ni+29e, i*ni+22, (i+1)*ni+22, (i+1)*ni+25, i*ni+25e, i*ni+25, (i+1)*ni+25, (i+1)*ni+26, i*ni+26e, i*ni+26, (i+1)*ni+26, (i+1)*ni+34, i*ni+34*enddo!生成横隔板单元type,1mat,1real,3*do, i, 1, 18, 1!生成横隔板单元*do, l, 0, 9, 1e, del(i)*ni+1+l, del(i)*ni+2+l, del(i)*ni+13+l, del(i)*ni+12+l,*enddo*do, l, 0, 7,1e, infl+i*nfl+l+1, infl+i*nfl+l+2, infl+i*nfl+l+11, infl+i*nfl+l+10,e, del(i)*ni+l+13,del(i)*ni+l+14, infl+i*nfl+l+2, infl+i*nfl+l+1,*enddo*do, l, 0, 5, 1e, infl+i*nfl+l+11, infl+i*nfl+l+12, infl+i*nfl+l+20, infl+i*nfl+l+19, *enddo*do, l, 0, 5, 1e, infl+i*nfl+l+19, infl+i*nfl+l+20, infl+i*nfl+l+27, infl+i*nfl+l+26 e, infl+i*nfl+l+33, infl+i*nfl+l+34, infl+i*nfl+l+20, infl+i*nfl+l+19 *enddoe, del(i)*ni+31, infl+i*nfl+19, infl+i*nfl+26, (del(i)-1)*ni+31,e, (del(i)+1)*ni+31, infl+i*nfl+33, infl+i*nfl+19, del(i)*ni+31,e, infl+i*nfl+25, del(i)*ni+32, (del(i)-1)*ni+32, infl+i*nfl+32,e, infl+i*nfl+nfl, (del(i)+1)*ni+32, del(i)*ni+32, infl+i*nfl+25,e, del(i)*ni+12, del(i)*ni+13, infl+i*nfl+1,del(i)*ni+23e, del(i)*ni+23, infl+i*nfl+1, infl+i*nfl+10, del(i)*ni+24,e, del(i)*ni+24, infl+i*nfl+10, del(i)*ni+31, del(i)*ni+29,e, infl+i*nfl+10, infl+i*nfl+11, infl+i*nfl+19, del(i)*ni+31,e, del(i)*ni+21, del(i)*ni+22, del(i)*ni+25, infl+i*nfl+9,e, infl+i*nfl+9, del(i)*ni+25, del(i)*ni+26, infl+i*nfl+18,e, infl+i*nfl+18, del(i)*ni+26, del(i)*ni+34, del(i)*ni+32,e, infl+i*nfl+17, infl+i*nfl+18, del(i)*ni+32, infl+i*nfl+25,*enddo!生成混凝土单元type,2mat,2real, 4*do, i, 1, 153, 1*do, j, 0, 130, 3*do, k, 2, 3, 1e,incr+nc*i+j+k,incr+nc*(i+1)+j+k,incr+nc*(i+1)+j+k+3,incr+nc*i+j+k+3,incr+nc*i +j+k-1,incr+nc*(i+1)+j+k-1,incr+nc*(i+1)+j+k+2, incr+nc*i+j+k+2,*enddo*enddo*enddo!生成Y1钢绞线单元type,3mat,1real, 5*do, i, 37, 126, 1e, insty1+i-36,insty1+i-ni+1e, insty1+i-36+127-37+1, insty1+i-ni+127-37+2e, insty1+(i-36)+2*(127-37)+2,insty1+(i-ni)+2*(127-37)+3*enddoe, instb1+4, insty1+1,e, instb1+8, insty1+1+127-37+1e, instb1+12, insty1+1+2*(127-37)+2e, instb1+13, insty1+(127-36)e, instb1+17, insty1+(127-36)+127-37+1,e, instb1+21, insty1+(127-36)+2*(127-37)+2*do, i, 1, 3*do, j ,0, 5e, instb1+i+j*4, instb1+i+j*4+1*enddo*enddofinish/SOL!进入求解器施加荷载ESEL,S,TYPE,,3 !选择现有单元类型为LINK8的单元EPLOTBFE,all,TEMP,1,tem1 !对目前显示的所有单元施加温度荷载，大小为-tem1ESEL,none !全不选择/PREP7!进入前处理器,生成Y2钢绞线的单元*set,insty2,insty1+3*(127-37+1)*do, i, 66,114,1n, insty2+(i-65), concrete(56,1),concrete(56,2),concrete(56,3)+z(i)n, insty2+(i-65)+141-66+1, concrete(80,1),concrete(80,2),concrete(80,3)+z(i) *enddo*do, i, 66, 113, 1e, insty2+i-65,insty2+i-65+1e, insty2+(i-65)+141-66+1, insty2+(i-65)+141-66+2*enddoe, instb1+28, insty2+1e, instb1+32, insty2+1+141-66+1e, instb1+33, insty2+(113-65)e, instb1+37,insty2+(113-65)+141-66+1*do, i, 1, 3*do, j ,0, 3e, instb1+i+j*4+24, instb1+i+j*4+1+24*enddo*enddoFINISH/SOL!进入求解器施加荷载BFE,all,TEMP,1,tem2 !对目前选择的单元即Y2施加温度荷载，大小为-tem2ESEL,all !选择所有的单元FINISH/PREP7!生成钢筋锚固区type,2mat,2real, 4*do, i, 1, 3e, instb2+i, instb2+4+i, instb2+5+i, instb2+1+i, incr+54+nc*(i+27),incr+69+nc*(i+27),incr+69+nc*(i+28), incr+54+nc*(i+28),e, instb2+4+i,instb2+8+i, instb2+9+i,instb2+5+i, incr+nc*(i+27)+69,incr+nc*(i+27)+84, incr+nc*(i+28)+84, incr+nc*(i+28)+69e, instb2+12+i, instb2+16+i, instb2+17+i, instb2+13+i, incr+54+nc*(i+134),incr+69+nc*(i+134),incr+69+nc*(i+135),incr+54+nc*(i+135)e, instb2+16+i, instb2+20+i, instb2+21+i, instb2+17+i, incr+69+nc*(i+134),incr+84+nc*(i+134),incr+84+nc*(i+135),incr+69+nc*(i+135)*enddo*do, i,1, 3e, instb2+24+i, instb2+28+i, instb2+29+i, instb2+25+i,incr+57+nc*(i+56),incr+81+nc*(i+56),incr+81+nc*(i+57),incr+57+nc*(i+57)e, instb2+32+i, instb2+36+i, instb2+37+i, instb2+33+i,incr+57+nc*(i+121),incr+81+nc*(i+121),incr+81+nc*(i+122),incr+57+nc*(i+122)*enddo!加载FINISH/SOL*do, i, 1, 2*do, j, 4,8d, ni*i+j, all,*enddo*enddo*do, i, 94, 96*do, j, 4,8d, ni*i+j, all,*enddo*enddo*do, i, 153, 154*do, j, 4,8d, ni*i+j, ux,*enddo*enddo*do, i, 153, 154*do, j, 4,8d, ni*i+j, uz,*enddo*enddoFINISH/PREP7cpintf, ux, 0.01cpintf, uy, 0.01cpintf, uz, 0.01nsel, all eplot/solu solve。
钢筋混凝土变截面梁有限元仿真分析
钢筋混凝土变截面梁有限元仿真分析打开文本图片集摘要：变截面构件在工程实践中得到了越来越广泛的应用，尤其是在大跨度桥主梁设计中，常采用变截面。
其优点在于构件随着弯矩及剪力的分布变化而合理地变化截面，从而达到了节约材料、降低成本的目的。
利用大型有限元软件ANSYS对变截面梁进行了静力分析，以期为相关工程提供参考。
关键词：桥梁；变截面连续箱梁；有限元分析Keyword：bridge；continuoubo某girderwithvariablecroection；finiteelementanalyi0引言桥梁作为线路的一个重要组成部分，随着时间的流逝，每当运输工具发生改变时，都会对桥梁的性能等各方面提出更高的要求，由此便推动了桥梁工程技术的发展。
在新桥设计、桥梁施工监控和成桥试验检测中，应力、应变，挠度计算都极其重要。
随着桥梁结构的日益复杂，传统计算方法由于耗时耗力，效率低下且精确度低，已经无法满足现代桥梁施工的需要，ANSYS具有计算功能强大、建模方便及计算结果直观等优点，近年来，得到越来越多桥梁建设者和研究者的重视[1]。
本文利用ANSYS10.0软件对某钢筋砼变截面连续箱梁桥进行有限元分析，将一个连续介质体的问题近似地离散为由适当选取的有限个单元、有限个结点连接起来的集合体加以分析，用主梁的一些重要截面的弯矩作为控制条件，首先简化工程问题，抽象模型，确定定解数据，然后根据待分析问题的数学力学模型和有限元软件的功能，输入待分析问题的定解数据，并进行有限元剖分，形成分析问题的有限元模型，使用选定的有限元软件进行分析。
主要包括单元分析，约束处理，求解，计算结果的后处理，得到其应力、变形分布规律。
1有限元模型某桥梁全长1161.8m，上部构造主桥为（60+3某96+60）m。
桥面宽度是13.7m保持不变，如图1。
纵梁的截面，如图2所示。
本文采用实体建模，自由分网。
采用单元类型是SOLID65。
为计算简便和建模方便，在整个建模及分析过程中，材料视为单一材料的混凝土，即将其简化为匀质，不考虑钢筋、预应力及预应力拱的作用。
ansys_变截面梁精品文档6页
1、选择菜单“Utility Menu>File>Change Title”。
2、输入“Lateral Torsional Buckling Analysis”并单击OK。
3、定义关键点。
选择“Main Menu>Preprocessor> -Modeling- Create>Keypoints> In Active4、在关键点1和2之间建立直线。
选择MainMenu>Preprocessor>-Modeling-Create>-Lines-Lines>Straight Line，出现“Create Straight Line picker”窗口，在图形窗口中拾取关键点1和2，然后按“OK”。
5、保存模型。
选择Utility Menu>File>Save As。
将文件名存为“buckle.DB”，并单击OK。
7.6.6 定义单元类型和横截面信息1、选择“Main Menu>Preferences”，单击“Structural”检查框。
按“OK”。
2、选择“Main Menu>Preprocessor>Element Type>Add/Edit/Delete”。
出现“Element Types”对话框。
3、单击“Add...”。
出现“Library of Element Types ”对话框。
4、在左列选择“Structural Beam”。
5、在右列选择“3D finite strain, 3 node 189”以选中 BEAM189。
7、定义梁的矩形截面。
选择“Main Menu>Preprocessor>Sections>-Beam-Common Sects”。
出现“BeamTool ”对话框。
缺省时 ANSYS 将截面号设置为1，将子类型设置为 RECT (在子类型处图示一个矩形)。
变截面箱梁 ANSYS 前处理程序开发与应用
变截面箱梁 ANSYS 前处理程序开发与应用王韬杰【摘要】A program is developed which can turn the geometric information of box girder with varying depth into ANSYS preprocess command flow,using AutoCAD and Visual Fortran program.This pro-gram can produce the high quality hexahedron mesh and adjust the grid density according to user's needs well.This program will save the time of finite element modeling and improve the design calcula-tion efficiency in the premise of the calculation accuracy.Finally,a detailed finite element analysis of a large continuous rigid frame arch bridge is presented with the help of this program.%基于通用有限元分析平台 ANSYS，利用计算机辅助设计软件AutoCAD，以及 Fortran 程序，设计了将变截面箱梁几何信息转换为 ANSYS 前处理命令流，从而进行快速网格剖分的程序。
该程序具有生成高质量六面体网格，以及根据用户需要随意调整网格密度的功能，在保证网格精度的前提下节约了有限元建模时间，提高了计算效率。
运用该程序对某大型连续刚构拱桥进行了局部有限元分析。
【期刊名称】《交通科技》【年(卷),期】2015(000)006【总页数】3页(P25-27)【关键词】程序设计;有限单元法;网格划分【作者】王韬杰【作者单位】天津大学建筑工程学院天津 300072【正文语种】中文混凝土箱梁广泛应用于简支梁、悬臂梁、连续梁、连续刚构、斜拉桥等各种桥梁结构中。
ANSYS定义变截面
ansys变截面梁的定义（ZZ）ansys专区2010-02-04 09:23:22 阅读198 评论0 字号：大中小（例1：转自/forum/viewthread.php?tid=31432&highlight=%B1%E4%BD%D8%C3%E6）变截面梁的定义定义变截面梁之前需定义两个截面：SECTYPE，1，BEAM，HREC，，0 ！空芯矩形截面1SECDATA,10,20,2,2,2,2,0,0,0,0 ！截面1尺寸定义SECTYPE，2，BEAM，HREC，，0 ！空芯矩形截面2SECDATA,8,10,2,2,2,2,0,0,0,0 ！截面2尺寸定义！定义两个截面后,再定义变截面SECTYPE,3,TAPER,, ！变截面3SECDATA,1,,,,SECDATA,2,10 ！截面由1变到2,位置由(0,0,0)变到(10,0,0,)（例2：转自/download/ANSYS_BEAM54.html）finish/clear/PREP7A_HYT1=0.4 !A端A_HYB1=0.1 !A端B_HYT1=0.2 !B端B_HYB1=0.1 !B端OFFSET=0.5 !偏移!*ET,1,BEAM54!*!**SET,_RC_SET,1,R,_RC_SET,0.08,0.0010666666666667,A_HYT1,A_HYB1,RMODIF,_RC_SET,9,0,-OFFSET,RMODIF,_RC_SET,14,0,RMODIF,_RC_SET,5,0.2*0.2,0.2*0.2**3/12,B_HYT1,B_HYB1, RMODIF,_RC_SET,11,0,-OFFSET,RMODIF,_RC_SET,15,0,RMODIF,_RC_SET,13,0,RMODIF,_RC_SET,16,0, , ,!**SET,_RC_SET,2,R,_RC_SET,0.08,0.0010666666666667,A_HYB1,A_HYT1, RMODIF,_RC_SET,9,0,OFFSET,RMODIF,_RC_SET,14,0,RMODIF,_RC_SET,5,0.2*0.2,0.2*0.2**3/12,B_HYB1,B_HYT1, RMODIF,_RC_SET,11,0,OFFSET,RMODIF,_RC_SET,15,0,RMODIF,_RC_SET,13,0,RMODIF,_RC_SET,16,0, , ,!*!*MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,200e3MPDATA,PRXY,1,,.2N,1,0,,,,,,N,2,5,,,,,,N,3,10,,,,,,TYPE, 1MAT, 1REAL, 1!*e,1,2REAL, 2e,3,2/eshap,1EPLOT（例3：转自/archiver/tid-64381.html）VM34是个变截面梁的例子，分别用SHELL63和BEAM44做的。
ansys workbench中变截面壳体的实现
Ansys Workbench中变截面壳体
的实现
目录
1问题描述： (1)
2实体单元（solid 186） (2)
2.1 网格 (2)
2.2 结果 (2)
3实体壳单元（solsh190） (4)
3.1 网格 (4)
3.2 结果 (4)
4壳体单元（shell181） (6)
4.1 网格 (6)
4.2 结果 (7)
5结果对比 (8)
1问题描述：
以下图简单模型为例，介绍如何在Ansys Workbench中实现变截面壳体：
截面尺寸如图1所示，梯形短边长10mm，长边长20mm，高50mm：
图1-1 截面尺寸
将截面拉伸100mm，得到所需模型，如图2所示。
图1-2 例子模型
对梯形长边面施加固定约束，顶端矩形面施加向下压力5Mpa，分别采用实体solid186单元，实体壳单元solsh190，壳体单元shell181进行模型。
2实体单元（solid 186）
2.1网格
图2-1
图2-2 2.2结果
图2-3
图2-4
3实体壳单元（solsh190）
3.1网格
图3-1
图3-2
图3-3 3.2结果
图3-3
图3-4
4壳体单元（shell181）
4.1网格
图4-1
图4-2
图4-3
图4-4 4.2结果
图4-5
图4-6
5结果对比
表5-1 各单元结果对比
通过结果对比可知，三种单元所得位移几乎一致，壳体单元并不能考虑应力集中，所得应力值偏低，实体单元最接近实际情况，实体壳单元其次，工程应用中在误差允许内可采用实体壳单元直接模拟，节省计算时间。
利用ANSYS生成变截面箱梁
利用ANSYS生成变截面箱梁Beam188/189 支持自定义的变截面1、首先建立变截面箱梁截面，并保存截面。
必须保证每个截面的关键点号相同，而且为保证生成的准确性，应尽量的使得关键点有足够多的数目。
2、要对截面进行面积分块，并指定各线段的段数，这样才能做出规整的箱梁截面网格划分，也就保证了变截面箱梁桥各截面的网格模式相同，建议对于变化急剧的两边截面应使得线足够的短。
并不是划分越密越可能成功，而是线越短越可能成功3。
通过梁节段两端对应的截面建立taper截面，定义好桥的线形后，指定每段线对应的taper截面。
4、可以用slist命令查看生成截面的性质，加深对secread 命令的理解以下例进行说明：Ⅰ-Ⅰ截面图全梁1/2图示Ⅰ-Ⅰ截面命令流finish/clear/prep7k,1k,2,2750k,3,3350,2400k,4,5450,2750k,5,6700,2850k,6,6700,3050k,7,0,3050k,8,0,280k,9,2359,280k,10,2886,2400k,11,1836,2750k,12,0,2750l,1,2l,2,3l,3,4l,4,5l,5,6l,6,7l,8,9l,9,10l,10,11l,11,12l,1,8l,7,12al,allarsym,x,all aadd,1,2 nummrg,all numcmp,all adele,1l,2,9l,13,20l,14,19l,10,3l,1,8al,22,17,25,19 al,25,1,21,7 lcomb,15,20 wpoffs,-5450,2750 wprota,,,90 lsbw,18wpcsys,-1 wpoffs,-3350,2400 wprota,,,90 lsbw,20wpcsys,-1 wpoffs,-2886,2400 wprota,,,90 lsbw,18wpcsys,-1l,15,18l,14,21l,19,22l,7,12al,14,13,18,26 al,18,12,29,27 al,29,23,30,28 al,30,15,31,20 lcomb,9,10 wpoffs,5450,2750 wprota,,,90 lsbw,6wpcsys,-1 wpoffs,3350,2400wprota,,,90 lsbw,32 wpcsys,-1 wpoffs,2886,2400 wprota,,,90 lsbw,6wpcsys,-1l,4,11l,3,23l,10,24al,5,4,6,10al,3,6,33,35al,24,35,34,36 al,9,36,32,31 ldiv,16,,,4 ldiv,11,,,4 ldiv,2,,,4ldiv,8,,,4l,30,27l,29,26l,28,25l,36,33l,35,32l,34,31al,23,42,49,39 al,49,41,50,38 al,50,40,51,37 al,51,11,22,16 al,24,48,52,45 al,52,47,53,44 al,53,46,54,43 al,54,8,21,2 nummrg,all numcmp,allet,1,82lesize,1,,,4 lesize,17,,,4 lesize,19,,,4 lesize,7,,,4 lesize,22,,,1 lesize,21,,,1 lesize,25,,,1 mshape,1,2d mshkey,2amesh,2 lesize,23,,,1 lesize,42,,,1 lesize,39,,,1 lesize,49,,,1 lesize,41,,,1 lesize,50,,,1 lesize,38,,,1 lesize,40,,,1 lesize,37,,,1 lesize,51,,,1 lesize,11,,,1 lesize,16,,,1 lesize,22,,,1 lesize,24,,,1 lesize,45,,,1 lesize,48,,,1 lesize,52,,,1 lesize,44,,,1 lesize,47,,,1 lesize,53,,,1 lesize,43,,,1 lesize,46,,,1 lesize,54,,,1 lesize,2,,,1 lesize,8,,,1 lesize,21,,,1 mshape,0,2d mshkey,1 amesh,11 amesh,12 amesh,13 amesh,14 amesh,15 amesh,16 amesh,17 amesh,18 lesize,13,,,2 lesize,26,,,2 lesize,14,,,1 lesize,18,,,1 mshape,1,2d mshkey,2lesize,12,,,3 lesize,27,,,3 lesize,18,,,1 lesize,29,,,1 mshape,1,2d mshkey,2 amesh,4 lesize,23,,,1 lesize,29,,,1 lesize,30,,,1 lesize,28,,,1 mshape,0,2d mshkey,2 amesh,5 lesize,30,,,1 lesize,20,,,4 lesize,31,,,1 lesize,15,,,4 mshape,1,2d mshkey,2 amesh,6 lesize,31,,,1 lesize,9,,,4 lesize,36,,,1 lesize,32,,,4 mshape,1,2d mshkey,2 amesh,10 lesize,35,,,1 lesize,34,,,1 lesize,24,,,1 lesize,36,,,1 mshape,0,2d mshkey,2 amesh,9 lesize,6,,,1 lesize,35,,,1 lesize,3,,,3 lesize,33,,,3 mshape,1,2d mshkey,2 amesh,8 lesize,6,,,1lesize,5,,,1lesize,10,,,2lesize,4,,,2 mshape,1,2dmshkey,2amesh,7secwrite,mybox1主命令流finish/clear/prep7k,5000,7800k,5001,8851k,5002,9899k,5003,10800k,5004,12300k,5005,0l,5000,5001l,5001,5002l,5002,5003l,5003,5004l,5000,5005k,100,8300,4000k,101,9500,4000k,102,10000,4000et,1,beam189mp,ex,1,3.0e10mp,prxy,1,0.167 sectype,1,beam,mesh secread,mybox1,,,mesh sectype,2,beam,mesh secread,mybox2,,,mesh sectype,3,taper secdata,1,7800 secdata,2,8851et,2,beam189mp,ex,2,3.0e10mp,prxy,2,0.167 sectype,4,beam,mesh secread,mybox2,,,mesh sectype,5,beam,mesh secread,mybox3,,,mesh sectype,6,taper secdata,4,8851secdata,5,9899et,3,beam189mp,ex,3,3.0e10mp,prxy,3,0.167 sectype,7,beam,mesh secread,mybox3,,,mesh sectype,8,beam,mesh secread,mybox4,,,mesh sectype,9,taper secdata,7,9899 secdata,8,10800et,4,beam189mp,ex,4,3.0e10mp,prxy,4,0.167 sectype,10,beam,mesh secread,mybox4,,,mesh et,5,beam189mp,ex,5,3.0e10mp,prxy,5,0.167 sectype,11,beam,mesh secread,mybox1,,,mesh lesize,1,,,2lesize,2,,,2lesize,3,,,2lesize,4,,,2lesize,5,,,8allsellsel,s,line,,1latt,1,,1,,100,,3 allsellsel,s,line,,2latt,2,,2,,101,,6 allsellsel,s,line,,3latt,3,,3,,102,,9 allsellsel,s,line,,4latt,4,,4,,102,,10 lsel,s,line,,5latt,5,,5,,102,,11 lmesh,1,5,1allsellmesh,1,5,1/eshape,1eplot。
在ANSYS中实现自定义梁截面
在ＡＮＳＹＳ中实现自定义梁截面ＡＮＳＹＳ软件是美国ＡＮＳＹＳ公司研制的一个功能强大的大型有限元分析软件,具有强大的前处理、求解和后处理功能,目前广泛应用于航空航天、核工业、铁道、石油化工、机械制造、水利水电、生物医学、土木工程、家用产品及科学研究等领域。
用有限元方法进行结构分析时,要将截面划分为若干个区格,在使用ＡＮＳＹＳ进行计算时,划分截面区格数目的多少,不仅影响计算的精度,也可能影响计算的收敛性。
在钢结构中梁、柱截面形状多为工字型、Ｔ型、Ｈ型等,当用ＡＮＳＹＳ软件进行钢结构分析时,ＡＮＳＹＳ提供了梁横截面库,在此截面库中包括了11种常用的截面形状。
但是选用ＡＮＳＹＳ截面库提供的截面无法对该截面的翼缘、腹板进行区格的划分。
为此采用一种自定义截面形状和尺寸的方法,解决了上述问题,实现了截面区格的划分。
具体原则及示例见下述.1.自我定义梁、柱截面1.1 单元类型的选择在ＡＮＳＹＳ中有许多单元,如梁单元、壳单元等,用户应根据计算类型、计算方法、材料等,选择单元类型。
1.2 数据输入并计算(1)输入材料属性:如材料的弹性模量、泊松比、屈服强度等。
(2)算出并依次输入截面各点坐标值,将输入的点依次用线连结起来并建立面。
1.2 对各个面进行布尔操作将各个面进行粘接以及分割操作,防止在截面连接处网格划分困难。
1.3 截面区格的划分其中包括定义网格尺寸和网格划分。
根据计算精度和计算工作量确定网格尺寸。
网格越密,计算精度越高,但是计算工作量越大。
1.4 模块的保存将自定义并且划分区格的截面存入事先建立的截面库中,以备日后使用。
1.5 残余应力问题先定义初应力文件,然后在结构分析中的第一个载荷部中施加残余应力。
通过上述过程就能实现截面区格的划分,使具体结构分析得以顺利进行。
2建模实例以工字形梁截面为例,该梁材料为Ｑ235钢,截面尺寸如下:翼缘宽ｂ=50ｍｍ、厚ｔ1=5ｍｍ;腹板高ｈ0=90ｍｍ、梁截面高ｈ=100ｍｍ、厚为ｔｗ=4.5ｍｍ;材料的弹性模量Ｅ=206000Ｎ/ｍｍ2,泊松比(即泊松系数)为μ=0.3。
自定义截面梁的ANSYS求解方法
自定义截面梁的ANSYS求解ANSYS有多种梁截面的形状可供选择，但是如果ANSYS里面没有你想要的梁截面该怎么办呢？那就需要自己定义自己喜欢的梁截面。
下面给出命令流，可以很好地定义自己喜欢的梁截面。
先把效果展示给大家，命令流只是用来学习的，放在文章的最后。
我要分析的一个问题是一根L=1m长的矩形截面梁，一端固定一端自由（悬臂梁）。
矩形截面的B×H=0.1m×0.1m，梁上表面受均布荷载q=1000N/m，梁的弹性模量为E= 210GPa。
这个是可以通过理论计算得到梁末端最大挠度：w=qL48EI=1000×148×210×109×112×0.1×0.13=57×10−4≈0.714×10−4（m）截面的定义可以直接在ANSYS中进行，先定义好面，然后用MESH200单元进行网格划分，最后保存为*.rect文件，然后在进行梁单元建模的时候，导入截面参数文件，最后赋予梁该截面形状，进行梁单元网格剖分，施加荷载和边界条件，进行求解。
如下图所示，即定义好的梁截面形状。
可以看到界面的各个系数。
然后建模求解。
静力学求解一个Solve命令就行了。
最后可以画出变形图和应力云图。
应力云图如下。
由变形图可以得到该悬臂梁模型用ANSYS计算出来的最大末端位移是0.724×10−4（m），应力云图可以获知最大拉应力和最大压应力的大小相等为2.72MPa，求解成功。
万冀2017.01.03 于珞珈山（附：ANSYS命令流于后两页）1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40finish/clear/prep7!利用MESH200单元对自定义梁截面进行网格剖分et,1,mesh200!*设置KeyOpt(1)=7即八节点四边形单元（截面划分网格必须这样）keyopt,1,1,7!*定义关键点k,1,0.0,0.0k,2,0.1,0.0k,3,0.1,0.1k,4,0.0,0.1!*连接关键点l,1,2l,2,3l,3,4l,4,1!*由线条构成截面al,1,2,3,4!*映射网格划分mshape,0,2dmshkey,1esize,0.01 !控制网格大小asel,,,,1amesh,all !划分网格!将自定义截面参数输出到文件mysec.sect文件中保存secwrite,mysec,sect,,1finish!**清空以往数据并开始新的计算/clear/prep7et,1,beam188mp,dens,1,7850mp,ex,1,2.1e11mp,prxy,1,0.3!***导入自定义截面sectype,1,beam,meshsecoffset,cent41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73secread,'mysec','sect',,meshsecplot,1,1 !**画出截面形状和网格r,1,k,1,0,0k,2,1,0k,3,0,1 ! 给梁单元定方向的3号关键点l,1,2,10 ! 连接关键点1和2并预划分10份!**LATT,-材料-实常数-单元类型-空-方向关键点号-空-梁截面编号latt,1,1,1,,3,,1esize,0.1 ! 划分单元大小lmesh,1sfbeam,all,1,pres,1000,,,,,,0 ! 施加均布压力(量纲：F/L) nsel,s,loc,x,0 ! 选取x=0处的节点施加固定约束(即模拟悬臂梁) d,all,allfinish!****静力求解/soluantype,0allsel,allsolvefinish/post1!*画变形图pldisp,1!*显示单元形状/eshape,1.0!*画x方向的应力云图plnsol,s,x/view,1,1,1,1/replot。
ANSYS论文：基于ANSYS的铝合金箱型截面悬臂梁模态分析
ANSYS论文：基于ANSYS的铝合金箱型截面悬臂梁模态分析摘要通过ansys对具有复杂约束条件的铝合金箱型截面悬臂梁进行了模态分析，并将结果与高精度的实验结果进行比较分析，验证了所使用的模态分析方法的正确性及可行性，为解决相似问题提供了一种新方法和新思路。
关键词ansys；模态分析；铝合金；悬臂梁；固有频率振动问题广泛存在于航空航天、机械动力、交通运输及军事国防工业等国民经济的各个领域。
模态分析是在振动测量中求解振动物体固有频率的重要方法。
通过模态分析，可以得到振动系统比较精确的固有频率、模态振型和模态刚度，从而为进一步解决振动问题打下重要基础。
但是在解决某些复杂约束情况下的模态分析问题时，由于无法较好地模拟真实的约束情况而使得求解结果误差很大，缺乏可信度。
本文通过使用ansys对一处于复杂约束情况下的实例进行数值模拟，得到了较精确的结果，为解决相似问题提供了新的思路和方法。
1ansys必要数据准备1.1试样类型及相关数据试样是某种型号的铝合金箱型截面梁，试样一端打有两个孔洞，通过螺丝安装在试验台上，使其成为悬臂梁。
试样的安装构造及横截面尺寸如图1所示。
由米尺测得试样的长度l=441.2mm，横截面上各尺寸及壁厚m由游标卡尺测得。
通过电子秤测得试样的质量m=73.29g。
1.2数据处理由试样长度和质量可求得试样的线密度，即ρ=m/l=0.166kg/m。
计算图1中所示试样横截面对x轴的惯性矩ix的值。
试样壁厚存在不均匀性，为计算简便，设横截面上的坐标原点位于外矩形的形心，上下左右四个小矩形的惯性矩分别为i上、i下、i左、i右，由惯性矩计算公式及移轴定理，可得横截面对x轴的惯性矩ix，即：ix=i上+i下+i左+i右=2579.7mm4。
2弹性模量的测量弹性模量是分析材料力学性能的一个极为重要的固有属性。
在使用ansys对试样进行模态分析时，弹性模量e是极其重要的，故下面来测定试样的弹性模量。
2.1测量方案由于没有对应的夹具，所以无法直接在拉伸试验机上进行该铝合金试样的拉伸试验，故给出以下两种测量其弹性模量的方案。
如何在ansys中添加用户所需的梁截面
1----创建截面的几何模型（二维面模型）或者通过第三方软件将autocad的dwf格式转化成通用格式iges(推荐algordos版),然后在ansys中导入。
2----对所有线设置单元份数或者单元最大尺寸(Main Menu>Preprocessor> -Meshing-Size Cntrls>-Lines-Picked Lines或使用MeshTool)。
记住：保证模型中的区格（cell）数目不能超过250个。
3----选择菜单Main Menu>Preprocessor>Sections>-Beam-Write Sec Mesh，弹出一个拾取窗口，单击Pick All拾取包含区格的所有面。
4---ANSYS自动在所有面上创建区格。
在划分网格时，ANSYS可能显示单元形状差的消息，也可以被忽略不显示，但是总能看到一条消息“Unable to mesh area....”。
如果已经完成上述工作，你清除所有面上的单元(Main Menu>Preprocessor>-Meshing-Clear> Areas)，并重复第2、3、4步，即必须重新控制网格密度。
5----如果第4步成功则弹出Write Section Library File对话框，File Name域填入一个未用过SECT文件名，Drives域指定一个截面文件存放驱动器，Directories域指定一个截面文件存放目录，然后单击按钮OK，完成用户截面文件建立。
如果在第3步中给线指定太多单元份数，区格和节点的数目可能超过限制范围，必须清除所有面上的单元2~4步，直到获得合适数量的区格和节点。
一旦完成上述工作，用户可以在以后分析中读取用户网格文件(Main Menu > Preprocessor > Sections > -Beam-Read Sect Mesh)，定义成适当的截面ID号，就与通用梁完全一致。
最新ansys_变截面梁汇总
a n s y s_变截面梁1、选择菜单“Utility Menu>File>Change Title”。
2、输入“Lateral Torsional Buckling Analysis”并单击OK。
3、定义关键点。
选择“Main Menu>Preprocessor> -Modeling- Create>Keypoints> In Active CS”，输入下表所示的关键点号和坐标：4、在关键点1和2之间建立直线。
选择Main Menu>Preprocessor>-Modeling-Create>-Lines-Lines>Straight Line，出现“Create Straight Line picker”窗口，在图形窗口中拾取关键点1和2，然后按“OK”。
5、保存模型。
选择Utility Menu>File>Save As。
将文件名存为“buckle.DB”，并单击OK。
7.6.6 定义单元类型和横截面信息1、选择“Main Menu>Preferences”，单击“Structural”检查框。
按“OK”。
2、选择“Main Menu>Preprocessor>Element Type>Add/Edit/Delete”。
出现“Element Types”对话框。
3、单击“Add...”。
出现“Library of Element Types ”对话框。
4、在左列选择“Structural Beam”。
5、在右列选择“3D finite strain, 3 node 189”以选中 BEAM189。
7、定义梁的矩形截面。
选择“Main Menu>Preprocessor>Sections>-Beam-Common Sects”。
出现“BeamTool ”对话框。
缺省时 ANSYS 将截面号设置为1，将子类型设置为 RECT (在子类型处图示一个矩形)。
Ansys上机例题-变截面梁
ANSYS 上机例题---变截面梁1 问题描述例：如图，一悬臂梁长1L m =，梁截面为实心圆截面,其中D 1200mm =,D2150mm =。
受均布载荷5q kN m =, q 方向向左。
梁材质为钢材，37850kg m ρ-=⋅， 112.110E Pa =⨯， 0.3μ=，用ANSYS 软件求解其应力分布情况（不计重力）。
图1 变截面悬臂梁2 求解步骤概述：前处理1）确定分析标题 2）设置菜单偏好 3）定义单元类型 4）定义材料类型 5）创建关键点6）创建关键点之间的梁，划分网格求解7）设置约束、施加载荷 8）求解后处理9）绘制梁的变形图10）绘制梁的应力图退出A N SYS3详细操作步骤：3.1进入ANSYS10.01 新建工作目录：在D盘新建一个文件夹xuanbiliang，作为ANSYS的工作目录，ANSYS 运行中产生的中间文件和结果文件将保存在这里。
文件名不要含有中文，以免产生乱码。
2开始>程序>ANSYS>ANSYS Product Launcher3 指定工作目录File Management(默认)>Working Directory>Browse选择工作目录为步骤1中新建的文件夹xuanbiliang>run3.2前处理4 确定分析标题：Utility Menu：File>Change Title<1>键入分析标题：“xuanbiliang”<2>OK5 设置单元偏好：Main Menu：Preference<1>仅仅打开“Structural”菜单过滤<2>OK2 输入工作名xuanbiliang 36 定义单元类型：Main Menu ：Preprocessor>Element Type>Add/Edit/Delete <1>Add<2>选择beam189 <3>close7 定义材料特性MainMenu ：Preprocessor>MaterialProps>MaterialModels>Structure>Linear>Elastic>Isotropic<2>在EX 项输入弹性模量：2.1e11 <3>在PRXY 项输入泊松比：0.31231 Add/Edit/Delete25<3>OK1234Main Menu：Preprocessor>Material Props>Material Models > Structure>Density<1>在DENS对画框键入：7850（定义钢材密度）<2>按“OK”完成1238 定义梁的截面Preprocessor>Sections>Beam>Common Sections<1>在BeamTool下拉菜单中选择实心圆截面，截面ID为1，输入对应尺寸R1=0.100，在name栏输入jm1，该截面取名为jm1<2>点击Apply后，定义截面ID为2，输入对应尺寸R2=0.075，在name 栏输入jm2，该截面取名为jm2<3>点“OK”2 取名jm11 Common sections349 建立梁的关键点Main Menu：Preprocessor>Modeling>Create>Key points>In Active CS<1>在关键点号菜单栏键入：1（关键点编号）<2>键入关键点1的X,Y,Z坐标：0,0,0（也可以用空着不填表示0）<3>按“Apply”完成关键点1的定义同理，按下表值建立所有关键点，在输入关键点3的坐标后，按OK键（而不是Apply）。
基于 ANSYS 的变截面异形连续梁桥裂缝和挠度验算
基于 ANSYS 的变截面异形连续梁桥裂缝和挠度验算张四国;闫旭【摘要】Based on the knowledge of FEM,3D Solid model of a special shaped continuous beam bridge was built using ANSYS considering the dead load,vehicle live load and temperature load conditions. Bridge deck stress distribution was explored under single operating condition and the most unfavorable conditions combination. Results were compared with the allowable values. Research shows that the bridge works in good condi-tion and longitudinal cracks,0. 08 mm in width,appeares at some parts of the bridge. The maximum crack width meet the specification requirements.%基于有限元法的相关知识，利用ANSYS软件对某异形连续梁桥进行了实体有限元建模及静力分析，考虑了恒载、车辆活载和温度荷载等工况，探索了桥梁在单一工况及最不利工况组合下桥面板的应力分布情况，并与规范容许值进行了对比和抗裂验算，结果表明，桥面受力状况良好，部分区域出现宽度为0.08 mm的纵向裂缝，但最大裂缝宽度满足规范要求。
【期刊名称】《山西建筑》【年(卷),期】2015(000)019【总页数】3页(P151-152,153)【关键词】有限元;变截面;连续梁;结构分析【作者】张四国;闫旭【作者单位】天津市市政工程设计研究院，天津 300051;天津市市政工程设计研究院，天津 300051【正文语种】中文【中图分类】U441公路桥涵的设计要满足技术先进、安全可靠、耐久适用和经济合理等要求[1]，其中安全性是设计的第一原则。
ansys变截面建模资料
梁壳体连接情况的讨论1 按“杆梁壳体”的顺序，只要后一种单元的自由度完全包容了前一种单元的自由度，则有公用结点即可，不需要约束方程。
例如：杆与梁、壳、体有公用结点即可，不需要写约束方程；梁与壳有公用结点即可，不需要写约束方程；梁与体则要同位置的不同结点，需要耦合自由度和约束方程；壳与体则要同位置的不同结点，需要耦合自由度和约束方程；2 壳梁自由度数目相同，自由度也相同，尽管壳的rotz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。
3 尽管可以采用耦合自由度和约束方程，但建议尽量不同时采用多种单元于一个结构中，除非你对结果的正确性有十足的把握。
4 当然，采用约束方程可能存在应力集中点，不必在意此点的应力。
一般有两个办法，其一是通过施加约束方程，即两单元某个自由度和板单元某个自由度之间的关系；其二是通过在二者连接部分施加刚域解决。
壳体用SHELL53单元建模，上下有加强的梁，用BEAM3建模，它们之间的在共用节点处一一耦合，通过CPINT命令实现。
壳单元的Z轴转角为零，而梁单元有完整的6个自由度，只需要耦合掉这个自由度应该就可以了。
复合材料的话，Hypermesh里面有专门的定义，可以设置铺层的厚度和方向~当然，你这样连接相当于仅在节点处壳体与加强筋连接，而有限元的原理之一是把节点几何尺寸设为不影响分析精度之小，那此处的应力集中就太大了。
解决办法：1。
壳体和梁连接处共线后mesh2。
用耦合方程连接壳体和梁1.!**********************2.!耦合自由度的尝试代码3.!范雨 2009-4-84.!国际单位制5.!**********************6.FINISH7./CLEAR8./TITLE,TEST APDL FOR COUPLING9.10./PREP711.!**********************12.!参数表13.!**********************14.*SET,R_SHELL,0.5!圆柱半径15.*SET,L_SHELL,2!圆柱长度16.*SET,LN_MESH,32!轴向分网17.*SET,RN_MESH,8!周向分网18.*SET,TH,0.01!壳单元厚度19.*SET,MEX,2.06E11!弹性模量20.*SET,MPRXY,0.3!弹性模量21.*SET,MDENS,7850!密度22.*SET,H_BEAM,0.1!梁界面参数：高23.*SET,B_BEAM,0.1!梁界面参数：宽24.*set,IZ_BEAM,(B_BEAM)*(H_BEAM**3)/12!绕Z的惯性矩25.!**********************26.!相关设置27.!**********************28.!设置材料29.MP,EX,1,MEX30.MP,PRXY,1,MPRXY31.MP,DENS,1,MENDS32.!设置壳单元33.ET,1,SHELL6334.R,1,TH35.!设置BEAM单元36.ET,2,BEAM337.R,2,H_BEAM*B_BEAM,IZ_BEAM,H_BEAM38.!显示SHELL厚度39./ESHAPE,140.!以箭头显示压力41./PSF,PRES,NORM,2,142.!**********************43.!柱面->几何->分网44.!**********************45.M,146.TYPE,147.REAL,148.!建立柱面49.K,1,0,0,050.K,2,0,L_SHELL,051.L,1,252.CIRCLE,1,R_SHELL,253.ADRAG,ALL,,,,,,154.LDELE,155.!设置分网数目56.ALLSEL,ALL57.LSEL,U,LOC,Y,058.LSEL,U,LOC,Y,L_SHELL59.LESIZE,ALL,,,LN_MESH60.ALLSEL,ALL61.LSEL,R,LOC,Y,062.LSEL,A,LOC,Y,L_SHELL63.LESIZE,ALL,,,RN_MESH64.!对圆柱侧面映射分网65.ALLSEL,ALL66.ASEL,U,LOC,Y,067.ASEL,U,LOC,Y,L_SHELL68.AMESH,ALL69.LESIZE,ALL,,,LN_MESH70.!**********************71.!梁->几何->分网72.!**********************73.!建立辅助直线74.K,100,R_SHELL,0,075.K,101,R_SHELL,L_SHELL,076.L,100,10177.K,102,-1*R_SHELL,0,078.K,103,-1*R_SHELL,L_SHELL,079.L,102,10380.!设置分网参数81.TYPE,282.REAL,283.ALLSEL,ALL84.LSLA,U85.LESIZE,ALL,,,LN_MESH86.!分网87.LMESH,ALL88.!**********************89.!自动耦合自由度90.!**********************91.ALLSEL,ALL92.CPINTF,ALL,0.0001,93.!**********************94.!两端固支：全约束95.!**********************96.NSEL,S,LOC,Y,097.NSEL,R,LOC,X,-R_SHELL98.D,ALL,ALL,099.NSEL,S,LOC,Y,L_SHELL100.NSEL,R,LOC,X,-R_SHELL 101.D,ALL,ALL,0102.!**********************103.!加载BEAM上的压力载荷104.!**********************105.ESEL,S,TYPE,,2106.ESEL,R,ELEM,,1025,1056107.SFBEAM,ALL,,PRES,-5E5108.ALLSEL,ALL109.!求解110./SOLU111.SOLVE112.!后处理113./POST1114.PLNSOL,S,EQV,0,,复制代码壳体用SHELL53单元建模，上下有加强的梁，用BEAM3建模，它们之间的在共用节点处一一耦合，通过CPINT命令实现。
ANSYS第8章 两单元的常用操作
ANSYS 入门教程 (10) –梁单元的常用操作第8章两单元的常用操作一、梁截面和实常数梁单元是用有限元法进行梁柱体系分析时最常采用的单元类型之一。
梁单元是个线单元，需要另外定义其截面特性。
在ANSYS 中，梁单元的截面特性通常通过实常数来定义，但是，对于 BEAM44、BEAM188 和 BEAM189 也可以通过截面 (Section) 来定义梁单元的截面形状，或者定义Section 后再由软件计算相应的截面特性，再作为实常数输入。
下面叙述定义Section 的具体操作。
1 定义标准截面-菜单：Main Menu > Preprocessor > Section > Beam > Common Sections > 出现标准截面定义对话框：在对话框中：(1) 给截面命名、选择标准截面形状；(2) 选择节点偏移方法；(3) 输入截面数据和选择截面网格粗细程度：如果需要，可以修改截面的横向剪切刚度和添加分布质量(单位是质量/ 长度)：然后可以显示截面的形状和截面上的网格划分情况，以及有关的截面参数。
截面上的网格是软件自动划分的，用于计算截面的几何参数。
首先选择要显示的截面和是否显示截面网格：然后显示所选择的截面形状和几何参数：也可以对截面列表：列表结果：右图显示了列表的部分内容：对于不能使用截面的梁单元，可以将列表数据作为实常数输入到相应的梁单元类型定义中。
需要时，可以删除已定义的截面：2. 定义用户截面–需要先绘制相应的截面形状，然后将相应的AREA 定义为截面。
操作过程如下：(1) 绘制截面形状 - AREAs(2) 对各线段设置网格尺寸：对本例设置所有线段的单元长度为 4设置网格尺寸之后：(3) 然后定义用户截面：ANSYS 并不马上生成新的截面，而是将新截面的数据写到一个“截面”文件中，如图，点击 OK 后将在工作目录中生成截面文件2-C18.sec：(4) 为了使用新定义的截面，需要读入该截面文件：然后可以绘制该截面：3. 定义变截面梁 (锥形梁)：定义锥形梁之前，必须先定义锥形梁两端的截面，比如我们先定义两个矩形截面beam-4 和 beam-5 如下：然后定义锥形 (Taper) 截面：采用 By XYZ Location 方法，定义截面 7：使用截面7 对线段划分梁单元：二、为划分网格设置单元属性在对实体划分网格之前，需要设置单元属性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
利用ANSYS生成变截面箱梁
Beam188/189 支持自定义的变截面
1、首先建立变截面箱梁截面，并保存截面。
必须保证每个截面的关键点号相同，而且为保证生成的准确性，应尽量的使得关键点有足够多的数目。
2、要对截面进行面积分块，并指定各线段的段数，这样才能做出规整的箱梁截面网格划分，也就保证了变截面箱梁桥各截面的网格模式相同，建议对于变化急剧的两边截面应使得线足够的短。
并不是划分越密越可能成功，而是线越短越可能成功
3。
通过梁节段两端对应的截面建立taper截面，定义好桥的线形后，指定每段线对应的taper截面。
4、可以用slist命令查看生成截面的性质，加深对secread 命令的理解
以下例进行说明：
Ⅰ-Ⅰ截面图
全梁1/2图示Ⅰ-Ⅰ截面命令流
finish
/clear
/prep7
k,1
k,2,2750
k,3,3350,2400
k,4,5450,2750
k,5,6700,2850
k,6,6700,3050
k,8,0,280
k,9,2359,280 k,10,2886,2400 k,11,1836,2750 k,12,0,2750 l,1,2
l,2,3
l,3,4
l,4,5
l,5,6
l,6,7
l,8,9
l,9,10
l,10,11
l,11,12
l,1,8
l,7,12
al,all arsym,x,all aadd,1,2 nummrg,all
adele,1
l,2,9
l,13,20
l,14,19
l,10,3
l,1,8
al,22,17,25,19 al,25,1,21,7 lcomb,15,20 wpoffs,-5450,2750 wprota,,,90 lsbw,18 wpcsys,-1 wpoffs,-3350,2400 wprota,,,90 lsbw,20 wpcsys,-1 wpoffs,-2886,2400 wprota,,,90 lsbw,18 wpcsys,-1
l,14,21
l,19,22
l,7,12
al,14,13,18,26 al,18,12,29,27 al,29,23,30,28 al,30,15,31,20 lcomb,9,10 wpoffs,5450,2750 wprota,,,90 lsbw,6 wpcsys,-1 wpoffs,3350,2400 wprota,,,90 lsbw,32 wpcsys,-1 wpoffs,2886,2400 wprota,,,90 lsbw,6 wpcsys,-1
l,4,11
l,10,24
al,5,4,6,10 al,3,6,33,35 al,24,35,34,36 al,9,36,32,31 ldiv,16,,,4 ldiv,11,,,4 ldiv,2,,,4 ldiv,8,,,4
l,30,27
l,29,26
l,28,25
l,36,33
l,35,32
l,34,31
al,23,42,49,39 al,49,41,50,38 al,50,40,51,37 al,51,11,22,16 al,24,48,52,45 al,52,47,53,44
al,53,46,54,43 al,54,8,21,2 nummrg,all numcmp,all et,1,82 lesize,1,,,4 lesize,17,,,4 lesize,19,,,4 lesize,7,,,4 lesize,22,,,1 lesize,21,,,1 lesize,25,,,1 mshape,1,2d mshkey,2 amesh,1 amesh,2 lesize,23,,,1 lesize,42,,,1 lesize,39,,,1 lesize,49,,,1 lesize,41,,,1 lesize,50,,,1
lesize,38,,,1 lesize,40,,,1 lesize,37,,,1 lesize,51,,,1 lesize,11,,,1 lesize,16,,,1 lesize,22,,,1 lesize,24,,,1 lesize,45,,,1 lesize,48,,,1 lesize,52,,,1 lesize,44,,,1 lesize,47,,,1 lesize,53,,,1 lesize,43,,,1 lesize,46,,,1 lesize,54,,,1 lesize,2,,,1 lesize,8,,,1 lesize,21,,,1 mshape,0,2d mshkey,1
amesh,12 amesh,13 amesh,14 amesh,15 amesh,16 amesh,17 amesh,18 lesize,13,,,2 lesize,26,,,2 lesize,14,,,1 lesize,18,,,1 mshape,1,2d mshkey,2 amesh,3 lesize,12,,,3 lesize,27,,,3 lesize,18,,,1 lesize,29,,,1 mshape,1,2d mshkey,2 amesh,4
lesize,29,,,1 lesize,30,,,1 lesize,28,,,1 mshape,0,2d mshkey,2 amesh,5 lesize,30,,,1 lesize,20,,,4 lesize,31,,,1 lesize,15,,,4 mshape,1,2d mshkey,2 amesh,6 lesize,31,,,1 lesize,9,,,4 lesize,36,,,1 lesize,32,,,4 mshape,1,2d mshkey,2 amesh,10 lesize,35,,,1
lesize,34,,,1 lesize,24,,,1 lesize,36,,,1 mshape,0,2d mshkey,2 amesh,9 lesize,6,,,1 lesize,35,,,1 lesize,3,,,3 lesize,33,,,3 mshape,1,2d mshkey,2 amesh,8 lesize,6,,,1 lesize,5,,,1 lesize,10,,,2 lesize,4,,,2 mshape,1,2d mshkey,2 amesh,7 secwrite,mybox1 主命令流
finish
/clear
/prep7
k,5000,7800
k,5001,8851
k,5002,9899
k,5003,10800
k,5004,12300
k,5005,0
l,5000,5001
l,5001,5002
l,5002,5003
l,5003,5004
l,5000,5005
k,100,8300,4000
k,101,9500,4000
k,102,10000,4000
et,1,beam189
mp,ex,1,3.0e10
mp,prxy,1,0.167 sectype,1,beam,mesh secread,mybox1,,,mesh
sectype,2,beam,mesh secread,mybox2,,,mesh sectype,3,taper secdata,1,7800 secdata,2,8851
et,2,beam189
mp,ex,2,3.0e10
mp,prxy,2,0.167 sectype,4,beam,mesh secread,mybox2,,,mesh sectype,5,beam,mesh secread,mybox3,,,mesh sectype,6,taper secdata,4,8851 secdata,5,9899
et,3,beam189
mp,ex,3,3.0e10
mp,prxy,3,0.167 sectype,7,beam,mesh secread,mybox3,,,mesh sectype,8,beam,mesh secread,mybox4,,,mesh
sectype,9,taper secdata,7,9899 secdata,8,10800
et,4,beam189
mp,ex,4,3.0e10
mp,prxy,4,0.167 sectype,10,beam,mesh secread,mybox4,,,mesh et,5,beam189
mp,ex,5,3.0e10
mp,prxy,5,0.167 sectype,11,beam,mesh secread,mybox1,,,mesh lesize,1,,,2
lesize,2,,,2
lesize,3,,,2
lesize,4,,,2
lesize,5,,,8
allsel
lsel,s,line,,1
latt,1,,1,,100,,3
allsel
lsel,s,line,,2 latt,2,,2,,101,,6 allsel
lsel,s,line,,3 latt,3,,3,,102,,9 allsel
lsel,s,line,,4 latt,4,,4,,102,,10 lsel,s,line,,5 latt,5,,5,,102,,11 lmesh,1,5,1 allsel lmesh,1,5,1
/eshape,1 eplot。