pred,off撸管会不会影响身高autots,on
点击联系发帖人
时间：2015-05-29 02:00
[转载]ANSYS经典例子---（2）
四层弹簧－质点模型的地震分析
包括频率分析，谱分析和时程分析
作者：陆新征，清华大学土木工程系
DAMPRATIO=0.02 !振型阻尼系数都为0.02
ET,1,BEAM188
KEYOPT,1,7,1
KEYOPT,1,8,1
ET,2,MASS21,,,4
MP,EX,1,210E9
MP,NUXY,1,0.3
MP,DENS,1,7800
MPDATA,DAMP,1,,DAMPRATIO
SECTYPE,3,BEAM,I
SECDATA,0.2,0.2,0.5,0.01,0.01,0.008
R,1,1.6E2,1.6E2
R,2,1.2E2,1.2E2
K,5,0,12,0
K,6,0,0,10000
LSEL,,,,1,4
LATT,1,0,1,,6,,3
LESIZE,ALL, , ,3 , , , , ,1
NUMMRG,ALL
NUMCMP,ALL
*SET,NT,1001
*SET,DT,0.02
*DIM,AC,,NT
*VREAD,AC(1),RECORD,TXT
/SOLU !模态分析
MODOPT,SUBSP,8
MXPAND,8, , ,1
!得到自振频率1
*GET,FREQ1,MODE,1,FREQ
/ESHAPE,1.0
ANTYPE,SPECTR !谱分析
SPOPT,SPRS
MDAMP,1,DAMPRATIO, , , , , ,
SVTYPE,2,1
FREQ,0.167,0.25,0.333,0.,0.667,1,1.25,1.667
SV,,0.........
FREQ,2.222,10,100000
SV,,1...3528
!按7度多遇地震，地震影响系数0.08 第一组III场地TG＝0.45
SRSS,0.01,DISP
/INP,,MCOM
NSORT,U,SUM
*GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形
!开始地震计算
/CONFIG,NRES,20000
ANTYPE,TRANS
TRNOPT,FULL
ALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926
BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)
*DO,I,1,1001
ACEL,AC(I),0,0
TIME,I*0.02
OUTRES,ALL,ALL
NSOL,2,20,U,X,
PLVAR,2, , , , , , , , , ,
STORE,MERGE
*GET,DMAX_TRANS,VARI,2,EXTREM,VMAX, , ! 时程分析最大位移
record.txt
悬臂梁地震分析
!定义单元类型
ET,1,BEAM3
R,1,0.09,0..3, !0.3*0.3m截面
MP,EX,1,2.06e11
MP,NUXY,1,0.3
MP,DAMP,1,0.02 !阻尼比0.02
MP,dens,1,7.85e3
!--------------------------------------------------
!定义节点单元和位移约束
*do,i,1,11
n,i,(i-1)*0.5,0,0 !5m长,分10段
*do,i,1,10
!--------------------------------------------------
!定义和读入时程曲线
NT=1000 !时程曲线有NT个点
DT=0.01 !时间间隔
*dim,ac,,NT
!读入数据,这个数据文件可以用excel等软件来写
/input,tianjin,txt !天津波东西向共10秒
!--------------------------------------------------
NSUBST,1, , ,1 !1个子步
OUTRES,ALL,1 !输出每个子步的结果
ANTYPE,TRANS !时程分析
*do,i,1,NT
ACEL,0,ac ,0
!--------------------------------------------------
NSOL,2,11,U,Y,
PLVAR,2,3, , ,
下面是文件tianjin.txt，这个文件可以用excel等软件制作或编写很短的程序就能写出来先
贴前500个点
ac( 1)= -0.
ac( 2)= -0.
ac( 3)= -0.
ac( 4)= -0.
ac( 5)= 0.
ac( 6)= 0.
ac( 7)= 0.
ac( 8)= 0.
ac( 9)= 0.
ac( 10)= 0.
ac( 11)= 0.
ac( 12)= 0.
ac( 13)= 0.
ac( 14)= 0.
ac( 15)= 0.
ac( 16)= 0.
ac( 17)= -0.
ac( 18)= -0.
ac( 19)= 0.
ac( 20)= 0.
ac( 21)= 0.
ac( 22)= 0.
ac( 23)= 0.
ac( 24)= 0.
ac( 25)= -0.
ac( 26)= -0.
ac( 27)= -0.
ac( 28)= -0.
ac( 29)= -0.
ac( 30)= -0.
ac( 31)= -0.
ac( 32)= -0.
ac( 33)= -0.
ac( 34)= -0.
ac( 35)= -0.
ac( 36)= -0.
ac( 37)= -0.
ac( 38)= -0.
ac( 39)= 0.
ac( 40)= 0.
ac( 41)= 0.
ac( 42)= 0.
ac( 43)= -0.
ac( 44)= -0.
ac( 45)= -0.
ac( 46)= -0.
ac( 47)= -0.
ac( 48)= -0.
ac( 49)= -0.
ac( 50)= -0.
ac( 51)= -0.
ac( 52)= 0.
ac( 53)= 0.
ac( 54)= 0.
ac( 55)= 0.
ac( 56)= 0.
ac( 57)= 0.
ac( 58)= -0.
ac( 59)= -0.
ac( 60)= -0.
ac( 61)= -0.
ac( 62)= 0.
ac( 63)= 0.
ac( 64)= 0.
ac( 65)= 0.
ac( 66)= 0.
ac( 67)= 0.
ac( 68)= 0.
ac( 69)= 0.
ac( 70)= 0.
ac( 71)= 0.
ac( 72)= 0.
ac( 73)= 0.
ac( 74)= 0.
ac( 75)= 0.
ac( 76)= -0.
ac( 77)= -0.
ac( 78)= -0.
ac( 79)= -0.
ac( 80)= -0.
ac( 81)= -0.
ac( 82)= -0.
ac( 83)= -0.
ac( 84)= -0.
ac( 85)= -0.
ac( 86)= -0.
ac( 87)= -0.
ac( 88)= -0.
ac( 89)= -0.
ac( 90)= -0.
ac( 91)= -0.
ac( 92)= -0.
ac( 93)= -0.
ac( 94)= -0.
ac( 95)= -0.
ac( 96)= -0.
ac( 97)= -0.
ac( 98)= -0.
ac( 99)= -0.
ac( 100)= -0.
ac( 101)= -0.
ac( 102)= -0.
ac( 103)= -0.
ac( 104)= -0.
ac( 105)= -0.
ac( 106)= -0.
ac( 107)= -0.
ac( 108)= -0.
ac( 109)= -0.
ac( 110)= -0.
ac( 111)= 0.
ac( 112)= 0.
ac( 113)= 0.
ac( 114)= 0.
ac( 115)= 0.
ac( 116)= 0.
ac( 117)= 0.
ac( 118)= 0.
ac( 119)= 0.
ac( 120)= 0.
ac( 121)= 0.
ac( 122)= 0.
ac( 123)= 0.
ac( 124)= -0.
ac( 125)= -0.
ac( 126)= -0.
ac( 127)= -0.
ac( 128)= -0.
ac( 129)= -0.
ac( 130)= -0.
ac( 131)= -0.
ac( 132)= 0.
ac( 133)= 0.
ac( 134)= 0.
ac( 135)= 0.
ac( 136)= 0.
ac( 137)= 0.
ac( 138)= 0.
ac( 139)= 0.
ac( 140)= 0.
ac( 141)= 0.
ac( 142)= 0.
ac( 143)= 0.
ac( 144)= 0.
ac( 145)= 0.
ac( 146)= 0.
ac( 147)= -0.
ac( 148)= -0.
ac( 149)= -0.
ac( 150)= -0.
ac( 151)= 0.
ac( 152)= 0.
ac( 153)= 0.
ac( 154)= 0.
ac( 155)= 0.
ac( 156)= 0.
ac( 157)= 0.
ac( 158)= 0.
ac( 159)= 0.
ac( 160)= 0.
ac( 161)= 0.
ac( 162)= -0.
ac( 163)= -0.
ac( 164)= -0.
ac( 165)= -0.
ac( 166)= -0.
ac( 167)= -0.
ac( 168)= -0.
ac( 169)= -0.
ac( 170)= -0.
ac( 171)= -0.
ac( 172)= -0.
ac( 173)= -0.
ac( 174)= -0.
ac( 175)= -0.
ac( 176)= -0.
ac( 177)= -0.
ac( 178)= -0.
ac( 179)= -0.
ac( 180)= -0.
ac( 181)= -0.
ac( 182)= 0.
ac( 183)= 0.
ac( 184)= 0.
ac( 185)= 0.
ac( 186)= 0.
ac( 187)= 0.
ac( 188)= 0.
ac( 189)= 0.
ac( 190)= -0.
ac( 191)= -0.
ac( 192)= -0.
ac( 193)= -0.
ac( 194)= -0.
ac( 195)= 0.
ac( 196)= 0.
ac( 197)= 0.
ac( 198)= 0.
ac( 199)= 0.
ac( 200)= 0.
ac( 201)= -0.
ac( 202)= -0.
ac( 203)= 0.
ac( 204)= 0.
ac( 205)= 0.
ac( 206)= 0.
ac( 207)= 0.
ac( 208)= -0.
ac( 209)= -0.
ac( 210)= -0.
ac( 211)= -0.
ac( 212)= -0.
ac( 213)= 0.
ac( 214)= 0.
ac( 215)= 0.
ac( 216)= 0.
ac( 217)= -0.
ac( 218)= -0.
ac( 219)= -0.
ac( 220)= -0.
ac( 221)= -0.
ac( 222)= -0.
ac( 223)= -0.
ac( 224)= -0.
ac( 225)= 0.
ac( 226)= 0.
ac( 227)= 0.
ac( 228)= 0.
ac( 229)= 0.
ac( 230)= 0.
ac( 231)= 0.
ac( 232)= 0.
ac( 233)= -0.
ac( 234)= -0.
ac( 235)= -0.
ac( 236)= -0.
ac( 237)= -0.
ac( 238)= -0.
ac( 239)= -0.
ac( 240)= 0.
ac( 241)= 0.
ac( 242)= 0.
ac( 243)= 0.
ac( 244)= -0.
ac( 245)= -0.
ac( 246)= -0.
ac( 247)= -0.
ac( 248)= 0.
ac( 249)= 0.
ac( 250)= 0.
ac( 251)= 0.
ac( 252)= 0.
ac( 253)= 0.
ac( 254)= 0.
ac( 255)= -0.
ac( 256)= -0.
ac( 257)= -0.
ac( 258)= -0.
ac( 259)= -0.
ac( 260)= -0.
ac( 261)= -0.
ac( 262)= -0.
ac( 263)= -0.
ac( 264)= -0.
ac( 265)= -0.
ac( 266)= -0.
ac( 267)= -0.
ac( 268)= -0.
ac( 269)= -0.
ac( 270)= -0.
ac( 271)= -0.
ac( 272)= -0.
ac( 273)= -0.
ac( 274)= -0.
ac( 275)= -0.
ac( 276)= -0.
ac( 277)= -0.
ac( 278)= 0.
ac( 279)= 0.
ac( 280)= 0.
ac( 281)= -0.
ac( 282)= -0.
ac( 283)= -0.
ac( 284)= -0.
ac( 285)= -0.
ac( 286)= -0.
ac( 287)= -0.
ac( 288)= -0.
ac( 289)= -0.
ac( 290)= 0.
ac( 291)= 0.
ac( 292)= 0.
ac( 293)= 0.
ac( 294)= 0.
ac( 295)= 0.
ac( 296)= 0.
ac( 297)= 0.
ac( 298)= 0.
ac( 299)= 0.
ac( 300)= 0.
ac( 301)= 0.
ac( 302)= 0.
ac( 303)= 0.
ac( 304)= 0.
ac( 305)= -0.
ac( 306)= -0.
ac( 307)= -0.
ac( 308)= -0.
ac( 309)= -0.
ac( 310)= -0.
ac( 311)= 0.
ac( 312)= 0.
ac( 313)= 0.
ac( 314)= 0.
ac( 315)= 0.
ac( 316)= 0.
ac( 317)= -0.
ac( 318)= -0.
ac( 319)= -0.
ac( 320)= -0.
ac( 321)= -0.
ac( 322)= 0.
ac( 323)= 0.
ac( 324)= 0.
ac( 325)= 0.
ac( 326)= 0.
ac( 327)= 0.
ac( 328)= 0.
ac( 329)= 0.
ac( 330)= -0.
ac( 331)= -0.
ac( 332)= -0.
ac( 333)= -0.
ac( 334)= -0.
ac( 335)= -0.
ac( 336)= 0.
ac( 337)= 0.
ac( 338)= 0.
ac( 339)= 0.
ac( 340)= 0.
ac( 341)= 0.
ac( 342)= 0.
ac( 343)= 0.
ac( 344)= 0.
ac( 345)= -0.
ac( 346)= -0.
ac( 347)= -0.
ac( 348)= -0.
ac( 349)= -0.
ac( 350)= -0.
ac( 351)= -0.
ac( 352)= -0.
ac( 353)= -0.
ac( 354)= 0.
ac( 355)= 0.
ac( 356)= 0.
ac( 357)= 0.
ac( 358)= 0.
ac( 359)= 0.
ac( 360)= 0.
ac( 361)= 0.
ac( 362)= 0.
ac( 363)= -0.
ac( 364)= -0.
ac( 365)= -0.
ac( 366)= -0.
ac( 367)= -0.
ac( 368)= 0.
ac( 369)= 0.
ac( 370)= 0.
ac( 371)= 0.
ac( 372)= 0.
ac( 373)= 0.
ac( 374)= 0.
ac( 375)= 0.
ac( 376)= 0.
ac( 377)= 0.
ac( 378)= 0.
ac( 379)= 0.
ac( 380)= 0.
ac( 381)= -0.
ac( 382)= -0.
ac( 383)= -0.
ac( 384)= -0.
ac( 385)= -0.
ac( 386)= 0.
ac( 387)= 0.
ac( 388)= 0.
ac( 389)= 0.
ac( 390)= 0.
ac( 391)= 0.
ac( 392)= 0.
ac( 393)= 0.
ac( 394)= 0.
ac( 395)= 0.
ac( 396)= -0.
ac( 397)= -0.
ac( 398)= -0.
ac( 399)= -0.
ac( 400)= -0.
ac( 401)= -0.
ac( 402)= -0.
ac( 403)= -0.
ac( 404)= 0.
ac( 405)= 0.
ac( 406)= 0.
ac( 407)= 0.
ac( 408)= 0.
ac( 409)= 0.
ac( 410)= 0.
ac( 411)= 0.
ac( 412)= 0.
ac( 413)= 0.
ac( 414)= 0.
ac( 415)= 0.
ac( 416)= 0.
ac( 417)= 0.
ac( 418)= -0.
ac( 419)= -0.
ac( 420)= -0.
ac( 421)= -0.
ac( 422)= -0.
ac( 423)= 0.
ac( 424)= 0.
ac( 425)= 0.
ac( 426)= 0.
ac( 427)= 0.
ac( 428)= 0.
ac( 429)= 0.
ac( 430)= -0.
ac( 431)= -0.
ac( 432)= -0.
ac( 433)= -0.
ac( 434)= -0.
ac( 435)= 0.
ac( 436)= 0.
ac( 437)= 0.
ac( 438)= 0.
ac( 439)= 0.
ac( 440)= 0.
ac( 441)= 0.
ac( 442)= -0.
ac( 443)= -0.
ac( 444)= -0.
ac( 445)= -0.
ac( 446)= -0.
ac( 447)= -0.
ac( 448)= -0.
ac( 449)= -0.
ac( 450)= 0.
ac( 451)= 0.
ac( 452)= 0.
ac( 453)= -0.
ac( 454)= -0.
ac( 455)= -0.
ac( 456)= -0.
ac( 457)= -0.
ac( 458)= 0.
ac( 459)= 0.
ac( 460)= 0.
ac( 461)= 0.
ac( 462)= 0.
ac( 463)= -0.
ac( 464)= -0.
ac( 465)= -0.
ac( 466)= 0.
ac( 467)= 0.
ac( 468)= 0.
ac( 469)= 0.
ac( 470)= 0.
ac( 471)= 0.
ac( 472)= 0.
ac( 473)= 0.
ac( 474)= 0.
ac( 475)= 0.
ac( 476)= 0.
ac( 477)= 0.
ac( 478)= -0.
ac( 479)= -0.
ac( 480)= -0.
ac( 481)= -0.
ac( 482)= 0.
ac( 483)= 0.
ac( 484)= 0.
ac( 485)= 0.
ac( 486)= 0.
ac( 487)= 0.
ac( 488)= 0.
ac( 489)= 0.
ac( 490)= 0.
ac( 491)= -0.
ac( 492)= -0.
ac( 493)= -0.
ac( 494)= -0.
ac( 495)= -0.
ac( 496)= -0.
ac( 497)= -0.
ac( 498)= -0.
ac( 499)= -0.
ac( 500)= -0.
再贴后500个点
ac( 501)= -0.
ac( 502)= -0.
ac( 503)= -0.
ac( 504)= -0.
ac( 505)= -0.
ac( 506)= -0.
ac( 507)= -0.
ac( 508)= -0.
ac( 509)= -0.
ac( 510)= -0.
ac( 511)= 0.
ac( 512)= 0.
ac( 513)= 0.
ac( 514)= 0.
ac( 515)= 0.
ac( 516)= 0.
ac( 517)= 0.
ac( 518)= -0.
ac( 519)= -0.
ac( 520)= -0.
ac( 521)= 0.
ac( 522)= 0.
ac( 523)= 0.
ac( 524)= 0.
ac( 525)= 0.
ac( 526)= -0.
ac( 527)= -0.
ac( 528)= -0.
ac( 529)= -0.
ac( 530)= -0.
ac( 531)= -0.
ac( 532)= -0.
ac( 533)= -0.
ac( 534)= -0.
ac( 535)= -0.
ac( 536)= -0.
ac( 537)= -0.
ac( 538)= -0.
ac( 539)= -0.
ac( 540)= -0.
ac( 541)= -0.
ac( 542)= -0.
ac( 543)= -0.
ac( 544)= -0.
ac( 545)= -0.
ac( 546)= -0.
ac( 547)= -0.
ac( 548)= -0.
ac( 549)= -0.
ac( 550)= -0.
ac( 551)= -0.
ac( 552)= 0.
ac( 553)= 0.
ac( 554)= 0.
ac( 555)= 0.
ac( 556)= 0.
ac( 557)= 0.
ac( 558)= 0.
ac( 559)= 0.
ac( 560)= 0.
ac( 561)= -0.
ac( 562)= -0.
ac( 563)= -0.
ac( 564)= -0.
ac( 565)= -0.
ac( 566)= -0.
ac( 567)= -0.
ac( 568)= -0.
ac( 569)= -0.
ac( 570)= -0.
ac( 571)= -0.
ac( 572)= 0.
ac( 573)= 0.
ac( 574)= 0.
ac( 575)= 0.
ac( 576)= -0.
ac( 577)= -0.
ac( 578)= -0.
ac( 579)= -0.
ac( 580)= -0.
ac( 581)= -0.
ac( 582)= -0.
ac( 583)= -0.
ac( 584)= -0.
ac( 585)= -0.
ac( 586)= -0.
ac( 587)= -0.
ac( 588)= -0.
ac( 589)= -0.
ac( 590)= -0.
ac( 591)= -0.
ac( 592)= 0.
ac( 593)= 0.
ac( 594)= 0.
ac( 595)= 0.
ac( 596)= 0.
ac( 597)= 0.
ac( 598)= -0.
ac( 599)= -0.
ac( 600)= -0.
ac( 601)= -0.
ac( 602)= -0.
ac( 603)= -0.
ac( 604)= -0.
ac( 605)= 0.
ac( 606)= 0.
ac( 607)= 0.
ac( 608)= 0.
ac( 609)= 0.
ac( 610)= 0.
ac( 611)= 0.
ac( 612)= 0.
ac( 613)= 0.
ac( 614)= 0.
ac( 615)= 0.
ac( 616)= 0.
ac( 617)= -0.
ac( 618)= -0.
ac( 619)= -0.
ac( 620)= -0.
ac( 621)= -0.
ac( 622)= 0.
ac( 623)= 0.
ac( 624)= 0.
ac( 625)= 0.
ac( 626)= 0.
ac( 627)= 0.
ac( 628)= -0.
ac( 629)= -0.
ac( 630)= -0.
ac( 631)= -0.
ac( 632)= -0.
ac( 633)= -0.
ac( 634)= -0.
ac( 635)= -0.
ac( 636)= -0.
ac( 637)= -0.
ac( 638)= -0.
ac( 639)= -0.
ac( 640)= -0.
ac( 641)= -0.
ac( 642)= -0.
ac( 643)= -0.
ac( 644)= -0.
ac( 645)= -0.
ac( 646)= -0.
ac( 647)= -0.
ac( 648)= -0.
ac( 649)= 0.
ac( 650)= 0.
ac( 651)= 0.
ac( 652)= 0.
ac( 653)= 0.
ac( 654)= 0.
ac( 655)= 0.
ac( 656)= 0.
ac( 657)= 0.
ac( 658)= 0.
ac( 659)= 0.
ac( 660)= 0.
ac( 661)= -0.
ac( 662)= -0.
ac( 663)= -0.
ac( 664)= -0.
ac( 665)= -0.
ac( 666)= -0.
ac( 667)= -0.
ac( 668)= -0.
ac( 669)= -0.
ac( 670)= -0.
ac( 671)= -0.
ac( 672)= -0.
ac( 673)= -0.
ac( 674)= 0.
ac( 675)= 0.
ac( 676)= 0.
ac( 677)= 0.
ac( 678)= 0.
ac( 679)= 0.
ac( 680)= 0.
ac( 681)= 0.
ac( 682)= 0.
ac( 683)= 0.
ac( 684)= 0.
ac( 685)= 0.
ac( 686)= 0.
ac( 687)= 0.
ac( 688)= 0.
ac( 689)= 0.
ac( 690)= 0.
ac( 691)= 0.
ac( 692)= 0.
ac( 693)= 0.
ac( 694)= 1.
ac( 695)= 1.
ac( 696)= 1.
ac( 697)= 1.
ac( 698)= 1.
ac( 699)= 1.
ac( 700)= 1.
ac( 701)= 1.
ac( 702)= 1.
ac( 703)= 1.
ac( 704)= 1.
ac( 705)= 1.
ac( 706)= 1.
ac( 707)= 1.
ac( 708)= 1.
ac( 709)= 1.
ac( 710)= 1.
ac( 711)= 1.
ac( 712)= 1.
ac( 713)= 1.
ac( 714)= 1.
ac( 715)= 1.
ac( 716)= 1.
ac( 717)= 1.
ac( 718)= 1.
ac( 719)= 1.
ac( 720)= 1.
ac( 721)= 1.
ac( 722)= 1.
ac( 723)= 1.
ac( 724)= 1.
ac( 725)= 1.
ac( 726)= 1.
ac( 727)= 1.
ac( 728)= 0.
ac( 729)= 0.
ac( 730)= 0.
ac( 731)= 0.
ac( 732)= 0.
ac( 733)= -0.
ac( 734)= -0.
ac( 735)= -0.
ac( 736)= -0.
ac( 737)= -0.
ac( 738)= -0.
ac( 739)= -0.
ac( 740)= -0.
ac( 741)= -0.
ac( 742)= -0.
ac( 743)= -1.
ac( 744)= -1.
ac( 745)= -1.
ac( 746)= -1.
ac( 747)= -1.
ac( 748)= -1.
ac( 749)= -2.
ac( 750)= -2.
ac( 751)= -2.
ac( 752)= -2.
ac( 753)= -2.
ac( 754)= -2.
ac( 755)= -3.
ac( 756)= -3.
ac( 757)= -3.
ac( 758)= -3.
ac( 759)= -3.
ac( 760)= -3.
ac( 761)= -3.
ac( 762)= -2.
ac( 763)= -2.
ac( 764)= -2.
ac( 765)= -2.
ac( 766)= -2.
ac( 767)= -2.
ac( 768)= -2.
ac( 769)= -2.
ac( 770)= -2.
ac( 771)= -2.
ac( 772)= -2.
ac( 773)= -2.
ac( 774)= -2.
ac( 775)= -2.
ac( 776)= -2.
ac( 777)= -2.
ac( 778)= -2.
ac( 779)= -2.
ac( 780)= -2.
ac( 781)= -2.
ac( 782)= -2.
ac( 783)= -2.
ac( 784)= -2.
ac( 785)= -1.
ac( 786)= -1.
ac( 787)= -0.
ac( 788)= -0.
ac( 789)= -0.
ac( 790)= 0.
ac( 791)= 0.
ac( 792)= 0.
ac( 793)= 0.
ac( 794)= 0.
ac( 795)= 0.
ac( 796)= 1.
ac( 797)= 1.
ac( 798)= 1.
ac( 799)= 1.
ac( 800)= 1.
ac( 801)= 1.
ac( 802)= 1.
ac( 803)= 0.
ac( 804)= 0.
ac( 805)= 0.
ac( 806)= 0.
ac( 807)= 0.
ac( 808)= 0.
ac( 809)= 0.
ac( 810)= 0.
ac( 811)= 0.
ac( 812)= 0.
ac( 813)= 0.
ac( 814)= 0.
ac( 815)= 1.
ac( 816)= 1.
ac( 817)= 1.
ac( 818)= 1.
ac( 819)= 1.
ac( 820)= 1.
ac( 821)= 2.
ac( 822)= 2.
ac( 823)= 2.
ac( 824)= 2.
ac( 825)= 1.
ac( 826)= 1.
ac( 827)= 1.
ac( 828)= 1.
ac( 829)= 1.
ac( 830)= 0.
ac( 831)= 0.
ac( 832)= 1.
ac( 833)= 1.
ac( 834)= 1.
ac( 835)= 1.
ac( 836)= 1.
ac( 837)= 1.
ac( 838)= 1.
ac( 839)= 0.
ac( 840)= 0.
ac( 841)= 0.
ac( 842)= 0.
ac( 843)= 1.
ac( 844)= 1.
ac( 845)= 1.
ac( 846)= 1.
ac( 847)= 1.
ac( 848)= 1.
ac( 849)= 1.
ac( 850)= 1.
ac( 851)= 0.
ac( 852)= 0.
ac( 853)= 0.
ac( 854)= 0.
ac( 855)= 0.
ac( 856)= -0.
ac( 857)= -0.
ac( 858)= -0.
ac( 859)= -1.
ac( 860)= -1.
ac( 861)= -1.
ac( 862)= -1.
ac( 863)= -1.
ac( 864)= -1.
ac( 865)= -0.
ac( 866)= -0.
ac( 867)= 0.
ac( 868)= 0.
ac( 869)= 0.
ac( 870)= 0.
ac( 871)= 0.
ac( 872)= 0.
ac( 873)= 0.
ac( 874)= 0.
ac( 875)= 0.
ac( 876)= -0.
ac( 877)= -0.
ac( 878)= -0.
ac( 879)= -1.
ac( 880)= -1.
ac( 881)= -1.
ac( 882)= -1.
ac( 883)= -0.
ac( 884)= -0.
ac( 885)= -0.
ac( 886)= 0.
ac( 887)= 0.
ac( 888)= 1.
ac( 889)= 1.
ac( 890)= 1.
ac( 891)= 1.
ac( 892)= 1.
ac( 893)= 1.
ac( 894)= 1.
ac( 895)= 0.
ac( 896)= 0.
ac( 897)= 0.
ac( 898)= -0.
ac( 899)= -0.
ac( 900)= -0.
ac( 901)= -0.
ac( 902)= -1.
ac( 903)= -1.
ac( 904)= -1.
ac( 905)= -1.
ac( 906)= -1.
ac( 907)= -1.
ac( 908)= -1.
ac( 909)= -1.
ac( 910)= -1.
ac( 911)= -1.
ac( 912)= -1.
ac( 913)= -1.
ac( 914)= -1.
ac( 915)= -1.
ac( 916)= -1.
ac( 917)= -0.
ac( 918)= -0.
ac( 919)= -0.
ac( 920)= -0.
ac( 921)= -0.
ac( 922)= 0.
ac( 923)= 0.
ac( 924)= 0.
ac( 925)= 0.
ac( 926)= 0.
ac( 927)= 0.
ac( 928)= 0.
ac( 929)= 0.
ac( 930)= 0.
ac( 931)= 0.
ac( 932)= 0.
ac( 933)= 0.
ac( 934)= 0.
ac( 935)= 0.
ac( 936)= 0.
ac( 937)= 0.
ac( 938)= 0.
ac( 939)= 0.
ac( 940)= 0.
ac( 941)= -0.
ac( 942)= -0.
ac( 943)= -0.
ac( 944)= -0.
ac( 945)= -0.
ac( 946)= 0.
ac( 947)= 0.
ac( 948)= 0.
ac( 949)= 0.
ac( 950)= 0.
ac( 951)= 0.
ac( 952)= 0.
ac( 953)= -0.
ac( 954)= -0.
ac( 955)= -0.
ac( 956)= -0.
ac( 957)= -0.
ac( 958)= -0.
ac( 959)= -0.
ac( 960)= -0.
ac( 961)= -0.
ac( 962)= -0.
ac( 963)= -0.
ac( 964)= 0.
ac( 965)= 0.
ac( 966)= 0.
ac( 967)= 0.
ac( 968)= 1.
ac( 969)= 1.
ac( 970)= 1.
ac( 971)= 1.
ac( 972)= 1.
ac( 973)= 1.
ac( 974)= 1.
ac( 975)= 1.
ac( 976)= 1.
ac( 977)= 1.
ac( 978)= 0.
ac( 979)= 0.
ac( 980)= 0.
ac( 981)= -0.
ac( 982)= -0.
ac( 983)= -0.
ac( 984)= -0.
ac( 985)= -0.
ac( 986)= -0.
ac( 987)= -0.
ac( 988)= -0.
ac( 989)= 0.
ac( 990)= 0.
ac( 991)= 0.
ac( 992)= 0.
ac( 993)= 0.
ac( 994)= 0.
ac( 995)= 0.
ac( 996)= 0.
ac( 997)= 0.
ac( 998)= 0.
ac( 999)= 0.
ac(1000)= 0.
54单元描述变截面梁的例子
作者：清华大学土木系，陆新征
A_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.6667,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.6667,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,0
MPDATA,EX,1,,200e3
MPDATA,PRXY,1,,.2
N,1,0,,,,,,
N,2,5,,,,,,
N,3,10,,,,,,
变截面梁实例
et,1,beam188
keyopt,1,8,2
keyopt,1,9,2
mp,ex,1,3.5e10
mp,dens,1,2549
mp,prxy,1,0.1667
sectype,1,beam,rect ! define cross section at first end point
sectype,2,beam,rect ! define cross section at far end
secdata,3,0.5
sectype,3,taper ! new Section ID for tapered beam analysis
secdata,1,0.0,0.0 ! section 1 at location (0,0,0)
secdata,2,0.0,,-20.0 ! section 2 at location (0,-20,0)
k, 1, -2.7 ,0, -27.85
k, 2, -2.7 ,0, -19
k, 3, -2.7 ,0, -13
lstr, 1, 2
lstr, 2, 3
lsel,s,,,1
latt,1,,1,,,,1
lesize,all,,,1
lsel,s,,,2
latt,1,,1,,,,3
lesize,all,,,1
/view,,1,2,3
拱桥浇筑过程分析-单元生死应用实例
et,1,BEAM3
mp,ex,1,2.1e11
mp,prxy,1,.3
mp,dens,1,7800
r,1,25.48e-4 ,,1.5
K,2,10,5,,
K,3,20,0,,
LARC, 1, 3, 2
LESIZE,ALL, , ,20, ,1, , ,1,
n,22,3.3986,6
n,23,7.695,6
n,24,12.305,6
n,25,16.601,6
nropt,full,,on
esel,s,,,25,27,1
antype,,rest,,,
f,12,fy,-1000
nSOL,2,12,u,y,
PLVAR,1, , , , , , , , , ,
简支梁实体与预应力钢筋分析实例
/COM, Structural
egjx=2e5 !Ey
agjx=140 !单根钢绞线面积
ehnt=4e4 !Eh
xzxs=1.0e-5 !线胀系数
yjl=200000 !定义预加力
et,1,link8 !定义link8单元
et,2,solid95 !定义solid95单元
r,1,agjx !定义link8单元的面积
r,2 !定义第2种实常数
mp,ex,1,egjx !定义link8单元的弹性模量
mp,prxy,1,0.3 !定义link8单元的泊松系数
mp,alpx,1,1.0e-5 !定义线膨胀系数
mp,ex,2,ehnt !定义solid95单元的弹性模量
mp,prxy,2,0.3 !定义solid95单元的泊松系数
blc4, , ,100,200,3000 !定义梁体
/view,1,1,1,1 !定义ISO查看
vplot !绘制梁体
kwpave,6 !工作平面移动到关键点6
wpoff,-30 !工作平面移动-30mm(X)
wprot,0,0,90 !工作平面旋转
vsbw,1 !分割梁体
wpoff,0,0,-40 !工作平面移动-40mm(Z)
vsbw,2 !分割梁体
wpoff,0,40 !工作平面移动40mm(Y)
wprot,0,90 !工作平面旋转
vsbw,all !分割梁体
wpstyl !关闭工作平面显示
nummrg,all,,,,low !整理
numcmp,all !压缩编号
esize,30 !定义网分时边长控制
lsel,s,,,28,38,10 !定义line28和38为新的选择集
latt,1,1,1 !定义选择集的属性
lmesh,all !对线划分单元
allsel,all !新的选择集为所有的实体
gplot !绘制所有的实体
vsel,s,,,all !定义所有体为选择集
vatt,2,2,2 !定义选择集的属性
mshape,0,3d !将体划分单元的形状定位HEX
mshkey,1 !采用MAPPED划分器
vmesh,all !对体进行划分单元
dl,3,,all !对线line7施加约束(UX,UY,UZ)
dl,16,,all !对线line31施加约束(UX,UY,UZ)
dl,23,,all !对线line23施加约束(UX,UY,UZ)
dl,2,,uy !对线line4施加约束(UY)
dl,15,,uy !对线line30施加约束(UY)
dl,22,,uy !对线line23施加约束(UY)
dk,2,,,,,ux,uy !对关键点2约束(UX,UY)
bfl,28,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度
bfl,38,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度
solve !求解
plnsol,s,z,0,1 !绘制Z方向的应力
etable,sigi,ls,1 !定义钢筋单元数据表
plls,sigi,sigi,1 !绘制上述应力
!/exit,nosav
简单的二维焊接分析－单元生死实例
下面的命令流进行的是一个简单的二维焊接分析, 利用ANSYS单元生死和热-结构耦合分析功能进行焊接过程仿真,
计算焊接过程中的温度分布和应力分布以及冷却后的焊缝残余应力。
/title,Weld Analysis by "Element Birth and Death"
/unit,si ! 采用国际单位制
et,1,13,4 ! 13号二维耦合单元, 同时具有温度和位移自由度
! 1号材料是钢
! 2号材料是铝
! 3号材料是铜
! 铝是本次分析中的焊料, 它将钢结构部分和铜结构部分焊接起来
! 下面是在几个温度点下, 各材料的弹性模量
mptemp,1,20,500,00
mpdata,ex,1,1,1.93e11,1.50e11,0.70e11,0.10e11,0.01e11
mpdata,ex,2,1,1.02e11,0.50e11,0.08e11,0.001e11,0.0001e11
mpdata,ex,3,1,1.17e11,0.90e11,0.30e11,0.05e11,0.005e11
! 假设各材料都是双线性随动硬化弹塑性本构关系
! 下面是各材料在各温度点下的屈服应力和屈服后的弹性模量
tb,bkin,1,5
tbtemp,20,1
tbdata,1,.193e11
tbtemp,500,2
tbdata,1, 933e6,0.150e11
tbtemp,1000,3
tbdata,1, 435e6,0.070e11
tbtemp,1500,4
tbdata,1, 70e6,0.010e11
tbtemp,2000,5
tbdata,1, 7e6,0.001e11
tb,bkin,2,5
tbtemp,20,1
tbdata,1,800e6,0.102e11
tbtemp,500,2
tbdata,1,400e6,0.050e11
tbtemp,1000,3
tbdata,1, 70e6,0.008e11
tbtemp,1500,4
tbdata,1, 1e6,0.0001e11
tbtemp,2000,5
tbdata,1,0.1e6,0.00001e11
tb,bkin,3,5
tbtemp,20,1
tbdata,1,900e6,0.117e11
tbtemp,500,2
tbdata,1,700e6,0.090e11
tbtemp,1000,3
tbdata,1,230e6,0.030e11
tbtemp,1500,4
tbdata,1, 40e6,0.005e11
tbtemp,2000,5
tbdata,1, 4e6,0.0005e11
! 材料密度(假设为常值)
mp,dens,1,8030
mp,dens,2,4850
mp,dens,3,8900
! 热膨胀系数(假设为常值)
mp,alpx,1,1.78e-5
mp,alpx,2,9.36e-6
mp,alpx,3,1.66e-5
! 泊松比(假设为常值)
mp,nuxy,1,0.29
mp,nuxy,2,0.30
mp,nuxy,3,0.30
! 热传导系数(假设为常值)
mp,kxx,1,16.3
mp,kxx,2,7.44
mp,kxx,3,393
! 比热(假设为常值)
mp,c,1,502
mp,c,2,544
mp,c,3,385
! 热膨胀系数(假设为常值)
! 由于该13号单元还有磁自由度, 此处假设一磁特性, 但并不影响我们所关心的结果
mp,murx,1,1
mp,murx,2,1
mp,murx,3,1
! 假设焊料(铝)焊上去后的初始温度是;
mp,reft,1,20
mp,reft,2,1500
mp,reft,3,20
! 下面建立几何模型
k,2,0.5,0,0
k,4,0,0.3,0
k,5,0.35,0.3,0
k,6,0.65,0.3,0
k,7,1,0.3,0
! 划分单元
esize,0.025
esize,0.05
antype,4 ! 瞬态分析
trnopt,full
! 在模型的左边界加位移约束
*get,minx,node,,mnloc,x
nsel,s,loc,x,minx
d,all,ux,0
*get,miny,node,,mnloc,y
nsel,r,loc,y,miny
d,all,uy,0
! 假设在模型的左右边界处温度始终保持在20℃, 其它边界条件如对流、辐射等都可
! 施加，此处因只是示意而已，故只施加恒温边界条件。
*get,minx,node,,mnloc,x
nsel,s,loc,x,minx
d,all,temp,20
*get,maxx,node,,mxloc,x
nsel,s,loc,x,maxx
d,all,temp,20
! 由于第二个面是焊料所在区域，因此，首先将该区域的单元"死"掉
*get,emax,elem,,num,max
asel,s,area,,nna
*get,nse,elem,,count
*dim,ne,,nse
*dim,nex,,nse
*dim,ney,,nse
*dim,neorder,,nse
! 下面的DO循环用于将焊料区的单元按其形心Y坐标排序，以便后面
! 模拟焊料由下向上逐步"生长"的过程
*do,i1,1,nse
esel,u,elem,,mine
*get,nse1,elem,,count
*do,i,1,emax
*if,esel(i),eq,1,then
*do,i,1,nse1
*get,ney(i),elem,ne(i),cent,y
*get,nex(i),elem,ne(i),cent,x
*do,i,1,nse1
*if,ney(i),lt,miny,then
miny=ney(i)
minx=nex(i)
mine=ne(i)
*if,ney(i),eq,miny,then
*if,nex(i),lt,minx,then
miny=ney(i)
minx=nex(i)
mine=ne(i)
neorder(i1)=mine
max_tem=1500 ! 按照前面假设，焊料的初始温度为;
dt1=1e-3 ! 用于建立初始条件的一个很小的时间段
dt=5 ! 焊接一个单元所需的时间
t=0 ! 起始时间
*do,i,1,nse
ekill,neorder(i)
esel,s,live
allsel,all
outres,all,all
ic,all,temp,20
timint,0,struct
timint,1,therm
timint,0,mag
tintp,0.005,,,1,0.5,0.2
*do,i,1,nse
ealive,neorder(i)
esel,s,live
! 下面的求解用于建立温度的初始条件
d,nelem(neorder(i),j),temp,max_tem
! 下面的求解用于保证初始的升温速度为零
! 下面的步骤用于求解温度分布
ddele,nelem(neorder(i),j),temp
t=t+dt-2*dt1
nsubst,nsub1
nsubst,nsub2
! 下面的一系列命令用于生成应力的动画文件
/cont,1,15,0,,1200e6
/dscale,1,1.0
avprin,0,0
/seg,multi,stress,0.1
*do,i,1,nse
esel,u,elem,,neorder(i)
*do,i,1,nse
esel,a,elem,,neorder(i)
set,(i-1)*3+1,1
plnsol,s,eqv
*do,j,1,nsub1
set,(i-1)*3+3,j
plnsol,s,eqv
*do,i,1,nsub2
set,(nse-1)*3+4,i
plnsol,s,eqv
/seg,off,stress,0.1
/anfile,save,stress,avi
! 下面的一系列命令用于生成温度的动画文件
/cont,1,15,0,0
/dscale,1,1.0
avprin,0,0
/seg,multi,temp,0.1
*do,i,1,nse
esel,u,elem,,neorder(i)
*do,i,1,nse
esel,a,elem,,neorder(i)
set,(i-1)*3+1,1
plnsol,temp
*do,j,1,nsub1
set,(i-1)*3+3,j
plnsol,temp
*do,i,1,nsub2
set,(nse-1)*3+4,i
plnsol,temp
/seg,off,temp,0.1
/anfile,save,temp,avi
隧道开挖（三维）的命令流
/FILENAME,'TUNNEL'
/CONFIG,NRES,1E6
!--------------------------------------------------------------------------------------------------------------------!
TM_YS=16*1E9&&&&&&
MU_YS=0.24&&&&&&&&
C_YS=1.80*1E6&&&&&&
FI_YS=43&&&&&&&&&&
MD_YS=0.027*1E6/10
TM_GXW=26*1E9&&&&&&
MU_GXW=0.21&&&&&&&&
C_GXW=2.4*1E6&&&&&&
FI_GXW=50&&&&&&&&&&
MD_GXW=0./10
TM_HNT=22*1E9&&&&&&
MU_HNT=0.20&&&&&&&&
C_HNT=1.73*1E6&&&&&&
FI_HNT=49&&&&&&&&&&
MD_HNT=0.024*1E6/10
TM_MG=TM_YS&&&&&&
MU_MG=MU_YS&&&&&&&&&&
MD_MG=MD_YS&&&&&&&&&&
NT_MG=90*1E3&&&&&&&&&&
C_MG=C_YS*1.10&&&&&&&&
FI_MG=FI_YS+5&&&&&&&&
ZH_MG=NT_MG/(1.2*1*1)
!----------------------------------------------------------------------------------------------------------------------!
R_HJ=6.230&&&&&&&&&&
H_ZQ=2.670&&&&&&&&&&
L_ZTZ=40&&&&&&&&&&&&
L_YTZ=40&&&&&&&&&&&&
Y_YJMIN=400&&&&&&&&&&
Y_YJMAX=415&&&&&&&&&&
L_YJ=Y_YJMAX-Y_YJMIN
L_MG=3.5&&&&&&&&&&&&
JJ_MG=5.0&&&&&&&&&&&&
PJ_MG=5.0&&&&&&&&&&
H_CENDB=1.47&&&&&&&&
PD=2/1000&&&&&&&&&&&&
W_GXW=0.2&&&&&&&&&&&&
H_JS=831&&&&&&&&&&&&
DMJL_SD=40&&&&&&&&&&
H_CTOB=1.47&&&&&&&&&&
*DIM,CEN_ZSD,ARRAY,2,1
*DIM,CEN_YSD,ARRAY,2,1
CEN_ZSD(1,1)=-DMJL_SD/2,880
CEN_YSD(1,1)=DMJL_SD/2,880 !+H_CTOB
H_ZQ=2.67&&&&&&&&&&&&&&&&
FLAG_RD=3&&&&&&&&&&
JL_RD=5&&&&&&&&&&&&&&
*DIM,Z_RDCC,ARRAY,3,2
Z_RDCC(1,1)=3,3,3
Z_RDCC(1,2)=CEN_ZSD(1,1),CEN_ZSD(2,1)+R_HJ+JL_RD+Z_RDCC(2,1)/2,L_ZTZ+7.5
*DIM,Y_RDCC,ARRAY,3,2
Y_RDCC(1,1)=3,3,3
Y_RDCC(1,2)=CEN_YSD(1,1),CEN_YSD(2,1)+R_HJ+JL_RD+Y_RDCC(2,1)/2,L_ZTZ+7.5
*DIM,DBKZDX,ARRAY,9,1
*DIM,DBKZDZ,ARRAY,(Y_YJMAX-Y_YJMIN)/PJ_MG+3,1
DBKZDX(1,1)=-150,-64,-38,-20,0,20,42,64,150&&
DBKZDZ(1,1)=0,L_ZTZ
DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1)=L_ZTZ+Y_YJMAX-Y_YJMIN+L_YTZ
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG
&DBKZDZ(I+2,1)=L_ZTZ+PJ_MG*I
*DIM,DBKZDGC,ARRAY,9,(Y_YJMAX-Y_YJMIN)/PJ_MG+3
DBKZDGC(1,1)=929,915,915,920,925,930,937,940,940
DBKZDGC(1,2)=929,915,915,920,925,930,937,940,940
DBKZDGC(1,3)=929,915,915,920,925,930,937,940,940
DBKZDGC(1,4)=929,915,915,920,925,930,937,940,940
DBKZDGC(1,5)=929,915,915,920,925,930,937,940,940
DBKZDGC(1,6)=929,915,915,920,925,930,937,940,940
!*************************************************************************************!
!----------------------------------------------------------------------------------------------------------------------!
KEYOPT,2,5,2
MP,EX,7,TM_YS
MP,PRXY,7,MU_YS
MP,DENS,7,MD_YS
MP,EX,1,TM_YS
MP,PRXY,1,MU_YS
MP,DENS,1,MD_YS
MP,EX,2,TM_GXW
MP,PRXY,2,MU_GXW
MP,DENS,2,MD_GXW
MP,EX,3,TM_MG
MP,PRXY,3,MU_MG
MP,DENS,3,MD_MG
MP,EX,4,TM_HNT
MP,PRXY,4,MU_HNT
MP,DENS,4,MD_HNT
MP,EX,5,TM_YS
MP,PRXY,5,MU_YS
MP,DENS,5,MD_YS
MP,EX,6,TM_YS
MP,PRXY,6,MU_YS
TB,DP,1,,,&&
TBMODIF,1,1,C_YS&&
TBMODIF,1,2,FI_YS&&&&
TB,DP,2,,,&&
TBMODIF,1,1,C_GXW&&
TBMODIF,1,2,FI_GXW&&&&
TB,DP,3,,,&&
TBMODIF,1,1,C_MG&&
TBMODIF,1,2,FI_MG&&&&&&
TB,DP,4,,,&&
TBMODIF,1,1,C_HNT&&
TBMODIF,1,2,FI_HNT&&
TB,DP,5,,,&&
TBMODIF,1,1,C_YS&&
TBMODIF,1,2,FI_YS&&
TB,DP,6,,,&&
TBMODIF,1,1,C_YS&&
TBMODIF,1,2,FI_YS
!**************************************************!
MAT_YS=1&&
MAT_MGKW=6
MAT_TXYS=1&&
*DO,I,1,9-1
&*DO,J,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1&
&&K,KBEGIN+1,DBKZDX(I,1),H_JS,DBKZDZ(J,1)
&&K,KBEGIN+2,DBKZDX(I,1),DBKZDGC(I,J),DBKZDZ(J,1)
&&K,KBEGIN+3,DBKZDX(I,1),DBKZDGC(I,J+1),DBKZDZ(J+1,1)
&&K,KBEGIN+4,DBKZDX(I,1),H_JS,DBKZDZ(J+1,1)
&&K,KBEGIN+5,DBKZDX(I+1,1),H_JS,DBKZDZ(J,1)
&&K,KBEGIN+6,DBKZDX(I+1,1),DBKZDGC(I+1,J),DBKZDZ(J,1)
&&K,KBEGIN+7,DBKZDX(I+1,1),DBKZDGC(I+1,J+1),DBKZDZ(J+1,1)
&&K,KBEGIN+8,DBKZDX(I+1,1),H_JS,DBKZDZ(J+1,1)
&&V,KBEGIN+1,KBEGIN+2,KBEGIN+3,KBEGIN+4,KBEGIN+5,KBEGIN+6,KBEGIN+7,KBEGIN+8
&&KBEGIN=KBEGIN+8
NUMCMP,ALL
VATT,MAT_YS&&
CYL4,CEN_ZSD(1,1),CEN_ZSD(2,1),R_HJ+L_MG,180,,,L_ZTZ+L_YTZ+Y_YJMAX-Y_YJMIN
CYL4,CEN_YSD(1,1),CEN_YSD(2,1),R_HJ+L_MG,180,,,L_ZTZ+L_YTZ+Y_YJMAX-Y_YJMIN
BLOCK,CEN_ZSD(1,1)-R_HJ-L_MG,CEN_ZSD(1,1)+R_HJ+L_MG,CEN_ZSD(2,1)-H_ZQ,CEN_ZSD(2,1),DBKZDZ(1,1),DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1)
BLOCK,CEN_YSD(1,1)-R_HJ-L_MG,CEN_YSD(1,1)+R_HJ+L_MG,CEN_YSD(2,1)-H_ZQ,CEN_YSD(2,1),DBKZDZ(1,1),DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1)
VSEL,S,MAT,,MAT_YS
CM,V_BSUB,VOLU&&&&&&
VSEL,INVE&&&&&&&&&&
CM,V_SUB,VOLU
VSBV,V_BSUB,V_SUB
CMDEL,V_BSUB&&&&&&&&
CMDEL,V_SUB
NUMCMP,ALL
VATT,MAT_YS&&
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&WPAVE,CEN_ZSD(1,1),CEN_ZSD(2,1),DBKZDZ(I,1)
&CYL4,0,0,R_HJ+L_MG,180,R_HJ,,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&
&CYL4,0,0,R_HJ,180,,,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&&&&&&&&&&&
&BLOCK,-R_HJ,R_HJ,-H_CTOB,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&
&BLOCK,-R_HJ,R_HJ,-H_ZQ,-H_CTOB,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&
&BLOCK,-R_HJ-L_MG,-R_HJ,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&
&BLOCK,R_HJ,R_HJ+L_MG,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&
&WPAVE,CEN_YSD(1,1),CEN_YSD(2,1),DBKZDZ(I,1)
&CYL4,0,0,R_HJ+L_MG,180,R_HJ,,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&
&CYL4,0,0,R_HJ,180,,,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&&&&&&&&&&&
&BLOCK,-R_HJ,R_HJ,-H_CTOB,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&
&BLOCK,-R_HJ,R_HJ,-H_ZQ,-H_CTOB,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&
&BLOCK,-R_HJ-L_MG,-R_HJ,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&
&BLOCK,R_HJ,R_HJ+L_MG,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&
VGLUE,ALL&&
NUMCMP,ALL
CM,ROCK,VOLU
WPAVE,0,0,0
*IF,FLAG_RD,EQ,-1,OR,FLAG_RD,EQ,2,THEN
&VSTART=1000
&NUMSTR,VOLU,VSTART
&BLOCK,Z_RDCC(1,2)-Z_RDCC(1,1)/2,Z_RDCC(1,2)+Z_RDCC(1,1)/2,Z_RDCC(2,2)-Z_RDCC(2,1)/2,Z_RDCC(2,2)+Z_RDCC(2,1)/2,Z_RDCC(3,2)-Z_RDCC(3,1)/2,Z_RDCC(3,2)+Z_RDCC(3,1)/2
&VSBV,ROCK,VSTART
&NUMCMP,VOLU
&CMDEL,ROCK
&CM,ROCK,VOLU
*IF,FLAG_RD,EQ,1,OR,FLAG_RD,EQ,2,THEN
&VSTART=1000
&NUMSTR,VOLU,VSTART
&BLOCK,Y_RDCC(1,2)-Y_RDCC(1,1)/2,Y_RDCC(1,2)+Y_RDCC(1,1)/2,Y_RDCC(2,2)-Y_RDCC(2,1)/2,Y_RDCC(2,2)+Y_RDCC(2,1)/2,Y_RDCC(3,2)-Y_RDCC(3,1)/2,Y_RDCC(3,2)+Y_RDCC(3,1)/2
&VSBV,ROCK,VSTART
&NUMCMP,VOLU
&CMDEL,ROCK
&CM,ROCK,VOLU
NUMCMP,ALL
!----------------------------------------------------------------------------------------------------------------------!
VATT,MAT_YS
VSEL,S,LOC,X,CEN_ZSD(1,1)-1E-7,CEN_ZSD(1,1)+1E-7&&&&&&&&
CM,CM_TEMPV,VOLU
ASEL,R,LOC,Y,CEN_ZSD(2,1)-H_ZQ-1E-7,CEN_ZSD(2,1)-H_ZQ+1E-7&&&&
CM,CM_ZDB,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_TEMPV
ASEL,R,LOC,Y,CEN_ZSD(2,1)-1E-7,CEN_ZSD(2,1)+1E-7&&&&&&&&&&
ASEL,U,LOC,X,CEN_ZSD(1,1)-R_HJ,CEN_ZSD(1,1)-R_HJ-L_MG&&&&
ASEL,U,LOC,X,CEN_ZSD(1,1)+R_HJ,CEN_ZSD(1,1)+R_HJ+L_MG
CM,CM_ZSD,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_ZSD
CMSEL,A,CM_ZDB
VATT,MAT_SD
CM,CM_ZSDKW,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_TEMPV
VSEL,U,MAT,,MAT_SD&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
ASEL,R,LOC,Y,CEN_ZSD(2,1)
VATT,MAT_MGKW
CM,CM_ZMG,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
VSEL,S,LOC,X,CEN_YSD(1,1)-1E-7,CEN_YSD(1,1)+1E-7&&&&&&&&
CM,CM_TEMPV,VOLU
ASEL,R,LOC,Y,CEN_YSD(2,1)-H_ZQ-1E-7,CEN_YSD(2,1)-H_ZQ+1E-7&&&&
CM,CM_YDB,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_TEMPV
ASEL,R,LOC,Y,CEN_YSD(2,1)-1E-7,CEN_YSD(2,1)+1E-7&&&&&&&&&&
ASEL,U,LOC,X,CEN_YSD(1,1)-R_HJ,CEN_YSD(1,1)-R_HJ-L_MG&&&&
ASEL,U,LOC,X,CEN_YSD(1,1)+R_HJ,CEN_YSD(1,1)+R_HJ+L_MG
CM,CM_YSD,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_YSD
CMSEL,A,CM_YDB
VATT,MAT_SD
CM,CM_YSDKW,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_TEMPV
VSEL,U,MAT,,MAT_SD&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
ASEL,R,LOC,Y,CEN_YSD(2,1)
VATT,MAT_MGKW
CM,CM_YMG,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMDEL,CM_TEMPV
ASEL,S,LOC,X,DBKZDX(1,1)
ASEL,A,LOC,X,DBKZDX(9,1)
DA,ALL,UX,0&&&&&&&&
ASEL,S,LOC,Z,DBKZDZ(1,1)
ASEL,A,LOC,Z,DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1)
DA,ALL,UZ,0&&&&&&
ASEL,S,LOC,Y,H_JS
DA,ALL,ALL&&&&&&&&
*DIM,ZLT_HFX,ARRAY,LTS_ZD,1
ZLT_HFX(1,1)=5,3 !4,2
*DIM,YLT_HFX,ARRAY,LTS_YD,1
YLT_HFX(1,1)=5,3 !4,2
*DIM,LT_HFZ,ARRAY,(Y_YJMAX-Y_YJMIN)/PJ_MG+2
LT_HFZ(1,1)=4&&&&&&&&&&&&&&&&&&
*DO,I,2,(Y_YJMAX-Y_YJMIN)/PJ_MG+1&&
<_HFZ(I,1)=2 !PJ_MG
LT_HFZ((Y_YJMAX-Y_YJMIN)/PJ_MG+2,1)=4 !10
*DO,J,1,LTS_ZD
&*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&&X_VCEN=(DBKZDX(J,1)+DBKZDX(J+1,1))/2
&&Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2
&&VSEL,S,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7
&&VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7
&&CM,CM_HFX,LINE
&&LSEL,R,LOC,X,X_VCEN&&&&&&&&&&
&&LESIZE,ALL,,,ZLT_HFX(J,1),,1
&&CMSEL,S,CM_HFX
&&LSEL,R,LOC,Z,Z_VCEN
&&LESIZE,ALL,,,LT_HFZ(I,1),,1
&&CMSEL,S,CM_HFX
&&LSEL,U,LOC,X,X_VCEN
&&LSEL,U,LOC,Z,Z_VCEN
&&LESIZE,ALL,,,9,,1&&
&&CMSEL,S,CM_HFX
&&MSHAPE,0,3D
&&MSHKEY,1
&&VMESH,ALL
&*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&&X_VCEN=(DBKZDX(9-J+1,1)+DBKZDX(9-J,1))/2
&&Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2
&&VSEL,S,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7
&&VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7
&&CM,CM_HFX,LINE
&&LSEL,R,LOC,X,X_VCEN&&&&&&&&&&
&&LESIZE,ALL,,,YLT_HFX(J,1),,1
&&CMSEL,S,CM_HFX
&&LSEL,R,LOC,Z,Z_VCEN
&&LESIZE,ALL,,,LT_HFZ(I,1),,1
&&CMSEL,S,CM_HFX
&&LSEL,U,LOC,X,X_VCEN
&&LSEL,U,LOC,Z,Z_VCEN
&&LESIZE,ALL,,,10,,1&&
&&CMSEL,S,CM_HFX
&&MSHAPE,0,3D
&&MSHKEY,1
&&VMESH,ALL
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&CMSEL,S,CM_ZSDKW
&X_VCEN=CEN_ZSD(1,1)
&Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2
&VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7
&CM,CM_VMESH,VOLU
&CM,CM_HFX,LINE
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7&&&&&&&&
&LESIZE,ALL,,,2,,1&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7&&&&&&&&&&
&LESIZE,ALL,,,2,,1&&&&&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7
&LSEL,U,LOC,Z,Z_VCEN
&LESIZE,ALL,2&&&&&&&&&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,Z,Z_VCEN&&&&
&LESIZE,ALL,,,LT_HFZ(I,1),,1
&CMSEL,S,CM_HFX
&LSEL,U,LOC,X,X_VCEN
&LSEL,U,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7
&LSEL,U,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7
&LSEL,U,LOC,Z,Z_VCEN
&LESIZE,ALL,,,1,,1&&&&
&CMSEL,S,CM_VMESH
&MSHAPE,0,3D
&VSWEEP,ALL
CMSEL,S,CM_ZMG
VSWEEP,ALL
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&CMSEL,S,CM_YSDKW
&X_VCEN=CEN_YSD(1,1)
&Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2
&VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7
&CM,CM_VMESH,VOLU
&CM,CM_HFX,LINE
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7&&&&&&&&
&LESIZE,ALL,,,2,,1&&&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7&&&&&&&&
&LESIZE,ALL,,,2,,1&&&&&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7
&LSEL,U,LOC,Z,Z_VCEN
&LSEL,U,LOC,Y,CEN_YSD(2,1)-H_CTOB&&&&&&&&&&&&&&&&&&&&
&LESIZE,ALL,2&&&&&&&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,Z,Z_VCEN&&
&LESIZE,ALL,,,LT_HFZ(I,1),,1
&CMSEL,S,CM_HFX
&LSEL,U,LOC,X,X_VCEN
&LSEL,U,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7
&LSEL,U,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7
&LSEL,U,LOC,Z,Z_VCEN
&LESIZE,ALL,,,1,,1&&&&
&CMSEL,S,CM_VMESH
&MSHAPE,0,3D
&VSWEEP,ALL
&CMSEL,S,CM_YMG
&ESIZE,2&&
&VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7
&VSWEEP,ALL
VSEL,U,MAT,,MAT_SD
VSEL,U,MAT,,MAT_MGKW
VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(9-LTS_YD,1)
CM,CM_VMESH,VOLU
VSEL,S,,,CM_VMESH
VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(LTS_ZD+2,1)&&
LSEL,R,LOC,X,(DBKZDX(LTS_ZD+1,1)+DBKZDX(LTS_ZD+2,1))/2
LESIZE,ALL,,,5&&
VSEL,S,,,CM_VMESH
VSEL,R,LOC,X,DBKZDX(LTS_ZD+2,1),DBKZDX(LTS_ZD+3,1)
LSEL,R,LOC,X,(DBKZDX(LTS_ZD+2,1)+DBKZDX(LTS_ZD+3,1))/2
LESIZE,ALL,,,3&&
VSEL,S,,,CM_VMESH
VSEL,R,LOC,X,DBKZDX(LTS_ZD+3,1),DBKZDX(LTS_ZD+4,1)&&
LSEL,R,LOC,X,(DBKZDX(LTS_ZD+3,1)+DBKZDX(LTS_ZD+4,1))/2
LESIZE,ALL,,,3&&
VSEL,S,,,CM_VMESH
VSEL,R,LOC,X,DBKZDX(LTS_ZD+4,1),DBKZDX(LTS_ZD+5,1)&&
LSEL,R,LOC,X,(DBKZDX(LTS_ZD+4,1)+DBKZDX(LTS_ZD+5,1))/2
LESIZE,ALL,,,5&&
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&!CMSEL,S,CM_VMESH
&VSEL,S,,,CM_VMESH
&VSEL,R,LOC,Z,DBKZDZ(I,1),DBKZDZ(I+1,1)
&VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(LTS_ZD+3,1)
&CM,CM_NEXTV,VOLU
&CMSEL,S,CM_NEXTV
&VSWEEP,ALL
LSEL,S,LOC,X,Y_RDCC(1,2)-Y_RDCC(1,1)/2-1E-7,Y_RDCC(1,2)+Y_RDCC(1,1)/2+1E-7
LSEL,R,LOC,Y,Y_RDCC(2,2)-Y_RDCC(2,1)/2-1E-7,Y_RDCC(2,2)+Y_RDCC(2,1)/2+1E-7
LSEL,R,LOC,Z,Y_RDCC(3,2)-Y_RDCC(3,1)/2-1E-7,Y_RDCC(3,2)+Y_RDCC(3,1)/2+1E-7
LESIZE,ALL,1&&
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&VSEL,S,,,CM_VMESH
&VSEL,R,LOC,Z,DBKZDZ(I,1),DBKZDZ(I+1,1)
&CM,CM_NEXTV,VOLU
&CMSEL,S,CM_NEXTV
&VSEL,R,LOC,X,DBKZDX(LTS_ZD+3,1),DBKZDX(LTS_ZD+3+2,1)
&VSWEEP,ALL
VSEL,S,,,CM_ZSDKW&&
VSEL,S,,,CM_ZMG
CM,CM_ZGXW,AREA
MAT,MAT_GXW
VSEL,S,,,CM_YSDKW&&
VSEL,S,,,CM_YMG
CM,CM_YGXW,AREA
MAT,MAT_GXW
NUMCMP,ALL
MPCHG,MAT_TXYS,ALL
SAVE,TNLMESH,DB&&
!------------------------------------------------------- 开始求解
---------------------------------------!
ANTYPE,STATIC,NEW
TIME,1&&&&
NROPT,FULL
AUTOTS,ON&&
DELTIM,0.5,0.01,0.5
LNSRCH,ON&&
PRED,ON&&&&
ESEL,S,TYPE,,2
/FILNAME,CALTUN
RESUME,TNLMESH,DB
ANTYPE,STATIC,NEW
TIME,1&&&&&&&&
NROPT,FULL
AUTOTS,ON&&
DELTIM,0.5,0.01,0.5
LNSRCH,ON&&
PRED,ON&&&&
ESEL,U,TYPE,,2
ISFILE,READ,TUNNEL,IST,,0
ESEL,S,TYPE,,2
EKILL,ALL!*
CNVTOL,F,0.05
!----------------------------------------------------------------------
开挖模拟 ---------------------------------------------------!
!-----------------------------------&&&&第一步开挖，这一步没有支护&&&&
----------------------------------!
VSEL,S,,,CM_ZSDKW
VSEL,A,,,CM_YSDKW
VSEL,R,LOC,Z,DBKZDZ(1,1),DBKZDZ(2,1)
!ESEL,A,TYPE,,2
EKILL,ALL&&
TIME,2&&&&
NROPT,FULL
AUTOTS,ON&&
DELTIM,0.3,0.001,0.5
SAVE,KW-1,DB&&
!---------------------------- 循环做第二步到第五步的开挖，第一步到第四步的支护
----------------------------!
ESEL,S,TYPE,,2
EKILL,ALL&&&&&&
!----------------------------&&第I步开挖,第I-1步支护
---------------------!
*DIM,KWD_TOU,ARRAY,4
*DIM,KWD_WEI,ARRAY,4
*DIM,ZHD_TOU,ARRAY,4
*DIM,ZHD_WEI,ARRAY,4
KWD_TOU(1,1)=3,4,5,6
KWD_WEI(1,1)=2,3,4,5
ZHD_TOU(1,1)=2,3,4,5
ZHD_WEI(1,1)=1,2,3,4
&TOU=KWD_TOU(I,1)
&WEI=KWD_WEI(I,1)
&VSEL,S,,,CM_ZSDKW
&VSEL,A,,,CM_YSDKW
&VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
&ESLV,R&&&&
&EKILL,ALL&&
&TOU=ZHD_TOU(I,1)
&WEI=ZHD_WEI(I,1)
&VSEL,S,,,CM_ZDB&&&&&&&&&&&&&&&&&&&&&&
&VSEL,A,,,CM_YDB
&VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
&EALIVE,ALL&&&&&&
&MPCHG,MAT_DB,ALL
&VSEL,S,,,CM_ZMG&&&&&&&&&&&&&&&&&&&&&&
&VSEL,A,,,CM_YMG
&VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
&MPCHG,MAT_MG,ALL
&ASEL,S,,,CM_ZGXW&&&&&&&&&&&&&&&&&&&&
&ASEL,A,,,CM_YGXW
&ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
&EALIVE,ALL&&&&&&&&
&MPCHG,MAT_GXW,ALL
&ASEL,S,,,CM_ZGXW&&
&ASEL,A,,,CM_YGXW
&ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
&SFE,ALL,,PRES,,ZH_MG
&TIME,I+2&&&&
&NROPT,FULL
&AUTOTS,ON&&
&DELTIM,0.3,0.001,0.5
&LNSRCH,ON&&
&PRED,ON&&
&!OUTRES,ALL,ALL
&SAVE,KW-%I+1%,DB&&
!-----------------------------------&&&&第5步支护，这一步没有开挖&&&&
----------------------------------!
VSEL,S,,,CM_ZDB&&&&&&&&&&&&&&&&&&&&&&
VSEL,A,,,CM_YDB
VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
EALIVE,ALL&&&&&&
MPCHG,MAT_DB,ALL
VSEL,S,,,CM_ZMG&&&&&&&&&&&&&&&&&&&&
VSEL,A,,,CM_YMG
VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
MPCHG,MAT_MG,ALL
ASEL,S,,,CM_ZGXW&&&&&&&&&&&&&&&&&&&&
ASEL,A,,,CM_YGXW
ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
EALIVE,ALL&&&&&&
MPCHG,MAT_GXW,ALL
ASEL,S,,,CM_ZGXW&&
ASEL,A,,,CM_YGXW
ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
SFA,ALL,,PRES,ZH_MG
TIME,7&&&&
NROPT,FULL
AUTOTS,ON&&
DELTIM,0.3,0.01,0.5
LNSRCH,ON&&
!OUTRES,ALL,ALL
SAVE,KW-6,DB&&
/EXIT,NOSAVE
岩土接触分析实例
例子1是为了获得初始应力，计算可以完成。命令流为：
!以便获得初始应力，收敛！！考虑土与岩石之间的摩擦，点面接触
fini !考虑土体的相互作用。自由划分网格，局部细化,不收敛！！！！！！！！！！！
/filname,diji11
/title,地基应力、应变的分布规律
!ET,1,PLANE2
ET,1,PLANE42
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
MP,EX,1,8e9 !岩石性质
MP,PRXY,1,0.29
MP,DENS,1,2600
MP,EX,2,5.0e7 !风化槽土体性质?
MP,PRXY,2,0.25
TBDATA,1,15e3,20,0
MP,DENS,2,2000
et,2,contac48,,,1
mp,mu,3,2.5
r,1,5.0e7,5.0e6,0.01,1.0,0.001,0.001
r,2,5.0e7,5.0e6,0.01,1.0,0.001,0.001
K,1,0,0,0,
K,2,2.5,-12,0,
K,3,6.5,-12,0,
K,4,13.5,-26.5,0,
K,5,20,-30,0,
K,6,24.5,-27,0,
K,7,31.5,-12,0
K,8,35.5,-12,0,
K,9,38,0,0,
k,10,-40,0,0,
k,11,-40,-60,0,
k,12,80,-60,0,
k,13,80,0,0,
K,14,6.5,-12,0,
K,15,13.5,-26.5,0,
K,16,20,-30,0,
K,17,24.5,-27,0,
K,18,31.5,-12,0
lstr,i,i+1
*do,i,10,12
lstr,i,i+1
*do,i,14,17
lstr,i,i+1
LARC,14,18,16,14.57
FLST,2,13,4
*do,i,1,13
FLST,2,5,4
*do,i,14,18
asel,s,area,,1
asel,s,area,,2
lsel,s,line,,3,4,1
CM,_TARGET1,node
lsel,s,line,,14,15,1
CM,_CONTACT1,node
lsel,s,line,,5,6,1
CM,_TARGET2,node
lsel,s,line,,16,17,1
CM,_CONTACT2,node
gcgen,_CONTACT1,_TARGET1
gcgen,_CONTACT2,_TARGET2
antype,static
nropt,full,,off
nsubst,100,1000,10
LSEL,s,line,,10,12,2
CM,_Y,LINE
DL,_Y, ,ux,0,0
LSEL,s,line,,11
CM,_Y,LINE
DL,_Y, ,uy,0,0
iswrite,on
outres,all,all
!读入初始应力，收敛！！考虑土和混凝土与岩石之间的摩擦，点面接触
fini !考虑土体和混凝土之间的相互作用。不收敛！！！！！！！！！！！
/filname,diji
/title,地基应力、应变的分布规律
ET,1,PLANE42
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
MP,EX,1,8e9 !岩石性质
MP,PRXY,1,0.29
MP,DENS,1,2600
MP,EX,2,5.0e7 !风化槽土体性质?
MP,PRXY,2,0.25
TBDATA,1,15e3,20,0
MP,DENS,2,2000
MP,EX,3,2.55e10 !混凝土性质
MP,PRXY,3,0.167
MP,DENS,3,2400
et,2,contac48,,,1
mp,mu,4,2.5
mp,mu,5,4.5
mp,mu,6,3.0
r,1,5.0e7,5.0e6,0.01,1.0,0.001,0.001
r,2,5.0e7,5.0e6,0.01,1.0,0.001,0.001
r,3,8.0e9,8.0e8,0.01,1.0,0.001,0.001
r,4,8.0e9,8.0e8,0.01,1.0,0.001,0.001
r,5,5.0e7,5.0e6,0.01,1.0,0.001,0.001
K,1,0,0,0,
K,2,2.5,-12,0,
K,3,6.5,-12,0,
K,4,13.5,-26.5,0,
K,5,20,-30,0,
K,6,24.5,-27,0,
K,7,31.5,-12,0
K,8,35.5,-12,0,
K,9,38,0,0,
k,10,-40,0,0,
k,11,-40,-60,0,
k,12,80,-60,0,
k,13,80,0,0,
K,14,6.5,-12,0,
K,15,13.5,-26.5,0,
K,16,20,-30,0,
K,17,24.5,-27,0,
K,18,31.5,-12,0
K,19,0,0,0,
K,20,2.5,-12,0,
K,21,6.5,-12,0,
K,22,31.5,-12,0
K,23,35.5,-12,0,
K,24,38,0,0,
K,25,20,0,0,
lstr,i,i+1
*do,i,10,12
lstr,i,i+1
*do,i,14,17
lstr,i,i+1
LARC,14,18,16,14.57
*do,i,19,20
lstr,i,i+1
LARC,21,22,16,14.57
*do,i,22,24
lstr,i,i+1
LSTR,25,19
FLST,2,13,4
*do,i,1,13
FLST,2,5,4
*do,i,14,18
FLST,2,7,4
*do,i,19,25
asel,s,area,,1
asel,s,area,,2
asel,s,area,,1,2
lsel,s,line,,3,4,1
CM,_TARGET1,node
lsel,s,line,,14,15,1
CM,_CONTACT1,node
lsel,s,line,,5,6,1
CM,_TARGET2,node
lsel,s,line,,16,17,1
CM,_CONTACT2,node
gcgen,_CONTACT1,_TARGET1
gcgen,_CONTACT2,_TARGET2
asel,s,area,,3
lsel,s,line,,1,8,7
CM,_TARGET3,node
lsel,s,line,,19,23,4
CM,_CONTACT3,node
lsel,s,line,,2,7,5
CM,_TARGET4,node
lsel,s,line,,20,22,2
CM,_CONTACT4,node
gcgen,_CONTACT3,_TARGET3
gcgen,_CONTACT4,_TARGET4
lsel,s,line,,18
CM,_TARGET5,node
lsel,s,line,,21
CM,_CONTACT5,node
gcgen,_CONTACT5,_TARGET5
antype,static
nropt,full,,off
esel,s,,,1,451
isfile,read,diji11,ist,,1
nsubst,100,1000,10
LSEL,s,line,,10,12,2
CM,_Y,LINE
DL,_Y, ,ux,0,0
LSEL,s,line,,11
CM,_Y,LINE
DL,_Y, ,uy,0,0
LSEL,s,line,,25
CM,_Y1,LINE
CMSEL,,_Y1
SFL,_Y1,PRES,4.5e5,6e5,
LSEL,s,line,,24
CM,_Y1,LINE
CMSEL,,_Y1
SFL,_Y1,PRES,4.5e5,4.5e5
iswrite,on
outres,all,all
钢筋混凝土管的动力响应特性分析实例
在管内壁激振后，采用瞬态模块，模拟半埋设钢筋混凝土管在内部震源作用下的动力特性。
请多提宝贵意见，如果给分就太好了！
/COM,ANSYS MEDIA REL. 7.1 (03-13-2003) REF. VERIF. MANUAL: REL.
/CONFIG,NRES,2000
/VERIFY,ACOUSTIC_WELL_RADIUS
/TITLE,ACOUSTIC, ACOUSTIC_WELL_RADIUS
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!! 按平面问题分析
!! 钢板＋钢筋混凝土＋砂浆/泥浆＋原位岩/土（瞬态）
!! 6种介质成分，存放在介质属性数组中。
!! 计算边界外延一定宽度，以抑制边界干扰*****************
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
JH_HSGD=10&&!计算高度
JH_GBHD=0.04&&!钢板厚度
JH_HNTHD=0.7&& !混凝土厚度
JH_JYHD=0.6&&!浆液厚度，高度方向上砂浆和泥浆各占一半，砂浆在下，泥浆在上
JH_YTTHD=0.6&&!岩土体厚度
JH_RADWELL=4
ZHD=JH_RADWELL+JH_GBHD+JH_HNTHD+JH_JYHD+JH_YTTHD
BJCR1=ZHD+0.5
BJCR2=ZHD+1.0
BJCR3=ZHD+1.5
BJCR4=ZHD+2.0
SCBCJG=250
ENDTIME=0.001 !总计算时间
DST=1E-06 !记录的时间步长
ZYPL=5E4 !震源频率
PI=2*ASIN(1)
TCCS=6 !介质分类数目
!钢&& 钢筋混凝土 泥浆 砂浆 原位岩/土
*DIM,TC,ARRAY,TCCS,5
!该数组为介质参数，第一列为介质密度&&&&
第二列为泊松比，
!&& 第三列为弹模
! 第四列为纵波波速，用于计算时间步长
! 第五列为横波波速，如为液体，横波波速为0
!土层自下而上排列，物理参数暂时按各土层相同考虑
TC(1,1)=00,00
TC(1,2)=0.285,0.32,0.48,0.35,0.35,0.35
!液体泊松比为0.5,计算出现错误，按提示改为0.499
TC(1,3)=210E9,36E9,2E6,6E9,20E9,20E9
!计算各种介质的P、S波速度,以计算时间步长
*DO,I,1,TCCS
&&LAME1=TC(I,3)*TC(I,2)/(1+TC(I,2))/(1-2*TC(I,2))
&&LAME2=TC(I,3)/2/(1+TC(I,2))
&&&&&&TC(I,4)=((LAME1+2*LAME2)/TC(I,1))**.5
&&&&&&TC(I,5)=(LAME2/TC(I,1))**.5
!《土中应力波及其量测》P26
!#########生成各介质区的边界控制点############
JZFQ=4&&&&&&
!介质分区数目
R1=JH_RADWELL
R2=JH_RADWELL+JH_GBHD
R3=JH_RADWELL+JH_GBHD+JH_HNTHD
R4=JH_RADWELL+JH_GBHD+JH_HNTHD+JH_JYHD
R5=JH_RADWELL+JH_GBHD+JH_HNTHD+JH_JYHD+JH_YTTHD
*DIM,JHRB,ARRAY,JZFQ+BJCSM+1,1
JHRB(1,1)=R1,R2,R3,R4,R5,R6,R7,R8,R9!
*DIM,JHRA,ARRAY,JZFQ+BJCSM+1,4&&&&
!定义数组，用于存放
*DIM,AMAT,ARRAY,JZFQ+BJCSM+1,1&&
*DIM,DYCC,ARRAY,JZFQ+BJCSM+1,1&&&&
JHRA(1,1)=R1,R2,R3,R3,R4,R5,R6,R7,R8!
JHRA(1,2)=R2,R3,R4,R4,R5,R6,R7,R8,R9!
JHRA(1,3)=0,0,0,180,0,0,0,0,0
JHRA(1,4)=360,360,180,360,360,360,360,360,360
AMAT(1,1)=1,2,3,4,5,5,5,5,5&&&&!介质编号&&&&&&
DYCC(1,1)=0.2,0.2,0.2,0.2,0.2,0.3,0.45,0.6,1!剖分尺寸
/REPLOT,RESIZE
*DO,I,1,JZFQ+BJCSM+1
&&CYL4,0,0,JHRA(I,1),JHRA(I,3),JHRA(I,2),JHRA(I,4)
NUMMRG,ALL
NUMCMP,ALL !压缩各定义项的编号,使之连续
!############计算时间步长##################
VMAX=TC(1,4) !钢板的纵波速度控制时间步长
DT=DYCC(1)/VMAX/2
*IF,DT,GE,DST,THEN !若时间步长大于设定的记录时间步长，采用记录时间步长记录
ZBC=ENDTIME/DT !时间总步数
!############震源处理##########
T0=1/ZYPL/2
!震源信号持续时间，采用正弦信号激震，也可采用RICKER信号，但要从数组文件中读取，可考虑设定一FLAG控制
XHCXSJ=NINT(1/ZYPL/DT)+1 !信号持续的时间步长
T0I=NINT(XHCXSJ/2)
F=ZYPL&&&&!经MATLAB实际输出后获得的信号实际频率
!######信号记录步长处理#########
*IF,DST/DT,LT,NINT(DST/DT),THEN
!计算信号记录的时间步长间隔，实际记录的时间步长与设定值可能不同
&&RECORDSTEP=NINT(DST/DT)-1
&&RECORDSTEP=NINT(DST/DT)
!#################时间步计算结束#########################
!##################单元剖分################
ET,1,PLANE183
*DO,I,1,JZFQ+BJCSM+1
&&LOCAL, 11, 1, 0,
0&&!定义局部坐标系，柱体
&&CSYS,11&&&&
&&LSEL,S,LOC,X,JHRB(I),&&&&
&&LESIZE,ALL,,R1ANGLE
&&ALLSEL,ALL
AREFINE,1,9,,1!精细剖分
AREFINE,1,4,,1!精细剖分
!#################赋予材料属性#################
*DO,I,1,TCCS
&&MP,DENS,I,TC(I,1)
&&MP,NUXY,I,TC(I,2)
&&MP,EX,I,TC(I,3)
LOCAL, 11, 1, 0,
0&&!定义局部坐标系，柱体
*DO,I,1,JZFQ+BJCSM+1
&&NSEL,S,LOC,X,JHRA(I,1),JHRA(I,2)
&&NSEL,R,LOC,Y,JHRA(I,3),JHRA(I,4)
&&ESLN,S&&&&&&&&
!边界节点单元
&&MAT,AMAT(I)
&&EMODIFY,ALL&&&&&&&&!确认材料参数修改
&&ALLSEL,ALL&&&&&&&&
!重新选择ALL，准备下一步操作
/PNUM,MAT,1
LOCAL, 11, 1, 0,
0&&!定义局部坐标系，柱体
CSYS,11&&&&
LSEL,S,LOC,X,R1
*GET, R1NUM, NODE, 0, COUNT
*GET, R1NMIN, NODE, 0,NUM,
MIN&&!获取最小节点编号
*DIM,R1NODE,ARRAY,R1NUM,1&&!定义内圈节点数组
*DIM,R1NODEXY,ARRAY,R1NUM,2&&!定义内圈节点坐标数组，存储极坐标
R1NODE(1,1)=R1NMIN&&
!赋内圈节点编号的最小值
*GET,RNX, NODE, R1NODE(1,1), LOC,
X&&!获取内圈的节点的坐标，
*GET,RNY, NODE, R1NODE(1,1), LOC,
Y&&!在此为极坐标
R1NODEXY(1,1)=RNX
R1NODEXY(1,2)=RNY
*DO,I,2,R1NUM
&&R1NODE(I,1)=NDNEXT(R1NODE(I-1))&&
!对内圈节点赋值，NDNEXT的输出为选定的节点中编号大于且最接近该节点的节点编号
&&*GET,RNX, NODE, R1NODE(I,1),
LOC, X&&!获取内圈的节点的坐标，
&&*GET,RNY, NODE, R1NODE(I,1),
LOC, Y&&!在此为极坐标
&&R1NODEXY(I,1)=RNX
&&R1NODEXY(I,2)=RNY&&
!###############按极坐标中角度，由小到大排列节点编号和对应坐标
*DO,I,1,R1NUM-1
*DO,J,I+1,R1NUM
&&*IF,R1NODEXY(I,2),GT,R1NODEXY(J,2),THEN
&&&&RY=R1NODEXY(I,2)
&&&&R1NODEXY(I,2)=R1NODEXY(J,2)
&&&&R1NODEXY(J,2)=RY
&&&&RX=R1NODEXY(I,1)
&&&&R1NODEXY(I,1)=R1NODEXY(J,1)
&&&&R1NODEXY(J,1)=RX
&&&&RNODE=R1NODE(I)
&&&&R1NODE(I)=R1NODE(J)
&&&&R1NODE(J)=RNODE
&&*ENDIF&&&&&&&&&&
!####################排序OVER#####################
*CFOPEN,R1NODEXY,DAT!输出井壁内表面节点编号和坐标
*DO,K,1,R1NUM
&&NUM_NODE=R1NODE(K)
&&RX=R1NODEXY(K,1)
&&RY=R1NODEXY(K,2)
&&*VWRITE,NUM_NODE,RX,RY
(F10.0,F13.2,F13.2)
ALLSEL,ALL
!#################求解过程##########################
!##与竖直方向采用直角坐标的不同：
!##1、震源位置变换方式
!##2、记录点位置
!##3、边界条件设定，不作专门的吸收边界，而以空间扩大来处理
*CFOPEN,USUM,DAT
!采用REIKER子波，在初始时段施加
NSUBST,1, , ,1 !1个子步
ANTYPE,TRANS !时程分析
OUTPR,BASIC,SCBCJG, !每隔一定时间步长输出基本值
NSEL,S,LOC,X,R9
D,ALL,UX,0
D,ALL,UY,0
ALLSEL,ALL
*DO,I,1,ZBC
&&*IF,I,GT,1,AND,I,LE,XHCXSJ,THEN&&!在信号持续时间段延续脉冲子波
&&&&*DO,J,1,R1NUM
&&&&&&*GET,RNX,
NODE,R1NODE(J,1), LOC, X
&&&&&&*GET,RNY,
NODE,R1NODE(J,1), LOC, Y
&&&&&&NN=R1NODE(J,1)
&&&&&&PIFT=PI*ZYPL*(I-T0I)*DT
&&&&&&PIFT=PIFT**2
&&&&&&AC=-(1-2*PIFT)*EXP(-PIFT)&&!RECKER子波，何兵寿，P48
&&&&&&F,NN,FX,AC*RNX/R1&&
!施加X方向作用力
&&&&&&F,NN,FY,AC*RNY/R1&&!施加Y方向作用力
&&&&*ENDDO
&&*IF,I,EQ,XHCXSJ+1,THEN&&!在XHCXSJ+1时刻，删除作用在节点上的位移&&
&&&&FDELE,ALL&&!在信号持续时间结束后，删除所有作用力
&&TIME,I*DT
*DO,K,1,R1NUM
&&NUM_NODE=R1NODE(K)
&&RX=R1NODEXY(K,1)
&&RY=R1NODEXY(K,2)
&&*GET, USUMR1, NODE,
NUM_NODE,U,SUM !输出节点编号，极坐标R，角度，USUM
*VWRITE,I,K,NUM_NODE,USUMR1
(F10.0,F10.0,F10.0,E13.3)
BNODEN=NODE(0,R1,0)
BNODES=NODE(0,-R1,0)
FILE,'file','rst','.'
NSOL,2,BNODES,U,X,uxs
NSOL,4,BNODES,U,Y,uys
NSOL,3,BNODEN,U,X,uxn
NSOL,5,BNODEN,U,Y,uyn
PLVAR,uxs,uxn
/IMAGE,SAVE,UXS_N,JPG
PLVAR,uys,uyn, , , , , , ,
/IMAGE,SAVE,UYS_N,JPG
LOCAL, 11, 1, 0,
0&&!定义局部坐标系，柱体
CSYS,11&&&&
NSEL,S,LOC,X,R1,R5
PLNS,U,SUM
ANDATA,0.5, ,2,1,zbc,10,0,1
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。}
四层弹簧－质点模型的地震分析
包括频率分析，谱分析和时程分析
作者：陆新征，清华大学土木工程系
DAMPRATIO=0.02 !振型阻尼系数都为0.02
ET,1,BEAM188
KEYOPT,1,7,1
KEYOPT,1,8,1
ET,2,MASS21,,,4
MP,EX,1,210E9
MP,NUXY,1,0.3
MP,DENS,1,7800
MPDATA,DAMP,1,,DAMPRATIO
SECTYPE,3,BEAM,I
SECDATA,0.2,0.2,0.5,0.01,0.01,0.008
R,1,1.6E2,1.6E2
R,2,1.2E2,1.2E2
K,5,0,12,0
K,6,0,0,10000
LSEL,,,,1,4
LATT,1,0,1,,6,,3
LESIZE,ALL, , ,3 , , , , ,1
NUMMRG,ALL
NUMCMP,ALL
*SET,NT,1001
*SET,DT,0.02
*DIM,AC,,NT
*VREAD,AC(1),RECORD,TXT
/SOLU !模态分析
MODOPT,SUBSP,8
MXPAND,8, , ,1
!得到自振频率1
*GET,FREQ1,MODE,1,FREQ
/ESHAPE,1.0
ANTYPE,SPECTR !谱分析
SPOPT,SPRS
MDAMP,1,DAMPRATIO, , , , , ,
SVTYPE,2,1
FREQ,0.167,0.25,0.333,0.,0.667,1,1.25,1.667
SV,,0.........
FREQ,2.222,10,100000
SV,,1...3528
!按7度多遇地震，地震影响系数0.08 第一组III场地TG＝0.45
SRSS,0.01,DISP
/INP,,MCOM
NSORT,U,SUM
*GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形
!开始地震计算
/CONFIG,NRES,20000
ANTYPE,TRANS
TRNOPT,FULL
ALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926
BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)
*DO,I,1,1001
ACEL,AC(I),0,0
TIME,I*0.02
OUTRES,ALL,ALL
NSOL,2,20,U,X,
PLVAR,2, , , , , , , , , ,
STORE,MERGE
*GET,DMAX_TRANS,VARI,2,EXTREM,VMAX, , ! 时程分析最大位移
record.txt
悬臂梁地震分析
!定义单元类型
ET,1,BEAM3
R,1,0.09,0..3, !0.3*0.3m截面
MP,EX,1,2.06e11
MP,NUXY,1,0.3
MP,DAMP,1,0.02 !阻尼比0.02
MP,dens,1,7.85e3
!--------------------------------------------------
!定义节点单元和位移约束
*do,i,1,11
n,i,(i-1)*0.5,0,0 !5m长,分10段
*do,i,1,10
!--------------------------------------------------
!定义和读入时程曲线
NT=1000 !时程曲线有NT个点
DT=0.01 !时间间隔
*dim,ac,,NT
!读入数据,这个数据文件可以用excel等软件来写
/input,tianjin,txt !天津波东西向共10秒
!--------------------------------------------------
NSUBST,1, , ,1 !1个子步
OUTRES,ALL,1 !输出每个子步的结果
ANTYPE,TRANS !时程分析
*do,i,1,NT
ACEL,0,ac ,0
!--------------------------------------------------
NSOL,2,11,U,Y,
PLVAR,2,3, , ,
下面是文件tianjin.txt，这个文件可以用excel等软件制作或编写很短的程序就能写出来先
贴前500个点
ac( 1)= -0.
ac( 2)= -0.
ac( 3)= -0.
ac( 4)= -0.
ac( 5)= 0.
ac( 6)= 0.
ac( 7)= 0.
ac( 8)= 0.
ac( 9)= 0.
ac( 10)= 0.
ac( 11)= 0.
ac( 12)= 0.
ac( 13)= 0.
ac( 14)= 0.
ac( 15)= 0.
ac( 16)= 0.
ac( 17)= -0.
ac( 18)= -0.
ac( 19)= 0.
ac( 20)= 0.
ac( 21)= 0.
ac( 22)= 0.
ac( 23)= 0.
ac( 24)= 0.
ac( 25)= -0.
ac( 26)= -0.
ac( 27)= -0.
ac( 28)= -0.
ac( 29)= -0.
ac( 30)= -0.
ac( 31)= -0.
ac( 32)= -0.
ac( 33)= -0.
ac( 34)= -0.
ac( 35)= -0.
ac( 36)= -0.
ac( 37)= -0.
ac( 38)= -0.
ac( 39)= 0.
ac( 40)= 0.
ac( 41)= 0.
ac( 42)= 0.
ac( 43)= -0.
ac( 44)= -0.
ac( 45)= -0.
ac( 46)= -0.
ac( 47)= -0.
ac( 48)= -0.
ac( 49)= -0.
ac( 50)= -0.
ac( 51)= -0.
ac( 52)= 0.
ac( 53)= 0.
ac( 54)= 0.
ac( 55)= 0.
ac( 56)= 0.
ac( 57)= 0.
ac( 58)= -0.
ac( 59)= -0.
ac( 60)= -0.
ac( 61)= -0.
ac( 62)= 0.
ac( 63)= 0.
ac( 64)= 0.
ac( 65)= 0.
ac( 66)= 0.
ac( 67)= 0.
ac( 68)= 0.
ac( 69)= 0.
ac( 70)= 0.
ac( 71)= 0.
ac( 72)= 0.
ac( 73)= 0.
ac( 74)= 0.
ac( 75)= 0.
ac( 76)= -0.
ac( 77)= -0.
ac( 78)= -0.
ac( 79)= -0.
ac( 80)= -0.
ac( 81)= -0.
ac( 82)= -0.
ac( 83)= -0.
ac( 84)= -0.
ac( 85)= -0.
ac( 86)= -0.
ac( 87)= -0.
ac( 88)= -0.
ac( 89)= -0.
ac( 90)= -0.
ac( 91)= -0.
ac( 92)= -0.
ac( 93)= -0.
ac( 94)= -0.
ac( 95)= -0.
ac( 96)= -0.
ac( 97)= -0.
ac( 98)= -0.
ac( 99)= -0.
ac( 100)= -0.
ac( 101)= -0.
ac( 102)= -0.
ac( 103)= -0.
ac( 104)= -0.
ac( 105)= -0.
ac( 106)= -0.
ac( 107)= -0.
ac( 108)= -0.
ac( 109)= -0.
ac( 110)= -0.
ac( 111)= 0.
ac( 112)= 0.
ac( 113)= 0.
ac( 114)= 0.
ac( 115)= 0.
ac( 116)= 0.
ac( 117)= 0.
ac( 118)= 0.
ac( 119)= 0.
ac( 120)= 0.
ac( 121)= 0.
ac( 122)= 0.
ac( 123)= 0.
ac( 124)= -0.
ac( 125)= -0.
ac( 126)= -0.
ac( 127)= -0.
ac( 128)= -0.
ac( 129)= -0.
ac( 130)= -0.
ac( 131)= -0.
ac( 132)= 0.
ac( 133)= 0.
ac( 134)= 0.
ac( 135)= 0.
ac( 136)= 0.
ac( 137)= 0.
ac( 138)= 0.
ac( 139)= 0.
ac( 140)= 0.
ac( 141)= 0.
ac( 142)= 0.
ac( 143)= 0.
ac( 144)= 0.
ac( 145)= 0.
ac( 146)= 0.
ac( 147)= -0.
ac( 148)= -0.
ac( 149)= -0.
ac( 150)= -0.
ac( 151)= 0.
ac( 152)= 0.
ac( 153)= 0.
ac( 154)= 0.
ac( 155)= 0.
ac( 156)= 0.
ac( 157)= 0.
ac( 158)= 0.
ac( 159)= 0.
ac( 160)= 0.
ac( 161)= 0.
ac( 162)= -0.
ac( 163)= -0.
ac( 164)= -0.
ac( 165)= -0.
ac( 166)= -0.
ac( 167)= -0.
ac( 168)= -0.
ac( 169)= -0.
ac( 170)= -0.
ac( 171)= -0.
ac( 172)= -0.
ac( 173)= -0.
ac( 174)= -0.
ac( 175)= -0.
ac( 176)= -0.
ac( 177)= -0.
ac( 178)= -0.
ac( 179)= -0.
ac( 180)= -0.
ac( 181)= -0.
ac( 182)= 0.
ac( 183)= 0.
ac( 184)= 0.
ac( 185)= 0.
ac( 186)= 0.
ac( 187)= 0.
ac( 188)= 0.
ac( 189)= 0.
ac( 190)= -0.
ac( 191)= -0.
ac( 192)= -0.
ac( 193)= -0.
ac( 194)= -0.
ac( 195)= 0.
ac( 196)= 0.
ac( 197)= 0.
ac( 198)= 0.
ac( 199)= 0.
ac( 200)= 0.
ac( 201)= -0.
ac( 202)= -0.
ac( 203)= 0.
ac( 204)= 0.
ac( 205)= 0.
ac( 206)= 0.
ac( 207)= 0.
ac( 208)= -0.
ac( 209)= -0.
ac( 210)= -0.
ac( 211)= -0.
ac( 212)= -0.
ac( 213)= 0.
ac( 214)= 0.
ac( 215)= 0.
ac( 216)= 0.
ac( 217)= -0.
ac( 218)= -0.
ac( 219)= -0.
ac( 220)= -0.
ac( 221)= -0.
ac( 222)= -0.
ac( 223)= -0.
ac( 224)= -0.
ac( 225)= 0.
ac( 226)= 0.
ac( 227)= 0.
ac( 228)= 0.
ac( 229)= 0.
ac( 230)= 0.
ac( 231)= 0.
ac( 232)= 0.
ac( 233)= -0.
ac( 234)= -0.
ac( 235)= -0.
ac( 236)= -0.
ac( 237)= -0.
ac( 238)= -0.
ac( 239)= -0.
ac( 240)= 0.
ac( 241)= 0.
ac( 242)= 0.
ac( 243)= 0.
ac( 244)= -0.
ac( 245)= -0.
ac( 246)= -0.
ac( 247)= -0.
ac( 248)= 0.
ac( 249)= 0.
ac( 250)= 0.
ac( 251)= 0.
ac( 252)= 0.
ac( 253)= 0.
ac( 254)= 0.
ac( 255)= -0.
ac( 256)= -0.
ac( 257)= -0.
ac( 258)= -0.
ac( 259)= -0.
ac( 260)= -0.
ac( 261)= -0.
ac( 262)= -0.
ac( 263)= -0.
ac( 264)= -0.
ac( 265)= -0.
ac( 266)= -0.
ac( 267)= -0.
ac( 268)= -0.
ac( 269)= -0.
ac( 270)= -0.
ac( 271)= -0.
ac( 272)= -0.
ac( 273)= -0.
ac( 274)= -0.
ac( 275)= -0.
ac( 276)= -0.
ac( 277)= -0.
ac( 278)= 0.
ac( 279)= 0.
ac( 280)= 0.
ac( 281)= -0.
ac( 282)= -0.
ac( 283)= -0.
ac( 284)= -0.
ac( 285)= -0.
ac( 286)= -0.
ac( 287)= -0.
ac( 288)= -0.
ac( 289)= -0.
ac( 290)= 0.
ac( 291)= 0.
ac( 292)= 0.
ac( 293)= 0.
ac( 294)= 0.
ac( 295)= 0.
ac( 296)= 0.
ac( 297)= 0.
ac( 298)= 0.
ac( 299)= 0.
ac( 300)= 0.
ac( 301)= 0.
ac( 302)= 0.
ac( 303)= 0.
ac( 304)= 0.
ac( 305)= -0.
ac( 306)= -0.
ac( 307)= -0.
ac( 308)= -0.
ac( 309)= -0.
ac( 310)= -0.
ac( 311)= 0.
ac( 312)= 0.
ac( 313)= 0.
ac( 314)= 0.
ac( 315)= 0.
ac( 316)= 0.
ac( 317)= -0.
ac( 318)= -0.
ac( 319)= -0.
ac( 320)= -0.
ac( 321)= -0.
ac( 322)= 0.
ac( 323)= 0.
ac( 324)= 0.
ac( 325)= 0.
ac( 326)= 0.
ac( 327)= 0.
ac( 328)= 0.
ac( 329)= 0.
ac( 330)= -0.
ac( 331)= -0.
ac( 332)= -0.
ac( 333)= -0.
ac( 334)= -0.
ac( 335)= -0.
ac( 336)= 0.
ac( 337)= 0.
ac( 338)= 0.
ac( 339)= 0.
ac( 340)= 0.
ac( 341)= 0.
ac( 342)= 0.
ac( 343)= 0.
ac( 344)= 0.
ac( 345)= -0.
ac( 346)= -0.
ac( 347)= -0.
ac( 348)= -0.
ac( 349)= -0.
ac( 350)= -0.
ac( 351)= -0.
ac( 352)= -0.
ac( 353)= -0.
ac( 354)= 0.
ac( 355)= 0.
ac( 356)= 0.
ac( 357)= 0.
ac( 358)= 0.
ac( 359)= 0.
ac( 360)= 0.
ac( 361)= 0.
ac( 362)= 0.
ac( 363)= -0.
ac( 364)= -0.
ac( 365)= -0.
ac( 366)= -0.
ac( 367)= -0.
ac( 368)= 0.
ac( 369)= 0.
ac( 370)= 0.
ac( 371)= 0.
ac( 372)= 0.
ac( 373)= 0.
ac( 374)= 0.
ac( 375)= 0.
ac( 376)= 0.
ac( 377)= 0.
ac( 378)= 0.
ac( 379)= 0.
ac( 380)= 0.
ac( 381)= -0.
ac( 382)= -0.
ac( 383)= -0.
ac( 384)= -0.
ac( 385)= -0.
ac( 386)= 0.
ac( 387)= 0.
ac( 388)= 0.
ac( 389)= 0.
ac( 390)= 0.
ac( 391)= 0.
ac( 392)= 0.
ac( 393)= 0.
ac( 394)= 0.
ac( 395)= 0.
ac( 396)= -0.
ac( 397)= -0.
ac( 398)= -0.
ac( 399)= -0.
ac( 400)= -0.
ac( 401)= -0.
ac( 402)= -0.
ac( 403)= -0.
ac( 404)= 0.
ac( 405)= 0.
ac( 406)= 0.
ac( 407)= 0.
ac( 408)= 0.
ac( 409)= 0.
ac( 410)= 0.
ac( 411)= 0.
ac( 412)= 0.
ac( 413)= 0.
ac( 414)= 0.
ac( 415)= 0.
ac( 416)= 0.
ac( 417)= 0.
ac( 418)= -0.
ac( 419)= -0.
ac( 420)= -0.
ac( 421)= -0.
ac( 422)= -0.
ac( 423)= 0.
ac( 424)= 0.
ac( 425)= 0.
ac( 426)= 0.
ac( 427)= 0.
ac( 428)= 0.
ac( 429)= 0.
ac( 430)= -0.
ac( 431)= -0.
ac( 432)= -0.
ac( 433)= -0.
ac( 434)= -0.
ac( 435)= 0.
ac( 436)= 0.
ac( 437)= 0.
ac( 438)= 0.
ac( 439)= 0.
ac( 440)= 0.
ac( 441)= 0.
ac( 442)= -0.
ac( 443)= -0.
ac( 444)= -0.
ac( 445)= -0.
ac( 446)= -0.
ac( 447)= -0.
ac( 448)= -0.
ac( 449)= -0.
ac( 450)= 0.
ac( 451)= 0.
ac( 452)= 0.
ac( 453)= -0.
ac( 454)= -0.
ac( 455)= -0.
ac( 456)= -0.
ac( 457)= -0.
ac( 458)= 0.
ac( 459)= 0.
ac( 460)= 0.
ac( 461)= 0.
ac( 462)= 0.
ac( 463)= -0.
ac( 464)= -0.
ac( 465)= -0.
ac( 466)= 0.
ac( 467)= 0.
ac( 468)= 0.
ac( 469)= 0.
ac( 470)= 0.
ac( 471)= 0.
ac( 472)= 0.
ac( 473)= 0.
ac( 474)= 0.
ac( 475)= 0.
ac( 476)= 0.
ac( 477)= 0.
ac( 478)= -0.
ac( 479)= -0.
ac( 480)= -0.
ac( 481)= -0.
ac( 482)= 0.
ac( 483)= 0.
ac( 484)= 0.
ac( 485)= 0.
ac( 486)= 0.
ac( 487)= 0.
ac( 488)= 0.
ac( 489)= 0.
ac( 490)= 0.
ac( 491)= -0.
ac( 492)= -0.
ac( 493)= -0.
ac( 494)= -0.
ac( 495)= -0.
ac( 496)= -0.
ac( 497)= -0.
ac( 498)= -0.
ac( 499)= -0.
ac( 500)= -0.
再贴后500个点
ac( 501)= -0.
ac( 502)= -0.
ac( 503)= -0.
ac( 504)= -0.
ac( 505)= -0.
ac( 506)= -0.
ac( 507)= -0.
ac( 508)= -0.
ac( 509)= -0.
ac( 510)= -0.
ac( 511)= 0.
ac( 512)= 0.
ac( 513)= 0.
ac( 514)= 0.
ac( 515)= 0.
ac( 516)= 0.
ac( 517)= 0.
ac( 518)= -0.
ac( 519)= -0.
ac( 520)= -0.
ac( 521)= 0.
ac( 522)= 0.
ac( 523)= 0.
ac( 524)= 0.
ac( 525)= 0.
ac( 526)= -0.
ac( 527)= -0.
ac( 528)= -0.
ac( 529)= -0.
ac( 530)= -0.
ac( 531)= -0.
ac( 532)= -0.
ac( 533)= -0.
ac( 534)= -0.
ac( 535)= -0.
ac( 536)= -0.
ac( 537)= -0.
ac( 538)= -0.
ac( 539)= -0.
ac( 540)= -0.
ac( 541)= -0.
ac( 542)= -0.
ac( 543)= -0.
ac( 544)= -0.
ac( 545)= -0.
ac( 546)= -0.
ac( 547)= -0.
ac( 548)= -0.
ac( 549)= -0.
ac( 550)= -0.
ac( 551)= -0.
ac( 552)= 0.
ac( 553)= 0.
ac( 554)= 0.
ac( 555)= 0.
ac( 556)= 0.
ac( 557)= 0.
ac( 558)= 0.
ac( 559)= 0.
ac( 560)= 0.
ac( 561)= -0.
ac( 562)= -0.
ac( 563)= -0.
ac( 564)= -0.
ac( 565)= -0.
ac( 566)= -0.
ac( 567)= -0.
ac( 568)= -0.
ac( 569)= -0.
ac( 570)= -0.
ac( 571)= -0.
ac( 572)= 0.
ac( 573)= 0.
ac( 574)= 0.
ac( 575)= 0.
ac( 576)= -0.
ac( 577)= -0.
ac( 578)= -0.
ac( 579)= -0.
ac( 580)= -0.
ac( 581)= -0.
ac( 582)= -0.
ac( 583)= -0.
ac( 584)= -0.
ac( 585)= -0.
ac( 586)= -0.
ac( 587)= -0.
ac( 588)= -0.
ac( 589)= -0.
ac( 590)= -0.
ac( 591)= -0.
ac( 592)= 0.
ac( 593)= 0.
ac( 594)= 0.
ac( 595)= 0.
ac( 596)= 0.
ac( 597)= 0.
ac( 598)= -0.
ac( 599)= -0.
ac( 600)= -0.
ac( 601)= -0.
ac( 602)= -0.
ac( 603)= -0.
ac( 604)= -0.
ac( 605)= 0.
ac( 606)= 0.
ac( 607)= 0.
ac( 608)= 0.
ac( 609)= 0.
ac( 610)= 0.
ac( 611)= 0.
ac( 612)= 0.
ac( 613)= 0.
ac( 614)= 0.
ac( 615)= 0.
ac( 616)= 0.
ac( 617)= -0.
ac( 618)= -0.
ac( 619)= -0.
ac( 620)= -0.
ac( 621)= -0.
ac( 622)= 0.
ac( 623)= 0.
ac( 624)= 0.
ac( 625)= 0.
ac( 626)= 0.
ac( 627)= 0.
ac( 628)= -0.
ac( 629)= -0.
ac( 630)= -0.
ac( 631)= -0.
ac( 632)= -0.
ac( 633)= -0.
ac( 634)= -0.
ac( 635)= -0.
ac( 636)= -0.
ac( 637)= -0.
ac( 638)= -0.
ac( 639)= -0.
ac( 640)= -0.
ac( 641)= -0.
ac( 642)= -0.
ac( 643)= -0.
ac( 644)= -0.
ac( 645)= -0.
ac( 646)= -0.
ac( 647)= -0.
ac( 648)= -0.
ac( 649)= 0.
ac( 650)= 0.
ac( 651)= 0.
ac( 652)= 0.
ac( 653)= 0.
ac( 654)= 0.
ac( 655)= 0.
ac( 656)= 0.
ac( 657)= 0.
ac( 658)= 0.
ac( 659)= 0.
ac( 660)= 0.
ac( 661)= -0.
ac( 662)= -0.
ac( 663)= -0.
ac( 664)= -0.
ac( 665)= -0.
ac( 666)= -0.
ac( 667)= -0.
ac( 668)= -0.
ac( 669)= -0.
ac( 670)= -0.
ac( 671)= -0.
ac( 672)= -0.
ac( 673)= -0.
ac( 674)= 0.
ac( 675)= 0.
ac( 676)= 0.
ac( 677)= 0.
ac( 678)= 0.
ac( 679)= 0.
ac( 680)= 0.
ac( 681)= 0.
ac( 682)= 0.
ac( 683)= 0.
ac( 684)= 0.
ac( 685)= 0.
ac( 686)= 0.
ac( 687)= 0.
ac( 688)= 0.
ac( 689)= 0.
ac( 690)= 0.
ac( 691)= 0.
ac( 692)= 0.
ac( 693)= 0.
ac( 694)= 1.
ac( 695)= 1.
ac( 696)= 1.
ac( 697)= 1.
ac( 698)= 1.
ac( 699)= 1.
ac( 700)= 1.
ac( 701)= 1.
ac( 702)= 1.
ac( 703)= 1.
ac( 704)= 1.
ac( 705)= 1.
ac( 706)= 1.
ac( 707)= 1.
ac( 708)= 1.
ac( 709)= 1.
ac( 710)= 1.
ac( 711)= 1.
ac( 712)= 1.
ac( 713)= 1.
ac( 714)= 1.
ac( 715)= 1.
ac( 716)= 1.
ac( 717)= 1.
ac( 718)= 1.
ac( 719)= 1.
ac( 720)= 1.
ac( 721)= 1.
ac( 722)= 1.
ac( 723)= 1.
ac( 724)= 1.
ac( 725)= 1.
ac( 726)= 1.
ac( 727)= 1.
ac( 728)= 0.
ac( 729)= 0.
ac( 730)= 0.
ac( 731)= 0.
ac( 732)= 0.
ac( 733)= -0.
ac( 734)= -0.
ac( 735)= -0.
ac( 736)= -0.
ac( 737)= -0.
ac( 738)= -0.
ac( 739)= -0.
ac( 740)= -0.
ac( 741)= -0.
ac( 742)= -0.
ac( 743)= -1.
ac( 744)= -1.
ac( 745)= -1.
ac( 746)= -1.
ac( 747)= -1.
ac( 748)= -1.
ac( 749)= -2.
ac( 750)= -2.
ac( 751)= -2.
ac( 752)= -2.
ac( 753)= -2.
ac( 754)= -2.
ac( 755)= -3.
ac( 756)= -3.
ac( 757)= -3.
ac( 758)= -3.
ac( 759)= -3.
ac( 760)= -3.
ac( 761)= -3.
ac( 762)= -2.
ac( 763)= -2.
ac( 764)= -2.
ac( 765)= -2.
ac( 766)= -2.
ac( 767)= -2.
ac( 768)= -2.
ac( 769)= -2.
ac( 770)= -2.
ac( 771)= -2.
ac( 772)= -2.
ac( 773)= -2.
ac( 774)= -2.
ac( 775)= -2.
ac( 776)= -2.
ac( 777)= -2.
ac( 778)= -2.
ac( 779)= -2.
ac( 780)= -2.
ac( 781)= -2.
ac( 782)= -2.
ac( 783)= -2.
ac( 784)= -2.
ac( 785)= -1.
ac( 786)= -1.
ac( 787)= -0.
ac( 788)= -0.
ac( 789)= -0.
ac( 790)= 0.
ac( 791)= 0.
ac( 792)= 0.
ac( 793)= 0.
ac( 794)= 0.
ac( 795)= 0.
ac( 796)= 1.
ac( 797)= 1.
ac( 798)= 1.
ac( 799)= 1.
ac( 800)= 1.
ac( 801)= 1.
ac( 802)= 1.
ac( 803)= 0.
ac( 804)= 0.
ac( 805)= 0.
ac( 806)= 0.
ac( 807)= 0.
ac( 808)= 0.
ac( 809)= 0.
ac( 810)= 0.
ac( 811)= 0.
ac( 812)= 0.
ac( 813)= 0.
ac( 814)= 0.
ac( 815)= 1.
ac( 816)= 1.
ac( 817)= 1.
ac( 818)= 1.
ac( 819)= 1.
ac( 820)= 1.
ac( 821)= 2.
ac( 822)= 2.
ac( 823)= 2.
ac( 824)= 2.
ac( 825)= 1.
ac( 826)= 1.
ac( 827)= 1.
ac( 828)= 1.
ac( 829)= 1.
ac( 830)= 0.
ac( 831)= 0.
ac( 832)= 1.
ac( 833)= 1.
ac( 834)= 1.
ac( 835)= 1.
ac( 836)= 1.
ac( 837)= 1.
ac( 838)= 1.
ac( 839)= 0.
ac( 840)= 0.
ac( 841)= 0.
ac( 842)= 0.
ac( 843)= 1.
ac( 844)= 1.
ac( 845)= 1.
ac( 846)= 1.
ac( 847)= 1.
ac( 848)= 1.
ac( 849)= 1.
ac( 850)= 1.
ac( 851)= 0.
ac( 852)= 0.
ac( 853)= 0.
ac( 854)= 0.
ac( 855)= 0.
ac( 856)= -0.
ac( 857)= -0.
ac( 858)= -0.
ac( 859)= -1.
ac( 860)= -1.
ac( 861)= -1.
ac( 862)= -1.
ac( 863)= -1.
ac( 864)= -1.
ac( 865)= -0.
ac( 866)= -0.
ac( 867)= 0.
ac( 868)= 0.
ac( 869)= 0.
ac( 870)= 0.
ac( 871)= 0.
ac( 872)= 0.
ac( 873)= 0.
ac( 874)= 0.
ac( 875)= 0.
ac( 876)= -0.
ac( 877)= -0.
ac( 878)= -0.
ac( 879)= -1.
ac( 880)= -1.
ac( 881)= -1.
ac( 882)= -1.
ac( 883)= -0.
ac( 884)= -0.
ac( 885)= -0.
ac( 886)= 0.
ac( 887)= 0.
ac( 888)= 1.
ac( 889)= 1.
ac( 890)= 1.
ac( 891)= 1.
ac( 892)= 1.
ac( 893)= 1.
ac( 894)= 1.
ac( 895)= 0.
ac( 896)= 0.
ac( 897)= 0.
ac( 898)= -0.
ac( 899)= -0.
ac( 900)= -0.
ac( 901)= -0.
ac( 902)= -1.
ac( 903)= -1.
ac( 904)= -1.
ac( 905)= -1.
ac( 906)= -1.
ac( 907)= -1.
ac( 908)= -1.
ac( 909)= -1.
ac( 910)= -1.
ac( 911)= -1.
ac( 912)= -1.
ac( 913)= -1.
ac( 914)= -1.
ac( 915)= -1.
ac( 916)= -1.
ac( 917)= -0.
ac( 918)= -0.
ac( 919)= -0.
ac( 920)= -0.
ac( 921)= -0.
ac( 922)= 0.
ac( 923)= 0.
ac( 924)= 0.
ac( 925)= 0.
ac( 926)= 0.
ac( 927)= 0.
ac( 928)= 0.
ac( 929)= 0.
ac( 930)= 0.
ac( 931)= 0.
ac( 932)= 0.
ac( 933)= 0.
ac( 934)= 0.
ac( 935)= 0.
ac( 936)= 0.
ac( 937)= 0.
ac( 938)= 0.
ac( 939)= 0.
ac( 940)= 0.
ac( 941)= -0.
ac( 942)= -0.
ac( 943)= -0.
ac( 944)= -0.
ac( 945)= -0.
ac( 946)= 0.
ac( 947)= 0.
ac( 948)= 0.
ac( 949)= 0.
ac( 950)= 0.
ac( 951)= 0.
ac( 952)= 0.
ac( 953)= -0.
ac( 954)= -0.
ac( 955)= -0.
ac( 956)= -0.
ac( 957)= -0.
ac( 958)= -0.
ac( 959)= -0.
ac( 960)= -0.
ac( 961)= -0.
ac( 962)= -0.
ac( 963)= -0.
ac( 964)= 0.
ac( 965)= 0.
ac( 966)= 0.
ac( 967)= 0.
ac( 968)= 1.
ac( 969)= 1.
ac( 970)= 1.
ac( 971)= 1.
ac( 972)= 1.
ac( 973)= 1.
ac( 974)= 1.
ac( 975)= 1.
ac( 976)= 1.
ac( 977)= 1.
ac( 978)= 0.
ac( 979)= 0.
ac( 980)= 0.
ac( 981)= -0.
ac( 982)= -0.
ac( 983)= -0.
ac( 984)= -0.
ac( 985)= -0.
ac( 986)= -0.
ac( 987)= -0.
ac( 988)= -0.
ac( 989)= 0.
ac( 990)= 0.
ac( 991)= 0.
ac( 992)= 0.
ac( 993)= 0.
ac( 994)= 0.
ac( 995)= 0.
ac( 996)= 0.
ac( 997)= 0.
ac( 998)= 0.
ac( 999)= 0.
ac(1000)= 0.
54单元描述变截面梁的例子
作者：清华大学土木系，陆新征
A_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.6667,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.6667,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,0
MPDATA,EX,1,,200e3
MPDATA,PRXY,1,,.2
N,1,0,,,,,,
N,2,5,,,,,,
N,3,10,,,,,,
变截面梁实例
et,1,beam188
keyopt,1,8,2
keyopt,1,9,2
mp,ex,1,3.5e10
mp,dens,1,2549
mp,prxy,1,0.1667
sectype,1,beam,rect ! define cross section at first end point
sectype,2,beam,rect ! define cross section at far end
secdata,3,0.5
sectype,3,taper ! new Section ID for tapered beam analysis
secdata,1,0.0,0.0 ! section 1 at location (0,0,0)
secdata,2,0.0,,-20.0 ! section 2 at location (0,-20,0)
k, 1, -2.7 ,0, -27.85
k, 2, -2.7 ,0, -19
k, 3, -2.7 ,0, -13
lstr, 1, 2
lstr, 2, 3
lsel,s,,,1
latt,1,,1,,,,1
lesize,all,,,1
lsel,s,,,2
latt,1,,1,,,,3
lesize,all,,,1
/view,,1,2,3
拱桥浇筑过程分析-单元生死应用实例
et,1,BEAM3
mp,ex,1,2.1e11
mp,prxy,1,.3
mp,dens,1,7800
r,1,25.48e-4 ,,1.5
K,2,10,5,,
K,3,20,0,,
LARC, 1, 3, 2
LESIZE,ALL, , ,20, ,1, , ,1,
n,22,3.3986,6
n,23,7.695,6
n,24,12.305,6
n,25,16.601,6
nropt,full,,on
esel,s,,,25,27,1
antype,,rest,,,
f,12,fy,-1000
nSOL,2,12,u,y,
PLVAR,1, , , , , , , , , ,
简支梁实体与预应力钢筋分析实例
/COM, Structural
egjx=2e5 !Ey
agjx=140 !单根钢绞线面积
ehnt=4e4 !Eh
xzxs=1.0e-5 !线胀系数
yjl=200000 !定义预加力
et,1,link8 !定义link8单元
et,2,solid95 !定义solid95单元
r,1,agjx !定义link8单元的面积
r,2 !定义第2种实常数
mp,ex,1,egjx !定义link8单元的弹性模量
mp,prxy,1,0.3 !定义link8单元的泊松系数
mp,alpx,1,1.0e-5 !定义线膨胀系数
mp,ex,2,ehnt !定义solid95单元的弹性模量
mp,prxy,2,0.3 !定义solid95单元的泊松系数
blc4, , ,100,200,3000 !定义梁体
/view,1,1,1,1 !定义ISO查看
vplot !绘制梁体
kwpave,6 !工作平面移动到关键点6
wpoff,-30 !工作平面移动-30mm(X)
wprot,0,0,90 !工作平面旋转
vsbw,1 !分割梁体
wpoff,0,0,-40 !工作平面移动-40mm(Z)
vsbw,2 !分割梁体
wpoff,0,40 !工作平面移动40mm(Y)
wprot,0,90 !工作平面旋转
vsbw,all !分割梁体
wpstyl !关闭工作平面显示
nummrg,all,,,,low !整理
numcmp,all !压缩编号
esize,30 !定义网分时边长控制
lsel,s,,,28,38,10 !定义line28和38为新的选择集
latt,1,1,1 !定义选择集的属性
lmesh,all !对线划分单元
allsel,all !新的选择集为所有的实体
gplot !绘制所有的实体
vsel,s,,,all !定义所有体为选择集
vatt,2,2,2 !定义选择集的属性
mshape,0,3d !将体划分单元的形状定位HEX
mshkey,1 !采用MAPPED划分器
vmesh,all !对体进行划分单元
dl,3,,all !对线line7施加约束(UX,UY,UZ)
dl,16,,all !对线line31施加约束(UX,UY,UZ)
dl,23,,all !对线line23施加约束(UX,UY,UZ)
dl,2,,uy !对线line4施加约束(UY)
dl,15,,uy !对线line30施加约束(UY)
dl,22,,uy !对线line23施加约束(UY)
dk,2,,,,,ux,uy !对关键点2约束(UX,UY)
bfl,28,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度
bfl,38,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度
solve !求解
plnsol,s,z,0,1 !绘制Z方向的应力
etable,sigi,ls,1 !定义钢筋单元数据表
plls,sigi,sigi,1 !绘制上述应力
!/exit,nosav
简单的二维焊接分析－单元生死实例
下面的命令流进行的是一个简单的二维焊接分析, 利用ANSYS单元生死和热-结构耦合分析功能进行焊接过程仿真,
计算焊接过程中的温度分布和应力分布以及冷却后的焊缝残余应力。
/title,Weld Analysis by "Element Birth and Death"
/unit,si ! 采用国际单位制
et,1,13,4 ! 13号二维耦合单元, 同时具有温度和位移自由度
! 1号材料是钢
! 2号材料是铝
! 3号材料是铜
! 铝是本次分析中的焊料, 它将钢结构部分和铜结构部分焊接起来
! 下面是在几个温度点下, 各材料的弹性模量
mptemp,1,20,500,00
mpdata,ex,1,1,1.93e11,1.50e11,0.70e11,0.10e11,0.01e11
mpdata,ex,2,1,1.02e11,0.50e11,0.08e11,0.001e11,0.0001e11
mpdata,ex,3,1,1.17e11,0.90e11,0.30e11,0.05e11,0.005e11
! 假设各材料都是双线性随动硬化弹塑性本构关系
! 下面是各材料在各温度点下的屈服应力和屈服后的弹性模量
tb,bkin,1,5
tbtemp,20,1
tbdata,1,.193e11
tbtemp,500,2
tbdata,1, 933e6,0.150e11
tbtemp,1000,3
tbdata,1, 435e6,0.070e11
tbtemp,1500,4
tbdata,1, 70e6,0.010e11
tbtemp,2000,5
tbdata,1, 7e6,0.001e11
tb,bkin,2,5
tbtemp,20,1
tbdata,1,800e6,0.102e11
tbtemp,500,2
tbdata,1,400e6,0.050e11
tbtemp,1000,3
tbdata,1, 70e6,0.008e11
tbtemp,1500,4
tbdata,1, 1e6,0.0001e11
tbtemp,2000,5
tbdata,1,0.1e6,0.00001e11
tb,bkin,3,5
tbtemp,20,1
tbdata,1,900e6,0.117e11
tbtemp,500,2
tbdata,1,700e6,0.090e11
tbtemp,1000,3
tbdata,1,230e6,0.030e11
tbtemp,1500,4
tbdata,1, 40e6,0.005e11
tbtemp,2000,5
tbdata,1, 4e6,0.0005e11
! 材料密度(假设为常值)
mp,dens,1,8030
mp,dens,2,4850
mp,dens,3,8900
! 热膨胀系数(假设为常值)
mp,alpx,1,1.78e-5
mp,alpx,2,9.36e-6
mp,alpx,3,1.66e-5
! 泊松比(假设为常值)
mp,nuxy,1,0.29
mp,nuxy,2,0.30
mp,nuxy,3,0.30
! 热传导系数(假设为常值)
mp,kxx,1,16.3
mp,kxx,2,7.44
mp,kxx,3,393
! 比热(假设为常值)
mp,c,1,502
mp,c,2,544
mp,c,3,385
! 热膨胀系数(假设为常值)
! 由于该13号单元还有磁自由度, 此处假设一磁特性, 但并不影响我们所关心的结果
mp,murx,1,1
mp,murx,2,1
mp,murx,3,1
! 假设焊料(铝)焊上去后的初始温度是;
mp,reft,1,20
mp,reft,2,1500
mp,reft,3,20
! 下面建立几何模型
k,2,0.5,0,0
k,4,0,0.3,0
k,5,0.35,0.3,0
k,6,0.65,0.3,0
k,7,1,0.3,0
! 划分单元
esize,0.025
esize,0.05
antype,4 ! 瞬态分析
trnopt,full
! 在模型的左边界加位移约束
*get,minx,node,,mnloc,x
nsel,s,loc,x,minx
d,all,ux,0
*get,miny,node,,mnloc,y
nsel,r,loc,y,miny
d,all,uy,0
! 假设在模型的左右边界处温度始终保持在20℃, 其它边界条件如对流、辐射等都可
! 施加，此处因只是示意而已，故只施加恒温边界条件。
*get,minx,node,,mnloc,x
nsel,s,loc,x,minx
d,all,temp,20
*get,maxx,node,,mxloc,x
nsel,s,loc,x,maxx
d,all,temp,20
! 由于第二个面是焊料所在区域，因此，首先将该区域的单元"死"掉
*get,emax,elem,,num,max
asel,s,area,,nna
*get,nse,elem,,count
*dim,ne,,nse
*dim,nex,,nse
*dim,ney,,nse
*dim,neorder,,nse
! 下面的DO循环用于将焊料区的单元按其形心Y坐标排序，以便后面
! 模拟焊料由下向上逐步"生长"的过程
*do,i1,1,nse
esel,u,elem,,mine
*get,nse1,elem,,count
*do,i,1,emax
*if,esel(i),eq,1,then
*do,i,1,nse1
*get,ney(i),elem,ne(i),cent,y
*get,nex(i),elem,ne(i),cent,x
*do,i,1,nse1
*if,ney(i),lt,miny,then
miny=ney(i)
minx=nex(i)
mine=ne(i)
*if,ney(i),eq,miny,then
*if,nex(i),lt,minx,then
miny=ney(i)
minx=nex(i)
mine=ne(i)
neorder(i1)=mine
max_tem=1500 ! 按照前面假设，焊料的初始温度为;
dt1=1e-3 ! 用于建立初始条件的一个很小的时间段
dt=5 ! 焊接一个单元所需的时间
t=0 ! 起始时间
*do,i,1,nse
ekill,neorder(i)
esel,s,live
allsel,all
outres,all,all
ic,all,temp,20
timint,0,struct
timint,1,therm
timint,0,mag
tintp,0.005,,,1,0.5,0.2
*do,i,1,nse
ealive,neorder(i)
esel,s,live
! 下面的求解用于建立温度的初始条件
d,nelem(neorder(i),j),temp,max_tem
! 下面的求解用于保证初始的升温速度为零
! 下面的步骤用于求解温度分布
ddele,nelem(neorder(i),j),temp
t=t+dt-2*dt1
nsubst,nsub1
nsubst,nsub2
! 下面的一系列命令用于生成应力的动画文件
/cont,1,15,0,,1200e6
/dscale,1,1.0
avprin,0,0
/seg,multi,stress,0.1
*do,i,1,nse
esel,u,elem,,neorder(i)
*do,i,1,nse
esel,a,elem,,neorder(i)
set,(i-1)*3+1,1
plnsol,s,eqv
*do,j,1,nsub1
set,(i-1)*3+3,j
plnsol,s,eqv
*do,i,1,nsub2
set,(nse-1)*3+4,i
plnsol,s,eqv
/seg,off,stress,0.1
/anfile,save,stress,avi
! 下面的一系列命令用于生成温度的动画文件
/cont,1,15,0,0
/dscale,1,1.0
avprin,0,0
/seg,multi,temp,0.1
*do,i,1,nse
esel,u,elem,,neorder(i)
*do,i,1,nse
esel,a,elem,,neorder(i)
set,(i-1)*3+1,1
plnsol,temp
*do,j,1,nsub1
set,(i-1)*3+3,j
plnsol,temp
*do,i,1,nsub2
set,(nse-1)*3+4,i
plnsol,temp
/seg,off,temp,0.1
/anfile,save,temp,avi
隧道开挖（三维）的命令流
/FILENAME,'TUNNEL'
/CONFIG,NRES,1E6
!--------------------------------------------------------------------------------------------------------------------!
TM_YS=16*1E9&&&&&&
MU_YS=0.24&&&&&&&&
C_YS=1.80*1E6&&&&&&
FI_YS=43&&&&&&&&&&
MD_YS=0.027*1E6/10
TM_GXW=26*1E9&&&&&&
MU_GXW=0.21&&&&&&&&
C_GXW=2.4*1E6&&&&&&
FI_GXW=50&&&&&&&&&&
MD_GXW=0./10
TM_HNT=22*1E9&&&&&&
MU_HNT=0.20&&&&&&&&
C_HNT=1.73*1E6&&&&&&
FI_HNT=49&&&&&&&&&&
MD_HNT=0.024*1E6/10
TM_MG=TM_YS&&&&&&
MU_MG=MU_YS&&&&&&&&&&
MD_MG=MD_YS&&&&&&&&&&
NT_MG=90*1E3&&&&&&&&&&
C_MG=C_YS*1.10&&&&&&&&
FI_MG=FI_YS+5&&&&&&&&
ZH_MG=NT_MG/(1.2*1*1)
!----------------------------------------------------------------------------------------------------------------------!
R_HJ=6.230&&&&&&&&&&
H_ZQ=2.670&&&&&&&&&&
L_ZTZ=40&&&&&&&&&&&&
L_YTZ=40&&&&&&&&&&&&
Y_YJMIN=400&&&&&&&&&&
Y_YJMAX=415&&&&&&&&&&
L_YJ=Y_YJMAX-Y_YJMIN
L_MG=3.5&&&&&&&&&&&&
JJ_MG=5.0&&&&&&&&&&&&
PJ_MG=5.0&&&&&&&&&&
H_CENDB=1.47&&&&&&&&
PD=2/1000&&&&&&&&&&&&
W_GXW=0.2&&&&&&&&&&&&
H_JS=831&&&&&&&&&&&&
DMJL_SD=40&&&&&&&&&&
H_CTOB=1.47&&&&&&&&&&
*DIM,CEN_ZSD,ARRAY,2,1
*DIM,CEN_YSD,ARRAY,2,1
CEN_ZSD(1,1)=-DMJL_SD/2,880
CEN_YSD(1,1)=DMJL_SD/2,880 !+H_CTOB
H_ZQ=2.67&&&&&&&&&&&&&&&&
FLAG_RD=3&&&&&&&&&&
JL_RD=5&&&&&&&&&&&&&&
*DIM,Z_RDCC,ARRAY,3,2
Z_RDCC(1,1)=3,3,3
Z_RDCC(1,2)=CEN_ZSD(1,1),CEN_ZSD(2,1)+R_HJ+JL_RD+Z_RDCC(2,1)/2,L_ZTZ+7.5
*DIM,Y_RDCC,ARRAY,3,2
Y_RDCC(1,1)=3,3,3
Y_RDCC(1,2)=CEN_YSD(1,1),CEN_YSD(2,1)+R_HJ+JL_RD+Y_RDCC(2,1)/2,L_ZTZ+7.5
*DIM,DBKZDX,ARRAY,9,1
*DIM,DBKZDZ,ARRAY,(Y_YJMAX-Y_YJMIN)/PJ_MG+3,1
DBKZDX(1,1)=-150,-64,-38,-20,0,20,42,64,150&&
DBKZDZ(1,1)=0,L_ZTZ
DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1)=L_ZTZ+Y_YJMAX-Y_YJMIN+L_YTZ
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG
&DBKZDZ(I+2,1)=L_ZTZ+PJ_MG*I
*DIM,DBKZDGC,ARRAY,9,(Y_YJMAX-Y_YJMIN)/PJ_MG+3
DBKZDGC(1,1)=929,915,915,920,925,930,937,940,940
DBKZDGC(1,2)=929,915,915,920,925,930,937,940,940
DBKZDGC(1,3)=929,915,915,920,925,930,937,940,940
DBKZDGC(1,4)=929,915,915,920,925,930,937,940,940
DBKZDGC(1,5)=929,915,915,920,925,930,937,940,940
DBKZDGC(1,6)=929,915,915,920,925,930,937,940,940
!*************************************************************************************!
!----------------------------------------------------------------------------------------------------------------------!
KEYOPT,2,5,2
MP,EX,7,TM_YS
MP,PRXY,7,MU_YS
MP,DENS,7,MD_YS
MP,EX,1,TM_YS
MP,PRXY,1,MU_YS
MP,DENS,1,MD_YS
MP,EX,2,TM_GXW
MP,PRXY,2,MU_GXW
MP,DENS,2,MD_GXW
MP,EX,3,TM_MG
MP,PRXY,3,MU_MG
MP,DENS,3,MD_MG
MP,EX,4,TM_HNT
MP,PRXY,4,MU_HNT
MP,DENS,4,MD_HNT
MP,EX,5,TM_YS
MP,PRXY,5,MU_YS
MP,DENS,5,MD_YS
MP,EX,6,TM_YS
MP,PRXY,6,MU_YS
TB,DP,1,,,&&
TBMODIF,1,1,C_YS&&
TBMODIF,1,2,FI_YS&&&&
TB,DP,2,,,&&
TBMODIF,1,1,C_GXW&&
TBMODIF,1,2,FI_GXW&&&&
TB,DP,3,,,&&
TBMODIF,1,1,C_MG&&
TBMODIF,1,2,FI_MG&&&&&&
TB,DP,4,,,&&
TBMODIF,1,1,C_HNT&&
TBMODIF,1,2,FI_HNT&&
TB,DP,5,,,&&
TBMODIF,1,1,C_YS&&
TBMODIF,1,2,FI_YS&&
TB,DP,6,,,&&
TBMODIF,1,1,C_YS&&
TBMODIF,1,2,FI_YS
!**************************************************!
MAT_YS=1&&
MAT_MGKW=6
MAT_TXYS=1&&
*DO,I,1,9-1
&*DO,J,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1&
&&K,KBEGIN+1,DBKZDX(I,1),H_JS,DBKZDZ(J,1)
&&K,KBEGIN+2,DBKZDX(I,1),DBKZDGC(I,J),DBKZDZ(J,1)
&&K,KBEGIN+3,DBKZDX(I,1),DBKZDGC(I,J+1),DBKZDZ(J+1,1)
&&K,KBEGIN+4,DBKZDX(I,1),H_JS,DBKZDZ(J+1,1)
&&K,KBEGIN+5,DBKZDX(I+1,1),H_JS,DBKZDZ(J,1)
&&K,KBEGIN+6,DBKZDX(I+1,1),DBKZDGC(I+1,J),DBKZDZ(J,1)
&&K,KBEGIN+7,DBKZDX(I+1,1),DBKZDGC(I+1,J+1),DBKZDZ(J+1,1)
&&K,KBEGIN+8,DBKZDX(I+1,1),H_JS,DBKZDZ(J+1,1)
&&V,KBEGIN+1,KBEGIN+2,KBEGIN+3,KBEGIN+4,KBEGIN+5,KBEGIN+6,KBEGIN+7,KBEGIN+8
&&KBEGIN=KBEGIN+8
NUMCMP,ALL
VATT,MAT_YS&&
CYL4,CEN_ZSD(1,1),CEN_ZSD(2,1),R_HJ+L_MG,180,,,L_ZTZ+L_YTZ+Y_YJMAX-Y_YJMIN
CYL4,CEN_YSD(1,1),CEN_YSD(2,1),R_HJ+L_MG,180,,,L_ZTZ+L_YTZ+Y_YJMAX-Y_YJMIN
BLOCK,CEN_ZSD(1,1)-R_HJ-L_MG,CEN_ZSD(1,1)+R_HJ+L_MG,CEN_ZSD(2,1)-H_ZQ,CEN_ZSD(2,1),DBKZDZ(1,1),DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1)
BLOCK,CEN_YSD(1,1)-R_HJ-L_MG,CEN_YSD(1,1)+R_HJ+L_MG,CEN_YSD(2,1)-H_ZQ,CEN_YSD(2,1),DBKZDZ(1,1),DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1)
VSEL,S,MAT,,MAT_YS
CM,V_BSUB,VOLU&&&&&&
VSEL,INVE&&&&&&&&&&
CM,V_SUB,VOLU
VSBV,V_BSUB,V_SUB
CMDEL,V_BSUB&&&&&&&&
CMDEL,V_SUB
NUMCMP,ALL
VATT,MAT_YS&&
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&WPAVE,CEN_ZSD(1,1),CEN_ZSD(2,1),DBKZDZ(I,1)
&CYL4,0,0,R_HJ+L_MG,180,R_HJ,,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&
&CYL4,0,0,R_HJ,180,,,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&&&&&&&&&&&
&BLOCK,-R_HJ,R_HJ,-H_CTOB,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&
&BLOCK,-R_HJ,R_HJ,-H_ZQ,-H_CTOB,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&
&BLOCK,-R_HJ-L_MG,-R_HJ,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&
&BLOCK,R_HJ,R_HJ+L_MG,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&
&WPAVE,CEN_YSD(1,1),CEN_YSD(2,1),DBKZDZ(I,1)
&CYL4,0,0,R_HJ+L_MG,180,R_HJ,,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&
&CYL4,0,0,R_HJ,180,,,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&&&&&&&&&&&
&BLOCK,-R_HJ,R_HJ,-H_CTOB,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&&&
&BLOCK,-R_HJ,R_HJ,-H_ZQ,-H_CTOB,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&
&BLOCK,-R_HJ-L_MG,-R_HJ,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&
&BLOCK,R_HJ,R_HJ+L_MG,-H_ZQ,0,0,DBKZDZ(I+1,1)-DBKZDZ(I,1)&&&&
VGLUE,ALL&&
NUMCMP,ALL
CM,ROCK,VOLU
WPAVE,0,0,0
*IF,FLAG_RD,EQ,-1,OR,FLAG_RD,EQ,2,THEN
&VSTART=1000
&NUMSTR,VOLU,VSTART
&BLOCK,Z_RDCC(1,2)-Z_RDCC(1,1)/2,Z_RDCC(1,2)+Z_RDCC(1,1)/2,Z_RDCC(2,2)-Z_RDCC(2,1)/2,Z_RDCC(2,2)+Z_RDCC(2,1)/2,Z_RDCC(3,2)-Z_RDCC(3,1)/2,Z_RDCC(3,2)+Z_RDCC(3,1)/2
&VSBV,ROCK,VSTART
&NUMCMP,VOLU
&CMDEL,ROCK
&CM,ROCK,VOLU
*IF,FLAG_RD,EQ,1,OR,FLAG_RD,EQ,2,THEN
&VSTART=1000
&NUMSTR,VOLU,VSTART
&BLOCK,Y_RDCC(1,2)-Y_RDCC(1,1)/2,Y_RDCC(1,2)+Y_RDCC(1,1)/2,Y_RDCC(2,2)-Y_RDCC(2,1)/2,Y_RDCC(2,2)+Y_RDCC(2,1)/2,Y_RDCC(3,2)-Y_RDCC(3,1)/2,Y_RDCC(3,2)+Y_RDCC(3,1)/2
&VSBV,ROCK,VSTART
&NUMCMP,VOLU
&CMDEL,ROCK
&CM,ROCK,VOLU
NUMCMP,ALL
!----------------------------------------------------------------------------------------------------------------------!
VATT,MAT_YS
VSEL,S,LOC,X,CEN_ZSD(1,1)-1E-7,CEN_ZSD(1,1)+1E-7&&&&&&&&
CM,CM_TEMPV,VOLU
ASEL,R,LOC,Y,CEN_ZSD(2,1)-H_ZQ-1E-7,CEN_ZSD(2,1)-H_ZQ+1E-7&&&&
CM,CM_ZDB,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_TEMPV
ASEL,R,LOC,Y,CEN_ZSD(2,1)-1E-7,CEN_ZSD(2,1)+1E-7&&&&&&&&&&
ASEL,U,LOC,X,CEN_ZSD(1,1)-R_HJ,CEN_ZSD(1,1)-R_HJ-L_MG&&&&
ASEL,U,LOC,X,CEN_ZSD(1,1)+R_HJ,CEN_ZSD(1,1)+R_HJ+L_MG
CM,CM_ZSD,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_ZSD
CMSEL,A,CM_ZDB
VATT,MAT_SD
CM,CM_ZSDKW,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_TEMPV
VSEL,U,MAT,,MAT_SD&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
ASEL,R,LOC,Y,CEN_ZSD(2,1)
VATT,MAT_MGKW
CM,CM_ZMG,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
VSEL,S,LOC,X,CEN_YSD(1,1)-1E-7,CEN_YSD(1,1)+1E-7&&&&&&&&
CM,CM_TEMPV,VOLU
ASEL,R,LOC,Y,CEN_YSD(2,1)-H_ZQ-1E-7,CEN_YSD(2,1)-H_ZQ+1E-7&&&&
CM,CM_YDB,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_TEMPV
ASEL,R,LOC,Y,CEN_YSD(2,1)-1E-7,CEN_YSD(2,1)+1E-7&&&&&&&&&&
ASEL,U,LOC,X,CEN_YSD(1,1)-R_HJ,CEN_YSD(1,1)-R_HJ-L_MG&&&&
ASEL,U,LOC,X,CEN_YSD(1,1)+R_HJ,CEN_YSD(1,1)+R_HJ+L_MG
CM,CM_YSD,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_YSD
CMSEL,A,CM_YDB
VATT,MAT_SD
CM,CM_YSDKW,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMSEL,S,CM_TEMPV
VSEL,U,MAT,,MAT_SD&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
ASEL,R,LOC,Y,CEN_YSD(2,1)
VATT,MAT_MGKW
CM,CM_YMG,VOLU&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
CMDEL,CM_TEMPV
ASEL,S,LOC,X,DBKZDX(1,1)
ASEL,A,LOC,X,DBKZDX(9,1)
DA,ALL,UX,0&&&&&&&&
ASEL,S,LOC,Z,DBKZDZ(1,1)
ASEL,A,LOC,Z,DBKZDZ((Y_YJMAX-Y_YJMIN)/PJ_MG+3,1)
DA,ALL,UZ,0&&&&&&
ASEL,S,LOC,Y,H_JS
DA,ALL,ALL&&&&&&&&
*DIM,ZLT_HFX,ARRAY,LTS_ZD,1
ZLT_HFX(1,1)=5,3 !4,2
*DIM,YLT_HFX,ARRAY,LTS_YD,1
YLT_HFX(1,1)=5,3 !4,2
*DIM,LT_HFZ,ARRAY,(Y_YJMAX-Y_YJMIN)/PJ_MG+2
LT_HFZ(1,1)=4&&&&&&&&&&&&&&&&&&
*DO,I,2,(Y_YJMAX-Y_YJMIN)/PJ_MG+1&&
<_HFZ(I,1)=2 !PJ_MG
LT_HFZ((Y_YJMAX-Y_YJMIN)/PJ_MG+2,1)=4 !10
*DO,J,1,LTS_ZD
&*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&&X_VCEN=(DBKZDX(J,1)+DBKZDX(J+1,1))/2
&&Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2
&&VSEL,S,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7
&&VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7
&&CM,CM_HFX,LINE
&&LSEL,R,LOC,X,X_VCEN&&&&&&&&&&
&&LESIZE,ALL,,,ZLT_HFX(J,1),,1
&&CMSEL,S,CM_HFX
&&LSEL,R,LOC,Z,Z_VCEN
&&LESIZE,ALL,,,LT_HFZ(I,1),,1
&&CMSEL,S,CM_HFX
&&LSEL,U,LOC,X,X_VCEN
&&LSEL,U,LOC,Z,Z_VCEN
&&LESIZE,ALL,,,9,,1&&
&&CMSEL,S,CM_HFX
&&MSHAPE,0,3D
&&MSHKEY,1
&&VMESH,ALL
&*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&&X_VCEN=(DBKZDX(9-J+1,1)+DBKZDX(9-J,1))/2
&&Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2
&&VSEL,S,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7
&&VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7
&&CM,CM_HFX,LINE
&&LSEL,R,LOC,X,X_VCEN&&&&&&&&&&
&&LESIZE,ALL,,,YLT_HFX(J,1),,1
&&CMSEL,S,CM_HFX
&&LSEL,R,LOC,Z,Z_VCEN
&&LESIZE,ALL,,,LT_HFZ(I,1),,1
&&CMSEL,S,CM_HFX
&&LSEL,U,LOC,X,X_VCEN
&&LSEL,U,LOC,Z,Z_VCEN
&&LESIZE,ALL,,,10,,1&&
&&CMSEL,S,CM_HFX
&&MSHAPE,0,3D
&&MSHKEY,1
&&VMESH,ALL
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&CMSEL,S,CM_ZSDKW
&X_VCEN=CEN_ZSD(1,1)
&Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2
&VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7
&CM,CM_VMESH,VOLU
&CM,CM_HFX,LINE
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7&&&&&&&&
&LESIZE,ALL,,,2,,1&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7&&&&&&&&&&
&LESIZE,ALL,,,2,,1&&&&&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7
&LSEL,U,LOC,Z,Z_VCEN
&LESIZE,ALL,2&&&&&&&&&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,Z,Z_VCEN&&&&
&LESIZE,ALL,,,LT_HFZ(I,1),,1
&CMSEL,S,CM_HFX
&LSEL,U,LOC,X,X_VCEN
&LSEL,U,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7
&LSEL,U,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7
&LSEL,U,LOC,Z,Z_VCEN
&LESIZE,ALL,,,1,,1&&&&
&CMSEL,S,CM_VMESH
&MSHAPE,0,3D
&VSWEEP,ALL
CMSEL,S,CM_ZMG
VSWEEP,ALL
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&CMSEL,S,CM_YSDKW
&X_VCEN=CEN_YSD(1,1)
&Z_VCEN=(DBKZDZ(I,1)+DBKZDZ(I+1,1))/2
&VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7
&CM,CM_VMESH,VOLU
&CM,CM_HFX,LINE
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7&&&&&&&&
&LESIZE,ALL,,,2,,1&&&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7&&&&&&&&
&LESIZE,ALL,,,2,,1&&&&&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,X,X_VCEN-1E-7,X_VCEN+1E-7
&LSEL,U,LOC,Z,Z_VCEN
&LSEL,U,LOC,Y,CEN_YSD(2,1)-H_CTOB&&&&&&&&&&&&&&&&&&&&
&LESIZE,ALL,2&&&&&&&&
&CMSEL,S,CM_HFX
&LSEL,R,LOC,Z,Z_VCEN&&
&LESIZE,ALL,,,LT_HFZ(I,1),,1
&CMSEL,S,CM_HFX
&LSEL,U,LOC,X,X_VCEN
&LSEL,U,LOC,X,X_VCEN-1E-7,X_VCEN-R_HJ+1E-7
&LSEL,U,LOC,X,X_VCEN+1E-7,X_VCEN+R_HJ-1E-7
&LSEL,U,LOC,Z,Z_VCEN
&LESIZE,ALL,,,1,,1&&&&
&CMSEL,S,CM_VMESH
&MSHAPE,0,3D
&VSWEEP,ALL
&CMSEL,S,CM_YMG
&ESIZE,2&&
&VSEL,R,LOC,Z,Z_VCEN-1E-7,Z_VCEN+1E-7
&VSWEEP,ALL
VSEL,U,MAT,,MAT_SD
VSEL,U,MAT,,MAT_MGKW
VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(9-LTS_YD,1)
CM,CM_VMESH,VOLU
VSEL,S,,,CM_VMESH
VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(LTS_ZD+2,1)&&
LSEL,R,LOC,X,(DBKZDX(LTS_ZD+1,1)+DBKZDX(LTS_ZD+2,1))/2
LESIZE,ALL,,,5&&
VSEL,S,,,CM_VMESH
VSEL,R,LOC,X,DBKZDX(LTS_ZD+2,1),DBKZDX(LTS_ZD+3,1)
LSEL,R,LOC,X,(DBKZDX(LTS_ZD+2,1)+DBKZDX(LTS_ZD+3,1))/2
LESIZE,ALL,,,3&&
VSEL,S,,,CM_VMESH
VSEL,R,LOC,X,DBKZDX(LTS_ZD+3,1),DBKZDX(LTS_ZD+4,1)&&
LSEL,R,LOC,X,(DBKZDX(LTS_ZD+3,1)+DBKZDX(LTS_ZD+4,1))/2
LESIZE,ALL,,,3&&
VSEL,S,,,CM_VMESH
VSEL,R,LOC,X,DBKZDX(LTS_ZD+4,1),DBKZDX(LTS_ZD+5,1)&&
LSEL,R,LOC,X,(DBKZDX(LTS_ZD+4,1)+DBKZDX(LTS_ZD+5,1))/2
LESIZE,ALL,,,5&&
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&!CMSEL,S,CM_VMESH
&VSEL,S,,,CM_VMESH
&VSEL,R,LOC,Z,DBKZDZ(I,1),DBKZDZ(I+1,1)
&VSEL,R,LOC,X,DBKZDX(LTS_ZD+1,1),DBKZDX(LTS_ZD+3,1)
&CM,CM_NEXTV,VOLU
&CMSEL,S,CM_NEXTV
&VSWEEP,ALL
LSEL,S,LOC,X,Y_RDCC(1,2)-Y_RDCC(1,1)/2-1E-7,Y_RDCC(1,2)+Y_RDCC(1,1)/2+1E-7
LSEL,R,LOC,Y,Y_RDCC(2,2)-Y_RDCC(2,1)/2-1E-7,Y_RDCC(2,2)+Y_RDCC(2,1)/2+1E-7
LSEL,R,LOC,Z,Y_RDCC(3,2)-Y_RDCC(3,1)/2-1E-7,Y_RDCC(3,2)+Y_RDCC(3,1)/2+1E-7
LESIZE,ALL,1&&
*DO,I,1,(Y_YJMAX-Y_YJMIN)/PJ_MG+3-1
&VSEL,S,,,CM_VMESH
&VSEL,R,LOC,Z,DBKZDZ(I,1),DBKZDZ(I+1,1)
&CM,CM_NEXTV,VOLU
&CMSEL,S,CM_NEXTV
&VSEL,R,LOC,X,DBKZDX(LTS_ZD+3,1),DBKZDX(LTS_ZD+3+2,1)
&VSWEEP,ALL
VSEL,S,,,CM_ZSDKW&&
VSEL,S,,,CM_ZMG
CM,CM_ZGXW,AREA
MAT,MAT_GXW
VSEL,S,,,CM_YSDKW&&
VSEL,S,,,CM_YMG
CM,CM_YGXW,AREA
MAT,MAT_GXW
NUMCMP,ALL
MPCHG,MAT_TXYS,ALL
SAVE,TNLMESH,DB&&
!------------------------------------------------------- 开始求解
---------------------------------------!
ANTYPE,STATIC,NEW
TIME,1&&&&
NROPT,FULL
AUTOTS,ON&&
DELTIM,0.5,0.01,0.5
LNSRCH,ON&&
PRED,ON&&&&
ESEL,S,TYPE,,2
/FILNAME,CALTUN
RESUME,TNLMESH,DB
ANTYPE,STATIC,NEW
TIME,1&&&&&&&&
NROPT,FULL
AUTOTS,ON&&
DELTIM,0.5,0.01,0.5
LNSRCH,ON&&
PRED,ON&&&&
ESEL,U,TYPE,,2
ISFILE,READ,TUNNEL,IST,,0
ESEL,S,TYPE,,2
EKILL,ALL!*
CNVTOL,F,0.05
!----------------------------------------------------------------------
开挖模拟 ---------------------------------------------------!
!-----------------------------------&&&&第一步开挖，这一步没有支护&&&&
----------------------------------!
VSEL,S,,,CM_ZSDKW
VSEL,A,,,CM_YSDKW
VSEL,R,LOC,Z,DBKZDZ(1,1),DBKZDZ(2,1)
!ESEL,A,TYPE,,2
EKILL,ALL&&
TIME,2&&&&
NROPT,FULL
AUTOTS,ON&&
DELTIM,0.3,0.001,0.5
SAVE,KW-1,DB&&
!---------------------------- 循环做第二步到第五步的开挖，第一步到第四步的支护
----------------------------!
ESEL,S,TYPE,,2
EKILL,ALL&&&&&&
!----------------------------&&第I步开挖,第I-1步支护
---------------------!
*DIM,KWD_TOU,ARRAY,4
*DIM,KWD_WEI,ARRAY,4
*DIM,ZHD_TOU,ARRAY,4
*DIM,ZHD_WEI,ARRAY,4
KWD_TOU(1,1)=3,4,5,6
KWD_WEI(1,1)=2,3,4,5
ZHD_TOU(1,1)=2,3,4,5
ZHD_WEI(1,1)=1,2,3,4
&TOU=KWD_TOU(I,1)
&WEI=KWD_WEI(I,1)
&VSEL,S,,,CM_ZSDKW
&VSEL,A,,,CM_YSDKW
&VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
&ESLV,R&&&&
&EKILL,ALL&&
&TOU=ZHD_TOU(I,1)
&WEI=ZHD_WEI(I,1)
&VSEL,S,,,CM_ZDB&&&&&&&&&&&&&&&&&&&&&&
&VSEL,A,,,CM_YDB
&VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
&EALIVE,ALL&&&&&&
&MPCHG,MAT_DB,ALL
&VSEL,S,,,CM_ZMG&&&&&&&&&&&&&&&&&&&&&&
&VSEL,A,,,CM_YMG
&VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
&MPCHG,MAT_MG,ALL
&ASEL,S,,,CM_ZGXW&&&&&&&&&&&&&&&&&&&&
&ASEL,A,,,CM_YGXW
&ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
&EALIVE,ALL&&&&&&&&
&MPCHG,MAT_GXW,ALL
&ASEL,S,,,CM_ZGXW&&
&ASEL,A,,,CM_YGXW
&ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
&SFE,ALL,,PRES,,ZH_MG
&TIME,I+2&&&&
&NROPT,FULL
&AUTOTS,ON&&
&DELTIM,0.3,0.001,0.5
&LNSRCH,ON&&
&PRED,ON&&
&!OUTRES,ALL,ALL
&SAVE,KW-%I+1%,DB&&
!-----------------------------------&&&&第5步支护，这一步没有开挖&&&&
----------------------------------!
VSEL,S,,,CM_ZDB&&&&&&&&&&&&&&&&&&&&&&
VSEL,A,,,CM_YDB
VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
EALIVE,ALL&&&&&&
MPCHG,MAT_DB,ALL
VSEL,S,,,CM_ZMG&&&&&&&&&&&&&&&&&&&&
VSEL,A,,,CM_YMG
VSEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
MPCHG,MAT_MG,ALL
ASEL,S,,,CM_ZGXW&&&&&&&&&&&&&&&&&&&&
ASEL,A,,,CM_YGXW
ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
EALIVE,ALL&&&&&&
MPCHG,MAT_GXW,ALL
ASEL,S,,,CM_ZGXW&&
ASEL,A,,,CM_YGXW
ASEL,R,LOC,Z,DBKZDZ(WEI,1),DBKZDZ(TOU,1)
SFA,ALL,,PRES,ZH_MG
TIME,7&&&&
NROPT,FULL
AUTOTS,ON&&
DELTIM,0.3,0.01,0.5
LNSRCH,ON&&
!OUTRES,ALL,ALL
SAVE,KW-6,DB&&
/EXIT,NOSAVE
岩土接触分析实例
例子1是为了获得初始应力，计算可以完成。命令流为：
!以便获得初始应力，收敛！！考虑土与岩石之间的摩擦，点面接触
fini !考虑土体的相互作用。自由划分网格，局部细化,不收敛！！！！！！！！！！！
/filname,diji11
/title,地基应力、应变的分布规律
!ET,1,PLANE2
ET,1,PLANE42
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
MP,EX,1,8e9 !岩石性质
MP,PRXY,1,0.29
MP,DENS,1,2600
MP,EX,2,5.0e7 !风化槽土体性质?
MP,PRXY,2,0.25
TBDATA,1,15e3,20,0
MP,DENS,2,2000
et,2,contac48,,,1
mp,mu,3,2.5
r,1,5.0e7,5.0e6,0.01,1.0,0.001,0.001
r,2,5.0e7,5.0e6,0.01,1.0,0.001,0.001
K,1,0,0,0,
K,2,2.5,-12,0,
K,3,6.5,-12,0,
K,4,13.5,-26.5,0,
K,5,20,-30,0,
K,6,24.5,-27,0,
K,7,31.5,-12,0
K,8,35.5,-12,0,
K,9,38,0,0,
k,10,-40,0,0,
k,11,-40,-60,0,
k,12,80,-60,0,
k,13,80,0,0,
K,14,6.5,-12,0,
K,15,13.5,-26.5,0,
K,16,20,-30,0,
K,17,24.5,-27,0,
K,18,31.5,-12,0
lstr,i,i+1
*do,i,10,12
lstr,i,i+1
*do,i,14,17
lstr,i,i+1
LARC,14,18,16,14.57
FLST,2,13,4
*do,i,1,13
FLST,2,5,4
*do,i,14,18
asel,s,area,,1
asel,s,area,,2
lsel,s,line,,3,4,1
CM,_TARGET1,node
lsel,s,line,,14,15,1
CM,_CONTACT1,node
lsel,s,line,,5,6,1
CM,_TARGET2,node
lsel,s,line,,16,17,1
CM,_CONTACT2,node
gcgen,_CONTACT1,_TARGET1
gcgen,_CONTACT2,_TARGET2
antype,static
nropt,full,,off
nsubst,100,1000,10
LSEL,s,line,,10,12,2
CM,_Y,LINE
DL,_Y, ,ux,0,0
LSEL,s,line,,11
CM,_Y,LINE
DL,_Y, ,uy,0,0
iswrite,on
outres,all,all
!读入初始应力，收敛！！考虑土和混凝土与岩石之间的摩擦，点面接触
fini !考虑土体和混凝土之间的相互作用。不收敛！！！！！！！！！！！
/filname,diji
/title,地基应力、应变的分布规律
ET,1,PLANE42
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
MP,EX,1,8e9 !岩石性质
MP,PRXY,1,0.29
MP,DENS,1,2600
MP,EX,2,5.0e7 !风化槽土体性质?
MP,PRXY,2,0.25
TBDATA,1,15e3,20,0
MP,DENS,2,2000
MP,EX,3,2.55e10 !混凝土性质
MP,PRXY,3,0.167
MP,DENS,3,2400
et,2,contac48,,,1
mp,mu,4,2.5
mp,mu,5,4.5
mp,mu,6,3.0
r,1,5.0e7,5.0e6,0.01,1.0,0.001,0.001
r,2,5.0e7,5.0e6,0.01,1.0,0.001,0.001
r,3,8.0e9,8.0e8,0.01,1.0,0.001,0.001
r,4,8.0e9,8.0e8,0.01,1.0,0.001,0.001
r,5,5.0e7,5.0e6,0.01,1.0,0.001,0.001
K,1,0,0,0,
K,2,2.5,-12,0,
K,3,6.5,-12,0,
K,4,13.5,-26.5,0,
K,5,20,-30,0,
K,6,24.5,-27,0,
K,7,31.5,-12,0
K,8,35.5,-12,0,
K,9,38,0,0,
k,10,-40,0,0,
k,11,-40,-60,0,
k,12,80,-60,0,
k,13,80,0,0,
K,14,6.5,-12,0,
K,15,13.5,-26.5,0,
K,16,20,-30,0,
K,17,24.5,-27,0,
K,18,31.5,-12,0
K,19,0,0,0,
K,20,2.5,-12,0,
K,21,6.5,-12,0,
K,22,31.5,-12,0
K,23,35.5,-12,0,
K,24,38,0,0,
K,25,20,0,0,
lstr,i,i+1
*do,i,10,12
lstr,i,i+1
*do,i,14,17
lstr,i,i+1
LARC,14,18,16,14.57
*do,i,19,20
lstr,i,i+1
LARC,21,22,16,14.57
*do,i,22,24
lstr,i,i+1
LSTR,25,19
FLST,2,13,4
*do,i,1,13
FLST,2,5,4
*do,i,14,18
FLST,2,7,4
*do,i,19,25
asel,s,area,,1
asel,s,area,,2
asel,s,area,,1,2
lsel,s,line,,3,4,1
CM,_TARGET1,node
lsel,s,line,,14,15,1
CM,_CONTACT1,node
lsel,s,line,,5,6,1
CM,_TARGET2,node
lsel,s,line,,16,17,1
CM,_CONTACT2,node
gcgen,_CONTACT1,_TARGET1
gcgen,_CONTACT2,_TARGET2
asel,s,area,,3
lsel,s,line,,1,8,7
CM,_TARGET3,node
lsel,s,line,,19,23,4
CM,_CONTACT3,node
lsel,s,line,,2,7,5
CM,_TARGET4,node
lsel,s,line,,20,22,2
CM,_CONTACT4,node
gcgen,_CONTACT3,_TARGET3
gcgen,_CONTACT4,_TARGET4
lsel,s,line,,18
CM,_TARGET5,node
lsel,s,line,,21
CM,_CONTACT5,node
gcgen,_CONTACT5,_TARGET5
antype,static
nropt,full,,off
esel,s,,,1,451
isfile,read,diji11,ist,,1
nsubst,100,1000,10
LSEL,s,line,,10,12,2
CM,_Y,LINE
DL,_Y, ,ux,0,0
LSEL,s,line,,11
CM,_Y,LINE
DL,_Y, ,uy,0,0
LSEL,s,line,,25
CM,_Y1,LINE
CMSEL,,_Y1
SFL,_Y1,PRES,4.5e5,6e5,
LSEL,s,line,,24
CM,_Y1,LINE
CMSEL,,_Y1
SFL,_Y1,PRES,4.5e5,4.5e5
iswrite,on
outres,all,all
钢筋混凝土管的动力响应特性分析实例
在管内壁激振后，采用瞬态模块，模拟半埋设钢筋混凝土管在内部震源作用下的动力特性。
请多提宝贵意见，如果给分就太好了！
/COM,ANSYS MEDIA REL. 7.1 (03-13-2003) REF. VERIF. MANUAL: REL.
/CONFIG,NRES,2000
/VERIFY,ACOUSTIC_WELL_RADIUS
/TITLE,ACOUSTIC, ACOUSTIC_WELL_RADIUS
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!! 按平面问题分析
!! 钢板＋钢筋混凝土＋砂浆/泥浆＋原位岩/土（瞬态）
!! 6种介质成分，存放在介质属性数组中。
!! 计算边界外延一定宽度，以抑制边界干扰*****************
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
JH_HSGD=10&&!计算高度
JH_GBHD=0.04&&!钢板厚度
JH_HNTHD=0.7&& !混凝土厚度
JH_JYHD=0.6&&!浆液厚度，高度方向上砂浆和泥浆各占一半，砂浆在下，泥浆在上
JH_YTTHD=0.6&&!岩土体厚度
JH_RADWELL=4
ZHD=JH_RADWELL+JH_GBHD+JH_HNTHD+JH_JYHD+JH_YTTHD
BJCR1=ZHD+0.5
BJCR2=ZHD+1.0
BJCR3=ZHD+1.5
BJCR4=ZHD+2.0
SCBCJG=250
ENDTIME=0.001 !总计算时间
DST=1E-06 !记录的时间步长
ZYPL=5E4 !震源频率
PI=2*ASIN(1)
TCCS=6 !介质分类数目
!钢&& 钢筋混凝土 泥浆 砂浆 原位岩/土
*DIM,TC,ARRAY,TCCS,5
!该数组为介质参数，第一列为介质密度&&&&
第二列为泊松比，
!&& 第三列为弹模
! 第四列为纵波波速，用于计算时间步长
! 第五列为横波波速，如为液体，横波波速为0
!土层自下而上排列，物理参数暂时按各土层相同考虑
TC(1,1)=00,00
TC(1,2)=0.285,0.32,0.48,0.35,0.35,0.35
!液体泊松比为0.5,计算出现错误，按提示改为0.499
TC(1,3)=210E9,36E9,2E6,6E9,20E9,20E9
!计算各种介质的P、S波速度,以计算时间步长
*DO,I,1,TCCS
&&LAME1=TC(I,3)*TC(I,2)/(1+TC(I,2))/(1-2*TC(I,2))
&&LAME2=TC(I,3)/2/(1+TC(I,2))
&&&&&&TC(I,4)=((LAME1+2*LAME2)/TC(I,1))**.5
&&&&&&TC(I,5)=(LAME2/TC(I,1))**.5
!《土中应力波及其量测》P26
!#########生成各介质区的边界控制点############
JZFQ=4&&&&&&
!介质分区数目
R1=JH_RADWELL
R2=JH_RADWELL+JH_GBHD
R3=JH_RADWELL+JH_GBHD+JH_HNTHD
R4=JH_RADWELL+JH_GBHD+JH_HNTHD+JH_JYHD
R5=JH_RADWELL+JH_GBHD+JH_HNTHD+JH_JYHD+JH_YTTHD
*DIM,JHRB,ARRAY,JZFQ+BJCSM+1,1
JHRB(1,1)=R1,R2,R3,R4,R5,R6,R7,R8,R9!
*DIM,JHRA,ARRAY,JZFQ+BJCSM+1,4&&&&
!定义数组，用于存放
*DIM,AMAT,ARRAY,JZFQ+BJCSM+1,1&&
*DIM,DYCC,ARRAY,JZFQ+BJCSM+1,1&&&&
JHRA(1,1)=R1,R2,R3,R3,R4,R5,R6,R7,R8!
JHRA(1,2)=R2,R3,R4,R4,R5,R6,R7,R8,R9!
JHRA(1,3)=0,0,0,180,0,0,0,0,0
JHRA(1,4)=360,360,180,360,360,360,360,360,360
AMAT(1,1)=1,2,3,4,5,5,5,5,5&&&&!介质编号&&&&&&
DYCC(1,1)=0.2,0.2,0.2,0.2,0.2,0.3,0.45,0.6,1!剖分尺寸
/REPLOT,RESIZE
*DO,I,1,JZFQ+BJCSM+1
&&CYL4,0,0,JHRA(I,1),JHRA(I,3),JHRA(I,2),JHRA(I,4)
NUMMRG,ALL
NUMCMP,ALL !压缩各定义项的编号,使之连续
!############计算时间步长##################
VMAX=TC(1,4) !钢板的纵波速度控制时间步长
DT=DYCC(1)/VMAX/2
*IF,DT,GE,DST,THEN !若时间步长大于设定的记录时间步长，采用记录时间步长记录
ZBC=ENDTIME/DT !时间总步数
!############震源处理##########
T0=1/ZYPL/2
!震源信号持续时间，采用正弦信号激震，也可采用RICKER信号，但要从数组文件中读取，可考虑设定一FLAG控制
XHCXSJ=NINT(1/ZYPL/DT)+1 !信号持续的时间步长
T0I=NINT(XHCXSJ/2)
F=ZYPL&&&&!经MATLAB实际输出后获得的信号实际频率
!######信号记录步长处理#########
*IF,DST/DT,LT,NINT(DST/DT),THEN
!计算信号记录的时间步长间隔，实际记录的时间步长与设定值可能不同
&&RECORDSTEP=NINT(DST/DT)-1
&&RECORDSTEP=NINT(DST/DT)
!#################时间步计算结束#########################
!##################单元剖分################
ET,1,PLANE183
*DO,I,1,JZFQ+BJCSM+1
&&LOCAL, 11, 1, 0,
0&&!定义局部坐标系，柱体
&&CSYS,11&&&&
&&LSEL,S,LOC,X,JHRB(I),&&&&
&&LESIZE,ALL,,R1ANGLE
&&ALLSEL,ALL
AREFINE,1,9,,1!精细剖分
AREFINE,1,4,,1!精细剖分
!#################赋予材料属性#################
*DO,I,1,TCCS
&&MP,DENS,I,TC(I,1)
&&MP,NUXY,I,TC(I,2)
&&MP,EX,I,TC(I,3)
LOCAL, 11, 1, 0,
0&&!定义局部坐标系，柱体
*DO,I,1,JZFQ+BJCSM+1
&&NSEL,S,LOC,X,JHRA(I,1),JHRA(I,2)
&&NSEL,R,LOC,Y,JHRA(I,3),JHRA(I,4)
&&ESLN,S&&&&&&&&
!边界节点单元
&&MAT,AMAT(I)
&&EMODIFY,ALL&&&&&&&&!确认材料参数修改
&&ALLSEL,ALL&&&&&&&&
!重新选择ALL，准备下一步操作
/PNUM,MAT,1
LOCAL, 11, 1, 0,
0&&!定义局部坐标系，柱体
CSYS,11&&&&
LSEL,S,LOC,X,R1
*GET, R1NUM, NODE, 0, COUNT
*GET, R1NMIN, NODE, 0,NUM,
MIN&&!获取最小节点编号
*DIM,R1NODE,ARRAY,R1NUM,1&&!定义内圈节点数组
*DIM,R1NODEXY,ARRAY,R1NUM,2&&!定义内圈节点坐标数组，存储极坐标
R1NODE(1,1)=R1NMIN&&
!赋内圈节点编号的最小值
*GET,RNX, NODE, R1NODE(1,1), LOC,
X&&!获取内圈的节点的坐标，
*GET,RNY, NODE, R1NODE(1,1), LOC,
Y&&!在此为极坐标
R1NODEXY(1,1)=RNX
R1NODEXY(1,2)=RNY
*DO,I,2,R1NUM
&&R1NODE(I,1)=NDNEXT(R1NODE(I-1))&&
!对内圈节点赋值，NDNEXT的输出为选定的节点中编号大于且最接近该节点的节点编号
&&*GET,RNX, NODE, R1NODE(I,1),
LOC, X&&!获取内圈的节点的坐标，
&&*GET,RNY, NODE, R1NODE(I,1),
LOC, Y&&!在此为极坐标
&&R1NODEXY(I,1)=RNX
&&R1NODEXY(I,2)=RNY&&
!###############按极坐标中角度，由小到大排列节点编号和对应坐标
*DO,I,1,R1NUM-1
*DO,J,I+1,R1NUM
&&*IF,R1NODEXY(I,2),GT,R1NODEXY(J,2),THEN
&&&&RY=R1NODEXY(I,2)
&&&&R1NODEXY(I,2)=R1NODEXY(J,2)
&&&&R1NODEXY(J,2)=RY
&&&&RX=R1NODEXY(I,1)
&&&&R1NODEXY(I,1)=R1NODEXY(J,1)
&&&&R1NODEXY(J,1)=RX
&&&&RNODE=R1NODE(I)
&&&&R1NODE(I)=R1NODE(J)
&&&&R1NODE(J)=RNODE
&&*ENDIF&&&&&&&&&&
!####################排序OVER#####################
*CFOPEN,R1NODEXY,DAT!输出井壁内表面节点编号和坐标
*DO,K,1,R1NUM
&&NUM_NODE=R1NODE(K)
&&RX=R1NODEXY(K,1)
&&RY=R1NODEXY(K,2)
&&*VWRITE,NUM_NODE,RX,RY
(F10.0,F13.2,F13.2)
ALLSEL,ALL
!#################求解过程##########################
!##与竖直方向采用直角坐标的不同：
!##1、震源位置变换方式
!##2、记录点位置
!##3、边界条件设定，不作专门的吸收边界，而以空间扩大来处理
*CFOPEN,USUM,DAT
!采用REIKER子波，在初始时段施加
NSUBST,1, , ,1 !1个子步
ANTYPE,TRANS !时程分析
OUTPR,BASIC,SCBCJG, !每隔一定时间步长输出基本值
NSEL,S,LOC,X,R9
D,ALL,UX,0
D,ALL,UY,0
ALLSEL,ALL
*DO,I,1,ZBC
&&*IF,I,GT,1,AND,I,LE,XHCXSJ,THEN&&!在信号持续时间段延续脉冲子波
&&&&*DO,J,1,R1NUM
&&&&&&*GET,RNX,
NODE,R1NODE(J,1), LOC, X
&&&&&&*GET,RNY,
NODE,R1NODE(J,1), LOC, Y
&&&&&&NN=R1NODE(J,1)
&&&&&&PIFT=PI*ZYPL*(I-T0I)*DT
&&&&&&PIFT=PIFT**2
&&&&&&AC=-(1-2*PIFT)*EXP(-PIFT)&&!RECKER子波，何兵寿，P48
&&&&&&F,NN,FX,AC*RNX/R1&&
!施加X方向作用力
&&&&&&F,NN,FY,AC*RNY/R1&&!施加Y方向作用力
&&&&*ENDDO
&&*IF,I,EQ,XHCXSJ+1,THEN&&!在XHCXSJ+1时刻，删除作用在节点上的位移&&
&&&&FDELE,ALL&&!在信号持续时间结束后，删除所有作用力
&&TIME,I*DT
*DO,K,1,R1NUM
&&NUM_NODE=R1NODE(K)
&&RX=R1NODEXY(K,1)
&&RY=R1NODEXY(K,2)
&&*GET, USUMR1, NODE,
NUM_NODE,U,SUM !输出节点编号，极坐标R，角度，USUM
*VWRITE,I,K,NUM_NODE,USUMR1
(F10.0,F10.0,F10.0,E13.3)
BNODEN=NODE(0,R1,0)
BNODES=NODE(0,-R1,0)
FILE,'file','rst','.'
NSOL,2,BNODES,U,X,uxs
NSOL,4,BNODES,U,Y,uys
NSOL,3,BNODEN,U,X,uxn
NSOL,5,BNODEN,U,Y,uyn
PLVAR,uxs,uxn
/IMAGE,SAVE,UXS_N,JPG
PLVAR,uys,uyn, , , , , , ,
/IMAGE,SAVE,UYS_N,JPG
LOCAL, 11, 1, 0,
0&&!定义局部坐标系，柱体
CSYS,11&&&&
NSEL,S,LOC,X,R1,R5
PLNS,U,SUM
ANDATA,0.5, ,2,1,zbc,10,0,1
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。}
我要回帖
更多推荐
- ·合成的帝王红普通朱砂和帝王砂的区别带上对身体有害吗？
- ·49岁在厂做内账3年，代账会计工资一般多少不到5000，月休才一天，太低了。想转行学技术，又没有什么基础，怎么办？
- ·请问这人这是谁叫什么名字？
- ·技校跟职校的区别在哪里和职业学校有什么区别吗
- ·技校是什么学历中专学历还是大专学历？
- ·为什么车轮是圆的，并且车轴都装在圆心到直线距离公式上
- ·求发一份完整的midas连续梁桥建模设计的midas文件
- ·《七层拔毒净斑霜方官 网是多少？？》
- ·一道除法下面算式的被除数相同商是24，余数是26，除数最小是多少？这时被除数是多少
- ·贵州省黔东南州凯里市凯里市嘎醉河水库何时开标
- ·现浇马路用的石子多少钱一立方一立方等于多少吨
- ·极近距离射击击形成的弹痕有什么特点
- ·武安西苑中学学录取出来了吗
- ·光敏传感器色素有哪2种存在形式，怎样调节种子萌发，开花的光周期诱导
- ·在三角形abc中 ab bcab等于ac若角一等于50度，则角b等于多少度
- ·milk thistlebuttermilk怎么吃吃
- ·水光潋滟晴方好，山色空蒙雨亦奇停顿.写的是哪个风景名胜
- ·国家规定第十层地狱电影下载楼的抗震级别应该是几级？
- ·可能补语前面可以加能源动词前面用什么的吗吗?
- ·pred,off撸管会不会影响身高autots,on
- ·查询不到怎么查加拿大签证进度旅游签证进度怎么回事
- ·2米2加1米的椭圆形中间点二点是多少
- ·电子脉冲电子围栏围栏有什么缺点啊，适合用于哪些地方
- ·水蜡洗车能擦铁具吗
- ·为什么候车室称大朗候机楼楼
- ·怎么宗黛坊清之贴足贴真假蜂蜜的鉴别方法鉴别
- ·请问今天农历几月几号？怎样快速减肥 急急急急拜托了各位
- ·如何破解mi语、最新脑筋急转弯弯呀
- ·这个是什么符号花？
- ·请问这是什么鱼，食人鱼生活在哪里里？
- ·高潮是好像是这个，听不太清，you want wayi just want youme是女生唱
- ·电脑看不了电影影
- ·你是哪里是你的拥抱的呀
- ·梦到他人兜很多很多钱从到外地买车 回来落户回来
