进行优化分析,设计变量有变化,但是状态变量没有变化,请帮忙看一下,优化分析部分的命令流有什么错误?
这是全部命令流,只看优化分析部分就行,问题就出在优化设计部分
finish
/clear
/filname,HERWDQ
/title,BL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!设计变量初值
E_x1=3.45E10
D_z1=2550
E_x2=3.45E10
D_z2=2550
E_x3=3.45E10
D_z3=2550
E_x4=3.45E10
D_z4=2550
E_x5=3.45E10
D_z5=2550
E_x6=3.45E10
D_z6=2550
E_x7=3.45E10
D_z7=2550
E_x8=3.45E10
D_z8=2550
E_x9=3.45E10
D_z9=2550
E_x10=3.45E10
D_z10=2550
E_x11=3.45E10
D_z11=2550
E_x12=3.45E10
D_z12=2550
E_x13=3.45E10
D_z13=2550
E_x14=3.45E10
D_z14=2550
E_x15=3.45E10
D_z15=2550
E_x16=3.45E10
D_z16=2550
E_x17=3.45E10
D_z17=2550
E_x18=3.45E10
D_z18=2550
E_x19=3.45E10
D_z19=2550
E_x20=3.45E10
D_z20=2550
E_x21=3.45E10
D_z21=2550
E_x22=3.45E10
D_z22=2550
E_x23=3.45E10
D_z23=2550
E_x24=3.45E10
D_z24=2550
/prep7
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!实测数据添加
*dim,sc_pinlv,,3,1,1
sc_pinlv(1)=2.8
sc_pinlv(2)=3.3
sc_pinlv(3)=5.1
f1=2.45e10 !!!!!!!!!!!!!!!!!弹模下限
f2=4.45e10 !!!!!!!!!!!!!!!!!弹模上限
h1=2450 !!!!!!!!!!!!!!!!!密度下限
h2=2650 !!!!!!!!!!!!!!!!!密度上限
q=5 !!!!!!!!!!!!!!!!!最大迭代次数
*do,I,1,24
et,i,44
mp,ex,i,E_x%i%
mp,prxy,i,0.2
mp,dens,i,D_z%i%
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义梁单元节点
*do,I,1,15,1
K,I,i-15,25.935-(i-1)*0.026,0
*enddo
K,17,1.709,25.526,0
*do,I,18,45,1
K,I,i-15,25.935-(i-1)*0.026,0
*enddo
*do,I,46,49,1
K,I,30+(i-45)*0.2,24.795-(i-45)*0.005,0
*enddo
*do,I,50,53,1
K,I,30+(i-44)*0.2,24.7698-(i-50)*0.0052,0
*enddo
*do,I,54,74,1
K,I,i-22 ,24.749-(i-54)*0.026,0
*enddo
K,75,52.5,24.223,0
K,76,85.5,23.233,0
*do,I,77,97,1
K,I,i+9 ,23.216-(i-77)*0.035,0
*enddo
*do,I,98,101,1
K,I,106.2+(i-98)*0.2,22.5222-(i-98)*0.0072,
*enddo
*do,I,102,106,1
K,I,107+(i-101)*0.2,22.488-(i-102)*0.007,0
*enddo
*do,I,107,134,1
K,I,i+2,22.426-(i-107)*0.034,0
*enddo
*do,I,135,161,1
K,I,i+3,21.43-(i-135)*0.035,0
*enddo
k,200,0,0,0
k,201,0.946,0.63,0
k,202,1,0.666,0
k,203,1.2,0.798,0
k,204,1.32,0.883,0
k,205,1.709,1.133,0
k,206,1.963,1.299,0
*do,i,207,257,1
k,i,i-205,0.671*(i-205)-0.0048*(i-205)* (i-205),0
*enddo
k,258,52.5,22.021 ,0
k,259,52.5,23.2155,0
k,260,53,23.248,0
k,261,54,23.313,0
k,262,55,23.3725,0
k,263,56,23.4275,0
k,264,57,23.478,0
k,265,58,23.523,0
k,266,59,23.5635,0
k,267,60,23.599,0
k,268,61,23.6295,0
k,269,62,23.6555,0
k,270,63,23.6765,0
k,271,64,23.6925,0
k,272,65,23.704,0
k,273,66,23.71,0
k,274,67,23.711,0
k,275,68,23.7085,0
k,276,68.75,23.7025,0
k,277,69,23.7,0
k,278,69.25,23.695,0
k,279,70,23.6785,0
k,280,71,23.6515,0
k,281,71.5,23.63575,0
k,282,72,23.62,0
k,283,73,23.583,0
k,284,74,23.5425,0
k,285,75,23.495,0
k,286,76,23.4455,0
k,287,77,23.3895,0
k,288,78,23.329,0
k,289,79,23.2635,0
k,290,80,23.193,0
k,291,81,23.118,0
k,292,82,23.0375,0
k,293,83,22.9525,0
k,294,84,22.863,0
k,295,85,22.768,0
k,296,85.5,21.913,0
k,297,85.5,22.7205,0
*do,i,298,348,1
k,i,i-212,-0.0048*(i-212)* (i-212)+0.6621*(i-212)+0.4126,0
*enddo
K,349,136.037,1.659,0
K,350,136.291,1.533,0
k,351,136.672,1.283,0
K,352,136.8,1.198,0
K,353,137,1.066,0
K,354,137.054,1.03,0
K,355,138,0.4,0
*do,I,500,520,1
K,I,0,(i-500)*1.2135+0.7,-2.975
*enddo
*do,I,521,530,1
K,I,(i-520)*1.6,23.493-(i-521)*1.477,-2.975
*enddo
*do,I,531,540,1
K,I,16,11.636+(i-531)*1.436,-2.975
*enddo
*do,I,541,549,1
K,I,17.44+(i-541)*1.44, 23.76734-(i-541)*0.79266,-2.975
*enddo
*do,I,550,554,1
K,I,30.4,16.6334+(i-550)*1.8854,-2.975
*enddo
*do,I,555,559,1
K,I,31.6,17.0792+(i-555)*1.77395,-2.975
*enddo
*do,I,560,563,1
K,I,34.2+(i-560)*2.6,23.2215-(i-560)*0.9535,-2.975
*enddo
K,563,42,20.361,-2.975
K,564,42,23.893,-2.975
K,565,96,20.361,-2.975
K,566,96,22.373,-2.975
K,567,98.6,20.7496,-2.975
K,568,101.2,21.1382,-2.975
K,569,103.8,21.5268,-2.975
*do,I,570,574,1
K,I,106.4,17.0792+(i-570)*1.20905,-2.975
*enddo
*do,I,575,579,1
K,I,107.6,16.6334+(i-575)*1.31015,-2.975
*enddo
K,580,111.2,17.82005,-2.975
K,581,114.8,19.0067,-2.975
K,582,118.4,20.19335,-2.975
*do,I,583,593,1
K,I,122,10.2+(i-583)*1.118,-2.975
*enddo
*do,I,594,602,1
K,I,123.6+(i-594)*1.6,11.263+(i-594)*1.063,-2.975
*enddo
*do,I,603,623,1
K,I,138,0.7+(i-603)*1.0065,-2.975
*enddo
*do,I,624,634,1
K,I,148,10+(i-624)*1.048066,-2.975
*enddo
*do,I,635,645,1
K,I,156,11+(i-635)*0.920118,-2.975
*enddo
*do,I,646,666,1
K,I,0, (i-646)*1.2135+0.7,2.975
*enddo
*do,I,667,676,1
K,I, (i-666)*1.6,23.493-(i-667)*1.477,2.975
*enddo
*do,I,677,686,1
K,I,16, 11.636+(i-677)*1.436,2.975
*enddo
*do,I,687,695,1
K,I, 17.44+(i-687)*1.44, 23.76734-(i-687)*0.79266,2.975
*enddo
*do,I,696,700,1
K,I, 30.4,16.6334+(i-696)*1.8854,2.975
*enddo
*do,I,701,705,1
K,I,31.6,17.0792+(i-701)*1.77395,2.975
*enddo
*do,I,706,709,1
K,I, 34.2+(i-706)*2.6,23.2215-(i-706)*0.9535,2.975
*enddo
K,709,42,20.361,2.975
K,710,42,23.893,2.975
K,711,96,20.361,2.975
K,712,96,22.373,2.975
K,713,98.6,20.7496,2.975
K,714,101.2,21.1382,2.975
K,715,103.8,21.5268,2.975
K,716,106.4,21.9154,2.975
*do,I,717,720,1
K,I,106.4,17.0792+(i-717)*1.20905,2.975
*enddo
*do,I,721,725,1
K,I,107.6,16.6334+(i-721)*1.31015,2.975
*enddo
K,726,111.2,17.82005,2.975
K,727,114.8,19.0067,2.975
K,728,118.4,20.19335,2.975
K,729,122,21.38,2.975
*do,I,730,739,1
K,I,122,10.2+(i-730)*1.118,2.975
*enddo
*do,I,740,749,1
K,I,123.6+(i-740)*1.6,11.263+(i-740)*1.063,2.975
*enddo
*do,I,750,769,1
K,I,138,0.7+(i-750)*1.0065,2.975
*enddo
*do,I,770,780,1
K,I,148,10+(i-770)*1.048066,2.975
*enddo
*do,I,781,791,1
K,I,156,11+(i-781)*0.920118,2.975
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义梁单元截面
K,1000, ,26.935
sectype,1,beam,mesh
secread,xiangliang,,,mesh
secoffset,cent
*do,I,1,14,1
L,I,i+1
Latt,1, ,1, ,1000,,1
Lsel,s,,,i
Lesize,I,,,2
Lmesh,i
Allsel,all
*enddo
/eshape,1
Eplot
L,15,17
Latt,1, ,1, ,1000,1
Lsel,s,,,15
Lesize,15,,,2
Lmesh,15
allsel,all
*do,I,17,48,1
L,I,i+1
Latt,1, ,1, ,1000,,1
Lsel,s,,,i-1
Lesize,I-1,,,2
Lmesh,i-1
Allsel,all
*enddo
*do,I,50,74,1
L,I,i+1
Latt,1, ,1, ,1000,,1
Lsel,s,,,i-2
Lesize,I-2,,,2
Lmesh,i-2
Allsel,all
*enddo
*do,I,76,100,1
L,I,i+1
Latt,1, ,1, ,1000,,1
Lsel,s,,,i-3
Lesize,I-3,,,2
Lmesh,i-3
Allsel,all
*enddo
*do,I,102,160,1
L,I,i+1
Latt,1, ,1, ,1000,,1
Lsel,s,,,i-4
Lesize,I-4,,,2
Lmesh,i-4
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义拱圈单元
sectype,2,beam,mesh
secread,gongquan,,,mesh
secoffset,cent
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义拱脚(弹模E2)
*do,I,200,206,1
L,I,i+1
Latt,2, ,2, ,1000,,2
Lsel,s,,,i-43
Lesize,I-43,,,2
Lmesh,i-43
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!拱圈定义
*do,I,207,257,1
L,I,i+1
Latt,1, ,1, ,1000,,2
Lsel,s,,,i-43
Lesize,I-43,,,2
Lmesh,i-43
Allsel,all
*enddo
K,1001,123,25
k,1002,85.5,25
L,296,298
Latt,1, ,1, ,1002,,2
Lsel,s,,,215
Lesize,215,,,2
Lmesh,215
allsel,all
*do,I,298,347,1
L,I,i+1
Latt,3, ,3, ,1001,,2
Lsel,s,,,i-82
Lesize,I-82,,,2
Lmesh,i-82
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义拱脚(弹模E3)
*do,I,348,354,1
L,I,i+1
Latt,3, ,3, ,1001,,2
Lsel,s,,,i-82
Lesize,I-82,,,2
Lmesh,i-82
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义立柱单元
sectype,3,beam,mesh
secread,i-i,,,mesh
secoffset,cent
*do,I,500,519,1
L,I,i+1
Latt,1, ,1, ,1000,,3
Lsel,s,,,i-227
Lesize,I-227,,,2
Lmesh,i-227
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义2A杆件(弹模E4)
sectype,4,beam,mesh
secread,j-j,,,mesh
secoffset,cent
*do,I,520,529,1
L,I,i+1
Latt,4, ,4, ,1001,,4
Lsel,s,,,i-227
Lesize,I-227,,,2
Lmesh,i-227
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义3A杆件(弹模E5)
*do,I,530,539,1
L,I,i+1
Latt,5, ,5, ,1001,,4
Lsel,s,,,i-227
Lesize,I-227,,,2
Lmesh,i-227
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义6A杆件(弹模E6)
*do,I,540,549,1
L,I,i+1
Latt,6, ,6, ,1001,,4
Lsel,s,,,i-227
Lesize,I-227,,,2
Lmesh,i-227
Allsel,all
*enddo
sectype,5,beam,mesh
secread,k-k,,,mesh
secoffset,cent
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义7A杆件(弹模E7)
*do,I,550,553,1
L,I,i+1
Latt,7, ,7, ,1001,,5
Lsel,s,,,i-227
Lesize,I-227,,,2
Lmesh,i-227
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义9A杆件(弹模E8)
*do,I,555,558,1
L,I,i+1
Latt,8, ,8, ,1001,,5
Lsel,s,,,i-228
Lesize,I-228,,,2
Lmesh,i-228
Allsel,all
*enddo
*do,I,559,562,1
L,I,i+1
Latt,1, ,1, ,1001,,5
Lsel,s,,,i-228
Lesize,I-228,,,2
Lmesh,i-228
Allsel,all
*enddo
L,563,564
Latt,1, ,1, ,1002,,5
Lsel,s,,,335
Lesize,335,,,2
Lmesh,335
allsel,all
L,565,566
Latt,1, ,1, ,1000,,5
Lsel,s,,,336
Lesize,336,,,2
Lmesh,336
allsel,all
L,565,567
Latt,1, ,1, ,1000,,4
Lsel,s,,,337
Lesize,337,,,2
Lmesh,337
allsel,all
*do,I,567,568,1
L,I,i+1
Latt,1, ,1, ,1000,,4
Lsel,s,,,i-229
Lesize,I-229,,,2
Lmesh,i-229
Allsel,all
*enddo
L,569,574
Latt,1, ,1, ,1000,,4
Lsel,s,,,340
Lesize,340,,,2
Lmesh,340
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义9B杆件(弹模E9)
*do,I,570,573,1
L,I,i+1
Latt,9, ,9, ,1000,,5
Lsel,s,,,i-229
Lesize,I-229,,,2
Lmesh,i-229
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义7B杆件(弹模E10)
*do,I,575,578,1
L,I,i+1
Latt,10, ,10, ,1000,,5
Lsel,s,,,i-230
Lesize,I-230,,,2
Lmesh,i-230
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义6B杆件(弹模E11)
L,575,580
Latt,11, ,11, ,1000,,4
Lsel,s,,,349
Lesize,349,,,2
Lmesh,349
allsel,all
*do,I,580,581,1
L,I,i+1
Latt,1, ,1, ,1000,,4
Lsel,s,,,i-230
Lesize,I-230,,,2
Lmesh,i-230
Allsel,all
*enddo
L,582,593
Latt,11, ,11, ,1000,,4
Lsel,s,,,352
Lesize,352,,,2
Lmesh,352
allsel,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!定义3B杆件(弹模E12)
*do,I,583,592,1
L,I,i+1
Latt,12, ,12, ,1000,,4
Lsel,s,,,i-230
Lesize,I-230,,,2
Lmesh,i-230
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义2B杆件(弹模E13)
L,583,594
Latt,13, ,13, ,1000,,4
Lsel,s,,,363
Lesize,363,,,2
Lmesh,363
allsel,all
*do,I,594,601,1
L,I,i+1
Latt,13, ,13, ,1000,,4
Lsel,s,,,i-230
Lesize,I-230,,,2
Lmesh,i-230
Allsel,all
*enddo
L,602,623
Latt,13, ,13, ,1000,,4
Lsel,s,,,372
Lesize,372,,,2
Lmesh,372
allsel,all
sectype,6,beam,mesh
secread,m-m,,,mesh
secoffset,cent
*do,I,603,622,1
L,I,i+1
Latt,1, ,1, ,1002,,6
Lsel,s,,,i-230
Lesize,I-230,,,2
Lmesh,i-230
Allsel,all
*enddo
*do,I,624,633,1
L,I,i+1
Latt,1, ,1, ,1002,,4
Lsel,s,,,i-231
Lesize,I-231,,,2
Lmesh,i-231
Allsel,all
*enddo
*do,I,646,665,1
L,I,i+1
Latt,1, ,1, ,1000,,3
Lsel,s,,,i-243
Lesize,I-243,,,2
Lmesh,i-243
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义2A杆件(弹模E14)
*do,I,666,675,1
L,I,i+1
Latt,14, ,14, ,1001,,4
Lsel,s,,,i-243
Lesize,I-243,,,2
Lmesh,i-243
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义3A杆件(弹模E15)
*do,I,676,685,1
L,I,i+1
Latt,15, ,15, ,1001,,4
Lsel,s,,,i-243
Lesize,I-243,,,2
Lmesh,i-243
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义6A杆件(弹模E16)
*do,I,686,695,1
L,I,i+1
Latt,16, ,16, ,1001,,4
Lsel,s,,,i-243
Lesize,I-243,,,2
Lmesh,i-243
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义7A杆件(弹模E17)
*do,I,696,699,1
L,I,i+1
Latt,17, ,17, ,1001,,5
Lsel,s,,,i-243
Lesize,I-243,,,2
Lmesh,i-243
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义9A杆件(弹模E18)
*do,I,701,704,1
L,I,i+1
Latt,18, ,18, ,1001,,5
Lsel,s,,,i-244
Lesize,I-244,,,2
Lmesh,i-244
Allsel,all
*enddo
*do,I,705,708,1
L,I,i+1
Latt,1, ,1, ,1001,,5
Lsel,s,,,i-244
Lesize,I-244,,,2
Lmesh,i-244
Allsel,all
*enddo
L,709,710
Latt,1, ,1, ,1001,,5
Lsel,s,,,465
Lesize,465,,,2
Lmesh,465
allsel,all
L,711,712
Latt,1, ,1, ,1001,,5
Lsel,s,,,466
Lesize,466,,,2
Lmesh,466
allsel,all
L,711,713
Latt,1, ,1, ,1000,,4
Lsel,s,,,467
Lesize,467,,,2
Lmesh,467
allsel,all
*do,I,713,715,1
L,I,i+1
Latt,1, ,1, ,1000,,4
Lsel,s,,,i-245
Lesize,I-245,,,2
Lmesh,i-245
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义9B杆件(弹模E19)
*do,I,717,719,1
L,I,i+1
Latt,19, ,19, ,1000,,5
Lsel,s,,,i-246
Lesize,I-246,,,2
Lmesh,i-246
Allsel,all
*enddo
L,720,716
Latt,19, ,19, ,1000,,5
Lsel,s,,,474
Lesize,474,,,2
Lmesh,474
allsel,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义7B杆件(弹模E20)
*do,I,721,724,1
L,I,i+1
Latt,20, ,20, ,1000,,5
Lsel,s,,,i-246
Lesize,I-246,,,2
Lmesh,i-246
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义6B杆件(弹模E21)
L,721,726
Latt,21, ,21, ,1000,,4
Lsel,s,,,479
Lesize,479,,,2
Lmesh,479
allsel,all
*do,I,726,728,1
L,I,i+1
Latt,21, ,21, ,1000,,4
Lsel,s,,,i-246
Lesize,I-246,,,2
Lmesh,i-246
Allsel,all
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义3B杆件(弹模E22)
*do,I,730,738,1
L,I,i+1
Latt,22, ,22, ,1000,,4
Lsel,s,,,i-247
Lesize,I-247,,,2
Lmesh,i-247
Allsel,all
*enddo
L,739,729
Latt,22, ,22, ,1000,,4
Lsel,s,,,492
Lesize,492,,,2
Lmesh,492
allsel,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义2B杆件(弹模E23)
L,730,740
Latt,23, ,23, ,1000,,4
Lsel,s,,,493
Lesize,493,,,2
Lmesh,493
allsel,all
*do,I,740,748,1
L,I,i+1
Latt,23, ,23, ,1000,,4
Lsel,s,,,i-246
Lesize,I-246,,,2
Lmesh,i-246
Allsel,all
*enddo
*do,I,750,768,1
L,I,i+1
Latt,1, ,1, ,1000,,6
Lsel,s,,,i-247
Lesize,I-247,,,2
Lmesh,i-247
Allsel,all
*enddo
L,769,749
Latt,1, ,1, ,1000,,6
Lsel,s,,,522
Lesize,522,,,2
Lmesh,522
allsel,all
*do,I,770,779,1
L,I,i+1
Latt,1, ,1, ,1000,,4
Lsel,s,,,i-247
Lesize,I-247,,,2
Lmesh,i-247
Allsel,all
*enddo
*do,I,781,790,1
L,I,i+1
Latt,1, ,1, ,1000,,4
Lsel,s,,,i-248
Lesize,I-248,,,2
Lmesh,i-248
Allsel,all
*enddo
*do,I,635,644,1
L,I,i+1
Latt,1, ,1, ,1000,,4
Lsel,s,,,i-92
Lesize,I-92,,,2
Lmesh,i-92
Allsel,all
*enddo
L,511,657
Latt,1, ,1, ,1000,,5
Lsel,s,,,553
Lesize,553,,,2
Lmesh,553
allsel,all
L,526,672
Latt,1, ,1, ,1001,,5
Lsel,s,,,554
Lesize,554,,,2
Lmesh,554
allsel,all
L,535,681
Latt,1, ,1, ,1001,,5
Lsel,s,,,555
Lesize,555,,,2
Lmesh,555
allsel,all
L,546,692
Latt,1, ,1, ,1001,,5
Lsel,s,,,556
Lesize,556,,,2
Lmesh,556
allsel,all
L,561,707
Latt,1, ,1, ,1001,,5
Lsel,s,,,557
Lesize,557,,,2
Lmesh,557
allsel,all
L,581,727
Latt,1, ,1, ,1001,,5
Lsel,s,,,558
Lesize,558,,,2
Lmesh,558
allsel,all
L,588,735
Latt,1, ,1, ,1001,,5
Lsel,s,,,559
Lesize,559,,,2
Lmesh,559
allsel,all
L,598,744
Latt,1, ,1, ,1000,,5
Lsel,s,,,560
Lesize,560,,,2
Lmesh,560
allsel,all
L,613,760
Latt,1, ,1, ,1000,,5
Lsel,s,,,561
Lesize,561,,,2
Lmesh,561
allsel,all
sectype,7,beam,mesh
secread,1-1,,,mesh
secoffset,cent
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义跨中单元(弹模E24)
*do,I,259,262,1
L,I,i+1
Latt,24, ,24, ,1000,,7
Lsel,s,,,i+303
Lesize,I+303,,,2
Lmesh,i+303
Allsel,all
*enddo
sectype,8,beam,mesh
secread,2-2,,,mesh
secoffset,cent
*do,I,263,266,1
L,I,i+1
Latt,24, ,24, ,1000,,8
Lsel,s,,,i+303
Lesize,I+303,,,2
Lmesh,i+303
Allsel,all
*enddo
sectype,9,beam,mesh
secread,3-3,,,mesh
secoffset,cent
*do,I,267,270,1
L,I,i+1
Latt,24, ,24, ,1000,,9
Lsel,s,,,i+303
Lesize,I+303,,,2
Lmesh,i+303
Allsel,all
*enddo
sectype,10,beam,mesh
secread,4-4,,,mesh
secoffset,cent
*do,I,271,274,1
L,I,i+1
Latt,24, ,24, ,1000,,10
Lsel,s,,,i+303
Lesize,I+303,,,2
Lmesh,i+303
Allsel,all
*enddo
sectype,11,beam,mesh
secread,5-5,,,mesh
secoffset,cent
*do,I,275,278,1
L,I,i+1
Latt,24, ,24, ,1000,,11
Lsel,s,,,i+303
Lesize,I+303,,,2
Lmesh,i+303
Allsel,all
*enddo
sectype,12,beam,mesh
secread,6-6,,,mesh
secoffset,cent
*do,I,279,282,1
L,I,i+1
Latt,24, ,24, ,1000,,12
Lsel,s,,,i+303
Lesize,I+303,,,2
Lmesh,i+303
Allsel,all
*enddo
sectype,13,beam,mesh
secread,7-7,,,mesh
secoffset,cent
*do,I,283,284,1
L,I,i+1
Latt,24, ,24, ,1001,,13
Lsel,s,,,i+303
Lesize,I+303,,,2
Lmesh,i+303
Allsel,all
*enddo
L,285,286
Latt,24, ,24, ,1001,,13
Lsel,s,,,588
Lesize,588,,,2
Lmesh,588
allsel,all
sectype,14,beam,mesh
secread,8-8,,,mesh
secoffset,cent
*do,I,286,289,1
L,I,i+1
Latt,24, ,24, ,1001,,14
Lsel,s,,,i+303
Lesize,I+303,,,2
Lmesh,i+303
Allsel,all
*enddo
sectype,15,beam,mesh
secread,9-9,,,mesh
secoffset,cent
*do,I,290,294,1
L,I,i+1
Latt,24, ,24, ,1001,,15
Lsel,s,,,i+303
Lesize,I+303,,,2
Lmesh,i+303
Allsel,all
*enddo
L,295,297
Latt,24, ,24, ,1001,,15
Lsel,s,,,598
Lesize,598,,,2
Lmesh,598
allsel,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!约束条件
dk,1,all
dk,161,all
dk,200,all
dk,355,all
dk,624,all
dk,635,all
dk,770,all
dk,781,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!刚性连接
CERIG,629,1620,all
CERIG,629,1095,all
CERIG,54,1697,all
CERIG,54,1172,all
CERIG,114,1777,all
CERIG,114,1252,all
CERIG,178,1833,all
CERIG,178,1308,all
CERIG,195,1325,all
CERIG,195,1850,all
CERIG,243,1870,all
CERIG,243,1345,all
CERIG,332,1350,all
CERIG,332,1875,all
CERIG,380,1891,all
CERIG,380,1366,all
CERIG,397,1924,all
CERIG,397,1399,all
CERIG,461,1940,all
CERIG,461,1415,all
CERIG,521,2020,all
CERIG,521,1495,all
CERIG,561,2141,all
CERIG,561,1616,all
CERIG,593,2182,all
CERIG,593,2223,all
CERIG,710,1737,all
CERIG,710,1212,all
CERIG,770,1292,all
CERIG,770,1312,all
CERIG,770,1817,all
CERIG,770,1837,all
CERIG,814,1866,all
CERIG,814,1341,all
CERIG,2254,858,all
CERIG,2254,287,all
CERIG,2399,291,all
CERIG,2399,862,all
CERIG,903,1349,all
CERIG,903,1874,all
CERIG,947,1370,all
CERIG,947,1386,all
CERIG,947,1895,all
CERIG,947,1911,all
CERIG,1007,1419,all
CERIG,1007,1944,all
CERIG,1091,2024,all
CERIG,1091,1499,all
finish
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!进入求解器
/solu
acel,0,9.8,0
antype,modal
msave,0
modopt,lanb,20,1,20,,on
eqslv,spar
mxpand,20,1,20,1,0.001
outpr,all
outres,all
/output,1,txt
solve
finish
/solu
expass,1
mxpand,20,1,20,1,0.001
outres,all,all
solve
finish
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!后处理
/post1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!提取振型参与系数并放入数组
*dim,pf,,20
*do,i,1,20,1
*GET,PF(i),mode,i,pfact
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!判断振型参与系数选取受弯振型频率与扭转振型频率
*dim,swpl,,20 !受弯频率
*dim,sjpl,,20 !扭转频率
*do,i,1,20,1
*if,pf(i),ABGT,10000,then !如果pf(1)绝对值大于10000
*get,swpl(i),mode,i,freq
*else
*get,sjpl(i),mode,i,freq
*endif
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!提取前3阶受弯频率数组
*dim,pin_lv,,3,1,1
m=1
*do,n,1,20,1
*if,swpl(n),ne,0,then
pin_lv(m)=swpl(n)
m=m+1
*elseif,m,gt,3,then
*exit
*endif
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!状态变量和目标函数
*do,i,1,3,1
*dim,obj_func%i%,,1,1
herwdq_OPT%i%=pin_lv(i)/sc_pinlv(i)
obj_func%i%(1)=(pin_lv(i)-sc_pinlv(i))*(pin_lv(i)-sc_pinlv(i))/sc_pinlv(i)
*enddo
obj_func=0
*do,i,1,3,1
obj_func=obj_func+obj_func%i%(1) !目标函数
*enddo
finish
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!一阶优化工具
LGWRITE,'herw Analysis','lgw'
finish
/opt
OPANL,'herw Analysis','lgw',' '
!设计变量
opvar,E_x2,dv,f1,f2
opvar,D_z2,dv,h1,h2
opvar,E_x3,dv,f1,f2
opvar,D_z3,dv,h1,h2
opvar,E_x4,dv,f1,f2
opvar,D_z4,dv,h1,h2
opvar,E_x5,dv,f1,f2
opvar,D_z5,dv,h1,h2
opvar,E_x6,dv,f1,f2
opvar,D_z6,dv,h1,h2
opvar,E_x7,dv,f1,f2
opvar,D_z7,dv,h1,h2
opvar,E_x8,dv,f1,f2
opvar,D_z8,dv,h1,h2
opvar,E_x9,dv,f1,f2
opvar,D_z9,dv,h1,h2
opvar,E_x10,dv,f1,f2
opvar,D_z10,dv,h1,h2
opvar,E_x11,dv,f1,f2
opvar,D_z11,dv,h1,h2
opvar,E_x12,dv,f1,f2
opvar,D_z12,dv,h1,h2
opvar,E_x13,dv,f1,f2
opvar,D_z13,dv,h1,h2
opvar,E_x14,dv,f1,f2
opvar,D_z14,dv,h1,h2
opvar,E_x15,dv,f1,f2
opvar,D_z15,dv,h1,h2
opvar,E_x16,dv,f1,f2
opvar,D_z16,dv,h1,h2
opvar,E_x17,dv,f1,f2
opvar,D_z17,dv,h1,h2
opvar,E_x18,dv,f1,f2
opvar,D_z18,dv,h1,h2
opvar,E_x19,dv,f1,f2
opvar,D_z19,dv,h1,h2
opvar,E_x20,dv,f1,f2
opvar,D_z20,dv,h1,h2
opvar,E_x21,dv,f1,f2
opvar,D_z21,dv,h1,h2
opvar,E_x22,dv,f1,f2
opvar,D_z22,dv,h1,h2
opvar,E_x23,dv,f1,f2
opvar,D_z23,dv,h1,h2
opvar,E_x24,dv,f1,f2
opvar,D_z24,dv,h1,h2
!
*do,i,1,3,1
opvar,herwdq_OPT%i%,sv,0.95,1.05 !状态变量
*enddo
!
opvar,obj_func,obj !目标函数
!
opdata,,,
oploop,prep,proc,,all
opprnt,full
opkeep,on
!
optype,first
opfrst,q !优化计算
opexe
*get,n,opt,0,best
/output,OPT_LIST,lis
oplist,n
/output