版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

shangguanlou

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 44.8
紅花: 1
帖子: 32
在線: 8.6小時
蟲號: 1921146
注冊: 2012-08-02
專業(yè): 防災(zāi)工程

[求助] 墩柱擬靜力分析，計算不收斂

剛剛學(xué)習(xí)ANSYS，建了一個墩柱的有限元模型，計算的時候不能收斂，請大家指導(dǎo)一下，謝謝！命令流上傳到附件。
                                       finish
/clear
/PREP7
!==定義程序中需要用到的參數(shù)變量
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 !===墩頂外包厚度
! ==定義材料常數(shù)==
!==混凝土材料模型==
ET,1,65                         !保護(hù)層混凝土
ft=2.01e6
ft1=2.70e6
!  保護(hù)層混凝土
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             !==墩身保護(hù)層             單元總長0.45

!cylind,0.175,0.2,2.25,2.55,0,360             !==墩身保護(hù)層             單元總長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                !  保護(hù)層混凝土 ,墩身及深入基座
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          !==鋼套、核心、保護(hù)層，共13條，單元長度0.15，劃分2分
LPLOT
lesize,all,,,1

LSEL,S,LENGTH,,2.4             !==墩身上半部約束混凝土、墩身保護(hù)層，共9條，單元長度0.05，劃分45份
LPLOT
lesize,all,,,48

LSEL,S,LENGTH,,0.6             !==基座（核心、保護(hù)層、外擴(kuò)），共13條，單元長度0.3，劃分2份
LPLOT
lesize,all,,,2

!==徑向劃分單元===
LSEL,S,LENGTH,,R-C             !==墩核心，共16條，劃分6份
LPLOT
lesize,all,,,6

LSEL,S,LENGTH,,C             !==墩保護(hù)層，共16條，，劃分1份
LPLOT
lesize,all,,,1

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

!==環(huán)向劃分單元=
LSEL,S,RADIUS,,R             !==保護(hù)層，共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

! 線劃分結(jié)束

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

回復(fù)此樓

» 猜你喜歡

323求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有4人回復(fù)
352求調(diào)劑已經(jīng)有3人回復(fù)
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑已經(jīng)有8人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有18人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
280分求調(diào)劑一志愿085802 已經(jīng)有7人回復(fù)
0854電子信息求調(diào)劑已經(jīng)有3人回復(fù)
263求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

1樓 2014-12-25 15:36:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

ansys 熱應(yīng)力分析加大載荷計算結(jié)果不收斂已經(jīng)有5人回復(fù)
解決非線性分析不收斂的技巧已經(jīng)有37人回復(fù)

點擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 shangguanlou 的主題更新