ansys表面效应单元的问题

lipard

如果我把一部分面加入表面效应单元，但是我不在表面效应单元上加载荷，最后对计算结构有影响吗？

lipard

我补充一下：如果我在一个面上建了一个表面效应单元，但是我加载加到原来的节点上，但是表面效应单元还在，这样算出来结果有差别的，这是为什么？

glise

应该是没有差别的
表面效应单元是为了施加载荷方便，只起到过度的作用

lipard

glise 发表于 2012-12-19 10:12

                               
登录/注册后可看大图

应该是没有差别的
表面效应单元是为了施加载荷方便，只起到过度的作用

可是我算出来的结果有差别 想不通为啥  

lipard

我附上我写的命令请高手帮我看看

/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

水面之下

我用简单例子做过对比，施加表面效应单元与不施加计算结果一样，施加在表面效应单元上的压强载荷。