Ansys碳纖維分析實(shí)例（CFRP）

2016-12-12 by:CAE仿真在線來(lái)源:互聯(lián)網(wǎng)

針對(duì)葉片中的梁、殼等復(fù)合材料層合結(jié)構(gòu),ANSYS提供了一系列的特殊單元——結(jié)構(gòu)多層復(fù)合材料單元,以模擬各種復(fù)合材料。鋪層單元中可以考慮復(fù)合材料特有的鋪層特性和各向異性特性。

本計(jì)算采用的是相對(duì)簡(jiǎn)單的線性鋪層單元Shell99。該單元是一種八節(jié)點(diǎn)3D殼單元,每個(gè)節(jié)點(diǎn)有六個(gè)自由度,主要適用于薄到中等厚度的板和殼結(jié)構(gòu),一般要求寬厚比應(yīng)大于10。Shell99可實(shí)現(xiàn)多達(dá)250層的等厚材料層,或者125層厚度在單元面內(nèi)呈現(xiàn)雙線性變化的不等厚材料層。如果材料層大于250,用戶可通過(guò)輸入自己的材料矩陣形式來(lái)建立模型,還可以通過(guò)一個(gè)選項(xiàng)將單元節(jié)點(diǎn)偏置到結(jié)構(gòu)的表層或底層。

關(guān)于ansys有限元軟件對(duì)復(fù)合材料的支持,可參考:ANSYS復(fù)合材料分析技術(shù)介紹

下面是ansy分析碳纖維FRP【纖維增強(qiáng)復(fù)合材料(Fiber Reinforced Polymer/Plastic,簡(jiǎn)稱FRP)】加固混凝土實(shí)例,APDL命令流如下:

*set,h,80
*set,B,500/2
*set,L,1800
*set,LF,100
*set,LS,50
*set,a,20
*set,PI,ACOS(-1)
*set,SR,PI*16
*set,RO,PI*(6/2)**2/200/h
*set,CB,300
*set,Ch,0.167
*set,CL,L-2*LS
*set,F1,0.002
*set,F2,30/LF
!定義完參數(shù),開(kāi)始建立模型
/prep7
et,1,solid65
et,2,link8
et,3,shell41
keyopt,3,1,2
et,4,solid45
mp,ex,1,18000
mp,prxy,1,0.2
!定義MISES本構(gòu)關(guān)系
tb,kinh,1,1,7
tbpt,,0.0001,1.8
tbpt,,0.0004,6.66
tbpt,,0.0008,11.84
tbpt,,0.0012,15.54
tbpt,,0.0016,17.76
tbpt,,0.002,18.5
tbpt,,0.0033,18.5
tb,conc,1,1,9
tbdata,,0.3,0.5,1.75,-1
!鋼筋雙線性隨動(dòng)強(qiáng)化
mp,ex,2,2.1e5
mp,prxy,2,0.3
tb,bkin,2,1,2,1
tbdata,,235,0
!碳纖維線彈性
mp,ex,3,2.35e5
mp,prxy,3,0
!支座及整體鋼筋
mp,ex,4,2.1e5
mp,prxy,4,0.3
r,1,2,RO,0,0,,,
rmore,,,,,,,
r,2,SR,,
r,3,CH,,
RMORE,,
!建立板模型
block,0,-B,0,h,-LF/2,-(L/2-LS)
lplot
/pnum,line,1
/replot
lgen,2,8,,,50,,,,0
lgen,2,13,,,100,,,,0
adrag,13,,,,,,9
adrag,14,,,,,,9
/pnum,line,0
/replot
v***a,1,7
v***a,3,8
!建立受力鋼筋模型
lsel,s,loc,z,-LF/2
Lsel,r,loc,Y,0
lgen,2,all,,,,a,,,0
lsel,s,loc,z,-50
lsel,r,loc,y,A
adrag,all,,,,,,9
allsel,all
v***a,2,15
v***a,4,14
v***a,1,13
lsel,s,loc,y,0
lsel,a,loc,y,A
lsel,a,loc,y,H
lesize,all,50
lsel,s,loc,z,-LF/2
lsel,a,loc,z,-(L/2-LS)
lsel,u,loc,y,0
lsel,u,loc,y,a
lsel,u,loc,y,h
lesize,all,20
allsel,all
type,1
mat,1
real,1
esys,0
mshape,0,3D
mshkey,1
vmesh,all
asel,s,loc,z,-LF/2
!完成另一半板模型
type,1
extopt,esize,LF/50,0
extopt,aclear,1
extopt,attr,0,0,0
mat,1
real,1
esys,0
vext,all,,,0,0,LF
asel,s,loc,z,LF/2
type,1
extopt,esize,(L/2-LF/2-LS)/50,0
extopt,aclear,1
extopt,attr,0,0,0
mat,1
real,1
esys,0
vext,all,,,0,0,L/2-LF/2-LS
allsel,all
!生成支座模型
asel,s,loc,z,L/2-LS
type,1
extopt,esize,1,o
extopt,aclear,1
extopt,attr,0,0,0
mat,1
real,1
esys,0
vext,all,,,0,0,25
asel,s,loc,z,L/2-LS+25
type,1
extopt,esize,2,o
extopt,aclear,1
extopt,attr,0,0,0
mat,1
real,1
esys,0
vext,all,,,0,0,LS
asel,s,loc,z,-(L/2-LS)
type,1
extopt,esize,1,o
extopt,aclear,1
extopt,attr,0,0,0
mat,1
real,1
esys,0
vext,all,,,0,0,-25
asel,s,loc,z,-(L/2-LS+25)
type,1
extopt,esize,2,o
extopt,aclear,1
extopt,attr,0,0,0
mat,1
real,1
esys,0
vext,all,,,0,0,-LS
allsel,all
eplot
!劃分鋼筋單元
lsel,s,loc,x,-200
lsel,a,loc,x,-100
lsel,r,loc,y,A
type,2
mat,2
real,2
esys,0
lmesh,all
lsel,s,loc,x,0
lsel,r,loc,y,A
type,2
mat,2
real,2
esys,0
lmesh,all
allsel,all
!劃分碳纖維單元
ksel,s,loc,z,850
ksel,r,loc,y,0
ksel,r,loc,x,0
kgen,2,all,,,-CB/2,,,,0
lstr,41,97
allsel,all
lstr,41,4
adrag,221,,,,,,222
asel,s,,,168
lsla,s
lesize,all,50,,,,,,,1
type,3
mat,3
real,3
esys,0
mshape,1,2D
amesh,all
allsel,all
nummrg,node,,,,low
numcmp,node
!劃分支座單元
asel,s,loc,y,0
asel,r,loc,z,L/2-25,L/2+25
type,4
extopt,esize,1,0
extopt,aclear,1
extopt,attr,0,0,0
mat,4
real,3
esys,0
vext,all,,,0,-20
asel,s,loc,y,0
asel,r,loc,z,-(L/2-25),-(L/2+25)
type,4
extopt,esize,1,0
extopt,aclear,1
extopt,attr,0,0,0
mat,4
real,3
esys,0
vext,all,,,0,-20
allsel,all
finish
!設(shè)置邊界
/solu
nsel,s,loc,x,0
/go
dsym,symm,x
nsel,s,loc,z,L/2
nsel,r,loc,y,-20
/go
D,all,,,,,,UY,
nsel,s,loc,z,-L/2
nsel,r,loc,y,-20
/go
D,all,,,,,,UY,UZ,
!加載1
antype,0
nlgeom,1
nropt,full
eqslv,spar,,0
allsel,all
outres,all,all
time,1
autots,1
nsubst,10,20,5,1
cnvtol,f,,0.05,2,,
neqit,30
pred,on,,on
esel,mat,3
ekill,all
allsel,all
solve
!加載2
asel,s,loc,y,h
asel,r,loc,z,-LF/2,LF/2
/go
sfa,all,1,pres,F2
outres,all,all
time,2
autots,0
nsubst,30,,,1
pred,-1
arclen,1,10,0,
esel,mat,3
ealive,all
allsel,all
solve
finish
!查看結(jié)果
/post1
esel,mat,1
set,1,last,1,
plnsol,s,z,0,1
set,2,last,1,
plnsol,s,z,0,1
finish
!查看碳纖維結(jié)果
/post26
allsel,all
nsel,s,loc,x,0
nsel,r,loc,y,0
nsel,r,loc,z,0
*get,N1,node,0,num,max
nsol,2,N1,U,Y
prod,3,2,,,,,,-1,1,1,
xvar,3
plvar,1
finish