继续昨天的问题,一个cfrp 加固梁,cfrp与梁是共用节点,模型能以收敛,能算出来结果,但是算出来结果不争取,CFRP应力很大,而主梁上应力很小。
/filename,gujinzhuliang4
/prep7
et,1,solid45
et,2,shell41
keyopt,2,1,2
r,1
r,2,34e-5
mp,ex,1,9819e6
mp,prxy,1,0.266
mp,dens,1,672.6
mp,ex,2,2.41e11
mp,prxy,2,0
mp,dens,2,0.3
tb,bkin,2,1,2,1
tbdata,,3558e6
k,1,0,0,0
k,2,0,0,0.13
k,3,0,0.1,0.13
k,4,0,0.1,0
k,5,0.14/3,0,0
k,6,0.14/3,0,0.13
k,7,0.14/3,0.1,0.13
k,8,0.14/3,0.1,0
l,1,5
l,5,6
l,2,6
l,2,1
l,5,8
l,4,8
l,4,1
l,8,7
l,7,3
l,3,4
l,7,6
l,2,3
al,1,2,3,4
al,1,5,6,7
al,6,8,9,10
al,3,11,9,12
FLST,3,4,5,ORDE,2
FITEM,3,1
FITEM,3,-4
AGEN, ,P51X, , ,0.14/3,0,0, , ,1
*do,ii,1,3
agen,2,1,4,1,ii*(0.56/3)
*SET,ii,ii+1
*enddo
*do,ii,1,1
agen,2,13,16,1,ii*(0.14)
*SET,ii,ii+1
*enddo
*do,ii,1,3
agen,2,17,20,1,ii*(0.56/3)
*SET,ii,ii+1
*enddo
type,1
real,1
mat,1
block,0,1.4,0,0.1,0,0.13
lsel,s,,,98
lsel,a,,,103
lesize,all,,,30
lsel,s,,,100
lsel,a,,,101
lesize,all,,,30
lsel,s,,,105
lsel,a,,,108
lesize,all,,,5
lsel,s,,,106
lsel,a,,,107
lesize,all,,,5
lsel,s,,,97
lsel,a,,,104
lesize,all,,,4
lsel,s,,,99
lsel,a,,,102
lesize,all,,,4
allsel
vmesh,1
/replot
asel,s,area,,1,32,1
aplot
cm,crp,area
aatt,2,2,2,0
lsel,s,line,,12,96,12
lsel,a,line,,7,91,12
lesize,all,,,4
lsel,s,line,,2,86,12
lsel,a,line,,8,92,12
lesize,all,,,5
lsel,s,line,,3,87,12
lsel,a,line,,1,85,12
lsel,a,line,,6,90,12
lsel,a,line,,9,93,12
lesize,all,,,1
mshape,1,2D
mshkey,1
amesh,all
allsel,all
eplot
NUMMRG,NODE,0.001,0, ,LOW
numcmp,all
nplot
/solu
antype,static
nsel,s,loc,x,0
nsel,r,loc,z,0
d,all,ux
d,all,uy
d,all,uz
nsel,s,loc,x,1.4
nsel,r,loc,z,0
d,all,uy
d,all,uz
allsel
f,252,fz,-35000
f,282,fz,-35000
allsel
solve
FINISH
这是命令流,结果应力图是这样的![毕业论文急需要ansys模型,算出来结果不对啊,]()
然后用了另外一种接触对的方法来模拟CFRP和主梁之间的接触。
命令流如下
/filename,gujinzhuliang4
/prep7
et,1,solid45
et,2,shell41
et,3,targe170
keyopt,3,5,1
ET, 4, CONTA175
KEYOPT, 4, 2, 2
KEYOPT, 4, 12, 5
r,1
r,2,34e-5
r,3
r,4
r,5
r,6
mp,ex,1,9819e6
mp,prxy,1,0.266
mp,dens,1,672.6
mp,ex,2,2.41e11
mp,prxy,2,0
mp,dens,2,0.3
tb,bkin,2,1,2,1
tbdata,,3558e6
k,1,0,0,0
k,2,0,0,0.13
k,3,0,0.1,0.13
k,4,0,0.1,0
k,5,0.14/3,0,0
k,6,0.14/3,0,0.13
k,7,0.14/3,0.1,0.13
k,8,0.14/3,0.1,0
l,1,5
l,5,6
l,2,6
l,2,1
l,5,8
l,4,8
l,4,1
l,8,7
l,7,3
l,3,4
l,7,6
l,2,3
al,1,2,3,4
al,1,5,6,7
al,6,8,9,10
al,3,11,9,12
FLST,3,4,5,ORDE,2
FITEM,3,1
FITEM,3,-4
AGEN, ,P51X, , ,0.14/3,0,0, , ,1
*do,ii,1,3
agen,2,1,4,1,ii*(0.56/3)
*SET,ii,ii+1
*enddo
*do,ii,1,1
agen,2,13,16,1,ii*(0.14)
*SET,ii,ii+1
*enddo
*do,ii,1,3
agen,2,17,20,1,ii*(0.56/3)
*SET,ii,ii+1
*enddo
type,1
real,1
mat,1
block,0,1.4,0,0.1,0,0.13
lsel,s,,,98
lsel,a,,,103
lesize,all,,,30
lsel,s,,,100
lsel,a,,,101
lesize,all,,,30
lsel,s,,,105
lsel,a,,,108
lesize,all,,,5
lsel,s,,,106
lsel,a,,,107
lesize,all,,,5
lsel,s,,,97
lsel,a,,,104
lesize,all,,,4
lsel,s,,,99
lsel,a,,,102
lesize,all,,,4
allsel
vmesh,1
/replot
asel,s,area,,1,32,1
aplot
cm,crp,area
aatt,2,2,2,0
lsel,s,line,,12,96,12
lsel,a,line,,7,91,12
lesize,all,,,4
lsel,s,line,,2,86,12
lsel,a,line,,8,92,12
lesize,all,,,5
lsel,s,line,,3,87,12
lsel,a,line,,1,85,12
lsel,a,line,,6,90,12
lsel,a,line,,9,93,12
lesize,all,,,1
mshape,1,2D
mshkey,1
amesh,all
allsel,all
eplot
ALLS
ASEL, S,,,33
NSLA, S, 1
TYPE, 3
REAL, 3
ESURF
ALLS
LSEL, S,5,89,12
LSEL, S,7,91,12
NSLL, S, 1
TYPE, 4
REAL,3
ESURF
ALLS
ASEL, S,,,34
NSLA, S, 1
TYPE, 3
REAL, 4
ESURF
ALLS
LSEL, S,11,95,12
LSEL, S,12,96,12
NSLL, S, 1
TYPE, 4
REAL,4
ESURF
ALLS
ASEL, S,,,35
NSLA, S, 1
TYPE, 3
REAL, 5
ESURF
ALLS
LSEL, S,2,86,12
LSEL, S,4,88,12
NSLL, S, 1
TYPE, 4
REAL,5
ESURF
ALLS
ASEL, S,,,36
NSLA, S, 1
TYPE, 3
REAL, 6
ESURF
ALLS
LSEL, S,8,92,12
LSEL, S,10,94,12
NSLL, S, 1
TYPE, 4
REAL,6
ESURF
ALLS
SHSD, 3, CREATE
SHSD, 4, CREATE
SHSD, 5, CREATE
SHSD, 6, CREATE
/solu
antype,static
nsel,s,loc,x,0
nsel,r,loc,z,0
d,all,ux
d,all,uy
d,all,uz
nsel,s,loc,x,1.4
nsel,r,loc,z,0
d,all,uy
d,all,uz
allsel
f,252,fz,-35000
f,282,fz,-35000
allsel
solve
FINISH
结果算出来也不对,求教大神啊, 在线等回复
![毕业论文急需要ansys模型,算出来结果不对啊,-1]() |