24小时热门版块排行榜

返回列表

nayunayunayu

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 672.3
红花: 1
帖子: 71
在线: 28.9小时
虫号: 2959745
注册: 2014-02-09
性别: GG
专业: 结构工程

[求助] 模拟与试验数据差距很大是怎么回事？？已有3人参与

我的命令流如下：
FINISH
/CLEAR
/filname,nayu
/title,DSTC-C40-3
/PREP7
! FRP 厚度, FRP thickness
Fthick=4
! steel 厚度, FRP thickness
Sthick=5
! 混凝土弹模, Concrete Elastic Modulus
Ec=3.0e4
! 混凝土泊松比, Concrete Poisson ratio
EMUc=0.2
! 混凝土抗压抗拉强度
! Concrete Compressive and Tensile Strength
fc=14.3
ft=1.43
! FRP弹模和泊松比
! Elastic Modulue and Poisson ratio for FRP
Ef=2.2e4
EMUf=0.3
!钢管弹性模量
Es=2.106e5
!钢管泊松比
Esmuf=0.25
!钢管抗压抗拉强度
sc=36.3e4
st=36.3e4
fy=353.2          !常温下屈服应力
exx=2e5          !常温下杨氏模量
! Height and Radius of Tube
! 圆柱高度和半径
H=750
Rf=121
Rs=72.5
! 选择混凝土单元, Concrete element
ET,1,SOLID65
keyopt,1,5,1
keyopt,1,6,3
keyopt,1,7,1
! 选择FRP单元
ET,2,solid45
! 输入混凝土材料参数
MPTEMP,,,,,,,,
MPTEMP,1,20
MPDATA,EX,1,,Ec
MPDATA,PRXY,1,,EMUc
! 输入应力应变曲线, Concrete stress-stain curve
TB,MISO,1,1,
TBTEMP,20
tbpt,,164e-6,4.92
tbpt,,328e-6,9.75
tbpt,,493e-6,14.36
tbpt,,657e-6,18.57
tbpt,,821e-6,22.23
tbpt,,985e-6,25.21
tbpt,,1149e-6,27.45
tbpt,,1314e-6,28.96
tbpt,,1478e-6,29.79
tbpt,,1642e-6,30.05
tbpt,,1970e-6,28.75
tbpt,,2463e-6,24.51
! 输入破坏曲面, Failure envelop surface
TB,CONC,1,1,9,
TBTEMP,20
TBDATA,,0.4,0.9,ft,-1,,
TBDATA,,,,1,,,
MPTEMP,,,,,,,,
MPTEMP,1,20
! 输入FRP材料参数
MPDATA,EX,2,,Ef
MPDATA,PRXY,2,,EMUf
TB,bkin,2,1,2
TBTEMP,20
TBDATA,,430,0
! 输入钢管材料参数
MPDATA,EX,3,,Es
MPDATA,PRXY,3,,Esmuf
TB,BISO,3,1,2, 1                                           !双线性各向同性强化模型
TBDATA,,235,2.106E4
！建立组合柱模型
pcirc,Rs-Sthick,Rs,0,360 !钢管
pcirc,Rf,Rf+Fthick,0,360 !frp管
pcirc,Rs,Rf,0,360       !混凝土
k,100,0,0,0
k,101,0,0,750
l,100,101
vdrag,all,,,,,,25
LDELE,25,,,1
numcmp,all
vglue,all
numcmp,all
!赋予材料属性
VSEL,S,,,1                !FRP管 2
VATT,2,0,2
VSEL,S,,,3                !混凝土 1
VATT,1,0,1
VSEL,S,,,2                !型钢    3
VATT,3,0,2
allsel,all
！划分网格
vsweep,1,1,11,1                !型钢
vsweep,2,2,20,1                !FRP管
vsweep,3,21,22,1             !混凝土
numcmp,all                      !压缩节点编号
nummrg,all
ALLSEL,ALL
!*************************
csys,1
et,4,184
keyopt,4,1,1
nsel,s,loc,z,H
nsel,r,loc,x,Rs-Sthick,Rf
n,350000,0,0,H+10
type,4
*get,nnum,node,0,count
*get,ND,node,0,num,min
*do,i,2,nnum
!!节点个数是nnum，只需要生成nnum个mpc单元
E,350000,ND
ND=NDNEXT(ND)
*enddo
allsel
!**************************
!施加荷载
!**************************
F,350000,FZ,-2702920
!**************************
!施加自由度约束
!**************************
nsel,s,loc,Z,0
d,all,ux,0,,,,uy,uz
allsel,all
FINISH
!***************************
!进入求解器
!***************************
/SOLU
!求解选项设置
ANTYPE,0          !选择分析类型
NLGEOM,1          !计入大变形效应
NSUBST,30,100,20,on !子步数
AUTOTS,ON          !自动时间步
NEQIT,50          !平衡迭代次数
CNVTOL,U,,0.015 !设置收敛准则
ncnv,2       !如果不收敛时结束而不退出
OUTRES,ALL,ALL    !输出所有子步的结果
LNSRCH,1                !用Newton-Raphson方法激活线性搜索
cutcontrol,nioterpred,0          !在非线性解析中控制步长的削减
solve
模拟完成！！！！所得结果如下
A1（实验数据） A2（模拟数据）
荷载位移荷载          位移
200 2.9 200    0.096
400 3.91 400          0.193
600 4.76    600    0.293
800 5.56          800       0.398
1000 6.39 1000 0.511
1100 6.81 1100 0.572
1200 7.37 1200 0.639
1360 8.57 1360 0.867
1460 9.75 1460 0.94
1560 11.21    1560 1.091
1660 13.42    1660 1.161
1920  25.74 1920 2.156
试验肯定是没有问题的，但是有限元模拟为什么跟试验差的这么多，问题出在哪里呢？我也检查不出来呀，求教各位大神！！！

A1是师兄试验所得荷载-位移曲线，A2是我模拟的结果，差的也太多了····

回复此楼