|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
[这个贴子最后由ziyifan1983在 2007/05/24 10:07am 第 1 次编辑]
大家好,我是PHOENICS的初学者,现在希望能在PLANT里编程序,有很多问题不明白,希望得到你们的帮助!
1 我想输入一种耐火材料的的比热和密度,但是单单输入比热时程序可以运行,但是再输入密度后就不能运行了(键入回车后屏幕闪一下而已)我传了Q1和ground文件麻烦大家看一下
2 听人说PLANT里面可以输入的句子比较少,可以在ground里面直接输入,可是我在ground里直接输入密度和比热时,运行不成功,而且之前在ground里写的句子也不见了,不知道为什么,是不是要修改什么设置才能保留在ground里写的句子
3 我想解四元一次方程组请问在PHOENICS里面有没有现成的程序调用
谢谢大家的帮助
奇怪为什么传不上呢?我复制一点大家帮我看看吧
GROUND
GROUP 9. Properties of the medium (or media)
C
C The sections in this group are arranged sequentially in their
C order of calling from EARTH. Thus, as can be seen from below,
C the temperature sections (10 and 11) precede the density
C sections (1 and 3); so, density formulae can refer to
C temperature stores already set.
9 GO TO (91,92,93,94,95,96,97,98,99,900,901,902,903,904,905),ISC
C*****************************************************************
900 CONTINUE
C * ------------------- SECTION 10 ---------------------------
C For TMP1.LE.GRND--------- phase-1 temperature Index TEMP1
RETURN
901 CONTINUE
C * ------------------- SECTION 11 ---------------------------
C For TMP2.LE.GRND--------- phase-2 temperature Index TEMP2
RETURN
902 CONTINUE
C * ------------------- SECTION 12 ---------------------------
C For EL1.LE.GRND--------- phase-1 length scale Index LEN1
RETURN
903 CONTINUE
C * ------------------- SECTION 13 ---------------------------
C For EL2.LE.GRND--------- phase-2 length scale Index LEN2
RETURN
904 CONTINUE
C * ------------------- SECTION 14 ---------------------------
C For SOLVE(TEM1)-------- phase-1 specific heat
C Property name: PRPT01
IF(ISTEP.GE.1 .AND.ISTEP.LE.LSTEP ) THEN
IF(IZ.GE.1 .AND.IZ.LE.NZ ) THEN
LFISPH =L0F(ISPH1 )
DO 91401 IX=IXF ,42
IADD=NY*(IX-1)
DO 91401 IY=14 ,IYL
I=IY+IADD
L0ISPH=LFISPH+I
91401 F(L0ISPH )=1
ENDIF
ENDIF
RETURN
905 CONTINUE
C * ------------------- SECTION 15 ---------------------------
C For SOLVE(TEM2)-------- phase-2 specific heat
RETURN
91 CONTINUE
C * ------------------- SECTION 1 ---------------------------
C For RHO1.LE.GRND--- density for phase 1 Index DEN1
C Property name: PRPT02
IF(ISTEP.GE.1 .AND.ISTEP.LE.LSTEP ) THEN
IF(IZ.GE.1 .AND.IZ.LE.NZ ) THEN
LFDEN1 =L0F(AUX(DEN1 ))
DO 90102 IX=IXF ,42
IADD=NY*(IX-1)
DO 90102 IY=14 ,IYL
I=IY+IADD
L0DEN1=LFDEN1+I
90102 F(L0DEN1 )=2.0
ENDIF
ENDIF
RETURN
92 CONTINUE
C * ------------------- SECTION 2 ---------------------------
C For DRH1DP.LE.GRND--- D(LN(DEN))/DP for phase 1
C Index D1DP
RETURN
Q1中
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
* Non-default variable names
NAME(149) =SPH1 ; NAME(150) =TEM1
* Solved variables list
SOLVE(TEM1)
* Stored variables list
STORE(SPH1)
* Additional solver options
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
TERMS (TEM1,Y,Y,Y,Y,Y,Y)
NEWRH1 = T
************************************************************
Group 9. Properties
PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02
SETPRPS(1,115)
RHO1 = GRND
ENUT = 0.000000E+00
DRH1DP = 1.000000E-09
PRNDTL(TEM1) = -8.400000E-01
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Echo PLANT settings
PLANTBEGIN
NAMSAT=MOSG
<RPT01> ISPH1=1
REGION(,42,14,,,,,)
<RPT02> DEN1=2.0
REGION(,42,14,,,,,)
PLANTEND
93 CONTINUE
在GROUP7中不出现密度DEN1不知道为什么
让大家费神了,刚开始学很头大,希望得到你们帮助谢谢
|
|