solid65 钢筋混凝土 非线性 整体式（高手请进）
浏览：67887
各位高手，先谢过了
命令流如下：
/PREP7
!*
! 混凝土弹模, Concrete Elastic Modulus
Ecc=2.9e4
! 混凝土泊松比, Concrete Poisson ratio
EMUc=0.2
! 混凝土抗压抗拉强度
! Concrete Compressive and Tensile Strength
fc=13.1
ft=1.35
! 钢筋弹模
Es=2.0e5
! 钢筋泊松比
EMUs=0.3
! 钢筋屈服强度
fy=310
!***************
!*
hw=6.28 !*输入水位
gh=4.4 !*输入钢筋率变化的高度
!ttd=-5.5,shw=0.78
!***************
! eec-sigp
NPT=26 !!!!!!!!!!NPT=nint(eu-0.0008)/STP+1
*Dim,eec,,NPT,,
*Dim,sigp,,NPT,,
Ca = 2.21
Cd = 0.41
ec = 1370e-6 !!!!!**
eu = 4.2*ec !!!!!**
*VFILL,EEC,RAMP,0.0008,0.0002
*Do,i,1,NPT,
tt=eec(i)/ec
*If,tt,le,1,then
sigp(i)=(Ca*tt+(3-2*Ca)*tt**2+(Ca-2)*tt**3)*fc
*else
sigp(i)=(tt/(Cd*(tt-1)**2+tt))*fc
*endif
*enddo
i=
tt=
!***************
ET,1,SOLID65
KEYOPT,1,1,0
KEYOPT,1,5,0
KEYOPT,1,6,0
KEYOPT,1,7,1
!*
! 设定混凝土单单元实参数
! Concrete element real constants
!* 配筋率
ROWa1=0.0075
ROWa2=0.057
ROWa3=0.057
ROWb1=0.0105
ROWb2=0.057
ROWb3=0.057
!*
R,1,3,ROWa1,90,0,3,ROWa2,
RMORE,0,0,3,ROWa3,90,90,
RMORE, ,
R,2,3,ROWb1,90,0,3,ROWb2,
RMORE,0,0,3,ROWb3,90,90,
RMORE, ,
! 输入混凝土材料, Concrete material parameters
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,Ecc !!!!!!!!!!!!!!!!!!!!!!!水上混凝土
MPDATA,PRXY,1,,EMUc
MPDATA,DENS,1,,2.45E-006
MPDATA,EX,2,,Ecc !!!!!!!!!!!!!!!!!!!!!!!水下混凝土
MPDATA,PRXY,2,,EMUc
MPDATA,DENS,2,,1.45E-006
!
! 输入应力应变曲线, Concrete stress-stain curve
TB,MISO,1,1,NPT,
TBTEMP,0
*Do,i,1,NPT,
TBPT,,eec(i),sigp(i)
*enddo
i=
TB,MISO,2,1,NPT,
TBTEMP,0
*Do,i,1,NPT,
TBPT,,eec(i),sigp(i)
*enddo
i=
!
! 输入破坏曲面, Failure envelop surface
TB,CONC,1,1,9,
TBTEMP,0
TBDATA,,.5,.9,ft,fc,,
TBDATA,,,,1,,,
TB,CONC,2,1,9,
TBTEMP,0
TBDATA,,.5,.9,ft,fc,,
TBDATA,,,,1,,,
!
!钢材基本属性
MPDATA,EX,3,,Es
MPDATA,PRXY,3,,EMUs
!屈服准则
TB,BISO,3,1,2,
TBTEMP,0
TBDATA,,fy,,,,,
!
!建立几何模型
K,1,0,0,-4.51 !!!!!!!!!!!!!!!!!!!!!建立圆筒断面
k,2,0,0,-4.75
k,3,0,8,-4.75
k,4,0,8,-4.51
a,1,2,3,4
K,1000,0,10,0
K,1001,0,0,0
VROTAT,1,,,,,,1000,1001,-180,2 !!!!!!!!!!!!!!!!!!!!!!!旋转成体
!!!!!!!!!!!!!!!!!!将筒体和土体水上水下分开
wpoff,0,hw,0
wprot,0,90,0
VSBW,all,,DELETE
NUMMRG,ALL
NUMCMP,ALL
WPCSYS,-1,0
wpoff,0,GH,0
wprot,0,90,0
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,3
VSBW,P51X,,DELETE
NUMMRG,ALL
NUMCMP,ALL
WPCSYS,-1,0
lsel,s,,,7 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!划分第一层网格
lsel,a,,,13
lesize,all,,,15,,,,,1
lsel,s,,,26
lesize,all,,,6,,,,,1
lsel,s,,,10
lesize,all,,,1,,,,,1
allsel
TYPE, 1
MAT, 1
REAL, 1
ESYS, 0
SECNUM,
FLST,5,2,6,ORDE,2
FITEM,5,1
FITEM,5,-2
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
lsel,s,,,46 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!划分第二层网格
lesize,all,,,6,,,,,1
allsel
TYPE, 1
MAT, 2
REAL, 1
ESYS, 0
SECNUM,
FLST,5,2,6,ORDE,2
FITEM,5,5
FITEM,5,-6
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
lsel,s,,,37 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!划分第三层网格
lesize,all,,,6,,,,,1
allsel
TYPE, 1
MAT, 2
REAL, 2
ESYS, 0
SECNUM,
FLST,5,2,6,ORDE,2
FITEM,5,3
FITEM,5,-4
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!!***********
/REPLOT
VPLOT
EPLOT
ALLSEL,ALL
!************
!施加位移约束
FLST,2,6,5,ORDE,6
FITEM,2,6
FITEM,2,11
FITEM,2,14
FITEM,2,19
FITEM,2,23
FITEM,2,27
DA,P51X,SYMM
DA,1,ALL,0
DA,3,ALL,0
!*******************
ACEL,0,9.8,0,
!施加均布荷载
LOADX=400
LOAD1=200
LOAD2=350
LOAD3=500
xLOAD1=100
xLOAD2=200
xLOAD3=300
pp1=LOAD1/(LOADX*LOADX*8)*10000
PP2=LOAD2/(LOADX*LOADX*8)*10000
PP3=LOAD3/(LOADX*LOADX*8)*10000
xpp1=xLOAD1/(LOADX*LOADX*8)*10000
xPP2=xLOAD2/(LOADX*LOADX*8)*10000
xPP3=xLOAD3/(LOADX*LOADX*8)*10000
VSEL,S, , ,1
VSEL,A, , ,2
ESLV,S
EPLOT
/solu
SFE,ALL,2,PRES, ,pp1, , ,
SFE,ALL,4,PRES, ,xpp1, , ,
ALLSEL,ALL
/REPLOT
pp1=
xpp1=
VSEL,S, , ,5
VSEL,A, , ,6
ESLV,S
EPLOT
/solu
SFE,ALL,5,PRES, ,pp2, , ,
SFE,ALL,3,PRES, ,xpp2, , ,
ALLSEL,ALL
/REPLOT
pp2=
xpp2=
VSEL,S, , ,3
VSEL,A, , ,4
ESLV,S
EPLOT
/solu
SFE,ALL,4,PRES, ,pp3, , ,
SFE,ALL,2,PRES, ,xpp3, , ,
ALLSEL,ALL
/REPLOT
pp3=
xpp3=
/PSF,PRES,NORM,2,0,1
/PBF,DEFA, ,1
/PIC,DEFA, ,1
/PSYMB,CS,0
/PSYMB,NDIR,0
/PSYMB,ESYS,0
/PSYMB,LDIV,0
/PSYMB,LDIR,0
/PSYMB,ADIR,0
/PSYMB,ECON,0
/PSYMB,XNODE,0
/PSYMB,DOT,1
/PSYMB,PCONV,
/PSYMB,LAYR,0
/PSYMB,FBCS,0
!*
/PBC,ALL, ,1
/REP
finish
!*
!--------------定义非线性分析选项--------------------
/solu
cnvtol,f,,0.01,2 !定义收敛条件,使用缺省的VALUE
nsubst,25 !定义子步数
outres,all,all !输出每一子步的结果
autots,1 !打开自动时间步控制
lnsrch,1 !打开线性搜索
ncnv,2 !如果不收敛时结束而不退出
neqit,50 !每一子步中方程的迭代次数限值
pred,on !打开预测器
time,50 !定义时间
solve !求解
*求解的时候要改一下混凝土初始的应力应变值，使其应力应变曲线的初始斜率等于混凝土的弹性模量*
*
在求解过程中出现问题，使求解不能继续：
/COM,ANSYS RELEASE 11.0 UP20070125 20:12:07 10/01/2009
*** WARNING *** CP = 7.516 TIME= 20:13:41
Both solid model and finite element model boundary conditions have been
applied to this model. As solid loads are transferred to the nodes or
elements, they can overwrite directly applied loads.
*** WARNING *** CP = 9.438 TIME= 20:14:34
The MINREF field of the CNVTOL command has been assigned a default
value of 1.0 for convergence on f. Check results carefully.
*** ERROR *** CP = 9.781 TIME= 20:14:34
For material 1 at temperature 0, the initial slope of the stress-strain
curve based on the TB,MISO table ( 13888.0114 ) is less than EX (
29000 ).
*** ERROR *** CP = 9.828 TIME= 20:14:37
For material 2 at temperature 0, the initial slope of the stress-strain
curve based on the TB,MISO table ( 13888.0114 ) is less than EX (
29000 ).
*** ERROR *** CP = 9.859 TIME= 20:14:37
The step data was checked and there were errors found.
Please check output or errors file ( E:\yy.err ) for messages.
*** WARNING *** CP = 11.125 TIME= 20:15:08
The MINREF field of the CNVTOL command has been assigned a default
value of 1.0 for convergence on f. Check results carefully.
*** FATAL *** CP = 70.828 TIME= 20:16:10
The supplied nonlinear material properties produced an incorrect
material tangent matrix. Please check your TB input carefully. A
very small yield stress may cause this failure.
************************************************************************
The above error is non-recoverable by ANSYS
ANSYS run terminated by the indicated error
Current data base saved if possible.
************************************************************************
再次谢谢，很急啊
命令流如下：
/PREP7
!*
! 混凝土弹模, Concrete Elastic Modulus
Ecc=2.9e4
! 混凝土泊松比, Concrete Poisson ratio
EMUc=0.2
! 混凝土抗压抗拉强度
! Concrete Compressive and Tensile Strength
fc=13.1
ft=1.35
! 钢筋弹模
Es=2.0e5
! 钢筋泊松比
EMUs=0.3
! 钢筋屈服强度
fy=310
!***************
!*
hw=6.28 !*输入水位
gh=4.4 !*输入钢筋率变化的高度
!ttd=-5.5,shw=0.78
!***************
! eec-sigp
NPT=26 !!!!!!!!!!NPT=nint(eu-0.0008)/STP+1
*Dim,eec,,NPT,,
*Dim,sigp,,NPT,,
Ca = 2.21
Cd = 0.41
ec = 1370e-6 !!!!!**
eu = 4.2*ec !!!!!**
*VFILL,EEC,RAMP,0.0008,0.0002
*Do,i,1,NPT,
tt=eec(i)/ec
*If,tt,le,1,then
sigp(i)=(Ca*tt+(3-2*Ca)*tt**2+(Ca-2)*tt**3)*fc
*else
sigp(i)=(tt/(Cd*(tt-1)**2+tt))*fc
*endif
*enddo
i=
tt=
!***************
ET,1,SOLID65
KEYOPT,1,1,0
KEYOPT,1,5,0
KEYOPT,1,6,0
KEYOPT,1,7,1
!*
! 设定混凝土单单元实参数
! Concrete element real constants
!* 配筋率
ROWa1=0.0075
ROWa2=0.057
ROWa3=0.057
ROWb1=0.0105
ROWb2=0.057
ROWb3=0.057
!*
R,1,3,ROWa1,90,0,3,ROWa2,
RMORE,0,0,3,ROWa3,90,90,
RMORE, ,
R,2,3,ROWb1,90,0,3,ROWb2,
RMORE,0,0,3,ROWb3,90,90,
RMORE, ,
! 输入混凝土材料, Concrete material parameters
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,Ecc !!!!!!!!!!!!!!!!!!!!!!!水上混凝土
MPDATA,PRXY,1,,EMUc
MPDATA,DENS,1,,2.45E-006
MPDATA,EX,2,,Ecc !!!!!!!!!!!!!!!!!!!!!!!水下混凝土
MPDATA,PRXY,2,,EMUc
MPDATA,DENS,2,,1.45E-006
!
! 输入应力应变曲线, Concrete stress-stain curve
TB,MISO,1,1,NPT,
TBTEMP,0
*Do,i,1,NPT,
TBPT,,eec(i),sigp(i)
*enddo
i=
TB,MISO,2,1,NPT,
TBTEMP,0
*Do,i,1,NPT,
TBPT,,eec(i),sigp(i)
*enddo
i=
!
! 输入破坏曲面, Failure envelop surface
TB,CONC,1,1,9,
TBTEMP,0
TBDATA,,.5,.9,ft,fc,,
TBDATA,,,,1,,,
TB,CONC,2,1,9,
TBTEMP,0
TBDATA,,.5,.9,ft,fc,,
TBDATA,,,,1,,,
!
!钢材基本属性
MPDATA,EX,3,,Es
MPDATA,PRXY,3,,EMUs
!屈服准则
TB,BISO,3,1,2,
TBTEMP,0
TBDATA,,fy,,,,,
!
!建立几何模型
K,1,0,0,-4.51 !!!!!!!!!!!!!!!!!!!!!建立圆筒断面
k,2,0,0,-4.75
k,3,0,8,-4.75
k,4,0,8,-4.51
a,1,2,3,4
K,1000,0,10,0
K,1001,0,0,0
VROTAT,1,,,,,,1000,1001,-180,2 !!!!!!!!!!!!!!!!!!!!!!!旋转成体
!!!!!!!!!!!!!!!!!!将筒体和土体水上水下分开
wpoff,0,hw,0
wprot,0,90,0
VSBW,all,,DELETE
NUMMRG,ALL
NUMCMP,ALL
WPCSYS,-1,0
wpoff,0,GH,0
wprot,0,90,0
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,3
VSBW,P51X,,DELETE
NUMMRG,ALL
NUMCMP,ALL
WPCSYS,-1,0
lsel,s,,,7 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!划分第一层网格
lsel,a,,,13
lesize,all,,,15,,,,,1
lsel,s,,,26
lesize,all,,,6,,,,,1
lsel,s,,,10
lesize,all,,,1,,,,,1
allsel
TYPE, 1
MAT, 1
REAL, 1
ESYS, 0
SECNUM,
FLST,5,2,6,ORDE,2
FITEM,5,1
FITEM,5,-2
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
lsel,s,,,46 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!划分第二层网格
lesize,all,,,6,,,,,1
allsel
TYPE, 1
MAT, 2
REAL, 1
ESYS, 0
SECNUM,
FLST,5,2,6,ORDE,2
FITEM,5,5
FITEM,5,-6
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
lsel,s,,,37 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!划分第三层网格
lesize,all,,,6,,,,,1
allsel
TYPE, 1
MAT, 2
REAL, 2
ESYS, 0
SECNUM,
FLST,5,2,6,ORDE,2
FITEM,5,3
FITEM,5,-4
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VSWEEP,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!!***********
/REPLOT
VPLOT
EPLOT
ALLSEL,ALL
!************
!施加位移约束
FLST,2,6,5,ORDE,6
FITEM,2,6
FITEM,2,11
FITEM,2,14
FITEM,2,19
FITEM,2,23
FITEM,2,27
DA,P51X,SYMM
DA,1,ALL,0
DA,3,ALL,0
!*******************
ACEL,0,9.8,0,
!施加均布荷载
LOADX=400
LOAD1=200
LOAD2=350
LOAD3=500
xLOAD1=100
xLOAD2=200
xLOAD3=300
pp1=LOAD1/(LOADX*LOADX*8)*10000
PP2=LOAD2/(LOADX*LOADX*8)*10000
PP3=LOAD3/(LOADX*LOADX*8)*10000
xpp1=xLOAD1/(LOADX*LOADX*8)*10000
xPP2=xLOAD2/(LOADX*LOADX*8)*10000
xPP3=xLOAD3/(LOADX*LOADX*8)*10000
VSEL,S, , ,1
VSEL,A, , ,2
ESLV,S
EPLOT
/solu
SFE,ALL,2,PRES, ,pp1, , ,
SFE,ALL,4,PRES, ,xpp1, , ,
ALLSEL,ALL
/REPLOT
pp1=
xpp1=
VSEL,S, , ,5
VSEL,A, , ,6
ESLV,S
EPLOT
/solu
SFE,ALL,5,PRES, ,pp2, , ,
SFE,ALL,3,PRES, ,xpp2, , ,
ALLSEL,ALL
/REPLOT
pp2=
xpp2=
VSEL,S, , ,3
VSEL,A, , ,4
ESLV,S
EPLOT
/solu
SFE,ALL,4,PRES, ,pp3, , ,
SFE,ALL,2,PRES, ,xpp3, , ,
ALLSEL,ALL
/REPLOT
pp3=
xpp3=
/PSF,PRES,NORM,2,0,1
/PBF,DEFA, ,1
/PIC,DEFA, ,1
/PSYMB,CS,0
/PSYMB,NDIR,0
/PSYMB,ESYS,0
/PSYMB,LDIV,0
/PSYMB,LDIR,0
/PSYMB,ADIR,0
/PSYMB,ECON,0
/PSYMB,XNODE,0
/PSYMB,DOT,1
/PSYMB,PCONV,
/PSYMB,LAYR,0
/PSYMB,FBCS,0
!*
/PBC,ALL, ,1
/REP
finish
!*
!--------------定义非线性分析选项--------------------
/solu
cnvtol,f,,0.01,2 !定义收敛条件,使用缺省的VALUE
nsubst,25 !定义子步数
outres,all,all !输出每一子步的结果
autots,1 !打开自动时间步控制
lnsrch,1 !打开线性搜索
ncnv,2 !如果不收敛时结束而不退出
neqit,50 !每一子步中方程的迭代次数限值
pred,on !打开预测器
time,50 !定义时间
solve !求解
*求解的时候要改一下混凝土初始的应力应变值，使其应力应变曲线的初始斜率等于混凝土的弹性模量*
*
在求解过程中出现问题，使求解不能继续：
/COM,ANSYS RELEASE 11.0 UP20070125 20:12:07 10/01/2009
*** WARNING *** CP = 7.516 TIME= 20:13:41
Both solid model and finite element model boundary conditions have been
applied to this model. As solid loads are transferred to the nodes or
elements, they can overwrite directly applied loads.
*** WARNING *** CP = 9.438 TIME= 20:14:34
The MINREF field of the CNVTOL command has been assigned a default
value of 1.0 for convergence on f. Check results carefully.
*** ERROR *** CP = 9.781 TIME= 20:14:34
For material 1 at temperature 0, the initial slope of the stress-strain
curve based on the TB,MISO table ( 13888.0114 ) is less than EX (
29000 ).
*** ERROR *** CP = 9.828 TIME= 20:14:37
For material 2 at temperature 0, the initial slope of the stress-strain
curve based on the TB,MISO table ( 13888.0114 ) is less than EX (
29000 ).
*** ERROR *** CP = 9.859 TIME= 20:14:37
The step data was checked and there were errors found.
Please check output or errors file ( E:\yy.err ) for messages.
*** WARNING *** CP = 11.125 TIME= 20:15:08
The MINREF field of the CNVTOL command has been assigned a default
value of 1.0 for convergence on f. Check results carefully.
*** FATAL *** CP = 70.828 TIME= 20:16:10
The supplied nonlinear material properties produced an incorrect
material tangent matrix. Please check your TB input carefully. A
very small yield stress may cause this failure.
************************************************************************
The above error is non-recoverable by ANSYS
ANSYS run terminated by the indicated error
Current data base saved if possible.
************************************************************************
再次谢谢，很急啊
