24小时热门版块排行榜

返回列表

shangguanlou

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 44.8
红花: 1
帖子: 32
在线: 8.6小时
虫号: 1921146
注册: 2012-08-02
专业: 防灾工程

[求助] 墩柱拟静力分析，计算不收敛

刚刚学习ANSYS，建了一个墩柱的有限元模型，计算的时候不能收敛，请大家指导一下，谢谢！命令流上传到附件。
                                       finish
/clear
/PREP7
!==定义程序中需要用到的参数变量
R=0.2
L=2.4       !==墩高
C=0.025
EL=0.05
NDIV=6
d1=0.012
d2=0.008
L1=-0.6 !==墩底基座模拟高度
L2=0.3 !==墩底基座外包厚度
L3=0.015 !===墩顶外包厚度
! ==定义材料常数==
!==混凝土材料模型==
ET,1,65                         !保护层混凝土
ft=2.01e6
ft1=2.70e6
!  保护层混凝土
MP,EX,1,2.949057E10 ! CO
MP,NUXY,1,0.2
MP,DENS,1,2500  !混凝土
R,1
TB,CONCR,1
TBDATA,,0.5,0.95,ft,-1
TB,MISO,1,,20
TBPT,,0.0001,2949057.079
TBPT,,0.0003,8103462.716
TBPT,,0.0005,12060000
TBPT,,0.0007,14924417.06
TBPT,,0.0009,16919842.28
TBPT,,0.0011,18264859.46
TBPT,,0.0013,19137309.83
TBPT,,0.0015,19671270.3
TBPT,,0.0017,19964378.78
TBPT,,0.0019,20086672.84
TBPT,,0.0021,20088140.59
TBPT,,0.0023,20004470.14
TBPT,,0.0025,19861215.38
TBPT,,0.0027,19676748.24
TBPT,,0.0029,19464333.23
TBPT,,0.0031,19233583.42
TBPT,,0.0033,18991485.23
TBPT,,0.0035,18743123.54
TBPT,,0.0037,18492198.05
TBPT,,0.0039,18241393.41
TBPLOT,MISO,1
!  核心混凝土
MP,EX,2,2.89058E10 ! CON
MP,NUXY,2,0.2
MP,DENS,2,2500  !混凝土
TB,CONCR,2
TBDATA,,0.5,0.95,ft1,-1
TB,MISO,2,,75
TBPT,,0.0001,2890584.442
TBPT,,0.0003,7815471.108
TBPT,,0.0005,11692290.5
TBPT,,0.0007,14732081.79
TBPT,,0.0009,17126361.68
TBPT,,0.0011,19025748.87
TBPT,,0.0013,20544197.8
TBPT,,0.0015,21767106.21
TBPT,,0.0017,22758586.37
TBPT,,0.0019,23567093.92
TBPT,,0.0021,24229575.49
TBPT,,0.0023,24774474.01
TBPT,,0.0025,25223899.61
TBPT,,0.0027,25595203.32
TBPT,,0.0029,25902125.61
TBPT,,0.0031,26155641.6
TBPT,,0.0033,26364588.76
TBPT,,0.0035,26536137.36
TBPT,,0.0037,26676146.37
TBPT,,0.0039,26789434.99
TBPT,,0.0041,26879991.78
TBPT,,0.0043,26951136.87
TBPT,,0.0045,27005648.92
TBPT,,0.0047,27045865.16
TBPT,,0.0049,27073760.74
TBPT,,0.0051,27091012.18
TBPT,,0.0053,27099048.43
TBPT,,0.0055,27099092.04
TBPT,,0.0057,27092192.82
TBPT,,0.0059,27079255.24
TBPT,,0.0061,27061061.03
TBPT,,0.0063,27038287.86
TBPT,,0.0065,27011524.77
TBPT,,0.0067,26981285.12
TBPT,,0.0069,26948017.41
TBPT,,0.0071,26912114.31
TBPT,,0.0073,26873920.37
TBPT,,0.0075,26833738.52
TBPT,,0.0077,26791835.55
TBPT,,0.0079,26748446.83
TBPT,,0.0081,26703780.33
TBPT,,0.0083,26658020.1
TBPT,,0.0085,26611329.18
TBPT,,0.0087,26563852.2
TBPT,,0.0089,26515717.61
TBPT,,0.0091,26467039.56
TBPT,,0.0093,26417919.59
TBPT,,0.0095,26368448.08
TBPT,,0.0097,26318705.55
TBPT,,0.0099,26268763.72
TBPT,,0.0101,26218686.53
TBPT,,0.0103,26168530.96
TBPT,,0.0105,26118347.82
TBPT,,0.0107,26068182.4
TBPT,,0.0109,26018075.04
TBPT,,0.0111,25968061.67
TBPT,,0.0113,25918174.29
TBPT,,0.0115,25868441.33
TBPT,,0.0117,25818888.04
TBPT,,0.0119,25769536.84
TBPT,,0.0121,25720407.55
TBPT,,0.0123,25671517.71
TBPT,,0.0125,25622882.77
TBPT,,0.0127,25574516.31
TBPT,,0.0129,25526430.21
TBPT,,0.0131,25478634.86
TBPT,,0.0133,25431139.23
TBPT,,0.0135,25383951.11
TBPT,,0.0137,25337077.13
TBPT,,0.0139,25290522.93
TBPT,,0.0141,25244293.22
TBPT,,0.0143,25198391.92
TBPT,,0.0145,25152822.18
TBPT,,0.0147,25107586.48
TBPT,,0.0149,25062686.72
TBPLOT,MISO,2
et,2,solid45          !==墩顶凸起
mp,EX,3,2.1E11
mp,PRXY,3,.27
tb,bkin,3
tbtemp,0
tbdata,,272.5e6,2.1e9

!==钢筋材料模型======

et,3,8
mp,ex,4,2.0e11
mp,prxy,4,0.28
tb,bkin,4
tbtemp,0
tbdata,,362e6,2e9
pi=acos(-1)
r,2,0.25*pi*d1*d1
r,3,0.25*pi*d2*d2

!======建立模型 =====
!=====钢筋模型=====
csys,1
*do,j,1,12
n,j,0.175,j*360/12,-0.6    !建立12个点，半径为0.175，角度间距为 360/12
*enddo
ngen,3,100,1,12,1,0,0,0.3
ngen,L/EL+1,100,201,212,1,0,0,EL

type,3
real,2
mat,4

*do,i,1,4901,100
e,i,i+100
*enddo
*do,i,2,4902,100
e,i,i+100
*enddo

*do,i,3,4903,100
e,i,i+100
*enddo
*do,i,4,4904,100
e,i,i+100
*enddo

*do,i,5,4905,100
e,i,i+100
*enddo

*do,i,6,4906,100
e,i,i+100
*enddo

*do,i,7,4907,100
e,i,i+100
*enddo

*do,i,8,4908,100
e,i,i+100
*enddo

*do,i,9,4909,100
e,i,i+100
*enddo

*do,i,10,4910,100
e,i,i+100
*enddo

*do,i,11,4911,100
e,i,i+100
*enddo

*do,i,12,4912,100
e,i,i+100
*enddo
numcmp,all
/eshape,1 !显示单元实际形状
eplot
!  建立箍筋单元
type,3                      ! 墩底箍筋
real,3
mat,4
*do,i,1,11,1
  e,i,i+1
*enddo
e,12,1
ESEL,S,REAL,,3
Egen,3,12,all,,,,,,,,0.0,0.0,0.3
ESEL,S,ELEM,,625,636,1
EPLOT
Egen,L/EL+1,12,all,,,,,,,,0.0,0.0,EL
nummrg,all
numcmp,all
save
!===建立混凝土模型===
cylind,0,0.175,-0.6,0,0,360          !==墩身在基座的部分单元总长0.6
cylind,0.175,0.2,-0.6,0,0,360          !==基座的部分单元总长0.6
cylind,0.2,0.3,-0.6,0,0,360          !==基座的部分单元总长0.6

cylind,0,0.175,0,2.4,0,360             !==墩身核心             单元总长0.45

!cylind,0,0.175,2.25,2.55,0,360             !==墩身核心             单元总长0.45

cylind,0.175,0.2,0,2.4,0,360             !==墩身保护层             单元总长0.45

!cylind,0.175,0.2,2.25,2.55,0,360             !==墩身保护层             单元总长0.45

cylind,0.0,0.175,2.4,2.45,0,360          !==墩顶凸起部分--SOLID45


WPROTA,,90
VSEL,ALL
VSBW,ALL
WPROTA,,,90
VSBW,ALL

allsel
vptn,all
numcmp,all
ALLSEL
!==体选择
/trlcy,volu,1,all          !  体透明度显示

VSEL,S,VOLU,, 1,8,1                !  核心混凝土 ,墩身及深入基座
VPLOT
vatt,2,1,1

VSEL,S,VOLU,, 13,24,1                !  保护层混凝土 ,墩身及深入基座
VPLOT
vatt,1,1,1

VSEL,S,VOLU,,9,12,1 !  墩顶支座
VPLOT
vatt,3,1,2                            !  SOLID45
ALLSEL
!/trlcy,defa
!==划分单元===
!==竖向单元划分

LSEL,S,LENGTH,,0.05          !==钢套、核心、保护层，共13条，单元长度0.15，划分2分
LPLOT
lesize,all,,,1

LSEL,S,LENGTH,,2.4             !==墩身上半部约束混凝土、墩身保护层，共9条，单元长度0.05，划分45份
LPLOT
lesize,all,,,48

LSEL,S,LENGTH,,0.6             !==基座（核心、保护层、外扩），共13条，单元长度0.3，划分2份
LPLOT
lesize,all,,,2

!==径向划分单元===
LSEL,S,LENGTH,,R-C             !==墩核心，共16条，划分6份
LPLOT
lesize,all,,,6

LSEL,S,LENGTH,,C             !==墩保护层，共16条，，划分1份
LPLOT
lesize,all,,,1

LSEL,S,LENGTH,,0.1             !==基座，共8条，，划分2份
LPLOT
lesize,all,,,2

!==环向划分单元=
LSEL,S,RADIUS,,R             !==保护层，共16条，，划分6份
LPLOT
lesize,all,,,6

LSEL,S,RADIUS,,R-C             !=约束，共16条，，划分6份
LPLOT
lesize,all,,,6

LSEL,S,RADIUS,,0.3             !==基座，共8条，，划分6份
LPLOT
lesize,all,,,6

! 线划分结束

mshkey,1
mshape,0,3d
vsel,s,,,1,24
vmesh,all

ALLSEL,ALL
nummrg,all
numcmp,all
ALLSEL,ALL
SAVE

/solu
nsel,s,loc,z,2.45
*GET,Nnod,NODE,0,COUNT　　　　　  　
*DIM,Xyz,ARRAY,NNOD,1
*GET,Nd,NODE,0,NUM,MIN
*DO,I,1,Nnod,1
Xyz(I,1)=Nd
Nd=NDNEXT(Nd)
*ENDDO
*do,I,1,Nnod
NSEL,S,NODE,,Xyz(I,1)
F,Xyz(I,1),FZ,-520e3/nnod
*enddo
CSYS,1
nsel,s,LOC,X,0.3
d,all,all

OUTRES,ALL,ALL
pred,on
cnvtol,f,,0.05,2,0.2
autots,1
lnsrch,1
ncnv,2
neqit,50

time,1
allsel
nsubst,2
acel,,,9.8
kbc,1
lswrite,1
time,10
nsel,s,,,NODE(,,L)
f,NODE(,,L),fx,90000
!d,NODE(,,L),ux,0.05
KBC,0
allsel,all
NSUBST,100,1000000,10
lswrite,2

lssolve,1,2,1

回复此楼