| 查看: 798 | 回复: 0 | ||
[求助]
墩柱拟静力分析，计算不收敛
|
|
刚刚学习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 |
回复此楼» 猜你喜欢
298求调剂
已经有5人回复
-
272求调剂
已经有4人回复
-
285求调剂
已经有6人回复
-
材料调剂
已经有4人回复
-
面上模板改不了页边距吧？
已经有6人回复
-
307求调剂
已经有6人回复
-
304求调剂
已经有5人回复
-
317一志愿华南理工电气工程求调剂
已经有8人回复
-
272求调剂
已经有3人回复
-
化工专硕348，一志愿985求调剂
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
点击这里搜索更多相关资源
