24小时热门版块排行榜

返回列表

jiushiwo0539

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 38.8
帖子: 18
在线: 5.1小时
虫号: 1505078
注册: 2011-11-22
专业: 机械结构强度学

[求助] ansys热分析高斯热源模型热流密度施加不上去温度场不发生变化

我用ansys做的热分析，是对一个钢板进行的热流密度施加，但是热流密度好像施加不上去，命令流能运行，但是用后处理器查看温度场没有变化，用时间历程处理器查看加热区某一节点的温度场也没有发生变化，不知道为什么，希望哪位高手能指点一下，具体命令流如下：
【我是刚注册的号，金币不多，望请见谅！】

我也附加传了一份txt文本附件

finish
/clear
/filname,Fname
/title,FEM model of heating
/unit,SI

!**********************modeling******************************
/prep7
et,1,solid70    !3D热实体
et,2,solid87    !三维10节点4面体热单元
et,3,surf152    !接触单元
keyopt,3,4,1
keyopt,3,8,1

!******************材料常数***********************************
MPTEMP,1,0,20,75,100,175,200
MPTEMP,7,225,300,325,375,400,450
MPTEMP,13,475,500,550,575,600,675
MPTEMP,19,700,720,725,775,800,900
MPTEMP,25,1000,1100,1200,1500

MPDATA,KXX,1,1,51.9,51.74,51.3,51.1,49.525,49
MPDATA,KXX,1,7,48.275,46.1,45.25,43.55,42.7,41.05
MPDATA,KXX,1,13,40.225,39.4,37.5,36.55,35.6,32.75
MPDATA,KXX,1,19,31.8,30.64,30.35,27.45,26,26.6
MPDATA,KXX,1,25,27.2,27.7,28.2,29.7

MPDATA,C,1,1,450,459.6,486,494.25,519,525.5
MPDATA,C,1,7,532,563.5,574,599,614.75,646.25
MPDATA,C,1,13,662,683.75,727.25,749,773.25,846
MPDATA,C,1,19,1139,1373.4,1432,950,931.03,855.17
MPDATA,C,1,25,779.31,703.45,627.59,400

MPDATA,DENS,1,1,7880,7880,7880,7880,7880,7800
MPDATA,DENS,1,7,7800,7800,7790,7790,7780,7750
MPDATA,DENS,1,13,7720,7710,7710,7710,7700,7700
MPDATA,DENS,1,19,7650,7650,7600,7600,7588,7500
MPDATA,DENS,1,25,7400,7320,7200,7258

MPDATA,EX,1,1,2.08e11,2.06e11,2.02e11,2.0e11,1.99e11,1.98e11
MPDATA,EX,1,7,1.92e11,1.9e11,1.87e11,1.86e11,1.85e11,1.8e11
MPDATA,EX,1,13,1.77e11,1.75e11,1.7e11,1.63e11,1.65e11,1.3e11
MPDATA,EX,1,19,1.2e11,1.38e11,1.35e11,0.9e11,0.8e11,0.6e11
MPDATA,EX,1,25,0.4e11,0.25e11,0.25e11,0.25e11

MPDATA,PRXY,1,1,0.29,0.29,0.3,0.31,0.325,0.33
MPDATA,PRXY,1,7,0.335,0.35,0.355,0.36,0.367,0.37
MPDATA,PRXY,1,13,0.38,0.39,0.395,0.4,0.405,0.42
MPDATA,PRXY,1,19,0.425,0.43,0.43,0.435,0.44,0.445
MPDATA,PRXY,1,25,0.45,0.475,0.475,0.475

MPDATA,ALPX,1,1,10e-6,10e-6,10e-6,11e-6,11e-6,11e-6
MPDATA,ALPX,1,7,11e-6,12e-6,12e-6,12e-6,12e-6,13e-6
MPDATA,ALPX,1,13e-6,13e-6,13e-6,14e-6,14e-6,14e-6,14e-6
MPDATA,ALPX,1,19,14e-6,14e-6,14e-6,14e-6,14e-6,14e-6
MPDATA,ALPX,1,25,14e-6,14e-6,15e-6,15e-6

block,0,500,0,400,0,-14    !plate
block,0,40,0,250,0,-14 !flame region 1
block,0,90,0,280,0,-14 !transition region 5
block,0,90,0,400,0,-14 !up left region 2
block,90,500,0,280,0,-14 !upper right rigion 3
boptn,keep,no    !布尔设置
vovl,all          !体分开
numcmp,kp          !压缩代号
numcmp,line
numcmp,area
numcmp,volu

!*********************meshing************************
!meshing the flame rigion
type,1
lsel,s,line,,21,22,1
lsel,a,line,,5,8,3
lesize,all,8,,,1
lsel,s,line,,6,7,1
lsel,a,line,,23,24,2
lesize,all,8,,,1
lsel,s,line,,1,11,10
lesize,all,,,3,2
lsel,a,line,,9,10,1
lesize,all,,,3,0.5
lsel,all
vsel,s,volu,,1
vmesh,all
vsel,all

!meshing the linear rigion
type,1
lsel,s,line,,25,26,1
lsel,a,line,,39,40,1
lesize,all,30,,,1
lsel,a,line,,29,30,1
lesize,all,30,,,1
lsel,s,line,,13,14,1
lsel,a,line,,33,34,1
lesize,all,30,,,1
lsel,s,line,,12,15,3
lsel,a,line,,27,28,1
lesize,all,40,,,1
lsel,s,line,,31,32,1
lsel,a,line,,37,38,1
lesize,all,40,,,1
lsel,s,line,,35,36,1
lesize,all,40,,,1
lsel,all
vsel,s,volu,,2,4,1
vmesh,all
vsel,all

!meshing the transition region    !过渡网格划分
type,2
lsel,s,line,,43,44,1
lesize,all,,,2,1
lsel,s,line,,41,42,1
lesize,all,,,3,1
mopt,expnd,1.5
mopt,tetexpnd,2
mopt,trans,1.5
lsel,all
vsel,s,volu,,5
vmesh,all
vsel,all
!表面有效单元的划分
type,3
asel,s,area,,8
amesh,all
asel,all
numcmp,node
numcmp,elem

save
finish

!****************施加热流 heating flux  *********************

finish
/clear
resume,Fname,db
/filname,Fname1
/title,termperature filed of heating

*set,vatl,817.3                            !flow rate of acetylene(l/h)
*set,ata,0.29                               !heat transfer efficiency
*set,v,4                                  !flame moving speed
*set,tmptint,23                            !initial temperature
*set,k,3/40/40                         !for heating source
*set,q0,k*ata*11470*4.184*vatl/3600/3.1416
*set,timestep,8/v                            !time step length(time stayed on each element)
*set,tend,210/v                         !the end time of heating

!************************** define components ******************
!define the node component of flame region.
asel,s,area,,8                         !upper flame area.
nsla,s,1                               !all nodes.
cm,flmrg,node                         !flame region component

!define the upper surface
asel,s,area,,12,17,5
asel,a,area,,22,26,4                   !upper surface except flame area.
nsla,s,1                               !all nodes
cm,upsfc,node
!define the whole upper surface
cmsel,s,flmrg
cmsel,a,upsfc
cm,whupsfc,node

!define the below surface
asel,s,area,,7,11,14
asel,a,area,,16,21,5
asel,a,area,,25                         !below surface.
nsla,s,1
cm,blsfc,node

!define the around surface
asel,s,area,,5,10,5
asel,a,area,,14,15,1
asel,a,area,,19,20,1
asel,a,area,,23                         !around surface.
nsla,s,1
cm,rdsfc,node

!define the 152 surface element component
esel,s,type,,3
cm,htflrg,elem                         !heat flux elements

!****************************  heating  ******************************
/solu
antype,TRANS

!solution setting
tunif,tmptint                         !uniform temperature for the fist load step
TREF,TMPTINT
nsubst,2,5,1                            !two substeps.
kbc,0                                  !ramped load.
outres,nsol,last                      !write node sesults at the last substep
outpr,all,none                         !none print out.
eqslv,sparse,le-5                      !pcg solver  精度设置

!**************************** ****************************
!load step loop during heating
*set,Tstart,0.01
*set,Tend,tend
*set,Tinc,timestep
*do,TM,Tstart,Tend,Tinc
   *set,xorg,0
   *set,yorg,v*(TM-Tstart)
   *set,zorg,0
   clocal,11,1,xorg,yorg,zorg
   *if,TM-Tstart,le,Tend,then       !flame heating
   csys,11
   !flame round
   *set,R0,40                         !the outside flame round radius
   *set,Rend,8                      !the inner flame round rsdius
   *set,Rinc,-8                   !flame round step.
   *do,r,R0,Rend,Rinc
   cmsel,s,flmrg                   !nodes in flame region surface.
   nsel,r,loc,x,0,r-0.05*8
   esln
   esel,r,type,,3
   *set,Rmid,r+Rinc/2
   sfe,all,,hflux,,q0*exp(-1*k*Rmid*Rmid)
   *enddo
   cmsel,s,flmrg
   nsel,r,loc,x,0,R0
   esln,s,1
   cm,flame,elem                   !current flame round
   cmsel,s,htflrg                   !surface 152 elelements
   cmsel,u,flame                   !non current flame region
   sfedel,all,,hflux                   !delete last flame load
*endif

!delete heat flux after the last flame heating step.
*if,TM-Tstrat,eq,Tend+timestep,then
   cmsel,s,flame
   sfedel,all,,hflux                !delete heat flux
*endif

   allsel
   time,TM
   csys,0
   solve
*enddo

回复此楼