关于ANSYS非线性屈曲分析,采用弧长法计算不收敛,大神帮我看看命令流
/filname,C3
/title,C3
/PREP7
! FRP 厚度, FRP thickness
Fthick=5
! 混凝土弹模, Concrete Elastic Modulus
Ec=33109.7651
! 混凝土泊松比, Concrete Poisson ratio
EMUc=0.2
! 混凝土抗压抗拉强度
! Concrete Compressive and Tensile Strength
fc=45
ft=2.00
! FRP弹模和泊松比
! Elastic Modulue and Poisson ratio for FRP
Ef=22000
EMUf=0.2
!型钢弹性模量
Es=2.06e5
!型钢泊松比
Esmuf=0.3
!型钢抗压抗拉强度
fy=235 !常温下屈服应力
exx=2.06e5 !常温下杨氏模量
! Height and Radius of Tube
! 圆柱高度和半径
H=1400
Rf=95
! 选择混凝土单元, Concrete element
ET,1,SOLID65
keyopt,1,5,1 !混凝土的线性解输出控制:给出每个积分点的解
keyopt,1,6,3 !混凝土的线非性解输出控制:给出每个积分点的解
keyopt,1,7,1 !考虑应力松弛
! 选择FRP单元
ET,2,solid45
! 输入混凝土材料
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,Ec
MPDATA,PRXY,1,,EMUc
! 输入应力应变曲线, Concrete stress-stain curve
TB,kinh,1,1,
TBTEMP,0
tbpt,,164e-6,5.43
tbpt,,493e-6,15.61
tbpt,,821e-6,24.84
tbpt,,1149e-6,33.15
tbpt,,1478e-6,40.54
tbpt,,1806e-6,46.99
tbpt,,2134e-6,52.5
tbpt,,2462e-6,57.09
tbpt,,2790e-6,60.75
tbpt,,3118e-6,63.49
tbpt,,3446e-6,65.35
tbpt,,3774e-6,66.7
tbpt,,4102e-6,68.06
tbpt,,4430e-6,69.42
tbpt,,4758e-6,70.78
! 输入破坏曲面, Failure envelop surface
TB,CONC,1,1,9,
TBTEMP,0
TBDATA,,0.5,1,ft,-1,,
TBDATA,,,,1,,,
MPTEMP,,,,,,,,
MPTEMP,1,0
! 输入FRP材料参数
MPDATA,EX,2,,Ef
MPDATA,PRXY,2,,EMUf
TB,bkin,2,1,2
TBTEMP,0
TBDATA,,430,0
! 输入型钢材料参数
MPDATA,EX,3,,Es
MPDATA,PRXY,3,,Esmuf
TB,BISO,3,1,2,
TBTEMP,0
TBDATA,,235,2.06E4,,,,
! 建立组合柱模型
k,1,-34,-50,0
k,2,34,-50,0
k,3,34,-42.4,
k,4,2.3,-42.4,0
k,5,2.3,42.4,0
k,6,34,42.4,0
k,7,34,50,0
k,8,-34,50,0
k,9,-34,42.4,0
k,10,-2.3,42.4,0
k,11,-2.3,-42.4,0
k,12,-34,-42.4,0
a,1,2,3,4,5,6,7,8,9,10,11,12
pcirc,Rf,Rf+Fthick,0,360
pcirc,0,Rf,0,360
asba,3,1,sepo,,keep,
k,100,100,100,0
k,101,100,100,H
l,100,101
vdrag,all,,,,,,21
LDELE,21,,,1
numcmp,all
vglue,all
numcmp,all
! 赋予材性
VSEL,S,,,2 !FRP管 2
VATT,2,0,2
VSEL,S,,,3 !混凝土 1
VATT,1,0,1
VSEL,S,,,1 !型钢 3
VATT,3,0,2
allsel,all
lesize,35,,,35$lesize,37,,,35$lesize,47,,,35$lesize,49,,,35$lesize,57,,,35
lesize,59,,,35$lesize,29,,,35$lesize,39,,,35$lesize,31,,,35$lesize,33,,,35
lesize,22,,,35$lesize,43,,,35$lesize,27,,,35$lesize,41,,,35$lesize,51,,,35
lesize,55,,,35$lesize,54,,,35$lesize,46,,,35$lesize,23,,,35$lesize,25,,,35
!划分环向
lesize,48,,,10$lesize,58,,,10$lesize,45,,,10$lesize,60,,,10$lesize,52,,,10
lesize,53,,,10$lesize,50,,,10$lesize,56,,,10$lesize,14,,,10$lesize,19,,,10
lesize,20,,,10$lesize,13,,,10$lesize,16,,,10$lesize,17,,,10$lesize,15,,,10
lesize,18,,,10
!划分工字钢截面
lesize,34,,,6$lesize,21,,,6$lesize,7,,,6$lesize,1,,,6
lesize,40,,,6$lesize,28,,,6$lesize,4,,,6$lesize,10,,,6
lesize,9,,,2.5$lesize,5,,,2.5$lesize,11,,,2.5$lesize,3,,,2.5
lesize,38,,,2.5$lesize,30,,,2.5$lesize,26,,,2.5$lesize,42,,,2.5
lesize,32,,,1$lesize,36,,,1$lesize,44,,,1$lesize,24,,,1
lesize,6,,,1$lesize,8,,,1$lesize,2,,,1$lesize,12,,,1
vsweep,1,1,15,1 !型钢
vsweep,3,25,26,1 !混凝土
vsweep,2,2,24,1 !FRP管
numcmp,all !压缩节点编号
nummrg,all
ALLSEL,ALL
nsel,s,loc,Z,H
Sf,all,pres,-0.0001
nsel,s,loc,Z,0
d,all,ux,0,,,,uy,uz
allsel,all
finish
/SOLU
eqslv,spar
pstres,on
Solve
finish
/SOLU
antype,1
bucopt,lanb,6,0,0
mxpand,6,0,0,1,0.001
Outres,all,all
solve
finish
/PREP7
Upgeom,0.01,1,1,’C3’,’rst’,’’
finish
/solu
antype,0
nlgeom,on
outres,all,all
nsel,s,loc,Z,H
Sf,all,pres,70
allsel,all
nsubst,200
arclen,on
NEQIT,20
solve