请求帮助！关于预应力施加问题！
大家好！本人最近刚刚学习ＡＮＳＹＳ不久，参照一个箱梁模型建立了命令流，但是不知道中间有没有那些地方不合适或是存在问题，请教论坛里德各位前辈，谢谢大家。还有就是我要对我的模型施加预应力，如何施加，用哪种方法比较好？
ＣＦＲＰ材料我选用哪种单元类型、材料模型比较合适？
FINISH
/CLEAR
/FILENAME,CFRP box beam under preload
/PREP7
/COM,Structural
/unit,si
!材料定义
ET,1,LINK8 !1号单元，LINK8
MP,EX,1,2.01e11 !底板纵向钢筋
MP,PRXY,1,0.3
MP,DENS,1,7.85e3
R,1,322.062e-6
TB,BKIN,1
TBDATA,,411.9e6,0
MP,EX,2,2.03e11 !顶板和腹板纵向钢筋
MP,PRXY,2,0.3
MP,DENS,2,7.85e3
R,2,116.90e-6
TB,BKIN,2
TBDATA,,420.3e6,0
MP,EX,3,2.14e11 !箍筋、横隔板式转向块钢筋
MP,PRXY,3,0.3
MP,DENS,3,7.85e3
R,3,51.53e-6
TB,BKIN,3
TBDATA,,325.0e6,0
ET,2,SOLID65 !2号单元，SOLID65
MP,EX,4,2.99e10
MP,PRXY,4,0.2
MP,DENS,4,2.45e3
ET,3,SHELL181
R,5,5E-3,5E-3,5E-3,5E-3 !定义碳纤维布的厚度5mm
MP,EX,5,2.3e11
MP,ALPX,5,1.0E-5
TB,CONCR,4,1,9
TBDATA,,0.5,1,4.937e6,-1
TB,MISO,4,,18 !多线性等向强化：定义混凝土的本构关系
tbpt,,200e-6,5.98e6
tbpt,,400e-6,11.96e6
tbpt,,600e-6,14.52e6
tbpt,,800e-6,18.16e6
tbpt,,1000e-6,22.02e6
tbpt,,1200e-6,25.34e6
tbpt,,1400e-6,27.94e6
tbpt,,1600e-6,29.81e6
tbpt,,1800e-6,30.92e6
tbpt,,2000E-6,31.56e6
tbpt,,2200e-6,30.03e6
tbpt,,2400e-6,29.38e6
tbpt,,2600e-6,28.10e6
tbpt,,2800e-6,26.97e6
tbpt,,3000e-6,25.64e6
tbpt,,3200E-6,23.30e6
tbpt,,3400e-6,20.64e6
tbpt,,3600E-6,17.35e6
!创建混凝土几何模型
K,1,310e-3,0
K,2,395e-3,250e-3
K,3,750e-3,270e-3
K,4,750e-3,320e-3
K,5,418.8e-3,320e-3
K,6,358.8e-3,320e-3
K,7,200e-3,320e-3
K,8,0,320e-3
K,9,0,270e-3
K,10,200e-3,270e-3
K,11,335e-3,250e-3
K,12,267e-3,50e-3
K,13,0,50e-3
K,14,0,0
K,15,-310e-3,0
K,16,-395e-3,250e-3
K,17,-750e-3,270e-3
K,18,-750e-3,320e-3
K,19,-418.8e-3,320e-3
K,20,-358.8e-3,320e-3
K,21,-200e-3,320e-3
K,22,-200e-3,270e-3
K,23,-335e-3,250e-3
K,24,-267e-3,50e-3
A,1,12,13,14
A,1,2,11,12
A,2,3,4,5
A,2,5,6,11
A,11,6,7,10
A,10,7,8,9
A,8,21,22,9
A,21,20,23,22
A,20,19,16,23
A,19,18,17,16
A,23,16,15,24
A,15,14,13,24
!切出底板开孔处位置
WPOFF,130E-3, !将工作平面沿其X轴相对偏移130mm
WPROTA,,,-90
ASBW,1
WPCSYS,-1
WPOFF,-130E-3 !将工作平面沿其Y轴相对偏移-130mm
WPROTA,,,-90
ASBW,12
WPCSYS,-1
NUMCMP,ALL !重新编号
!复制梁体截面
AGEN,2,1,14,1,,,400E-3,,,
AGEN,2,1,14,1,,,1400E-3,,,
AGEN,2,1,14,1,,,1525E-3,,,
AGEN,2,1,14,1,,,1725E-3,,,
AGEN,2,1,14,1,,,1875E-3,,,
AGEN,2,1,14,1,,,2925E-3,,,
AGEN,2,1,14,1,,,3075E-3,,,
AGEN,2,1,14,1,,,3275E-3,,,
AGEN,2,1,14,1,,,3400E-3,,,
AGEN,2,1,14,1,,,4400E-3,,,
!创建梁体混凝土几何模型
VEXT,1,14,,,,400E-3 !选中面单元沿Z方向生成长度为400
VEXT,15,28,,,,1000E-3
VEXT,29,42,,,,125E-3
VEXT,56,,,,,200E-3
VEXT,44,54,,,,200E-3
VEXT,57,70,,,,150E-3
VEXT,71,84,,,,1050E-3
VEXT,85,98,,,,150E-3
VEXT,100,110,,,,200E-3
VEXT,112,,,,,200E-3
VEXT,113,126,,,,125E-3
VEXT,127,140,,,,1000E-3
VEXT,141,154,,,,400E-3
!创建横隔板截面
A,13,27,22,9
A,27,24,23,22
A,13,26,10,9
A,26,10,11,12
!复制横隔板截面
AGEN,2,759,762,1,,,1725E-3,,,
AGEN,2,759,762,1,,,2925E-3,,,
AGEN,2,759,762,1,,,4400E-3,,,
!创建横隔板混凝土几何模型
VEXT,759,762,,,,400E-3 !选中面单元沿Z方向生成长度为400
VEXT,763,766,,,,150E-3
VEXT,767,770,,,,150E-3
VEXT,771,774,,,,400E-3
VGLUE,ALL !粘贴体
!划分混凝土模型
VATT,4,4,2 !设置属性
MSHKEY,1 !映射网格划分
ESIZE,50E-3 !单元边长50E-3
VMESH,ALL !划分体
ALLSEL,ALL
!创建底板纵向钢筋
lsel,none
k,,0E-3,25E-3,0E-3
k,,0E-3,25E-3,4800E-3
l,29,30
lgen,3,43,,,140E-3,,
lgen,3,43,,,-140E-3,,
cm,dbzxj0,line !把底板纵向钢筋线定义为元件dbzxj0
lsel,all
cmsel,s,dbzxj0
latt,1,1,1
lesize,all,50E-3,,,
lmesh,all
allsel,all
!定义底板纵筋约束方程
cmsel,s,dbzxj0
nsll,s,1 !选择力筋线所属的所有力筋节点
ceintf,,ux,uy,uz
allsel,all
!创建顶板纵向钢筋
lsel,none
k,,0E-3,295E-3,0E-3
k,,0E-3,295E-3,4800E-3
l,39,40
lgen,4,48,,,127E-3,,
lgen,4,48,,,-127E-3,,
lgen,4,51,,,115E-3,,
lgen,4,54,,,-115E-3,,
cm,dbzxj1,line !把顶板纵向钢筋线定义为元件dbzxj1
lsel,all
cmsel,s,dbzxj1
latt,2,2,1
lesize,all,50E-3,,,
lmesh,all
allsel,all
!定义顶板纵筋约束方程
cmsel,s,dbzxj1
nsll,s,1 !选择力筋线所属的所有力筋节点
ceintf,,ux,uy,uz
allsel,all
!创建腹板纵向钢筋
lsel,none
k,,312.75E-3,92.5E-3,0E-3
k,,312.75E-3,92.5E-3,4800E-3
l,65,66
k,,-312.75E-3,92.5E-3,0E-3
k,,-312.75E-3,92.5E-3,4800E-3
l,67,68
k,,335.5E-3,160E-3,0E-3
k,,335.5E-3,160E-3,4800E-3
l,69,70
k,,-335.5E-3,160E-3,0E-3
k,,-335.5E-3,160E-3,4800E-3
l,71,72
k,,358.25E-3,227.5E-3,0E-3
k,,358.25E-3,227.5E-3,4800E-3
l,73,74
k,,-358.25E-3,227.5E-3,0E-3
k,,-358.25E-3,227.5E-3,4800E-3
l,75,76
cm,fbzxj,line !把腹板纵向钢筋线定义为元件fbzxj
lsel,all
cmsel,s,fbzxj
latt,2,2,1
lesize,all,50E-3,,,
lmesh,all
allsel,all
!定义腹板纵筋约束方程
cmsel,s,fbzxj
nsll,s,1 !选择力筋线所属的所有力筋节点
ceintf,,ux,uy,uz
allsel,all
!创建底板箍筋
lsel,none
k,,280E-3,25E-3,35E-3
k,,-280E-3,25E-3,35E-3
l,77,78
lgen,44,67,,,,,110E-3
!创建顶板箍筋
k,,726E-3,295E-3,35E-3
k,,-726E-3,295E-3,35E-3
l,173,174
lgen,44,111,,,,,110E-3
!创建腹板箍筋
k,,381E-3,295E-3,35E-3
k,,280E-3,25E-3,35E-3
l,269,270
lgen,44,159,,,,,110E-3
k,,-381E-3,295E-3,35E-3
k,,-280E-3,25E-3,35E-3
l,661,662
lgen,44,207,,,,,110E-3
cm,gj,line !把箍筋线定义为元件gj
lsel,all
cmsel,s,gj
latt,3,3,1
lesize,all,50E-3,,,
lmesh,all
allsel,all
!定义箍筋约束方程
cmsel,s,gj
nsll,s,1 !选择力筋线所属的所有力筋节点
ceintf,,ux,uy,uz
allsel,all
!创建横隔板式转向块钢筋
lsel,none
k,,0E-3,50E-3,1800E-3
k,,0E-3,270E-3,1800E-3
l,749,750
lgen,3,251,,,120E-3,,
lgen,3,251,,,-120E-3,,
k,,0E-3,50E-3,3000E-3
k,,0E-3,270E-3,3000E-3
l,759,760
lgen,3,256,,,120E-3,,
lgen,3,256,,,-120E-3,,
cm,zxkj,line !把横隔板式转向块钢筋线定义为元件zxkj
lsel,all
cmsel,s,zxkj
latt,3,3,1
lesize,all,50E-3,,,
lmesh,all
allsel,all
!定义隔板式转向块钢筋约束方程
cmsel,s,zxkj
nsll,s,1 !选择力筋线所属的所有力筋节点
ceintf,,ux,uy,uz
allsel,all
!创建嵌入混凝土部分CFRP模型(考虑到弯折角)
Asel,none
K,,60E-3,193E-3,125E-3
K,,110E-3,193E-3,125E-3
K,,60E-3,208E-3,0E-3
K,,110E-3,208E-3,0E-3
A,769,770,772,771
cm,Carea0,area !把锚固端CFRP嵌入部分定义为元件Carea0
asel,all
cmsel,s,Carea0
aatt,5,5,3
Aesize,all,50E-3,,,
Amesh,all
allsel,all
!定义锚固端CFRP嵌入部分约束方程
cmsel,s,Carea0
nsla,s,1 !选择CFRP面上节点
ceintf,,ux,uy,uz
allsel,all
!创建倾斜部分CFRP模型
Asel,none
K,,60E-3,0E-3,1725E-3
K,,110E-3,0E-3,1725E-3
A,769,770,774,773
cm,Carea1,area !把锚固端至底板弯折处CFRP部分定义为元件Carea1
asel,all
cmsel,s,Carea1
aatt,5,5,3
Aesize,all,50E-3,,,
Amesh,all
allsel,all
!创建两开孔中间部分CFRP模型
Asel,none
K,,60E-3,0E-3,3075E-3
K,,110E-3,0E-3,3075E-3
A,773,774,776,775
cm,Carea2,area !把两开孔部分CFRP部分定义为元件Carea2
asel,all
cmsel,s,Carea2
aatt,5,5,3
Aesize,all,50E-3,,,
Amesh,all
allsel,all
!定义两开孔部分CFRP部分约束方程
cmsel,s,Carea2
nsla,s,1 !选择CFRP面上节点
ceintf,,ux,uy,uz
allsel,all
!创建另一倾斜部分CFRP模型
Asel,none
K,,60E-3,193E-3,4675E-3
K,,110E-3,193E-3,4675E-3
A,775,776,778,777
cm,Carea3,area !把另一锚固端至底板弯折处CFRP部分定义为元件Carea3
asel,all
cmsel,s,Carea3
aatt,5,5,3
Aesize,all,50E-3,,,
Amesh,all
allsel,all
!创建另一嵌入混凝土部分CFRP模型
Asel,none
K,,60E-3,208E-3,4800E-3
K,,110E-3,208E-3,4800E-3
A,777,778,780,779
cm,Carea4,area !把另一锚固端CFRP嵌入部分定义为元件Carea4
asel,all
cmsel,s,Carea4
aatt,5,5,3
Aesize,all,50E-3,,,
Amesh,all
allsel,all
!定义另一锚固端CFRP嵌入部分约束方程
cmsel,s,Carea4
nsla,s,1 !选择CFRP面上节点
ceintf,,ux,uy,uz
allsel,all
!创建第二块嵌入混凝土部分CFRP模型(考虑到弯折角)
Asel,none
K,,-60E-3,193E-3,125E-3
K,,-110E-3,193E-3,125E-3
K,,-60E-3,208E-3,0E-3
K,,-110E-3,208E-3,0E-3
A,781,782,784,783
cm,Carea5,area !把第二块锚固端CFRP嵌入部分定义为元件Carea5
asel,all
cmsel,s,Carea5
aatt,5,5,3
Aesize,all,50E-3,,,
Amesh,all
allsel,all
!定义第二块锚固端CFRP嵌入部分约束方程
cmsel,s,Carea5
nsla,s,1 !选择CFRP面上节点
ceintf,,ux,uy,uz
allsel,all
!创建第二块倾斜部分CFRP模型
Asel,none
K,,-60E-3,0E-3,1725E-3
K,,-110E-3,0E-3,1725E-3
A,783,784,786,785
cm,Carea6,area !把第二块锚固端至底板弯折处CFRP部分定义为元件Carea6
asel,all
cmsel,s,Carea6
aatt,5,5,3
Aesize,all,50E-3,,,
Amesh,all
allsel,all
!创建第二块两开孔中间部分CFRP模型
Asel,none
K,,-60E-3,0E-3,3075E-3
K,,-110E-3,0E-3,3075E-3
A,785,786,788,787
cm,Carea7,area !把第二块两开孔部分CFRP部分定义为元件Carea7
asel,all
cmsel,s,Carea7
aatt,5,5,3
Aesize,all,50E-3,,,
Amesh,all
allsel,all
!定义第二块两开孔部分CFRP部分约束方程
cmsel,s,Carea7
nsla,s,1 !选择CFRP面上节点
ceintf,,ux,uy,uz
allsel,all
!创建第二块另一倾斜部分CFRP模型
Asel,none
K,,-60E-3,193E-3,4675E-3
K,,-110E-3,193E-3,4675E-3
A,787,788,790,789
cm,Carea8,area !把第二块另一锚固端至底板弯折处CFRP部分定义为元件Carea8
asel,all
cmsel,s,Carea8
aatt,5,5,3
Aesize,all,50E-3,,,
Amesh,all
allsel,all
!创建第二块另一嵌入混凝土部分CFRP模型
Asel,none
K,,-60E-3,208E-3,4800E-3
K,,-110E-3,208E-3,4800E-3
A,789,790,792,791
cm,Carea9,area !把第二块另一锚固端CFRP嵌入部分定义为元件Carea9
asel,all
cmsel,s,Carea9
aatt,5,5,3
Aesize,all,50E-3,,,
Amesh,all
allsel,all
!定义第二块另一锚固端CFRP嵌入部分约束方程
cmsel,s,Carea9
nsla,s,1 !选择CFRP面上节点
ceintf,,ux,uy,uz
allsel,all
!对CFRP施加初始应力
/SOLU
ISTRESS,2,2,2,1,1,1,Carea1,Carea2,Carea3,Carea4,Carea5,Carea6,Carea7,Carea8,Carea9,
SOLVE 求助大家了啊
希望大家能帮我看看
页:
[1]