1年前
好的,我使用link33模拟热传导,使用link34模拟对流换热,由于link34单元节点可以重合,所以我建立link34与link33单元节点重合,设立结构初始温度为500度,环境温度为40度,对流系数为12.5,代码如下:
finish
/clear
/prep7
k,1,0,-2.5,1.66
k,2,0,0,1.315
k,3,0,2.525,1.756
k,4,0,4.975,2.125
k,5,0,7.425,2.436
k,6,0,9.875,2.69
k,7,0,12.325,2.888
k,8,0,14.575,3.02
k,9,0,16.675,3.105
k,10,0,18.775,3.141
k,11,0,19.75,3.145
k,12,-2,-2.5,1.66
k,13,-2,0,1.315
k,14,-2,2.525,1.756
k,15,-2,4.975,2.125
k,16,-2,7.425,2.436
k,17,-2,9.875,2.69
k,18,-2,12.325,2.888
k,19,-2,14.575,3.02
k,20,-2,16.675,3.105
k,21,-2,18.775,3.141
k,22,-2,19.75,3.145
!雨棚下弦杆
k,23,0,-1.1,0.615
k,24,-2,-1.1,0.615
k,25,0,0,-0.205
k,26,-2,0,0
k,27,-1,1.263,0.141
k,28,-1,3.75, 0.752
k,29,-1,6.2, 1.264
k,30,-1,8.65,1.688
k,31,-1,11.1,2.025
k,32,-1,13.45,2.27
k,33,-1,15.625,2.427
k,34,-1,17.725,2.517
k,35,-1,19.75,2.545
!柱
k,36,0,0,-2.705
k,37,0,0,-5.205
k,38,0,0,-7.205
k,39,-2,0,-2.705
k,40,-2,0,-5.205
k,41,-2,0,-7.205
k,42,-1,13.45,-7.205
k,101,0,-2.4
k,102,0,-0.5
k,103,0,0.5
k,104,0,2.3
k,105,-2,4.1
k,106,-2,5.9
k,107,-2,7.7
k,108,-2,9.7
k,109,-2,11.5
!=============================建立几何模型=============================
!================董瑞琪雨棚桁架部分==================
*do,i,36,37,1
l,i,i+1
*enddo
*do,i,39,40,1
l,i,i+1
*enddo
l,25,36
l,26,39
*do,i,36,37,1
l,i,i+3
*enddo
l,25,2
l,13,26
!中间柱
l,42,32
!上(右)弦杆连接
*do,i,1,10,1
l,i,i+1
*enddo
!上(左)弦杆连接
*do,i,12,21,1
l,i,i+1
*enddo
!下弦杆连接
*do,i,27,34,1
l,i,i+1
*enddo
!上线横杆
*do,i,1,10,1
l,i,i+11
*enddo
l,23,24
l,25,26
!斜腹杆
*do,i,2,10,1
l,i,i+25
*enddo
*do,i,3,10,1
l,i,i+24
*enddo
*do,i,13,21,1
l,i,i+14
*enddo
*do,i,14,21,1
l,i,i+13
*enddo
l,1,23
l,12,24
l,23,25
l,24,26
l,25,27
l,26,27
l,13,24
l,2,23
!上弦斜腹杆
*do,i,1,9,2
l,i,i+12
*enddo
*do,i,3,9,2
l,i,i+10
*enddo
wprota,,90
*do,i,101,109,1
wpave,,ky(i)
lsel,s,loc,x,-2
lsbw,all
lsel,s,loc,x,0
lsbw,all
*enddo
allsel,all
wpcsys,-1
nummrg,all
!============定义单元类型================
!********定义并赋予杆件材料属性,实常数、单元类型、截面信息、
et,1,link33!*******单元类型
MP,KXX,1,58.2
MP,dens,1,7850
MP,c,1,480
et,2,link34!*******单元类型
keyopt,2,3,2
w1=0.168/2
n1=(0.168-0.008)/2
w2=0.180/2
n2=(0.180-0.012)/2
w3=0.180/2
n3=(0.180-0.008)/2
w4=0.121/2
n4=(0.121-0.006)/2
w5=0.07/2
n5=(0.07-0.004)/2
w6=0.402/2
n6=(0.402-0.012)/2
pi=3.1415
*do,i,1,6
a%i%=pi*(w%i%*w%i%-n%i%*n%i%)
*enddo
*do,i,1,6
r,i,a%i%
*enddo
!================董瑞琪==================
!*************赋截面************
!*************下弦杆1左四根************
lsel,s,loc,x,-2,0
lsel,r,loc,y,1.263+0.1,11.1-0.1
lsel,r,loc,z,0.141+0.1,2.025-0.1
lsel,r,loc,x,-1
lsel,u,loc,y,2.525
latt,1,3,1
lesize,all,,,2
lmesh,all
allsel,all
!*************下弦杆1中2根************
lsel,s,loc,x,-2,0
lsel,r,loc,y,15.625,19.75
lsel,r,loc,z,2.427,2.545
latt,1,3,1
lesize,all,,,2
lmesh,all
allsel,all
!*************下弦杆3中左1×2根************
lsel,s,loc,x,-2,0
lsel,r,loc,y,0,1.263
lsel,r,loc,z,0,0.141
lsel,u,loc,x,-1
latt,1,2,1
lesize,all,,,2
lmesh,all
allsel,all
lsel,s,loc,x,-1+0.001,0
lsel,r,loc,y,0,1.263
lsel,r,loc,z,-0.205,0.141
latt,1,2,1
lesize,all,,,2
lmesh,all
allsel,all
!*************下弦杆3中3根************
lsel,s,loc,x,-2,0
lsel,r,loc,y,11.1,15.625
lsel,r,loc,z,2.025,2.427
latt,1,2,1
lesize,all,,,2
lmesh,all
allsel,all
!*************上弦杆3左1×2根************
lsel,s,loc,x,-2,0
lsel,r,loc,y,0,2.525
lsel,r,loc,z,1.315,1.756
lsel,u,loc,x,-1
latt,1,2,1
lesize,all,,,2
lmesh,all
allsel,all
!*************上弦杆3中左两根************
lsel,s,loc,x,-2,0
lsel,r,loc,y,12.325,14.575
lsel,r,loc,z,2.888,3.02
lsel,u,loc,x,-1
latt,1,2,1
lesize,all,,,2
lmesh,all
allsel,all
!*************上弦杆2左4×2根************
lsel,s,loc,x,-2,0
lsel,r,loc,y,2.525,12.325
lsel,r,loc,z,1.756,2.888
lsel,u,loc,x,-0.5,-1.5
latt,1,1,1
lesize,all,,,2
lmesh,all
allsel,all
!*************上弦杆1中2.5×2根************
lsel,s,loc,x,-2,0
lsel,r,loc,y,14.575,19.75
lsel,r,loc,z,3.02,3.145
lsel,u,loc,x,-1
latt,1,1,1
lesize,all,,,2
lmesh,all
allsel,all
!*************上弦与下弦之间的腹杆2左侧************
lsel,s,loc,x,-2,0
lsel,r,loc,x,-0.5
lsel,u,loc,z,0,0.141
latt,1,4,1
lesize,all,,,2
lmesh,all
allsel,all
!*************上弦与下弦之间的腹杆2右侧************
lsel,s,loc,x,-2,0
lsel,r,loc,x,-1.5
lsel,u,loc,z,0,0.141
latt,1,4,1
lesize,all,,,2
lmesh,all
allsel,all
!*************雨棚柱腹杆ZFG************
lsel,s,loc,x,-2,0
lsel,r,loc,y,0
lsel,r,loc,x,-1
lsel,r,loc,z,0,-7.205
latt,1,1,1
lesize,all,,,2
lmesh,all
allsel,all
!*************雨棚柱YPZ2左边柱************
lsel,s,loc,x,-2,0
lsel,r,loc,y,0
lsel,u,loc,x,-1
latt,1,6,1
lesize,all,,,2
lmesh,all
allsel,all
!*************雨棚柱YPZ2中间柱************
lsel,s,loc,x,-2,0
lsel,r,loc,y,13.45
lsel,r,loc,z,-7.205,2.27
latt,1,6,1
lesize,all,,,2
lmesh,all
allsel,all
!*************上弦杆3左角1×2根************
lsel,s,loc,x,-2,0
lsel,r,loc,y,0,-2.5
lsel,r,loc,z,1.315,1.66
lsel,u,loc,x,-1
latt,1,2,1
lesize,all,,,2
lmesh,all
allsel,all
!*************下弦杆1左角1×2根上************
lsel,s,loc,x,-2,0
lsel,r,loc,y,-1.1,-2.5
lsel,r,loc,z,0.615+0.1,1.66-0.2
lsel,u,loc,x,-1
latt,1,3,1
lesize,all,,,2
lmesh,all
allsel,all
!*************下弦杆1左角1×2根下************
lsel,s,loc,x,-2,0
lsel,r,loc,y,0,-1.1
lsel,r,loc,z,0,0.615
lsel,u,loc,x,-1
latt,1,3,1
lesize,all,,,2
lmesh,all
allsel,all
!*************上弦斜杆和横杆************
lsel,s,loc,x,-2,0
lsel,r,loc,x,-1
lsel,r,loc,z,1.315,2.436
lsel,u,loc,y,6.2+0.3,15.625
latt,1,5,1
lesize,all,,,2
lmesh,all
allsel,all
!*************上弦斜杆和横杆************
lsel,s,loc,x,-2,0
lsel,r,loc,x,-1
lsel,r,loc,z,2.436,3.141
lsel,u,loc,y,15.625,19.75
latt,1,5,1
lesize,all,,,2
lmesh,all
allsel,all
!*************上弦斜杆和横杆************
lsel,s,loc,x,-2,0
lsel,r,loc,x,-1
lsel,r,loc,z,3.02,3.141
lsel,u,loc,y,14.575
latt,1,5,1
lesize,all,,,2
lmesh,all
allsel,all
!*************左边腹杆1×2根FG2************
lsel,s,loc,x,-2,0
lsel,r,loc,y,-1.1,0
lsel,u,loc,x,-1
lsel,r,loc,z,0.615,1.315
lsel,u,loc,y,0
latt,1,4,1
lesize,all,,,2
lmesh,all
allsel,all
!*************边上腹杆(补充)************
lsel,s,loc,y,-1.1
latt,1,4,1
lesize,all,,,2
lmesh,all
allsel,all
!*************(补充)************
lsel,s,loc,y,-0.5+0.001,0-0.001
lsel,r,loc,x,0
lsel,r,loc,z,-0.205,0.1677273
latt,1,3,1
lesize,all,,,2
lmesh,all
allsel,all
esel,s,type,,1
nsle,s,1
bfunif,temp,500
*GET,Nmax,NODE,,NUM,Max
*GET,Lmax,line,,NUM,Max
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!===按实常数选择线====
*do,i,1,6
lsel,s,real,,i
!====得到所选模型线总数====
*GET,LTOTAL%i%,LINE,,COUNT
!====定义数组====
*DIM,LNUM%i%,ARRAY,LTOTAL%i%
!====得到所选模型最小线号====
*GET,LNUM%i%(1),LINE,,NUM,MIN
!====将模型中线号从小到大存储到数组LNUM%i%中====
*DO,J,2,LTOTAL%i%
!====将比线号LNUM%i%(j-1)大的最小线号赋予给lnum%i%(j)
LNUM%i%(J)=LSNEXT(LNUM%i%(J-1))
*ENDDO
!====按线号从小到大顺序计算长度,并建立实常数
!====定义数组====
*DIM,LLENGTH%i%,ARRAY,LTOTAL%i%
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!====i是实常数编号
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!====ltotal%i%:线总数数组,ltotal1是实常数为1时,线的总数变量
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!====m是为ltotal1中线个数设立的变量,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!====LLENGTH%i%:线长度数组,LLENGTH1实常数为1时线的长度数组
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!====三层循环嵌套从m开始=========
*DO,m,1,LTOTAL%i%
*GET,LLENGTH%i%(m),LINE,LNUM%i%(m),LENG
r,LNUM%i%(m)+LMAX,2*pi*w%i%*LLENGTH%i%(m)/2
!====LNUM%i%(m):线号数组LNUM%i%第m个元素
lsel,s,line,,LNUM%i%(m)
nsll,s,1
!====得到所选模型节点总数====
*GET,NTOTAL%m%,NODE,,COUNT
!====定义数组====
*DIM,NNUM%i%%m%,ARRAY,NTOTAL%m%
!====得到所选模型最小节点号====
*GET,NNUM%i%%m%(1),NODE,,NUM,MIN
!====将模型中节点号从小到大存储到数组NNUM%im%中====
*DO,w,2,NTOTAL%m%
!====将比节点NNUM%im%(w-1)大的最小节点号赋予给NNUM%im%(w)
NNUM%i%%m%(w)=NDNEXT(NNUM%i%%m%(w-1))
*ENDDO
!====建立节点====
*DO,q,1,NTOTAL%m%
N,NNUM%i%%m%(q)+nmax,nx(NNUM%i%%m%(q)),ny(NNUM%i%%m%(q)),nz(NNUM%i%%m%(q))
allsel,all
!====建立对流单元====
Type,2
real,LNUM%i%(m)+Lmax
E,NNUM%i%%m%(q),NNUM%i%%m%(q)+nmax
*ENDDO
*ENDDO
*ENDDO
esel,s,type,,2
sfe,all,1,conv,1,12.5
sfe,all,1,conv,2,40
allsel,all
/solu
antype,trans
timinf,on
time,600
autots,on
outres,all
nsubst,20
allsel,all
kbc,1
allsel,all
eqslv,jcg
solve
@水哥ANSYS评论于:link34对流单元使用