请教单元生死施加约束问题

发布于：2011-08-01 10:54:01 来自：建筑结构/ANSYS [复制转发]

建模没有问题，在求解设置方面我怀疑出现了问题。因为借助杀死单元模拟基坑开挖，前三个求解步中，被杀死的单元没有施加约束，但求解不收敛。
!建立模型
/COM, Structural
!*
/PREP7
!*
ET,1,PLANE82
!*
ET,2,PLANE42
!*
KEYOPT,1,6,0
KEYOPT,1,9,0
KEYOPT,1,10,0
!*
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,2
KEYOPT,2,5,0
KEYOPT,2,6,0
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,3e10
MPDATA,PRXY,1,,0.2
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,2.5e3

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,1.8e7
MPDATA,PRXY,2,,0.33
TB,DP,2,,,
TBMODIF,1,1,20000
TBMODIF,1,2,20
TBMODIF,1,3,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,2e3

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,3,,4.8e7
MPDATA,PRXY,3,,0.27
TB,DP,3,,,
TBMODIF,1,1,45000
TBMODIF,1,2,31
TBMODIF,1,3,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,2.3e3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,4,,5e9
MPDATA,PRXY,4,,0.25
TB,DP,4,,,
TBMODIF,1,1,3e6
TBMODIF,1,2,35
TBMODIF,1,3,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,4,,2.45e3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,5,,1.8e7
MPDATA,PRXY,5,,0.33
TB,DP,5,,,
TBMODIF,1,1,20000
TBMODIF,1,2,20
TBMODIF,1,3,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,5,,2e3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,6,,4.8e7
MPDATA,PRXY,6,,0.27
TB,DP,6,,,
TBMODIF,1,1,45000
TBMODIF,1,2,31
TBMODIF,1,3,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,6,,2.3e3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,7,,5e9
MPDATA,PRXY,7,,0.25
TB,DP,7,,,
TBMODIF,1,1,3e6
TBMODIF,1,2,35
TBMODIF,1,3,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,7,,2.45e3
n,1,,,,,,
ngen,2,1,1,,,,-1.73,,
ngen,5,1,2,,,,-2,,
ngen,2,1,6,,,,-2.22,,
ngen,2,1,7,,,,-2.28,,
ngen,2,1,8,,,,-0.72,,
ngen,18,1,9,,,,-2,,
ngen,51,26,1,26,1,2,,,
ngen,2,26,1301,1326,1,1,,,
ngen,2,20,1333,1352,1,1.865,,,
ngen,2,20,1353,1372,1,2,,,
ngen,2,20,1373,1392,1,2.11,,,
ngen,2,20,1393,1412,1,3,,,
ngen,2,106,1327,1352,1,11.95,,,
ngen,2,26,1433,1458,1,1.05,,,
ngen,7,26,1459,1484,1,2,,,

ngen,2,313,1328,,,0.865,,,
ngen,2,1,1641,,,0.435,,,
ngen,2,1,1642,,,0.43,,,
ngen,2,315,1329,,,1.865,,,
ngen,4,1,1644,,,,-2,,
ngen,2,4,1644,,,0.93,,,
ngen,2,1,1648,,,0.935,,,
ngen,2,5,1645,,,1,,,
ngen,2,1,1650,,,1.43,,,
ngen,2,1,1651,,,1.435,,,
ngen,2,7,1646,,,2,,,
ngen,2,1,1653,,,,-2,,
ngen,2,2,1653,,,1.93,,,
ngen,2,1,1655,,,1.935,,,
ngen,2,3,1654,,,1,,,
ngen,2,1,1657,,,2.43,,,
ngen,2,1,1658,,,2.435,,,
ngen,2,17,1643,,,,1.73,,
ngen,5,1,1660,,,2,,,
ngen,2,4,1661,1664,1,,-1.73,,
ngen,2,3,1666,1668,1,,-2,,
ngen,2,2,1670,1671,1,,-2,,
ngen,2,1,1673,,,,-2,,
type,1
mat,1
*do,j,1301,1308,1
e,j,j+26,j+27,j+1
*enddo
type,2
mat,5
*do,j,1301,1305,1
e,j,j+26,j+27,j+1
*enddo

type,2
mat,6
*do,j,1306,1307,1
e,j,j+26,j+27,j+1
*enddo
type,2
mat,7
e,1308,1334,1335,1309
type,2
mat,2
*do,j,1,1275,26
e,j,j+26,j+27,j+1
*enddo
*do,j,2,1276,26
e,j,j+26,j+27,j+1
*enddo
*do,j,3,1277,26
e,j,j+26,j+27,j+1
*enddo
*do,j,4,1278,26
e,j,j+26,j+27,j+1
*enddo
*do,j,5,1279,26
e,j,j+26,j+27,j+1
*enddo
*do,j,1661,1663,1
e,j,j+1,j+5,j+4
*enddo
*do,j,1666,1667,1
e,j,j+1,j+4,j+3
*enddo
*do,j,1433,1589,26
e,j,j+26,j+27,j+1
*enddo
*do,j,1434,1590,26
e,j,j+26,j+27,j+1
*enddo
*do,j,1435,1591,26
e,j,j+26,j+27,j+1
*enddo
*do,j,1436,1592,26
e,j,j+26,j+27,j+1
*enddo
*do,j,1437,1593,26
e,j,j+26,j+27,j+1
*enddo
type,2
mat,2
e,1327,1641,1328
e,1327,1642,1641
e,1327,1643,1642
e,1328,1641,1644,1329
e,1641,1642,1648,1644
e,1642,1643,1649,1648
e,1329,1644,1645,1330
e,1644,1650,1645
e,1644,1648,1651,1650
e,1648,1649,1652,1651
e,1330,1645,1646,1331
e,1645,1650,1653,1646
e,1650,1651,1655,1653
e,1651,1652,1656,1655
e,1331,1646,1647,1332
e,1646,1653,1654,1647
e,1653,1657,1654
e,1653,1655,1658,1657
e,1655,1656,1659,1658
e,1327,1660,1643
e,1660,1661,1665,1643
e,1643,1665,1649
e,1665,1666,1669,1649
e,1649,1669,1652
e,1669,1670,1672,1652
e,1652,1672,1656
e,1672,1673,1674,1656
e,1673,1436,1437,1674
e,1656,1674,1659
e,1674,1437,1438,1659
e,1664,1433,1434,1668
e,1668,1434,1435,1671
e,1671,1435,1436,1673
e,1670,1671,1673,1672
type,2
mat,3
*do,j,6,1280,26
e,j,j+26,j+27,j+1
*enddo
*do,j,7,1281,26
e,j,j+26,j+27,j+1
*enddo
*do,j,1333,1393,20
e,j,j+20,j+21,j+1
*enddo
*do,j,1413,1595,26
e,j,j+26,j+27,j+1
*enddo
*do,j,1438,1594,26
e,j,j+26,j+27,j+1
*enddo
type,2
mat,3
e,1332,1647,1353,1333
e,1647,1654,1373,1353
e,1654,1657,1393,1373
e,1657,1658,1413,1393
e,1658,1659,1439,1413
e,1659,1438,1439
type,2
mat,4
*do,j,8,1282,26
e,j,j+26,j+27,j+1
*enddo
*do,j,9,1283,26
e,j,j+26,j+27,j+1
*enddo
*do,j,10,1284,26
e,j,j+26,j+27,j+1
*enddo
*do,j,11,1285,26
e,j,j+26,j+27,j+1
*enddo
*do,j,12,1286,26
e,j,j+26,j+27,j+1
*enddo
*do,j,13,1287,26
e,j,j+26,j+27,j+1
*enddo
*do,j,14,1288,26
e,j,j+26,j+27,j+1
*enddo
*do,j,15,1289,26
e,j,j+26,j+27,j+1
*enddo
*do,j,16,1290,26
e,j,j+26,j+27,j+1
*enddo
*do,j,17,1291,26
e,j,j+26,j+27,j+1
*enddo
*do,j,18,1292,26
e,j,j+26,j+27,j+1
*enddo
*do,j,19,1293,26
e,j,j+26,j+27,j+1
*enddo
*do,j,20,1294,26
e,j,j+26,j+27,j+1
*enddo
*do,j,21,1295,26
e,j,j+26,j+27,j+1
*enddo
*do,j,22,1296,26
e,j,j+26,j+27,j+1
*enddo
*do,j,23,1297,26
e,j,j+26,j+27,j+1
*enddo
*do,j,24,1298,26
e,j,j+26,j+27,j+1
*enddo
*do,j,25,1299,26
e,j,j+26,j+27,j+1
*enddo
*do,j,1309,1325,1
e,j,j+26,j+27,j+1
*enddo

*do,j,1334,1351,1
e,j,j+20,j+21,j+1
*enddo
*do,j,1354,1371,1
e,j,j+20,j+21,j+1
*enddo
*do,j,1374,1391,1
e,j,j+20,j+21,j+1
*enddo
*do,j,1394,1411,1
e,j,j+20,j+21,j+1
*enddo
*do,j,1414,1431,1
e,j,j+26,j+27,j+1
*enddo
*do,j,1440,1457,1
e,j,j+26,j+27,j+1
*enddo
*do,j,1466,1483,1
e,j,j+26,j+27,j+1
*enddo
*do,j,1492,1509,1
e,j,j+26,j+27,j+1
*enddo
*do,j,1518,1535,1
e,j,j+26,j+27,j+1
*enddo
*do,j,1544,1561,1
e,j,j+26,j+27,j+1
*enddo
*do,j,1570,1587,1
e,j,j+26,j+27,j+1
*enddo
*do,j,1596,1613,1
e,j,j+26,j+27,j+1
*enddo
!施加约束和荷载
/solu
antype,static
nropt,full !设定牛顿-拉普森选项
nlgeom,on !打开大位移效果
autots,1
nsubst,1000
KBC,0
outres,basic,all
AUTOTS,0

SOLCONTROL,ON
autots,ON !使用自动时间步
lnsrch,ON !打开线性搜索
ncnv,2
neqit,50
pred,on !打开时间步长预测器
time,200

nsel,s,loc,x,0
nsel,a,loc,x,126
d,all,ux
allsel
nsel,s,loc,y,-48.95
d,all,uy
ACEL,0,9.8,0,

NLGEOM,1
NROPT,MODI, ,
LUMPM,0
EQSLV,SPAR, ,0,
PRECISION,0
MSAVE,0
PIVCHECK,1
SSTIF
PSTRES
TOFFST,0,
ALLSEL,ALL
!*
!*
CNVTOL,U, , ,2,,
!*
CNVTOL,F, , ,2, ,
!*
ANTYPE,4
NLGEOM,0
!*
NLGEOM,0
SSTIF,0
NROPT,MODI, ,
EQSLV,SPAR, ,0,
PRECISION,0
MSAVE,0
TOFFST,0,
!*
ANTYPE,0
FINISH
/PREP7
!*
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,2
KEYOPT,1,5,0
KEYOPT,1,6,0
!*
FINISH
/SOL
FLST,2,65,1,ORDE,62
FITEM,2,1
FITEM,2,27
FITEM,2,53
FITEM,2,79
FITEM,2,105
FITEM,2,131
FITEM,2,157
FITEM,2,183
FITEM,2,209
FITEM,2,235
FITEM,2,261
FITEM,2,287
FITEM,2,313
FITEM,2,339
FITEM,2,365
FITEM,2,391
FITEM,2,417
FITEM,2,443
FITEM,2,469
FITEM,2,495
FITEM,2,521
FITEM,2,547
FITEM,2,573
FITEM,2,599
FITEM,2,625
FITEM,2,651
FITEM,2,677
FITEM,2,703
FITEM,2,729
FITEM,2,755
FITEM,2,781
FITEM,2,807
FITEM,2,833
FITEM,2,859
FITEM,2,885
FITEM,2,911
FITEM,2,937
FITEM,2,963
FITEM,2,989
FITEM,2,1015
FITEM,2,1041
FITEM,2,1067
FITEM,2,1093
FITEM,2,1119
FITEM,2,1145
FITEM,2,1171
FITEM,2,1197
FITEM,2,1223
FITEM,2,1249
FITEM,2,1275
FITEM,2,1301
FITEM,2,1327
FITEM,2,1433
FITEM,2,1459
FITEM,2,1485
FITEM,2,1511
FITEM,2,1537
FITEM,2,1563
FITEM,2,1589
FITEM,2,1615
FITEM,2,1660
FITEM,2,-1664
!*
/GO
F,P51X,FY,-58000

!初始应力场模拟
!杀死支护单元
FINISH
/SOL
!*
!*
TIME,1
AUTOTS,-1
DELTIM, , , ,1
KBC,0
!*
TSRES,ERASE
esel,s,mat,,1
ekill,all

ALLSEL,ALL
/STATUS,SOLU
SOLVE

SAVE,'step1','db','D:\ansys\2\'
!激活支护单元，杀死支护桩部位土体单元
allsel
FINISH
/SOL
TIME,2
AUTOTS,-1
DELTIM, , , ,1
KBC,0
!*
TSRES,ERASE
esel,s,mat,,1
ealive,all
esel,s,mat,,5,6,1
ekill,all
esel,s,mat,,7
ekill,all
ALLSEL,ALL
/STATUS,SOLU
SOLVE
esel,s,elem,,1,14
SAVE,'step2','db','D:\ansys\2\'
!7月31日进度
!开挖模拟
allsel
FINISH
/SOL
TIME,3
AUTOTS,-1
DELTIM, , , ,1
KBC,0
!*
TSRES,ERASE
!杀死开挖单元，高度为0，坡度1:1
FLST,2,63,2
FITEM,2,278
FITEM,2,285
FITEM,2,284
FITEM,2,277
FITEM,2,276
FITEM,2,283
FITEM,2,275
FITEM,2,282
FITEM,2,274
FITEM,2,281
FITEM,2,273
FITEM,2,280
FITEM,2,272
FITEM,2,279
FITEM,2,337
FITEM,2,338
FITEM,2,269
FITEM,2,271
FITEM,2,268
FITEM,2,270
FITEM,2,267
FITEM,2,329
FITEM,2,327
FITEM,2,328
FITEM,2,326
FITEM,2,330
FITEM,2,331
FITEM,2,332
FITEM,2,340
FITEM,2,333
FITEM,2,339
FITEM,2,334
FITEM,2,286
FITEM,2,293
FITEM,2,287
FITEM,2,294
FITEM,2,301
FITEM,2,300
FITEM,2,336
FITEM,2,335
FITEM,2,465
FITEM,2,453
FITEM,2,454
FITEM,2,455
FITEM,2,295
FITEM,2,302
FITEM,2,288
FITEM,2,289
FITEM,2,296
FITEM,2,303
FITEM,2,456
FITEM,2,457
FITEM,2,297
FITEM,2,304
FITEM,2,290
FITEM,2,291
FITEM,2,298
FITEM,2,305
FITEM,2,458
FITEM,2,459
FITEM,2,306
FITEM,2,299
FITEM,2,292
EKILL,P51X
FLST,2,77,1,ORDE,23
FITEM,2,1327
FITEM,2,1433
FITEM,2,-1439
FITEM,2,1459
FITEM,2,-1465
FITEM,2,1485
FITEM,2,-1491
FITEM,2,1511
FITEM,2,-1517
FITEM,2,1537
FITEM,2,-1543
FITEM,2,1563
FITEM,2,-1569
FITEM,2,1589
FITEM,2,-1595
FITEM,2,1615
FITEM,2,-1621
FITEM,2,1643
FITEM,2,1649
FITEM,2,1652
FITEM,2,1656
FITEM,2,1659
FITEM,2,-1674
!*
ALLSEL,ALL
/STATUS,SOLU
SOLVE