!!连续梁桥
/prep7
et,1,4 !!!!定义梁单元
et,2,21 !!!!定义mass21单元
!!定义粱材料!!泊松比!!密度
mp,ex,2,3.45e10 !!直线段梁材料和1M段梁材料mp,nuxy,2,0.2
mp,dens,2,3302.153125
mp,ex,3,3.45e10
mp,nuxy,3,0.2
mp,dens,3,3301.658695
mp,ex,4,3.45e10
mp,nuxy,4,0.2
mp,dens,4,3299.906778
mp,ex,5,3.45e10
mp,nuxy,5,0.2
mp,dens,5,3298.327219
mp,ex,6,3.45e10
mp,nuxy,6,0.2
mp,dens,6,3292.351605
mp,ex,7,3.45e10
mp,nuxy,7,0.2
mp,dens,7,3284.137255
mp,ex,8,3.45e10
mp,nuxy,8,0.2
mp,dens,8,3271.802136
mp,ex,9,3.45e10
mp,nuxy,9,0.2
mp,dens,9,3260.41903
mp,ex,10,3.45e10
mp,nuxy,10,0.2
mp,dens,10,3248.193657
mp,nuxy,11,0.2
mp,dens,11,3235.117644
mp,ex,12,3.45e10
mp,nuxy,12,0.2
mp,dens,12,3221.585664
mp,ex,13,3.45e10
mp,nuxy,13,0.2
mp,dens,13,3208.826871
mp,ex,14,3.45e10
mp,nuxy,14,0.2
mp,dens,14,3194.279207
mp,ex,15,3.45e10
mp,nuxy,15,0.2
mp,dens,15,3179.924673
mp,ex,16,3.45e10
mp,nuxy,16,0.2
mp,dens,16,3166.445716
mp,ex,17,3.45e10
mp,nuxy,17,0.2
mp,dens,17,3152.555731
mp,ex,18,3.45e10
mp,nuxy,18,0.2
mp,dens,18,3138.312105
mp,ex,19,3.45e10
mp,nuxy,19,0.2
mp,dens,19,3124.795334
mp,ex,20,3.45e10
mp,nuxy,20,0.2
mp,dens,20,3110.7135
mp,ex,21,3.45e10
mp,nuxy,21,0.2
mp,dens,21,3097.080875
mp,nuxy,22,0.2
mp,dens,22,3083.186268
mp,ex,23,3.45e10
mp,nuxy,23,0.2
mp,dens,23,3068.968824
mp,ex,24,3.45e10
mp,nuxy,24,0.2
mp,dens,24,3055.612436
mp,ex,25,3.45e10
mp,nuxy,25,0.2
mp,dens,25,3045.857147
mp,ex,26,3.45e10
mp,nuxy,26,0.2
mp,dens,26,3035.174287
mp,ex,27,3.45e10
mp,nuxy,27,0.2
mp,dens,27,3026.696551
mp,ex,28,3.45e10
mp,nuxy,28,0.2
mp,dens,28,3015.795365
mp,ex,29,3.45e10
mp,nuxy,29,0.2
mp,dens,29,3007.710181
mp,ex,30,3.45e10
mp,nuxy,30,0.2
mp,dens,30,3000.513837
mp,ex,31,3.45e10
mp,nuxy,31,0.2
mp,dens,31,2978.611375
mp,ex,32,3.45e10
mp,nuxy,32,0.2
mp,dens,32,2958.618861
mp,nuxy,33,0.2
mp,dens,33,2937.888072
mp,ex,34,3.45e10
mp,nuxy,34,0.2
mp,dens,34,2919.475751
mp,ex,35,3.45e10
mp,nuxy,35,0.2
mp,dens,35,2903.359983 !!6700处
mp,ex,36,3.45e10
mp,nuxy,36,0.2
mp,dens,36,3302.153125 !!可以不用
mp,ex,37,3.45e10
mp,nuxy,37,0.2
mp,dens,37,3302.153125
mp,ex,38,3.45e10 !!-700处梁
mp,nuxy,38,0.2
mp,dens,38,3180.578901
mp,ex,39,3.45e10 !!边支点横隔板
mp,nuxy,39,0.2
mp,dens,39,2868.674818
mp,ex,40,3.45e10 !!合拢段横隔板
mp,nuxy,40,0.2
mp,dens,40,2868.674818
mp,ex,41,3.45e10 !!中支点横隔板
mp,nuxy,41,0.2
mp,dens,41,2757.470588
mp,ex,1,3.25e10 !!!!定义墩材料属性
mp,nuxy,1,0.2
mp,dens,1,2650
!!定义实常数编号，面积，IYY，IZZ，宽，高，，RMORE，，抗扭惯距!!主梁截面
r,1,11.851,28.52,215.151,16.95,4,, !!直线段rmore,,56.8972
r,2,11.862,28.57,215.458,16.95,4.001,, !!截面100处rmore,,56.905
r,3,11.871,28.91,215.622,16.95,4.009,,
rmore,,57.345
r,4,11.903,29.03,215.78,16.95,4.025,,
rmore,,57.938
r,5,11.932,29.626,215.88,16.95,4.049,,
rmore,,58.502
r,6,12.043,30.402,217.093,16.95,4.081,, rmore,,59.768
r,7,12.199,31.407,218.899,16.95,4.121,, rmore,,61.445
r,8,12.441,32.716,221.806,16.95,4.169,, rmore,,63.652
r,9,12.673,34.21,224.418,16.95,4.226,,
rmore,,66.058
r,10,12.932,35.94,227.279,16.95,4.290,, rmore,,68.801
r,11,13.221,37.941,230.512,16.95,4.362,, rmore,,71.948
r,12,13.534,40.233,233.872,16.95,4.442,, rmore,,75.432
r,13,13.843,42.747,237.133,16.95,4.530,, rmore,,79.173
r,14,14.213,45.688,241.054,16.95,4.627,, rmore,,83.484
r,15,14.598,48.964,245.004,16.95,4.731,, rmore,,88.155
r,16,14.979,52.547,248.817,16.95,4.843,, rmore,,93.111
r,17,15.393,56.562,252.927,16.95,4.963,, rmore,,98.528
r,18,15.842,61.069,257.347,16.95,5.092,, rmore,,104.47
r,19,16.293,65.992,261.673,16.95,5.228,, rmore,,110.742
r,20,16.791,71.484,266.426,16.95,5.372,, rmore,,117.648
r,21,17.303,77.568,271.225,16.95,5.525,, rmore,,124.861
r,22,17.858,84.282,276.363,16.95,5.685,, rmore,,132.631
r,23,18.464,91.768,281.923,16.95,5.854,, rmore,,141.068
r,24,19.072,99.936,287.364,16.95,6.03,, rmore,,149.875
r,25,19.524,108.482,291.159,16.95,6.214,, rmore,,158.46
r,26,20.084,118.046,295.597,16.95,6.407,, rmore,,167.857
r,27,20.536,128.202,298.9,16.95,6.607,, rmore,,177.137
r,28,21.148,139.669,303.86,16.95,6.816,, rmore,,187.495
r,29,21.626,151.737,307.186,16.95,7.032,, rmore,,197.507
r,30,22.07,164.734,310.057,16.95,7.257,, rmore,,207.681
r,31,23.541,182.225,323.433,16.95,7.490,, rmore,,223.527
r,32,25.066,201.932,336.467,16.95,7.730,, rmore,,241.046
r,33,26.871,223.700,351.448,16.95,7.979,, rmore,,257.765
r,34,28.707,247.928,365.663,16.95,8.235,, rmore,,275.922
r,35,30.533,274.36,378.934,16.95,8.500,, rmore,,294.461
r,46,14.58,35.599,238.702,16.95,4,, rmore,,238.702
!!横隔板截面
!!边支点
r,36,35.376,54.393,339.84,16.95,4.00,, rmore,,128.747
!!合拢段
r,37,35.376,53.405,339.84,16.95,4.00,, rmore,,127.015
!!中支点
r,38,71.981,454.871,543.871,16.95,8.5,, rmore,,684.565
!!主墩截面
!!截面1（实心）
r,39,39.932,82.876,212.138,8,5,, rmore,,203.709
!!截面（1/2空心）
r,40,23.863,72.571,168.024,8,5,, rmore,,168.491
!!截面3(空心）
r,41,14.236,52.366,114.529,8,5,, rmore,,116.391
!!边墩截面
!!截面1（实心1）
r,42,36.113,66.217,177.703,7.7,4.7,, rmore,,165.119
!!截面2（空心1/2）
r,43,20.753,35.743,99.711,7,4,, rmore,,88.486
!!截面3（空心）
r,44,9.103,16.082,45.911,6.5,3.5,, rmore,,39.39
!!截面4（实心2）
r,45,22.673,23.005,79.361,6.5,3.5,, rmore,,61.527
!!截面5（1/2）空心2
r,49,16,21.958,68.708,6.5,3.5,, rmore,,55.682
!!建立第一跨梁节点
xl=atan(1.029/80)
n,1,0,
n,2,2.25,2.25*xl
n,3,6,6*xl
n,4,9,9*xl
n,5,11,11*xl
n,6,13,13*xl-0.0045
n,7,15,15*xl-0.0135
n,8,17,17*xl-0.023
n,9,19,19*xl-0.038
n,10,21,21*xl-0.062
n,11,23,23*xl-0.092
n,12,25,25*xl-0.1275
n,13,27,27*xl-0.167
n,14,29,29*xl-0.2115
n,15,31,31*xl-0.261
n,16,33,33*xl-0.314
n,17,35,35*xl-0.3715
n,18,37,37*xl-0.435
n,19,39,39*xl-0.502
n,20,41,41*xl-0.5725
n,21,43,43*xl-0.648
n,22,45,45*xl-0.728
n,23,47,47*xl-0.8115
n,24,49,49*xl-0.8995
n,25,51,51*xl-0.9925
n,26,53,53*xl-1.0895
n,27,55,55*xl-1.1905
n,28,57,57*xl-1.325
n,29,59,59*xl-1.433
n,30,61,61*xl-1.5155
n,31,63,63*xl-1.6325
n,32,65,65*xl-1.7535
n,33,67,67*xl-1.879
n,34,69,69*xl-2.007
n,35,71,71*xl-2.1325
n,36,73,73*xl-2.2635
n,37,75,75*xl-2.399
n,38,77,77*xl-2.533
n,39,79,79*xl-2.631
xl2=atan(0.7/140)
local,11,0,80,80*xl-2.631,0,xl2
nsym,x,35,5,39,1 !!复制粱结点从5到39结点编号增加35 csys,0
n,75,149.6,0.15852
local,12,0,150,0.15852,0,xl2
nsym,x,36,40,75,1
local,13,0,220,-1.6027,0
nsym,x,111,1,111,1
csys,0
n,500,80,80*xl-2.631
n,501,220,-1.6027
n,502,360,-1.6148
csys,0 !!返回普通坐标
!!建立墩结点
!!PM112边墩截面
n,300,0,-2.165 !!(实心）2米长
n,301,0,-4.365
n,302,0,-9.365 !!（1/2空心）5M长
n,303,0,-13.365
n,304,0,-17.365
n,305,0,-21.365
n,306,0,-24.93
n,307,0,-28.93
n,308,0,-30.63
!!PM113主敦截面
n,309,80,0.5-6.665-0.8
n,310,80,0.5-3.5-6.665 !!(实心)3.5m
n,311,80,0.5-8.5-6.665 !!(1/2空心）5m n,312,80,0.5-12.5-6.665 !!(空心）4M
n,313,80,0.5-16.5-6.665 !!
n,314,80,0.5-18.594-6.665
n,315,80,0.5-22.594-6.665
n,316,80,0.5-24.594-6.665+0.4!(实心）2m
!!pm114主墩截面
n,350,220,0.5-6.665-0.8 !
n,351,220,0.5-3.5-6.665 !!(1/2空心)5m n,352,220,0.5-8.5-6.665 !!(空心）5m
n,353,220,0.5-12.5-6.665 !!(空心）6.064M n,354,220,0.5-16.5-6.665 !!(1/2空心）4M n,355,220,0.5-19.294-6.665 !!(实心）2M n,356,220,0.5-23.294-6.665
n,357,220,0.5-25.294-6.665+0.4
local,13,0,220,-1.6027,0
nsym,x,17,300,316,1
csys,0
!!生成边跨1
mat,39
real,36
e,1,2
mat,38
real,46
e,2,3
mat,2
real,1
e,3,4
e,4,5
mat,2
real,2
e,5,6
*do,i,3,35,1
mat,i
real,i
e,3+i,4+i
*enddo
*do,i,3,35,1
mat,i
real,i
e,38+i,39+i
*enddo
*do,i,3,35,1
mat,i
real,i
e,74+i,75+i
*enddo
*do,i,3,35,1
real,i
e,114+i,115+i *enddo
*do,i,3,35,1 mat,i
real,i
e,149+i,150+i *enddo
*do,i,3,35,1 mat,i
real,i
e,185+i,186+i *enddo
mat,41 real,38
e,39,500
mat,41 real,38
e,500,74
mat,2
real,2
e,41,40
mat,2
real,1
e,40,75
mat,40 real,37
e,75,111
mat,2
real,1
e,111,76
mat,2
real,2
mat,41 real,38 e,110,501
mat,41 real,38 e,501,221
mat,2 real,2 e,188,187
mat,2 real,1 e,187,222
mat,41 real,37 e,222,186
mat,2 real,1 e,186,151
mat,2 real,2 e,151,152
mat,41 real,38 e,185,502
mat,41 real,38 e,502,150
mat,39 real,36 e,112,113
real,46 e,113,114
mat,2 real,1
e,114,115 e,115,116
mat,2 real,2
e,116,117
!!建立桥墩!!边墩1 mat,1 real,42 e,300,301
mat,1 real,43 e,301,302
mat,1 real,44 e,302,303
mat,1 real,44 e,303,304
mat,1 real,44 e,304,305
mat,1 real,44 e,305,306
mat,1 real,49 e,306,307
mat,1 real,45 e,307,308
!!主墩1 mat,1 real,39 e,309,310
mat,1 real,40 e,310,311
mat,1 real,41 e,311,312
mat,1 real,41 e,312,313
mat,1 real,41 e,313,314
mat,1 real,40 e,314,315
mat,1 real,39 e,315,316
!!主墩2 mat,1 real,39 e,350,351
mat,1 real,40 e,351,352
mat,1
real,41 e,352,353
mat,1 real,41 e,353,354
mat,1 real,41 e,354,355
mat,1 real,40 e,355,356
mat,1 real,39 e,356,357
!!主墩3 mat,1 real,39 e,326,327
mat,1 real,40 e,327,328
mat,1 real,41 e,328,329
mat,1 real,41 e,329,330
mat,1 real,41 e,330,331
mat,1 real,40 e,331,332
mat,1
real,39
e,332,333
!!边墩2
mat,1
real,42
e,317,318
mat,1
real,43
e,318,319
mat,1
real,44
e,319,320
mat,1
real,44
e,320,321
mat,1
real,44
e,321,322
mat,1
real,44
e,322,323
mat,1
real,49
e,323,324
mat,1
real,45
e,324,325
!!添加质量惯距
!!两米段质量惯距
r,100,,,,1293348.4 !!100米处和直线段的质量惯距r,101,,,,1296019.6
r,102,,,,1297493
r,103,,,,1301181.8
r,104,,,,1311723
r,105,,,,1326622
r,106,,,,1348967
r,107,,,,1370728
r,108,,,,1395061
r,109,,,,1422801
r,110,,,,1452757
r,111,,,,1483258
r,112,,,,1519733
r,113,,,,1558030
r,114,,,,1597229
r,115,,,,1640292
r,116,,,,1687605
r,117,,,,1736625
r,118,,,,1790923
r,119,,,,1848602
r,120,,,,1911418
r,121,,,,1980562
r,122,,,,2052690
r,123,,,,2118097
r,124,,,,2192307
r,125,,,,2263641
r,126,,,,2350704
r,127,,,,2432292
r,128,,,,2516392
r,129,,,,2679987
r,130,,,,2853514
r,131,,,,3048284
r,132,,,,3252032
r,133,,,,3462506
!!边跨直线3米段
r,148,,,,1293348.4*3/2
!!-700处3.75米段的质量惯距r,134,,,,2725866.2
!!边支点横隔
r,135,,,,2246142.5
!!合龙段横隔
r,136,,,,833679/2
!!中支点横隔
r,137,,,,5292760
!!添加边墩质量惯距从上到下
!!变截面实心2米段
r,138,,,,1163498
!!1/2空心5米段
r,139,,,,1794766
!!空心段4米
r,140,,,,657211
!!1/2空心4米段
r,141,,,,1085079
!!实心2米段
r,142,,,,1064567
!!添加主墩质量惯性距从上到下
!!实心3米段
r,143,,,,2345361
!!1/2空心5米段
r,144,,,,2869095
!!空心4米段
r,145,,,,1769087
!!1/2空心4米段
r,146,,,,2869095*4/5
!!实心2米段
r,147,,,,2345361*2/3
!!在结点处加入质量惯性距
!!边跨合龙段
!左边边跨
type,2
real,135
e,2
type,2
real,134
e,3
type,2
real,148
e,4
type,2
real,100
e,5
!!右边边跨
type,2
real,135
e,113
type,2
real,134
e,114
type,2
real,148
e,115
type,2
real,100
e,116
!!对第一个主跨的循环
*do,i,6,39,1
type,2
real,94+i
e,i
*enddo
!!对第二个主跨的循环*do,i,41,74,1
type,2
real,59+i
e,i
*enddo
!!对第3个循环
*do,i,77,110,1
type,2
real,33+i
e,i
*enddo
!!对第4个循环
*do,i,188,221,1 type,2
real,(-88)+i
e,i
*enddo
!!对第5个循环
*do,i,152,185,1 type,2
real,i-52
e,i
*enddo
!!对第6个循环
*do,i,117,150,1 type,2
real,(-17)+i
e,i
*enddo
!!添加质量惯距到桥墩
!!左边墩
type,2
real,138
e,301
type,2
real,139
. !!连续梁桥 /prep7 et,1,4 !!!!定义梁单元 et,2,21 !!!!定义mass21单元 !!定义粱材料!!泊松比!!密度 mp,ex,2,3.45e10 !!直线段梁材料和1M段梁材料mp,nuxy,2,0.2 mp,dens,2,3302.153125 mp,ex,3,3.45e10 mp,nuxy,3,0.2 mp,dens,3,3301.658695 mp,ex,4,3.45e10 mp,nuxy,4,0.2 mp,dens,4,3299.906778 mp,ex,5,3.45e10 mp,nuxy,5,0.2 mp,dens,5,3298.327219 mp,ex,6,3.45e10 mp,nuxy,6,0.2
. mp,dens,6,3292.351605 mp,ex,7,3.45e10 mp,nuxy,7,0.2 mp,dens,7,3284.137255 mp,ex,8,3.45e10 mp,nuxy,8,0.2 mp,dens,8,3271.802136 mp,ex,9,3.45e10 mp,nuxy,9,0.2 mp,dens,9,3260.41903 mp,ex,10,3.45e10 mp,nuxy,10,0.2 mp,dens,10,3248.193657 mp,ex,11,3.45e10 mp,nuxy,11,0.2 mp,dens,11,3235.117644 mp,ex,12,3.45e10 mp,nuxy,12,0.2 mp,dens,12,3221.585664
. mp,ex,13,3.45e10 mp,nuxy,13,0.2 mp,dens,13,3208.826871 mp,ex,14,3.45e10 mp,nuxy,14,0.2 mp,dens,14,3194.279207 mp,ex,15,3.45e10 mp,nuxy,15,0.2 mp,dens,15,3179.924673 mp,ex,16,3.45e10 mp,nuxy,16,0.2 mp,dens,16,3166.445716 mp,ex,17,3.45e10 mp,nuxy,17,0.2 mp,dens,17,3152.555731 mp,ex,18,3.45e10 mp,nuxy,18,0.2 mp,dens,18,3138.312105 mp,ex,19,3.45e10
Finish(退出四大模块，回到BEGIN层) /clear (清空内存，开始新的计算) 1．定义参数、数组，并赋值. 2．/prep7(进入前处理) 定义几何图形：关键点、线、面、体 定义几个所关心的节点，以备后处理时调用节点号。 设材料线弹性、非线性特性 设置单元类型及相应KEYOPT 设置实常数 设置网格划分，划分网格 根据需要耦合某些节点自由度 定义单元表 3．/solu 加边界条件 设置求解选项 定义载荷步 求解载荷步 4./post1（通用后处理） 5./post26 （时间历程后处理） 6.PLOTCONTROL菜单命令 7.参数化设计语言 8.理论手册 Finish(退出四大模块，回到BEGIN层) /clear (清空内存，开始新的计算) 1.定义参数、数组，并赋值. dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组 par: 数组名 type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符） table imax,jmax, kmax 各维的最大下标号 var1,var2,var3 各维变量名，缺省为row,column,plane(当type为table时) 2./prep7(进入前处理) 2.1 设置单元类型及相应KEYOPT ET, itype, ename, kop1……kop6, inopr 设定当前单元类型 Itype：单元号 Ename：单元名设置实常数 Keyopt, itype, knum, value itype: 已定义的单元类型号 knum: 单元的关键字号
1 旋转摩擦 (1) 2. 电磁三d命令流实例（论坛看到） (11) 3. 帮助感应加热例子induction heating of a solid cylinder billet (15) 4. 感应加热温度场的数值模拟（论文）inducheat30命令流 (19) 5. 如何施加恒定的角速度？Simwe仿真论坛 (24) 6. 旋转一个已经生成好的物体 (27) 7. 产生这样的磁力线 (28) 8. 旋转摩擦生热简单例子（二维旋转） (32) 8.1. 原版 (32) 8.2. 部分gui操作 (35) 9. VM229 Input Listing (39) 10 轴承---耦合+接触分析 (47) 11. 板的冲压仿真 (52) 1 旋转摩擦 FINISH /FILNAME,Exercise24 !定义隐式热分析文件名 /PREP7 !进入前处理器 ET,1,SOLID5 !选择单元类型 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7800 !定义材料1的密度 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,460 !定义材料1的比热 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,KXX,1,,66.6 !定义材料1的热传导系数 MPTEMP,,,,,,,, MPTEMP,1,0 UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度 MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,206e9 !定义材料1的弹性模量 MPDATA,PRXY,1,,0.3 !定义材料1的泊松比 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,8900 !定义材料2的密度 MPTEMP,,,,,,,, MPTEMP,1,0
目录 一、工程背景 (1) 二、工程模型 (1) 三、ANSYS分析 (2) (一)前处理 (2) (1)定义单元类型 (2) (2)定义材料属性 (3) (3)建立工程简化模型 (3) (4)有限元网格划分 (5) (二)模态分析 (5) (1)选择求解类型 (5) (2)建立边界条件 (6) (3)输出设置 (6) (4)求解 (6) (5)读取结果 (6) (6)结果分析 (8) (三)结构试验载荷分析 (8) (1)第二跨跨中模拟车载分析 (8) (2)边跨跨中模拟车载分析 (9) 四、结果分析与强度校核 (10) (一)结果分析 (10)
(二)简单强度校核 (10) 参考文献 (11)
连续刚构桥分析 一、工程背景： 随着我国经济的发展，对交通运输的要求也不断提高；高速路，高铁线等遍布全国，这就免不了要架桥修路。截至2014年年底，我国公路桥梁总数已达万座，万延米i。进百万的桥梁屹立在我国交通线上，其安全便是头等大事。随着交通运输线的再扩大，连续刚构桥跨越能力大，施工难度小，行车舒顺，养护简便，造价较低等优点将被广泛应用。 二、工程模型： 现有某预应力混凝土连续刚构桥，桥梁全长为184m，宽13m，其中车行道宽，两侧防撞栏杆各主梁采用C50混凝土。桥梁设计载荷为公路—— 级。 图2-1桥梁侧立面图 上部结构为48m+88m+48m三跨预应力混凝土边界面连续箱梁。箱梁为单箱双室箱形截面，箱梁根部高5m，中跨梁高，边跨梁端高。箱梁顶板宽,底板宽,翼缘板悬臂长，箱梁高度从距墩中心处到跨中合龙段处按二次抛物线变化。0号至3号块长3m（4x3m)，4、5号块长，6号块到合龙段长4m(6x4m)，合龙段长2m。边跨端部设横隔板，墩顶0号块设两道厚横隔板。0号块范围内箱梁底板厚度为，1号块范围内底板厚度由线性变化到，2号块到合龙段范围内底板厚度由线性变化到。全桥顶板厚度为。0到5号块范围内腹板厚度为，6至7号块范围内腹板厚度由线性变化到，8号块到合龙段范围内夫板厚度为。 下部结构桥采用C50混凝土双薄壁墩，横向宽，厚，高25m双壁间设系梁，下设10mX10m矩形承台，厚。ii 图2-2主梁纵抛面图 图2-3 箱梁截面图 三、ANSYS分析： （一）前处理
ansys连续梁桥（ansys连续梁桥） The key steps of modeling and analysis are as follows: The center line 1, box beam to simulate the plate thickness is the sideline, box girder floor, roof, web and flange plate. 2, determine the location of each key point. The thickness of flange and floor 3, correct simulation of chamfering and the thickness of the gradient. 4. After entering into the process, analyze the stress and deformation. /prep7 /title, three, span, continus, Grider K, 1,0,0 K, 2, -2.1,0 K, 3, -2.6, -0.125 K, 4, -2.8, -0.125 K, 5, -3, -0.125 K, 6, -3.4857, -0.1036 K, 7, -3.9714, -0.0821 K, 8, -4.4571, -0.0607 K, 9, -4.9429, -0.0393 K, 10, -5.4286, -0.0179 K, 11, -5.9143,0.0036 K, 12, -6.4,0.025 K, 13, -2.800, -1.85 K, 14,0.0000, -1.85 Kgen, 9,1,12,1,0,0,49/8100 Kgen, 2,1,12,1,0,0,50900 Kgen, 9901912,1,0,0,34.5/8100 Kgen, 2901912,1,0,0,35900 C1=0.000843399 C2=0.001701323
腾讯朋友 ?首页 ?好友 ?社交 ?应用 ?消息 ?我的主页 ?设置 ?换肤 ?建议 ?退出 搜索搜索 郑军涛 ?主页 ?说说new ?日志 ?相册 ?分享 ?留言板 ?投票 ?礼物 ?好友 利用ANSYS建立变截面箱梁 分享 利用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,7,0,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 numcmp,all adele,1 l,2,9 l,13,20
!!连续梁桥 /prep7 et,1,4 !!!!定义梁单元 et,2,21 !!!!定义mass21单元 !!定义粱材料!!泊松比!!密度 mp,ex,2,3.45e10 !!直线段梁材料和1M段梁材料mp,nuxy,2,0.2 mp,dens,2,3302.153125 mp,ex,3,3.45e10 mp,nuxy,3,0.2 mp,dens,3,3301.658695 mp,ex,4,3.45e10 mp,nuxy,4,0.2 mp,dens,4,3299.906778 mp,ex,5,3.45e10 mp,nuxy,5,0.2 mp,dens,5,3298.327219 mp,ex,6,3.45e10 mp,nuxy,6,0.2 mp,dens,6,3292.351605 mp,ex,7,3.45e10 mp,nuxy,7,0.2 mp,dens,7,3284.137255 mp,ex,8,3.45e10 mp,nuxy,8,0.2 mp,dens,8,3271.802136 mp,ex,9,3.45e10 mp,nuxy,9,0.2 mp,dens,9,3260.41903 mp,ex,10,3.45e10 mp,nuxy,10,0.2 mp,dens,10,3248.193657
mp,nuxy,11,0.2 mp,dens,11,3235.117644 mp,ex,12,3.45e10 mp,nuxy,12,0.2 mp,dens,12,3221.585664 mp,ex,13,3.45e10 mp,nuxy,13,0.2 mp,dens,13,3208.826871 mp,ex,14,3.45e10 mp,nuxy,14,0.2 mp,dens,14,3194.279207 mp,ex,15,3.45e10 mp,nuxy,15,0.2 mp,dens,15,3179.924673 mp,ex,16,3.45e10 mp,nuxy,16,0.2 mp,dens,16,3166.445716 mp,ex,17,3.45e10 mp,nuxy,17,0.2 mp,dens,17,3152.555731 mp,ex,18,3.45e10 mp,nuxy,18,0.2 mp,dens,18,3138.312105 mp,ex,19,3.45e10 mp,nuxy,19,0.2 mp,dens,19,3124.795334 mp,ex,20,3.45e10 mp,nuxy,20,0.2 mp,dens,20,3110.7135 mp,ex,21,3.45e10 mp,nuxy,21,0.2 mp,dens,21,3097.080875
第一天目标： 熟悉ANSYS基本关键字的含义k --> Keypoints关键点l --> Lines线a --> Area 面v --> Volumes体e --> Elements单元n --> Nodes节点cm --> component组元et --> element type单元类型mp --> material property材料属性r --> real constant实常数d --> DOF constraint约束f --> Force Load集中力sf --> Surface load on nodes 表面载荷bf --> Body Force on Nodes体载荷ic --> Initial Conditions初始条件第二天目标： 了解命令流的整体结构,掌握每个模块的标识!文件说明段/BATCH/TILE,test analysis!定义工作标题/FILENAME,test!定义工作文件名/PREP7!进入前处理模块标识!定义单元,材料属性,实常数段ET,1,SHELL63!指定单元类型ET,2,SOLID45!指定体单元MP,EX,1,2E8!指定弹性模量MP,PRXY,1, 0.3!输入泊松比MP,DENS,1, 7.8E3!输入材料密度R,1, 0.001!指定壳单元实常数-厚度......!建立模型K,1,0,0,,!定义关键点 K,2,50,0,,K,3,50,10,,K,4,10,10,,K,5,10,50,,K,6,0,50,,A,1,2,3,4,5,6,!由关键点生成面......!划分网格ESIZE,1,0,AMESH, 1......FINISH!前处理结束标识/SOLU!进入求解模块标识!施加约束和载荷DL,5,,ALLSFL,3,PRES,1000SFL,2,PRES, 1000......SOLVE!求解标识FINISH!求解模块结束标识/POST1!进入通用后处理器标识....../POST26!进入时间历程后处理器……/EXIT,SAVE!退出并存盘以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助/ANGLE!指定绕轴旋转视图/DIST!说明对视图进行缩放/DEVICE!设置图例的显示,如： 风格,字体等/REPLOT!重新显示当前图例/RESET!恢复缺省的图形设置/VIEW!设置观察方向/ZOOM!对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分 1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例：
桥梁计算（常用的计算方法） 在Ansys单元库中，有近200种单元类型，在本章中将讨论一些在桥梁 工程中常用到的单元，包括一些单元的输人参数，如单元名称、节点、自由度、实常数、材料特性、表面荷载、体荷载、专用特性、关键选项KEYOPl等。＊＊＊关于单元选择问题 这是一个大问题，方方面面很多，主要是掌握有限元的理论知识。首先 当然是由问题类型选择不同单元，二维还是三维，梁，板壳，体，细梁，粗梁，薄壳，厚壳，膜等等，再定义你的材料：各向同性或各向异性，混凝土的各项?参数，粘弹性等等。接下来是单元的划分与网格、精度与求解时间的要求等 选择，要对各种单元的专有特性有个大概了解。 使用Ansys，还要了解Ansys的一个特点是笼统与通用，因此很多东西 被掩盖到背后去了。比如单元类型，在Solid里面看到十几种选择，Solid45，Solidl85，Solid95等，看来区别只是节点数目上。但是实际上每种类型里还 有Keyopt分成多种类型，比如最常用的线性单元Solid45，其Keyopt(1)：in●cludeorexclude extradisplacement shapes，就分为非协调元和协调元，Keyopt (2)：fullintegration。rreducedintegration其实又是两种不同的单元，这样不同 组合一下这个Solid45实际上是包含了6种不同单元，各有各的不同特点和 用处。因此使用Ansys要注意各单元的Keyopt选项。不同的选项会产生不 同的结果。· 举例来说：对线性元例如Solid45，要想把弯曲问题计算得比较精确，必 须要采用非协调模式。采用完全积分会产生剪切锁死，减缩积分又会产生 零能模式(ZEM)，非协调的线性元可以达到很高的精度，并且计算量比高阶 刷、很多，在变形较大时，用Enhanced Strain比非协调位移模式(Enhaced Displacement)更好(Solidl85)。但是这些非协调元都要求网格比较规则才 行，网格不规则的话，精度会大大下降，所以如何划分网格也是一门实践性 很强的学问。 采用高阶单元是提高精度的好办法，拿不定主意时采用高阶元是个比 较保险的选择，但是高阶单元在某些情况下也会出现剪切锁死，并且很难发 现，因此用减缩积分的高阶元通常是最保险的选择，但是在大位移时，网格 扭曲较大，减缩积分就不适用。 不同结构形式的桥梁具有不同的力学行为，必须针对性地创建其模型，? 选择维数最低的单元去获得预期的效果(尽量做到能选择点而不选择线，能 选择线而不选择平面，能选择平面而不选择壳，能选择壳而不选择三维实 体)。下面的几节介绍一下桥梁工程计算中经常会用到的单元。 ＊＊＊桥梁仿真单元类型
利用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
!****************************************************************************** *********************** ! case2:无偏载(以跨径布置30m+40m+30m,桥宽8.5为例） ! 上海城市设计研究院L1+L2+L3预应力混凝土曲线连续梁桥结构分析 ! 两端为抗扭支座，中间支座为点铰支座 ! 每次要记得修改横隔梁的参数，即Mass21单元的实常数 !****************************************************************************** *********************** FINI /CLE /prep7 !DEFINE THE ELEMENTARY PARAMETERS *DIM,L,ARRAY,10 *DIM,H,ARRAY,10 *DIM,CITA,ARRAY,10 !*****以下参数均可修改*************** N=3 !跨数 L(1)=30 !第一跨 L(2)=40 !第二跨 L(3)=30 !第三跨 e1=1.25 !1#墩处内支座到中心线的间距 e2=1.25 !1#墩处外支座到中心线的间距 e3=0 !2#墩处的支座偏心距（正的表示外偏） e4=0 !3#墩处的支座偏心距 e5=1.25 !4#墩处内支座到中心线的间距 e6=1.25 !4#墩处外支座到中心线的间距 R=10000 !曲线桥半径 H0=1.0 !梁底到截面形心处的高度 M=16146 !mass21单元质量 J=27246.38 !mass21单元转动惯量 !************************************* LL=0.0 *DO,I,1,N LL=LL+L(I) CITA(I)=L(I)/R/3.1415925*180 *ENDDO CITA0=LL/R/3.1415925*180
/prep7 !DEFINE THE ELEMENTARY PARAMETERS *DIM,L,ARRAY,10 *DIM,DISTC,ARRAY,10 N=4 NBOX=2 L(1)=20 L(2)=30 L(3)=30 L(4)=20 NUML=NBOX+3 NUMC=(L(1)+L(2)+L(3)+L(4))+1 CITA=25/180*3.1415926 DISTC(1)=2.9584 DISTC(2)=4.0802 DISTC(3)=0.00 DISTC(4)=4.0802 DISTC(5)=2.9584 !DEFINE THE NODES OF BRIDGE N, 1, 0, 0, 0.00 N, NUMC, NUMC-1, 0, 0.00 FILL,1,NUMC N, 200+1, (DISTC(1)+DISTC(2))*SIN(CITA)/COS(CITA), 0.0,-(DISTC(1)+DISTC(2)) N, 200+NUMC,NUMC-1+(DISTC(1)+DISTC(2))* SIN(CITA)/COS(CITA),0.0,-(DISTC(1)+DISTC(2)) FILL,201,200+NUMC N, 400+1, DISTC(2)*SIN(CITA)/COS(CITA), 0.0, -DISTC(2) N, 400+NUMC,NUMC-1+DISTC(2)*SIN(CITA)/ COS(CITA), 0.0,-DISTC(2) FILL,401,400+NUMC N, 600+1, -DISTC(2)*SIN(CITA)/COS(CITA), 0.0, DISTC(2) N, 600+NUMC,NUMC-1-DISTC(4)*SIN(CITA)/C OS(CITA),0.0, DISTC(4) FILL,601,600+NUMC N, 800+1, -(DISTC(4)+DISTC(5))*SIN(CITA)/COS(CITA ),0.0, DISTC(4)+DISTC(5) N, 800+NUMC,NUMC-1-(DISTC(4)+DISTC(5))* SIN(CITA)/COS(CITA),0.0, DISTC(4)+DISTC(5) FILL,801,800+NUMC ! 定义纵梁单元材料、几何参数 ET,1,BEAM4 MP,EX , 1, 3.45E+10 MP,NUXY, 1, 0.2000 MP,DENS, 1, 2600 MP,ALPX, 1, 1.00E-05 R,1,0.3133,0.0871,0.0959,1.0,0.25,, RMORE ,,0.00265,,,,, ET,2,BEAM4 MP,EX , 2, 3.45E+10 MP,NUXY, 2, 0.2000 MP,DENS, 2, 2600 MP,ALPX, 2, 1.00E-05 R,2,1.6048,0.5486,1.5400,2.0,1.0,, RMORE ,,0.036,,,,, ET,3,BEAM4 MP,EX , 3, 3.45E+10 MP,NUXY, 3, 0.2000 MP,DENS, 3, 2600 MP,ALPX, 3, 1.00E-05 R,3,2.3625,0.9885,2.2600,2.0,1.0,, RMORE ,,0.0956,,,, ! 定义横梁单元材料、几何参数 ET,4,BEAM4 MP,EX , 4, 3.45E+10 MP,NUXY, 4, 0.2000 MP,DENS, 4, 0 MP,ALPX, 4, 1.00E-05 R,4,0.500,0.4582,0.0416,1.0,1.0,, RMORE ,,0.00879,,,,, ET,5,BEAM4 MP,EX , 5, 3.45E+10 MP,NUXY, 5, 0.2000 MP,DENS, 5, 0 MP,ALPX, 5, 1.00E-05 R,5,1.0250,0.2770,0.0495,1.0,1.0,, RMORE ,,0.0609,,,,, ! 生成纵向梁格单元 TYPE,1 MAT,1
ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣，建立 的时候一针一线，小心而漫长，拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生 一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT
1模型简介 模型以壳单元shell63为主，横梁采用beam44单元，吊杆采用link10单元，横梁、吊杆与拱肋的连接采用节点耦合的方式。不考虑纵坡和横披。 11 1213 14 15 16X Y Z ROT CP ACEL 图1 三维图
111213 14 15 16X Y Z 图2 正立面图 2建模命令流 finish /clear /config,nproc,2 /prep7 !双核运行 !材料 !钢材 EX,1,2.06e11 PRXY ,1,0.3 DENS,1,7800 MPDA TA,ALPX,1,,1.2e-5 !吊杆钢材 EX,2,2.06e11 PRXY ,2,0.3 DENS,2,7800 MPDA TA,ALPX,2,,1.2e-5 !桥面钢材包括铺装层重量 EX,3,2.06e11
PRXY,3,0.3 DENS,3,7800*3 MPDA TA,ALPX,3,,1.2e-5 !单元 ET,1,plane82 ET,2,shell63 !设置shell63关键字 KEYOPT,2,1,0 KEYOPT,2,2,0 KEYOPT,2,3,0 KEYOPT,2,5,2 KEYOPT,2,6,0 KEYOPT,2,7,0 KEYOPT,2,8,0 KEYOPT,2,9,0 KEYOPT,2,11,0 et,3,beam189 !设置beam189关键字 KEYOPT,3,1,1 KEYOPT,3,2,0 KEYOPT,3,4,2 KEYOPT,3,6,3 KEYOPT,3,7,2 KEYOPT,3,8,3 KEYOPT,3,9,3 KEYOPT,3,10,0 KEYOPT,3,11,0 KEYOPT,3,12,0 ET,4,solid45 et,5,link10 r,1,0.010 r,2,0.016 r,3,0.02 r,4,0.04 !建立各种截面 !钢横梁 !HL1/HL3 a=1.0 b=0.680 th=0.02 ts=0.016 RECTNG,-a/2,a/2,-b/2,b/2, RECTNG,-a/2+ts,a/2-ts,-b/2+th,b/2-th,
平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end
/TITLE, BRIDGE XX /PREP7 ! 定义单元类型 ET,1,SHELL63 ET,2,SOLID45 ET,3,LINK8 !定义1,钢，2，混凝土，3，无限刚度MP,EX,1, MP,DENS,1, MP,PRXY,1, mp,alpx,1, MP,EX,2, MP,DENS,2, MP,PRXY,2, MP,EX,3, MP,DENS,3, MP,PRXY,3, *set,nz,154 !定义箱梁N1 R,1,7,7,7,7 !定义箱梁N2 R,2,8,8,8,8
!定义箱梁N3 R,3,4,4,4,4 !定义混凝土solid65的实常数 R, 4 !定义钢绞线的实常数 R, 5, 140 !降温=预加力/线膨胀系数*钢筋截面积*弹性模量 tem1=* tem2=* !定义纵向梯度参数表 *dim, z, array, nz *SET,z(1),130,260,371,482,593,702,811,920,1029,1138,1247,1356,146 5,1574,1683,,, *SET,z(19),,,,,,,,,,,,, *SET,z(31),,,,,, *SET,z(37),,,,, *SET,z(41),,,,,,,,,,,,,, *SET,z(55),,,,,,,,,,, *SET,z(65),,,,,,,, *SET,z(73),,,,,,,, *SET,z(80),,,,,,,,,,, *SET,z(91),,,,,,,,,,
*SET,z(100),,,,,,,,, *SET,z(109),,,,,,, *SET,z(115),,,,,,,,,,,, *SET,z(127),,,,,,,,, *SET,z(135),,,,,,,,,, *SET,z(145),,,,,,,,,, ! !定义横隔梁的梯度参数表 *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,,,-133,0,,,400,,633,-633,-533,,,-133, *SET,initial(17,1),0,,,400,,633,-633,-633,633,633,-766,,-633,,-50 0,500,,633, *SET,initial(35,1),,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,
/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT ! Line numbers off !