有关ground的问题

ziyifan1983

[这个贴子最后由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不知道为什么   

让大家费神了，刚开始学很头大，希望得到你们帮助谢谢