提问link34对流单元使用

好的，我使用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