finish
/clear
/config,nproc,2 !双核处理
/config,nres,10000000
/filnam,solid segment creep model of ansys
/title,solid segment creep model of ansys
/prep7
Ec=3.0e10 !定义弹性模量N/m^2
Pc=2550 !定义密度kg/m^3
Mu=0.1667 !定义泊松比
GMu=-0.9999 !不考虑剪切变形时泊松比
alpha=1.0E-5 !材料膨胀率
rftemp=100 !参考温度℃
P=1000e3 !外荷载N
BD=1.0 !梁高
BB=1.0 !梁宽
BL=5.0 !梁长
Ld=0.50 !梁纵向划分长度
BDN=6 !梁高度方向划分份数
BBN=6 !梁宽度方向划分份数
!---------------------------------------------------------------------------施工阶段数目及每个施工阶段的持续时间
CSstep=3 !施工步骤数目
*dim,CSperiod,,CSstep,1 !各施工阶段持续时间
CSperiod(1)=500,500,500 !每个=可以给10个数组元素赋值
!---------------------------------------------------------------------------为了不同施工阶段输入"tb,creep,n",在此定义与施工阶段数目相同的材料类型
et,1,solid45
keyopt,1,2,1 !缩减积分
mp,ex,1,Ec
mp,nuxy,1,Mu
mp,dens,1,Pc
mp,alpx,1,alpha !材料膨胀率
mp,reft,1,rftemp !参考温度
mp,ex,2,Ec
mp,nuxy,2,Mu
mp,dens,2,Pc
mp,alpx,2,alpha !材料膨胀率
mp,reft,2,rftemp !参考温度
mp,ex,3,Ec
mp,nuxy,3,Mu
mp,dens,3,Pc
mp,alpx,3,alpha !材料膨胀率
mp,reft,3,rftemp !参考温度
!---------------------------------------------------------------------------
!mp,ex,1,Ec
!mp,ey,1,Ec
!mp,ez,1,Ec
!mp,nuxy,1,Mu !固端梁设泊松比为零,不考虑剪切变形,弹性分析结果与Midas\Dr.Bridge结果接近.收缩徐变分析结果偏小
!mp,nuyz,1,Mu !
!mp,nuxz,1,Mu !固端梁设泊松比正常,考虑剪切变形,弹性分析结果较Midas\Dr.Bridge结果偏大.收缩徐变分析结果偏小,最终结果接近.
!mp,gxy,1,Ec/2/(1+Mu)!Ec/2/(1+GMu)! !
!mp,gyz,1,Ec/2/(1+Mu)!Ec/2/(1+GMu)! !悬臂梁设泊松比正常,考虑剪切变形
!mp,gxz,1,Ec/2/(1+Mu)!Ec/2/(1+GMu)! !
!mp,dens,1,Pc
!mp,alpx,1,alpha !材料膨胀率
!mp,reft,1,rftemp !参考温度
r,1
k,101, 0, BD/2, BB/2
k,102, 0, BD/2,-BB/2
k,103, 0,-BD/2,-BB/2
k,104, 0,-BD/2, BB/2
k,201,BL, BD/2, BB/2
k,202,BL, BD/2,-BB/2
k,203,BL,-BD/2,-BB/2
k,204,BL,-BD/2, BB/2
v,101,102,103,104,201,202,203,204
k,301,2*BL, BD/2, BB/2
k,302,2*BL, BD/2,-BB/2
k,303,2*BL,-BD/2,-BB/2
k,304,2*BL,-BD/2, BB/2
v,201,202,203,204,301,302,303,304
k,401,3*BL, BD/2, BB/2
k,402,3*BL, BD/2,-BB/2
k,403,3*BL,-BD/2,-BB/2
k,404,3*BL,-BD/2, BB/2
v,301,302,303,304,401,402,403,404
allsel
ksel,s,kp,,101
ksel,a,kp,,102
lslk,s,1
lesize,all,,,BBN
ksel,s,kp,,102
ksel,a,kp,,103
lslk,s,1
lesize,all,,,BDN
ksel,s,kp,,101
ksel,a,kp,,201
lslk,s,1
lesize,all,Ld
ksel,s,kp,,201
ksel,a,kp,,202
lslk,s,1
lesize,all,,,BBN
ksel,s,kp,,202
ksel,a,kp,,203
lslk,s,1
lesize,all,,,BDN
ksel,s,kp,,201
ksel,a,kp,,301
lslk,s,1
lesize,all,Ld
ksel,s,kp,,301
ksel,a,kp,,302
lslk,s,1
lesize,all,,,BBN
ksel,s,kp,,302
ksel,a,kp,,303
lslk,s,1
lesize,all,,,BDN
ksel,s,kp,,301
ksel,a,kp,,401
lslk,s,1
lesize,all,Ld
!---------------------------------------------------------------------------不同施工阶段对应不同材料属性,需要分段网格划分
allsel
vsel,s,loc,x,0,BL
vatt,1,1,1
mshkey,1 !映射网格划分
mshape,0,3d !平面四边形网格划分MSHAPE, KEY, Dimension
! KEY
! Key indicating the element shape to be used:
! 0 — Mesh with quadrilateral-shaped elements when Dimension = 2-D mesh with hexahedral-shaped elements when Dimension = 3-D.
! 1 — Mesh with triangle-shaped elements when Dimension = 2-D mesh with tetrahedral-shaped elements when Dimension = 3-D.
! Dimension
! Specifies the dimension of the model to be meshed:
! 2D — 2-D model (area mesh).
! 3D — 3-D model (volume mesh).
vmesh,all
allsel
vsel,s,loc,x,BL,2*BL
vatt,2,1,1
mshkey,1 !映射网格划分
mshape,0,3d !平面四边形网格划分MSHAPE, KEY, Dimension
! KEY
! Key indicating the element shape to be used:
! 0 — Mesh with quadrilateral-shaped elements when Dimension = 2-D mesh with hexahedral-shaped elements when Dimension = 3-D.
! 1 — Mesh with triangle-shaped elements when Dimension = 2-D mesh with tetrahedral-shaped elements when Dimension = 3-D.
! Dimension
! Specifies the dimension of the model to be meshed:
! 2D — 2-D model (area mesh).
! 3D — 3-D model (volume mesh).
vmesh,all
allsel
vsel,s,loc,x,2*BL,3*BL
vatt,3,1,1
mshkey,1 !映射网格划分
mshape,0,3d !平面四边形网格划分MSHAPE, KEY, Dimension
! KEY
! Key indicating the element shape to be used:
! 0 — Mesh with quadrilateral-shaped elements when Dimension = 2-D mesh with hexahedral-shaped elements when Dimension = 3-D.
! 1 — Mesh with triangle-shaped elements when Dimension = 2-D mesh with tetrahedral-shaped elements when Dimension = 3-D.
! Dimension
! Specifies the dimension of the model to be meshed:
! 2D — 2-D model (area mesh).
! 3D — 3-D model (volume mesh).
vmesh,all
allsel
!nsel,s,loc,x,0,BL
!esln,s,1,all
esel,s,mat,,1
cm,seg_1,elem
allsel
!nsel,s,loc,x,BL,2*BL
!esln,s,1,all
esel,s,mat,,2
cm,seg_2,elem
allsel
!nsel,s,loc,x,2*BL,3*BL
!esln,s,1,all
esel,s,mat,,3
cm,seg_3,elem
allsel
!---------------------------------------------------------------------------施加边界条件
allsel
nsel,s,loc,x,0
d,all,all
allsel
!---------------------------------------------------------------------------荷载及求解选项
/solu
antype,0
eqslv,pcg,,
solcontrol,on
autots,on !Use automatic time stepping.自动步进法
kbc,1 !1=stepped=Useful for rate-dependent behavior(e.g., creep, viscoplasticity, etc.)or transient load steps only
!nsubst,2000,2000,5, !AUTOTS,ON carry=on AUTOTS,off carry=off
bfunif,temp,rftemp !Assigns a uniform body force load to all nodes
rate,on !Specifies whether the effect of creep strain rate will be used in the solution of a load step
cutcontrol,crplimit,0.00,1 !1=隐式蠕变分析(隐式蠕变分析方法更强大\更快\更精确,一般推荐使用隐式蠕变分析);0=显式蠕变分析(对于需要很小时间步的情况,显式蠕变分析就非常有用)
tref,rftemp !Defines the reference temperature for the thermal strain calculations
allsel
sfcum,all,repl !面荷载覆盖
fcum,repl !force loads 覆盖
bfcum,all,add !体荷载叠加
!---------------------------------------------------------------------------JTG D62-2004
RH=80 !环境年平均相对湿度%
RH0=100 !相对湿度%
t0=28 !加载时砼龄期d
t1=1 !d
BA=BD*BB !砼梁截面面积m^2
BU=2*(BD+BB) !砼梁截面与大气接触长度m
H=2*BA/BU !构件理论厚度m
H0=100/1000 !m
Fcuk=30 !龄期28d,砼立方体抗压强度标准值MPa
Fcm=0.8*Fcuk+8 !龄期28d,砼平均立方体抗压强度MPa
Fcm0=10 !MPa
betaH=((150*(1+(1.2*RH/RH0)**18)*H/H0+250)<1500)
betaT0=1/(0.1+(t0/t1)**0.2)
betaFcm=5.3/(Fcm/Fcm0)**0.5
phiRH=1+(1-RH/RH0)/0.46/(H/H0)**(1/3)
phi0=phiRH*betaFcm*betaT0
!---------------------------------------------------------------------------施工阶段数目及每个施工阶段的持续时间
CSstep=3 !施工步骤数目
*dim,CSperiod,,CSstep,1 !各施工阶段持续时间
CSperiod(1)=500,500,500 !每个=可以给10个数组元素赋值
!---------------------------------------------------------------------------
t=t0
*do,i,1,CSstep
t=t+CSperiod(i) !计算终值时刻(Midas与Dr.Bridge的施工阶段持续时间均为:以t0为零点基础上的时间,因此软件中输入1000d,此处应输入1000+t0d)
*enddo
!nday=t-t0+1
nday=(t-t0)/500+1
*dim,day,,nday,1 !时间序列
*dim,CC,,nday,1 !蠕变方程系数C1序列
*dim,betaC,,nday,1 !加载后徐变随时间发展系数βc序列
*dim,phi,,nday,1 !加载龄期为t0,计算考虑龄期为t时的混凝土徐变系数φ序列
*dim,betaS,,nday,1 !收缩随时间发展系数βs序列
*dim,epsilonCS,,nday,1 !收缩开始时的龄期为ts,计算考虑龄期为t时的混凝土收缩应变εcs序列
*do,i,1,nday
!day(i)=t0+i-1
day(i)=t0+(i-1)*500
*enddo
betaRH=1.55*(1-(RH/RH0)**3)
betaSC=5.0
epsilonS=(160+10*betaSC*(9-Fcm/Fcm0))*10**(-6)
epsilonCS0=epsilonS*betaRH
ts=3 !d
!---------------------------------------------------------------------------步骤1:激活1号节段,定义1号节段初始徐变及收缩系数
i=1
time,day(i)
nlgeom,off !大变形
nropt,full !Use full Newton-Raphson
estif,1e-6
allsel
ekill,all
esel,s,,,seg_1 !激活第1节段
ealive,all
!--------------------------------------徐变应变与应力成正比
tb,creep,1,1,,1, !1 through 13 =Implicit creep option隐式蠕变分析;0=Explicit creep option显式蠕变分析
tbdata,1,0,1,0,0,0,0
!--------------------------------------徐变应变与应变成正比
!tb,creep,1,1,,1, !1 through 13 =Implicit creep option隐式蠕变分析;0=Explicit creep option显式蠕变分析
!tbdata,1,0,0,1,0,0,0
!--------------------------------------收缩换算温度
!nsle,s
!bf,all,temp,rftemp !收缩换算温度
!--------------------------------------约束未激活单元(约束后位移非累加;不约束正确,但收敛需谨慎)
allsel
esel,s,live
nsle,s
nsel,inve
d,all,all
allsel
nsel,s,loc,x,BL !施加荷载
*get,nd,node,0,count
f,all,fy,-P/nd
allsel
nsel,s,loc,x,0
d,all,all
allsel
solve
!---------------------------------------------------------------------------步骤2:重新选择1号节段,定义1号节段徐变及收缩系数
i=2
time,day(i)
allsel
esel,s,,,seg_1 !重新选择第1节段
betaC(i)=((day(i)-t0)/t1/(betaH+(day(i)-t0)/t1))**0.3
betaC(i-1)=((day(i-1)-t0)/t1/(betaH+(day(i-1)-t0)/t1))**0.3
phi(i)=phi0*betaC(i)
phi(i-1)=phi0*betaC(i-1)
CC(i)=(phi(i)-phi(i-1))/Ec/(day(i)-day(i-1))
!--------------------------------------徐变应变与应力成正比
tb,creep,1,1,,1, !1 through 13 =Implicit creep option隐式蠕变分析;0=Explicit creep option显式蠕变分析
tbdata,1,CC(i),1,0,0,0,0
!--------------------------------------徐变应变与应变成正比
!tb,creep,1,1,,1, !1 through 13 =Implicit creep option隐式蠕变分析;0=Explicit creep option显式蠕变分析
!tbdata,1,CC(i),0,1,0,0,0
!--------------------------------------收缩换算温度
!betaS(i-1)=((day(i-1)-ts)/t1/(350*(H/H0)**2+(day(i-1)-ts)/t1))**0.5
!betaS(i)=((day(i)-ts)/t1/(350*(H/H0)**2+(day(i)-ts)/t1))**0.5
!epsilonCS(i)=epsilonCS0*(betaS(i)-betaS(i-1))
!nsle,s
!bf,all,temp,-epsilonCS(i)/alpha !收缩换算温度
!--------------------------------------约束未激活单元(约束后位移非累加;不约束正确,但收敛需谨慎)
allsel
esel,s,live
nsle,s
nsel,inve
d,all,all
allsel
nsel,s,loc,x,BL !施加荷载
*get,nd,node,0,count
f,all,fy,-P/nd
allsel
nsel,s,loc,x,0
d,all,all
allsel
solve
!---------------------------------------------------------------------------步骤3:根据1号节段施工周期,修改1号节段徐变及收缩系数,并施加外荷载;激活2号节段,定义2号节段初始徐变及收缩系数
i=3
time,day(i)
esel,s,,,seg_1 !激活第1节段
betaC(i)=((day(i)-t0)/t1/(betaH+(day(i)-t0)/t1))**0.3
betaC(i-1)=((day(i-1)-t0)/t1/(betaH+(day(i-1)-t0)/t1))**0.3
phi(i)=phi0*betaC(i)
phi(i-1)=phi0*betaC(i-1)
CC(i)=(phi(i)-phi(i-1))/Ec/(day(i)-day(i-1))
!--------------------------------------徐变应变与应力成正比
tb,creep,1,1,,1, !1 through 13 =Implicit creep option隐式蠕变分析;0=Explicit creep option显式蠕变分析
tbdata,1,CC(i),1,0,0,0,0
!--------------------------------------徐变应变与应变成正比
!tb,creep,1,1,,1, !1 through 13 =Implicit creep option隐式蠕变分析;0=Explicit creep option显式蠕变分析
!tbdata,1,CC(i),0,1,0,0,0
!--------------------------------------收缩换算温度
!betaS(i-1)=((day(i-1)-ts)/t1/(350*(H/H0)**2+(day(i-1)-ts)/t1))**0.5
!betaS(i)=((day(i)-ts)/t1/(350*(H/H0)**2+(day(i)-ts)/t1))**0.5
!epsilonCS(i)=epsilonCS0*(betaS(i)-betaS(i-1))
!nsle,s
!bf,all,temp,-epsilonCS(i)/alpha !收缩换算温度
i=2
allsel
esel,s,,,seg_2 !激活第2节段
ealive,all
betaC(i)=((day(i)-t0)/t1/(betaH+(day(i)-t0)/t1))**0.3
betaC(i-1)=((day(i-1)-t0)/t1/(betaH+(day(i-1)-t0)/t1))**0.3
phi(i)=phi0*betaC(i)
phi(i-1)=phi0*betaC(i-1)
CC(i)=(phi(i)-phi(i-1))/Ec/(day(i)-day(i-1))
!--------------------------------------徐变应变与应力成正比
tb,creep,2,1,,1, !1 through 13 =Implicit creep option隐式蠕变分析;0=Explicit creep option显式蠕变分析
tbdata,1,CC(i),1,0,0,0,0
!--------------------------------------徐变应变与应变成正比
!tb,creep,2,1,,1, !1 through 13 =Implicit creep option隐式蠕变分析;0=Explicit creep option显式蠕变分析
!tbdata,1,CC(i),0,1,0,0,0
!--------------------------------------收缩换算温度
!betaS(i-1)=((day(i-1)-ts)/t1/(350*(H/H0)**2+(day(i-1)-ts)/t1))**0.5
!betaS(i)=((day(i)-ts)/t1/(350*(H/H0)**2+(day(i)-ts)/t1))**0.5
!epsilonCS(i)=epsilonCS0*(betaS(i)-betaS(i-1))
!nsle,s
!bf,all,temp,-epsilonCS(i)/alpha !收缩换算温度
!--------------------------------------约束未激活单元(约束后位移非累加;不约束正确,但收敛需谨慎)
allsel
ddele,all,all
esel,s,live
nsle,s
nsel,inve
d,all,all
allsel
nsel,s,loc,x,BL !施加荷载
*get,nd,node,0,count
f,all,fy,-P/nd
allsel
nsel,s,loc,x,2*BL !施加荷载
*get,nd,node,0,count
f,all,fy,-P/nd
allsel
nsel,s,loc,x,0
d,all,all
allsel
solve |
IP卡
狗仔卡
发表于 2011-3-16 20:00:34
提升卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
