ansys表面效应单元的问题
如果我把一部分面加入表面效应单元,但是我不在表面效应单元上加载荷,最后对计算结构有影响吗? 我补充一下:如果我在一个面上建了一个表面效应单元,但是我加载加到原来的节点上,但是表面效应单元还在,这样算出来结果有差别的,这是为什么? 应该是没有差别的表面效应单元是为了施加载荷方便,只起到过度的作用 glise 发表于 2012-12-19 10:12 static/image/common/back.gif
应该是没有差别的
表面效应单元是为了施加载荷方便,只起到过度的作用
可是我算出来的结果有差别 想不通为啥 我附上我写的命令请高手帮我看看
/FILNAME,FILE121201
/TITLE,TEST1212
/PREP7
SMRT,OFF
ET,1,PLANE55
ET,2,SURF151
!KEYOPT,2,2,0
KEYOPT,2,4,1
KEYOPT,2,5,0
KEYOPT,2,8,2
MP,KXX,1,0.0045
MP,DENS,1,1.093
MP,C,1,3.37
N=0
!通过关键点建立模型,模型大小为0.1*0.1,在x=0.025 及0.05出生成血管,施加对流
K,N+1,0,0
K,N+2,0.1,0 !ccccc
K,N+3,0,0.1 !ccccc
K,N+4,0.1,0.1 !ccccc
K,N+5,0.025,0 !ccccc
K,N+6,0.025,0.1 !ccccc
K,N+7,0.05,0 !ccccc
K,N+8,0.05,0.1 !ccccc
L,1,5
L,5,7
L,7,2
L,2,4
L,4,8
L,8,6
L,6,3
L,3,1
L,5,6
L,7,8
AL,8,1,9,7
AL,2,10,6,9
AL,3,4,5,10
!划分血管网格
ASEL,S,AREA,,2
AATT,1,,1,0
MSHAPE,0,2D
MSHKEY,1
ESIZE,0.00125 !cccccccc
AMESH,ALL
!血管壁上覆盖表面单元,生成表面效应单元,line 9 and 10 分别为x=0.025,x=0.05
LSEL,S,LINE,,9,10
!LSEL,A,LINE,,2
NSLL,S,1
TYPE,2
ESURF
!划分其他area网格
ASEL,S,AREA,,1,3,2
AATT,1,,1,0
MSHAPE,0,2D
MSHKEY,1
ESIZE,0.00125 !ccccccccc
AMESH,ALL
!下边界和左边界生成表面单元
!LSEL,S,LINE,,1,3,2
!LSEL,A,LINE,,8
!NSLL,S,1
!TYPE,2
!ESURF
!给血管表面单元施加对流载荷
!/PSYMB,ESYS,1
!ESEL,S,TYPE,,2
!LSEL,S,,,9,10
!NSLL,S,1
!ESLN,S,1
!SFE,ALL,,CONV,1,0.0684
!SFE,ALL,,CONV,2,37
!/psf,conv,,2,0,1
!对下边界施加常温对流载荷
!/PSYMB,ESYS,1
!LSEL,S,,,1,3
!NSLL,S,1
!ESLN,S,1
!SFE,ALL,,CONV,1,0.00062
!SFE,ALL,,CONV,2,25
!/psf,conv,,2,0,1
!对轴上施加绝热边界条件
!LSEL,S,,,8
!NSLL,S,1
!ESLN,S,1
!SFE,ALL,,CONV,1,0
!SFE,ALL,,CONV,2,37
!D,ALL,TEMP,37
!对右边界和上边界施加恒温条件
!LSEL,S,,,4,7
!NSLL,S,1
!D,ALL,TEMP,37
/SOLU
!定义分析类型
ANTYPE,4 !瞬态分析
!定义初始条件
TUNIF,37, !初始温度37
!对节点施加热生成率
*DIM,uhgen,,81,81
*VREAD,uhgen(1,1),B105T4,txt,,JIK,81,81
(81F8.3)
*STATUS, uhgen
YVALUE=0
XVALUE=0
XSTEP=0.00125 !ccccc
YSTEP=0.00125 !ccccc
*DO,II,1,81 !X Direction
*DO,JJ,1,81 !Y Direction
!选定节点时候需要给定一个较小的范围
Nn=NODE(XVALUE,YVALUE,0)
BF,Nn,HGEN,UHGEN(II,JJ)
ALLSEL
YVALUE=YVALUE+YSTEP
*ENDDO
XVALUE=XVALUE+XSTEP
YVALUE=0
*ENDDO
allsel
*DO,T1F1,1,300
time, T1F1
kbc,1
autots,on
deltim,0.5,0.5,1
timint,on
tintp,,,,1
outres,nsol,2,
LSWRITE,T1F1,
*ENDDO
*DEL,UHGEN
*DEL,YVALUE
*DEL,XVALUE
*DEL,ZVALUE
*DEL,XSTEP
*DEL,YSTEP
*DEL,ZSTEP
*DEL,II
*DEL,JJ
*DEL,KK
BFDELE,ALL,HGEN
ALLSEL
*DO,T1F2,301,600
time, T1F2
kbc,1
autots,off
deltim,1,1,2
timint,on
tintp,,,,1
outres,nsol,2,
LSWRITE,T1F2,
*ENDDO
Lssolve,T1F1,T1F2
我用简单例子做过对比,施加表面效应单元与不施加计算结果一样,施加在表面效应单元上的压强载荷。
页:
[1]