ansys建钢管混凝土模型

港控/mmm° 2021-11-23 11:16 214阅读 0赞

FINISH  
/CLEAR  
!(1)工作环境设置  
/FILNAME,CLO,1 ！指定文件名  
/TITLE,CLO ！指定图形标题  
/PREP7  
!(2)定义单元类型  
ET,1,LINK8 !定义钢筋的单元类型  
ET,2,SOLID65 !定义混凝土的单元类型  
KEYOPT,2,1,0  
KEYOPT,2,3,0  
KEYOPT,2,5,1  
KEYOPT,2,6,3  
KEYOPT,2,7,1  
KEYOPT,2,8,0  
ET,3,SOLID45 !定义钢管的单元类型  
KEYOPT,3,1,0  
KEYOPT,3,2,1  
KEYOPT,3,4,0  
KEYOPT,3,5,0  
KEYOPT,3,6,0  
ET,4,SHELL181  
KEYOPT,4,1,0  
KEYOPT,4,3,0  
KEYOPT,4,8,2  
KEYOPT,4,9,0  
KEYOPT,4,10,0  
!(3)定义实常数  
R,1,50.24E6, ， !定义钢筋的截面面积  
R,2, , , , , , ,  
R,3, ,  
R,4,0.06,0.06,0.06,0.06,0,0,  
!(4)定义钢材的材料模型及参数  
MP,EX,1,2.06E11  
MP,PRXY,1,0.25  
MP,DENS,1,7850  
TB,BISO,1,1,2, 1 !双线性各向同性强化模型  
TBDATA,235E6,2.06E10 !理想弹塑性模型  
!(5)定义混凝土材料  
MP,EX,2,32500E6  
MP,PRXY,2,0.173  
MP,DENS,2,2450  
TB,MISO,2,2,14,0 !多线性各向同性强化模型  
TBTEMP,0  
TBPT,0.0000593,1.9286E6 !应力应变数据表  
TBPT,0.0003,5.5782E6  
TBPT,0.0006,10.5336E6  
TBPT,0.0009,14.8662E6  
TBPT,0.0012,18.576E6  
TBPT,0.0016,22.5536E6  
TBPT,0.002,25.4239E6  
TBPT,0.0033,27.739E6  
TBPT,0.005,22.267E6  
TBPT,0.01,10.7514E6  
TBPT,0.015,7.911E6  
TBPT,0.02,6.658E6  
TBPT,0.025,5.915E6  
TBPT,0.03,5.40255E6  
TB,CONC,2,1,9,  
TBDATA,0.3,0.9,1.71E6,-1,0.9 !定义混凝土的破坏参数  
!(6)建立几何模型  
K,1,0.244,0.775,0.244  
K,2,-0.244,0.775,0.244  
K,3,0.244,0.775,-0.244  
K,4,-0.244 ,0.775,-0.244  
K,5,0.244,-0.775,0.244  
K,6,-0.244,-0.775,0.244  
K,7,0.244,-0.775,-0.244  
K,8,-0.244 ,-0.775,-0.244  
K,9,0.244,0.8,0.244  
K,10,-0.244,0.8,0.244  
K,11,0.244,0.8,-0.244  
K,12,-0.244 ,0.8,-0.244  
K,13,0.244,-0.8,0.244  
K,14,-0.244,-0.8,0.244  
K,15,0.244,-0.8,-0.244  
K,16,-0.244 ,-0.8,-0.244  
V,5,7,8,6,1,3,4,2 !建立混凝土的体1  
V,1,3,4,2,9,11,12,10 !建立盖板的体2和3  
V,13,15,16,14,5,7,8,6,  
VGLUE,ALL !粘结在一起  
/VIEW,1,1,1,1  
/ANG,1  
VSEL,1 !混凝土的体1赋予属性  
VATT,2,2,2,0  
VSEL,S,2,3,1 !盖板的体2和3赋予属性  
VATT,1,3,3,0  
ASEL,S,2,5,1 !混凝土的体四个面和盖板的四个面赋予shell单元及其他的属性ASEL,A,7,10,1  
ASEL,A,13,16,1  
AATT,1,4,4,0  
ALLSEL,ALL  
LSEL,S,LINE,5,11,2 !混凝土的竖向线划分为0.05  
LESIZE,ALL,0.05  
LSEL,A,LINE,13,19,2 !两个盖板的竖向线划分为0.025(和其长度相等)  
LESIZE,ALL,0.025  
LSEL,A,LINE,25,28,1  
LESIZE,ALL,0.025  
LSEL,INVE !横向截面的线划分为0.048(可整除)  
LESIZE,ALL,0.0488  
ASEL,ALL  
MSHAPE,0,3D  
MSHKEY,2 !采用映射网格划分  
AMESH,2,5,1 !对具有壳单元属性的面进行划分  
AMESH,7,10,1  
AMESH,13,16,1  
ALLSEL,ALL  
VMESH,1,3,1 !划分三个体  
NUMMRG,NODE  
NUMMRG,ELEM  
NUMMRG,KP  
FINISH  
/SOL  
ANTYPE,0 !静力分析类型  
NLGEOM,1 !打开大变形  
NROPT,AUTO, !自动选择求解方法  
EQSLV,PCG,O, !PCG求解器  
ALLSEL,ALL  
!(9)施加位移约束  
NSEL,S,LOC,Y,0.8  
D,ALL,UX,UZ, !盖板顶加x和z向的水平约束  
ALLSEL,ALL  
NSEL,S,LOC,Y,-0.8  
D,ALL,UX,UY,UZ, !盖板底加x和y和z向的水平约束  
ALLSEL,ALL  
!(10)施加均布荷载(分7个荷载步)  
ALLSEL,ALL  
NSEL,S,LOC,Y,0.8  
SF ,ALL,PRES,-10000E3  
TIME,1  
KBC,0  
NSUBST,10  
ALLSEL,ALL  
LSWRITE,1,  
NSEL,S,LOC,Y,0.8  
SF ,ALL,PRES,-20000E3  
TIME,2  
KBC,0  
NSUBST,10  
ALLSEL,ALL  
LSWRITE,2,  
NSEL,S,LOC,Y,0.8  
SF ,ALL,PRES,-30000E3  
TIME,3  
KBC,0  
NSUBST,10  
ALLSEL,ALL  
LSWRITE,3,  
NSEL,S,LOC,Y,0.8  
SF ,ALL,PRES,-40000E3  
TIME,4  
KBC,0  
NSUBST,10  
ALLSEL,ALL  
LSWRITE,4,  
NSEL,S,LOC,Y,0.8  
SF ,ALL,PRES,-45000E3  
TIME,5  
KBC,0  
NSUBST,10  
ALLSEL,ALL  
LSWRITE,5,  
NSEL,S,LOC,Y,0.8  
SF ,ALL,PRES,-48000E3  
TIME,6  
KBC,0  
NSUBST,10  
ALLSEL,ALL  
LSWRITE,6,  
NSEL,S,LOC,Y,0.8  
SF ,ALL,PRES,-51212.76E3  
TIME,7  
KBC,0  
NSUBST,10  
ALLSEL,ALL  
LSWRITE,7,  
AUTOTS,1 !打开自动时间步长  
CNVTOL,F,5000,0.05,2, !定义收敛准则  
CNVTOL,U,0.03,0,  
NEQIT,50, !平衡迭代次数  
PRED,ON,ON !打开预测器  
OUTRES,ALL,ALL !输出所有子步结果  
LSSOLVE,1,7,1