CadCaeCam.com是目前CAD/CAE/CAM/PLM类专业网站中,用户最多,技术含量最高的网站之一, 涵盖目前所有常用的C3P类软件技术讨论。

C3P专业门户

 找回密码
 注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

查看: 1601|回复: 0

[求助] 请教高手:温度为什么不向外扩散?

[复制链接]
发表于 2010-5-13 21:39:32 | 显示全部楼层 |阅读模式
本人最近做了一个例子,模拟冻结管冻结土体的过程。问题是不论冻结管的温度有多低或时间有多长,其冻结效果就是冻结管周边一点点,冻结影响范围不向外扩散,百思不得其解。请教高手帮我看看。期待您的回应!
; R. C: q' s' V源程序如下:' `7 P0 c" S  ]+ Y' P, l0 E  K: c
: ]' R& i' |% P3 R8 s
finish( {. w" \6 P& T. }
/clear
9 y8 I0 l+ B. P& o9 p0 G5 x/filname,0511,1; T! b' g! ^0 E) j# P3 r  O( ?
/title,rec
% K& m) Z# y3 t: O- `- ?/prep7
5 X* [" U4 R) a+ aet,1,plane559 f6 w8 p5 N* ?: k: D1 t" K$ S
!mptemp,1,-27,23          !定义材料1---土体属性, E, Q5 h4 {5 P( ~- q, c  b/ L) n* R3 A" d
!mpdata,c,1,1,2.54,3.3459 h: w/ o6 X9 T
MP,C,1,2.36 _8 |6 Q  n; n

* a  h" `. i1 Tmptemp,1,-27,23) R% ~' |0 u% h4 N+ B, u
mpdata,kxx,1,1,2.12,1.42
1 T% h% D9 @1 U; Z% s/ Fmptemp,1,-27,23# u$ p( p9 m! t6 q! ~; I
mpdata,dens,1,1,1528,1910
0 q* B1 ^5 o0 U7 B! O; n+ W5 m- A- @mptemp,1,-27,-1.4,0,23
8 L: {: }( Y& M4 x! H+ fmpdata,enth,1,1,0,5.525e4,1.782e5,2.50e5
1 W( I  u$ c( E" t9 g4 x8 r! w# Noutres,,all" C8 `+ q$ S3 N
mp,dens,2,1200         !定义材料2--冻结管属性
( P2 D; M/ A; l4 R1 Nmp,c,2,2.500/ k" U+ Y2 E' m
mp,kxx,2,2.5009 P+ ]; l5 R7 o* j& G. I

# y2 O. [6 r$ b6 o2 C' f- e3 d$ bk,1,0,5.748,0          !定义46根冻结管   
% B$ v( G  ^7 g( a, vk,2,-0.699,5.693,0
2 m' \. s* c  P6 \! K# ^# }k,3,-1.391,5.589,0
6 I4 @1 [6 w* p+ ]# {7 m7 Zk,4,-2.070,5.417,0
6 w3 Z9 [, m- Nk,5,-2.728,5.179,0
6 T0 R) Z% E6 ~k,6,-3.359,4.877,0
$ W4 l& F- o1 U; l4 R7 x& ek,7,-4.040,4.457,00 O& P2 U- C9 r0 ?9 q
k,8,-4.710,4.062,0( p$ \6 O) s. F! s- h4 Z& Q
k,9,-5.238,3.400,0
* ~3 ~0 D$ \$ s3 M8 sk,10,-5.78,2.72,04 c) j$ E4 s3 t' R; I
k,11,-6.231,1.997,0
3 v4 B& _! M" B6 l! a8 Xk,12,-6.543,1.165,0
. I$ r% Q( F% k# Hk,13,-6.762,0.339,04 l# |9 V& _* P, B. r9 T  f% g" g
k,14,-6.669,-0.661,0
3 O/ c8 E+ i- v: ?* Lk,15,-6.522,-1.569,01 s- u  _, k5 M8 L8 j; `
k,16,-6.296,-2.45,01 K7 _7 o: J) F2 P$ k
k,17,-5.849,-3.275,0* y  m* ]  B0 J7 T4 C) o
k,18,-5.239,-3.963,08 P% H/ z7 H2 U, f
k,19,-4.466,-4.463,0
) K5 Q/ r4 Z4 Vk,20,-3.612,-4.804,0
5 x9 L2 P3 M+ b0 @k,21,-2.732,-5.071,0
/ B8 Q7 u6 C6 c8 O5 d% Tk,22,-1.832,-5.264,0
" x% q2 ]6 l" X7 X' H, x1 z8 xk,23,-0.919,-5.379,0
, }& @: o3 m8 _k,24,0,-5.418,0
7 m% T  D; M: }k,25,0.919,-5.279,0
. H; F3 ^' t6 ]% e0 T& @' B- Kk,26,1.832,-5.264,0( E* L- ?# C: w8 o" w
k,27,2.732,-5.071,0
& ?) `9 o8 \1 w3 }9 W6 rk,28,3.612,-4.804,0, L0 g& ?* f& n% e- j$ x5 ]
k,29,4.466,-4.463,0
2 m5 B+ `% G8 h  c" x5 W$ T: sk,30,5.239,-3.963,0
  E8 p  s2 x0 u5 T; Kk,31,5.849,-3.275,0! X' T8 D" B) ?
k,32,6.269,-2.45,06 q$ h/ @0 D6 m( R# R
k,33,6.522,-1.569,0
1 U8 ~& M" i6 M# R+ D1 dk,34,6.669,-0.661,04 M! @, E; L2 y8 d1 j- {1 j1 I8 b
k,35,6.692,0.259,0! ?) J9 y/ \) k- S! g1 I. b% l
k,36,6.543,1.165,0. m9 ~$ x5 c* S) {4 e/ H3 S/ [) @
k,37,6.251,1.977,0$ l5 ^+ f# o/ Y
k,38,5.780,2.72,0/ U! u( B: v* a" t* G9 g9 V
k,39,5.238,3.4,0! O  x6 P/ ~8 l+ j
k,40,4.669,3.962,0
3 f% B- n' i* Zk,41,4.04,4.457,0, w1 [* u! k# Z- ~6 p) Y( ^
k,42,3.359,4.877,08 m% W+ Q6 c% f7 ], ?' M9 J9 F- u$ ?' S
k,43,2.728,5.179,0* }  X3 x- E6 v# [
k,44,2.07,5.417,0
' v# g9 W; J; w& n0 Z; Qk,45,1.391,5.589,0
  d" e. t5 N" C( n6 X( g. D4 w' uk,46,0.699,5.693,0
" |, c% g% w2 ^9 c, F& c% t7 z  t, w& j' X1 o7 D5 Y/ ^2 k% a& d
circle,1,0.1, u8 [+ P* S1 r+ l' m! I# ~+ F
circle,2,0.19 R$ D7 E8 N' m' V- A: ^
circle,3,0.10 i# v) V+ e/ U; ]1 g( p, ?) u
circle,4,0.1( i8 g* Y( u* ]- `6 T
circle,5,0.1
" t/ w7 V7 j% Z/ c$ j5 w' kcircle,6,0.1
9 W1 E- I! b6 K( h4 A7 J# `circle,7,0.1
6 K7 ?: r; c9 v0 ?5 u4 h3 S4 Tcircle,8,0.1
2 }9 {, g2 a/ P2 Gcircle,9,0.14 I2 E7 a% S, E# r/ p, z# n, S+ V1 ~
circle,10,0.1
4 S% H# E) {9 T  k; k0 X1 ]% Lcircle,11,0.1
7 s( A0 N3 C! n* c" }" B7 ecircle,12,0.1: k% v* b% [0 q6 H
circle,13,0.1
5 \9 I; B0 {0 c, O% ?circle,14,0.1
9 _* }0 C9 x7 Kcircle,15,0.13 H$ S/ O# \. z) l% s' I9 z
circle,16,0.11 S% X, s- ^6 q7 x4 ~1 e) h
circle,17,0.1
# y2 F) p0 R6 z. D& fcircle,18,0.1* c2 ^6 b' i- {$ p9 n9 `8 p, u* e
circle,19,0.11 Z8 T7 V0 O: s
circle,20,0.1
7 J( K. q7 X) Z! M* H! Pcircle,21,0.1: ^3 F8 t* A' p1 i2 ]3 H
circle,22,0.1
* x3 [6 i3 @( K5 S  E" ccircle,23,0.19 d$ P$ X1 }5 R9 q- \) i) o) q* v
circle,24,0.18 E! u4 \" A) d- b1 ~: v1 Y
circle,25,0.1
% T" S4 v) H+ b( S: Lcircle,26,0.1
$ B! @5 m9 F2 Q; Z. U' h- fcircle,27,0.1
% N: t$ h) u* s- Z1 rcircle,28,0.1# g8 h! M2 @, d3 S8 Y: h
circle,29,0.14 W0 D* P& E" j: F
circle,30,0.1. N; ^1 [# t4 B0 T5 j
circle,31,0.11 g2 q1 i6 K- j& V& l' [' [7 v/ V
circle,32,0.1- n( u- q! u. y$ V, Y8 Q( I( w
circle,33,0.10 c4 r; j0 `' |, _5 s9 ]- Z8 R+ I
circle,34,0.1+ M4 t( z# N- _, O1 q! w  Q+ d: N
circle,35,0.1! H7 p/ [; i( F) r
circle,36,0.1/ T, N- r# N  j/ h+ d% l  E
circle,37,0.1
: G3 }8 x# J) W) Ecircle,38,0.1
4 |0 \+ F) K5 R) u6 |' S. Icircle,39,0.1$ w+ X7 ^! F4 s# @, }
circle,40,0.1* q6 w2 V/ O( v2 V( ]
circle,41,0.1" Y! P3 W. c3 y) H
circle,42,0.1
$ H$ e( N2 u6 U' s; n  ^* Ecircle,43,0.1
3 o) K4 Z6 l! ]$ B, T* E; ecircle,44,0.1* ?" F( S3 v. N
circle,45,0.17 T; G7 }/ p+ k4 ~
circle,46,0.1
6 f8 Z4 I8 R( |' Arectng,-9.92,9.92,-8.66,13.75   !建立矩形6 G+ @7 j1 ^3 P: B1 G1 w+ ]# U
al,1,2,3,4                      !建46根冻结管( B1 F, X5 k! j( [
al,5,6,7,8
2 y3 u1 t) q5 oal,9,10,11,120 I; U3 V# o" L+ C. N
al,13,14,15,16
  X. ~5 i( n4 q& e6 R8 f4 ^' Yal,17,18,19,20  i8 L" b& r+ w2 V5 a4 R2 J, a
al,21,22,23,24/ p3 c0 ~: o' {5 T) T3 G& Z
al,25,26,27,28+ i0 F$ {8 V3 P3 C+ p
al,29,30,31,32
! m' b0 n0 x) Z9 H0 Y- [al,33,34,35,36% Z; c3 C  ^% H5 M1 t: [
al,37,38,39,401 @( c; g* m9 J6 t3 P5 f
al,41,42,43,44
* L7 O( z& r" [; n6 ~8 I5 `al,45,46,47,486 [: j& a9 H: W, b
al,49,50,51,52
  d' u- R8 ~4 K! o  {. X! eal,53,54,55,56! M2 H* }2 W) |) {$ s% N
al,57,58,59,601 M! ?5 j! P, J
al,61,62,63,64
) _7 d$ R! J9 Ial,65,66,67,68
( k& j. k! h5 J% [- R# A: u* Eal,69,70,71,72
2 W6 m  w+ ]9 {( h; d; Cal,73,74,75,76( |7 i9 p7 j! M) u
al,77,78,79,80+ c. w6 i9 A, v& j5 L
al,81,82,83,84
% ~! U! H+ g# O  j% Q/ T. Gal,85,86,87,88
8 l: N- v8 g' n) D% w  e. z. dal,89,90,91,92- m& `1 B1 v0 @" u% S
al,93,94,95,96+ Q- L5 l. p  S- l
al,97,98,99,100
, Q% k1 q! ^& O$ ?4 nal,101,102,103,1042 H5 ~, Y9 g$ m9 I0 L
al,105,106,107,108
! d+ K% p* V: v; z2 Gal,109,110,111,112
$ F5 y( ]1 v7 u; Y4 i' Z4 xal,113,114,115,116% x$ O$ i/ m7 m% E5 R! \
al,117,118,119,120
; e* D# y! L  I5 m/ p. hal,121,122,123,124. r3 p) r; C, A1 K5 U
al,125,126,127,1285 D; j$ z) l3 K8 }, g+ [" d
al,129,130,131,1327 }9 s# s2 N  }/ u" z
al,133,134,135,136
1 ~3 H( p5 {- B4 W& q1 hal,137,138,139,1406 M& |" j# e% u1 \4 k& r
al,141,142,143,144
+ P( q0 c1 S& W+ w, }  @+ C  w& _* g1 Kal,145,146,147,148: x+ t7 o5 M3 f" O1 p# O" y
al,149,150,151,1527 N3 T4 R- `; ^1 p2 A
al,153,154,155,156( i9 K: S% w; ?9 k8 G& l0 w
al,157,158,159,160" k4 Y# H& X+ N
al,161,162,163,1645 R2 w, r' |3 Y2 ^
al,165,166,167,168
# l" k0 Z" v2 h! ral,169,170,171,1724 c8 g& B2 q- U" y2 u
al,173,174,175,176. N. u5 C, P; U9 {2 l! ?
al,177,178,179,180. j9 ~5 r" S4 ^( o! H
al,181,182,183,184
! d, Z& G% l8 ]6 _5 g' u: eaoverlap,all            !冻结管与土体进行布尔运算
! b0 I( X5 i% p) |; i* ?numcmp,area# {  e; y: p. o
numcmp,line3 x  o( J) v9 o" k- M
allsel1 m1 G/ ?+ r# l! e( ^9 ~+ X! `( v
/pnum,area,18 i  O: |% R" P' G) Q# S* _0 T
/pnum,line,0% \% X! G/ ~- t5 P# K, h
aplot
2 L+ I) m8 ?/ I. F3 iasel,s,area,,1,46,1     !选择要设置属性的面5 y: y+ W' `8 Y1 S! q# _
aatt,2,,1
* R2 }- J# x8 i9 y! I# M# F# Tasel,s,area,,47/ U0 g' J5 U1 T5 _1 v0 p1 C9 v
aatt,1,,1
6 d- a4 i) E9 O# Mlsel,s,line,,1,184,1        !选择冻结管面上的线+ s5 _6 ~) s# z' w" M
lesize,all,,,2                   !设置该线等分2份
6 f/ O& B7 v9 e1 x; P" D) Elsel,s,line,,185,187,2      !选择土体边框上的线      
$ _, x& Y3 |% J. A7 |6 ^  W4 [- wlesize,all,,,50+ J# O4 V3 k; G/ `5 D
lsel,s,line,,186,188,2          ) H+ j8 d& u6 M; f6 n
lesize,all,,,80, [8 J9 I* b7 M6 A
allsel
- v+ D2 R- k  Q/ K4 }MSHAPE,0,2D
5 n. {  |0 \( k) S9 kMSHKEY,1
1 z9 e" C6 `/ Y  w2 J4 H!*  
( a& l6 `7 f2 h' V7 QFLST,5,46,5,ORDE,2  
  G( K" z, }0 \2 f% G9 XFITEM,5,1   
+ e  e3 s) A1 ^FITEM,5,-46 - ^: P* {" h+ d5 `' _) A2 e. C/ ~
CM,_Y,AREA  
" }. b6 t% \2 ~ASEL, , , ,P51X
+ F% ^* s  m0 K' q9 V7 N  tCM,_Y1,AREA 5 b/ P6 R" K. P1 D4 p$ s
CHKMSH,'AREA'   $ k! N. e8 @" ]1 U- i
CMSEL,S,_Y  
4 A8 s9 y; h! C) A1 K  d) A, L2 ^- n!*  
+ Z% p7 j, Y  _9 j! E( c; @+ aAMESH,_Y1   
  _: v2 `& a2 M) P, T!*  ) A0 d2 k# K: @) s
CMDELE,_Y   
' D' a2 T9 n2 T0 Q3 JCMDELE,_Y1  " }$ _1 N9 K, E$ n
CMDELE,_Y2  
0 r+ f5 o6 X( e  C# B! PMSHKEY,0
  \9 P5 m1 s+ S+ n/ _, v: s8 B8 D!*  
  s: K' j8 {! W; `" E  sCM,_Y,AREA  
) \& S* |" a5 k" ^$ `ASEL, , , ,      47 4 J( E) V9 Q/ x2 v5 J2 u* v" }
CM,_Y1,AREA 1 N1 e! G  o2 C+ F7 n
CHKMSH,'AREA'   
3 V. p5 L( @% L3 SCMSEL,S,_Y  ! a* i' j* c( I
!*  
5 ~  r# N* M# K1 G  W4 r! X2 ]AMESH,_Y1   
% ^' N2 @5 S& F% Y!*  
7 a+ c+ R% V, R3 z) GCMDELE,_Y   4 w; |- o  c# m- G+ E. ^$ r
CMDELE,_Y1  
5 Z) j& j2 i+ Y2 }1 L$ e, @6 s9 }CMDELE,_Y2  * f0 O4 _, s# i7 F0 W% |
!*  
& D+ \! B8 N# M) X/ m, O1 g# p!!!加载求解阶段
, a" T$ w+ j: e. b1 H/solu8 y6 _5 y; ?, h8 B4 R) x# S
antype,trans
3 N0 h" J/ p# d4 t  x) Ztimint,off       !先瞬态分析: C& H7 G6 e; P4 h; _
time,0.01
. [) @) _& Q: O- U; b5 Mdeltim,0.01
$ f, z& r, p) J5 W2 ^esel,s,mat,,1
3 `# e$ A# }2 r! |/ Insle,s
0 T  y: d4 q$ {$ A) {, x" F& jd,all,temp,23     !设置土体背景值为23°4 S- u% F9 \( W( |+ o  Q
allsel: c  @  [: ]& U# O$ e. t5 R: s1 ^
solve
2 S/ J4 Z0 v  r6 u: z. ~/post1               # \" V/ V. G& _% E1 u1 W; v
set,last             3 m0 J7 B2 G/ e1 t& z" V" G9 l
plnsol,temp4 S2 [6 B0 R4 t
finish
, f  |$ _4 l  z) x
; {5 p  H" F* X& P) B4 Y- I8 G* j/solu                !再稳态分析8 Z8 p$ F  l  |. u+ M& _7 i& w9 ]+ t, T
autots,on. O/ ~% d) h& T- g
outres,,all         
; H) [" P& G' |, ]) X) M1 stimint,on         
2 ]$ r8 |$ b( h: ttime,8640000         !设置终止计算时间
1 T: k8 \$ D/ K; k6 Ddeltim,432000        !设置时间步长
2 O& L% q5 `; y1 i7 b; Ikbc,1                !设置为阶跃加载形式
( k1 ^& W) i- Eesel,s,mat,,1
- m7 D9 X0 V( Knsle,s
/ p3 M6 W: W9 o3 u, [ddelet,all,temp                 
! g) C, O2 B! ]. j- q( r7 J. ?7 [toffst,273   
" h1 q  u; ^0 I4 R' TESEL,S,MAT,,2        ! 加温度载荷为-27°
! R  U( c& I8 \; [2 YNSLE,S  
# z  C5 K( ^+ |( O! w  H9 WD,ALL,TEMP,-276 K% d1 I$ O3 k- h' G. ]$ a
solve7 P% [1 N# f8 o5 r0 P
/post1               !  后处理# Y7 M# K0 d+ W0 L
set,last             !  读取最后一步结果
$ A7 A! [. Y! ~0 C' S! x0 X# kallsel+ v# o% g5 z+ }3 s: b
plnsol,temp
  z" ?0 ^4 O5 L8 s" D( X/ H/ u, w3 wfinish/ c- o) U; G  S- l- H
8 b9 ?. z, r' l+ ^. d6 u
  r9 f8 ~) I( J" ], e) r% b
谢谢了
您需要登录后才可以回帖 登录 | 注册

本版积分规则

QQ|手机版|小黑屋|CadCaeCamPlm专业门户网站 ( 蜀ICP备05008798号 ) 网站速度测试

GMT+8, 2019-8-20 10:21 , Processed in 0.088858 second(s), 20 queries .

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表