这是全部命令流，只看优化分析部分就行，问题就出在优化设计部分

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