| 查看: 869 | 回复: 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 |
回复此楼» 猜你喜欢
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
点击这里搜索更多相关资源
