| ²é¿´: 1157 | »Ø¸´: 4 | ||
lijunjiexmÌú³æ (СÓÐÃûÆø)
|
[ÇóÖú]
ÇóÖúÒ»¸ö±àÒë³ö´íµÄ¼òµ¥³ÌÐò£¬Çë¶à¶àÖ¸½Ì£¬Ð»Ð»
|
|
ÔÚ½¨Á¢µÄÒ»¸öprojectµÄSource filesµÄÆäÖÐÃüÃûΪSTRESSAVGµÄfÎļþÈçÏ£º C C********************************************************************** C********************************************************************** C********************************************************************** C SUBROUTINE STRESSAVG(NUMEL, STORE, I1, I2, STR, NKT, > ORDWHERE, FIRST, LAST, > COUNTPLUS, COUNTMINS, COUNTFLAG) IMPLICIT NONE INCLUDE 'PARAMETERS.INC' INCLUDE 'NUMBERS.INC' INTEGER*4 ORDWHERE(NDMB) INTEGER*4 I INTEGER*4 FIRST (NUMEL + 1), LAST (NUMEL + 1) REAL*8 STORE (NDMB) REAL*8 SUM1, SUM2, NKT REAL*8 STR(2, 3, NUMEL) INTEGER*4 II, J, I1, I2, INDEX, NUMEL INTEGER*4 COUNTPLUS(NUMEL), COUNTMINS(NUMEL) LOGICAL*1 COUNTFLAG C C---------------------------------------------------------------------- C STR (i , j , k) C C ^ ^ ^ C | | |_________ element number C | |_____________ components in order xx, yx, yy C |_________________ triangle number (1 or 2) C C---------------------------------------------------------------------- C IF (COUNTFLAG) THEN DO I = 1 , NUMEL COUNTPLUS(I) = 0 COUNTMINS(I) = 0 ENDDO ENDIF IF ((I1 .EQ. 1) .AND. (I2 .EQ. 1)) THEN INDEX = 1 ELSEIF ((I1 .EQ. 1) .AND. (I2 .EQ. 2)) THEN INDEX = 2 ELSEIF ((I1 .EQ. 2) .AND. (I2 .EQ. 2)) THEN INDEX = 3 ELSE PRINT *,'Wrong indices for stress tensor' ENDIF DO J = 1 , NUMEL SUM1 = ZERO SUM2 = ZERO *vdir nodep CDIR$ IVDEP DO II = FIRST(J) , LAST(J) I = ORDWHERE(II) IF (I .GT. 0) THEN SUM1 = SUM1 + STORE(I) IF (COUNTFLAG) THEN COUNTPLUS(J) = COUNTPLUS(J) + 1 ENDIF ELSE SUM2 = SUM2 + STORE(-I) IF (COUNTFLAG) THEN COUNTMINS(J) = COUNTMINS(J) + 1 ENDIF ENDIF ENDDO STR(1 , INDEX , J) = SUM1 STR(2 , INDEX , J) = SUM2 ENDDO DO J = 1 , NUMEL STR(1 , INDEX , J) = STR(1 , INDEX , J) * NKT / > MAX(1 , COUNTPLUS(J)) STR(2 , INDEX , J) = STR(2 , INDEX , J) * NKT / > MAX(1 , COUNTMINS(J)) ENDDO IF ((INDEX .EQ. 1) .OR. (INDEX .EQ. 3)) THEN DO J = 1 , NUMEL STR(1 , INDEX , J) = STR(1 , INDEX , J) - NKT STR(2 , INDEX , J) = STR(2 , INDEX , J) - NKT ENDDO ENDIF RETURN END ¶øPARAMETERS.INCºÍNUMBERS.INCͬʱ±£´æÔÚͬһ¸öprojectϵÄExtenal DependenciecÎļþ¼ÐÖУº PARAMETERS.INCÎļþΪ£º INTEGER*4 NSD, NEN, NDOF INTEGER*4 NNPMAX, NELMAX, NEQMAX INTEGER*4 NSBDIV, NBWMAX INTEGER*4 NBUF , MAXTRACE INTEGER*4 NDMB C PARAMETER ( NSD = 2, > NEN = 4, > NDOF = 2, > NNPMAX = 3000, > NEQMAX = NDOF * NNPMAX, > NSBDIV = 110, > NBWMAX = NDOF * NSBDIV, > NDMB = 524288, C > NDMB = 65536, C > NDMB = 2048, > NBUF = 1024, > NELMAX = 3000, > MAXTRACE = 10000 ) NUMBERS.INC ÎļþΪ£º REAL*8 ZERO, ONE, TWO, THREE, FOUR, FIVE, SIX, > HALF, THIRD, SIXTH, > FOURTH, TWSVTH, NINTH, TWFOUR, > PI, TWOPI, FTPI, > RSMALL PARAMETER (ZERO = 0.0 , > ONE = 1.0 , > TWO = 2.0 , > THREE = 3.0 , > FOUR = 4.0 , > FIVE = 5.0 , > SIX = 6.0 , > HALF = 1.0 / 2.0 , > THIRD = 1.0 / 3.0 , > SIXTH = 1.0 / 6.0 , > FOURTH = 1.0 / 4.0 , > TWSVTH = 1.0 / 27.0 , > NINTH = 1.0 / 9.0 , > TWFOUR = 24. , > PI = 3.14159265358979323846 , > TWOPI = 2. * 3.14159265358979323846 , > FTPI = 4. * 3.14159265358979323846 / 3. , > RSMALL = 1.0E-30 ) |
» ²ÂÄãϲ»¶
Ò»Ö¾Ô¸Ö£ÖÝ´óѧ085600Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
Çóµ÷¼Á Ò»Ö¾Ô¸Î÷ÄϽ»Í¨´óѧ085701»·¾³¹¤³Ì 282·Ö
ÒѾÓÐ7È˻ظ´
Çóµ÷¼Á
ÒѾÓÐ15È˻ظ´
²ÄÁÏר˶283Çóµ÷¼Á
ÒѾÓÐ15È˻ظ´
£¨µ÷¼Á£©Ò»Ö¾Ô¸±¨¿¼¹þ¶û±õ¹¤Òµ´óѧ0857×ÊÔ´Óë»·¾³×¨Òµ378·Ö¿¼Éú
ÒѾÓÐ3È˻ظ´
±¾¿Æ211£¬293·ÖÇëÇóµ÷¼Á
ÒѾÓÐ9È˻ظ´
315Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
308Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
301Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
304Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
ÇóÖúÒ»¼òµ¥matlab³ÌÐòµÄ³ö´íÎÊÌ⣬лл
ÒѾÓÐ3È˻ظ´
Çë½ÌÒ»¸öfortranС³ÌÐò±àÒë³ö´íµÄÎÊÌ⣬лл
ÒѾÓÐ9È˻ظ´
¡¾ÇóÖú¡¿Çë½ÌÒ»¸ömatlab³ÌÐò£¬Ð»Ð»°ï棡£¡
ÒѾÓÐ15È˻ظ´
¡¾ÇóÖú¡¿VASP5.2 ±àÒë³ö´ímake: *** [vasp] ´íÎó 1,Çë´ó¼Ò°ïæ·ÖÎö·ÖÎö,лл
ÒѾÓÐ7È˻ظ´
¡¾ÇóÖú¡¿f90×Ó³ÌÐòÔõô±àÒë´íÎó£¿Çë´ó¼Ò°ï棬лл
ÒѾÓÐ3È˻ظ´
¡¾ÇóÖú¡¿Çë½Ì±àÒë¹ý³ÌÖеÄÁ½¸öÎÊÌ⡾ÒÑÍê½á¡¿
ÒѾÓÐ7È˻ظ´
¡¾ÇóÖú¡¿ÎªÊ²Ã´ÎÒ±àÒëÁËÒ»¸öÊé±¾ÉϵÄÎÞ´íÎóµÄC++³ÌÐò£¬È´ÓÐÕâô¶àµÄ´íÎó£¿
ÒѾÓÐ8È˻ظ´
¡¾ÇóÖú¡¿fortran³ö´í ±àÒëͨ¹ý ÔËÐÐÒç³ö
ÒѾÓÐ10È˻ظ´
¡¾ÇóÖú¡¿VASP±àÒë³ö´í
ÒѾÓÐ5È˻ظ´
lijunjiexm
Ìú³æ (СÓÐÃûÆø)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 22.8
- ºì»¨: 1
- Ìû×Ó: 182
- ÔÚÏß: 87.1Сʱ
- ³æºÅ: 1198843
- ×¢²á: 2011-02-03
- רҵ: ³£Î¢·Ö·½³ÌÓ붯Á¦ÏµÍ³
¡ï
dubo(½ð±Ò+1): »¶Ó³£À´³ÌÐòÓïÑÔ°æÌÖÂÛ 2011-07-31 13:31:08
dubo(½ð±Ò+1): »¶Ó³£À´³ÌÐòÓïÑÔ°æÌÖÂÛ 2011-07-31 13:31:08
|
ÔÚwindowsÏÂCompaq Visual Fotran±àÒë³öÏÖÈçÏ´íÎóÐÅÏ¢£¨Çë¶à¶àÖ¸½Ì£¬Ð»Ð»£©£º -------------------Configuration: sucontr - Win32 Debug-------------------- Compiling Fortran... G:\Program Files\sucontr\src\stressavg.f G:\Program Files\sucontr\src\stressavg.f(14) : Error: This name does not have a type, and must have an explicit type. [NUMEL] INTEGER*4 FIRST (NUMEL + 1), LAST (NUMEL + 1) ------------------------^ G:\Program Files\sucontr\src\stressavg.f(14) : Error: A specification expression is invalid. [NUMEL] INTEGER*4 FIRST (NUMEL + 1), LAST (NUMEL + 1) ------------------------^ G:\Program Files\sucontr\src\stressavg.f(14) : Error: A specification expression is invalid. [NUMEL] INTEGER*4 FIRST (NUMEL + 1), LAST (NUMEL + 1) ------------------------------------------^ G:\Program Files\sucontr\src\stressavg.f(17) : Error: A specification expression is invalid. [NUMEL] REAL*8 STR(2, 3, NUMEL) ---------------------------^ Error executing df.exe. stressavg.obj - 4 error(s), 0 warning(s) |
2Â¥2011-07-17 21:58:30
snoopyzhao
ÖÁ×ðľ³æ (Ö°Òµ×÷¼Ò)
- ³ÌÐòÇ¿Ìû: 16
- Ó¦Öú: 157 (¸ßÖÐÉú)
- ¹ó±ö: 0.02
- ½ð±Ò: 18844.7
- ºì»¨: 29
- Ìû×Ó: 3803
- ÔÚÏß: 1422.4Сʱ
- ³æºÅ: 183750
- ×¢²á: 2006-02-13
- רҵ: ÎÛȾÉú̬»¯Ñ§
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
¡ï
jjdg(½ð±Ò+1): ¸ÐлӦÖú 2011-07-17 22:18:08
lijunjiexm(½ð±Ò+2): ¸ÄÁË£¬»¹ÊDz»ÐС£¡£ 2011-07-17 23:14:39
jjdg(½ð±Ò+1): ¸ÐлӦÖú 2011-07-17 22:18:08
lijunjiexm(½ð±Ò+2): ¸ÄÁË£¬»¹ÊDz»ÐС£¡£ 2011-07-17 23:14:39
|
ÄãÔÚINCLUDE 'NUMBERS.INC'ºó¼ÓÒ»¾ä£º INTEGER*4 NUMEL Ó¦¸Ã¾Í¿ÉÒÔÁË¡¡ |
3Â¥2011-07-17 22:12:46
snoopyzhao
ÖÁ×ðľ³æ (Ö°Òµ×÷¼Ò)
- ³ÌÐòÇ¿Ìû: 16
- Ó¦Öú: 157 (¸ßÖÐÉú)
- ¹ó±ö: 0.02
- ½ð±Ò: 18844.7
- ºì»¨: 29
- Ìû×Ó: 3803
- ÔÚÏß: 1422.4Сʱ
- ³æºÅ: 183750
- ×¢²á: 2006-02-13
- רҵ: ÎÛȾÉú̬»¯Ñ§
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
¡ï
jjdg(½ð±Ò+1): ¸Ðл²ÎÓë 2011-07-18 14:51:34
jjdg(½ð±Ò+1): ¸Ðл²ÎÓë 2011-07-18 14:51:34
|
¿´³ö´íÐÅÏ¢£º G:\Program Files\sucontr\src\stressavg.f(14) : Error: This name does not have a type, and must have an explicit type. [NUMEL] ˵µÄ¾ÍÊÇÎÒÇ°Ãæ¸ø³öµÄÄǸö½â¾ö·½°¸¡£ ²»¹ý£¬ÎÒ¸Õ²ÅÓÃgfortranÊÔÁËһϣ¬¼´Ê¹²»¼ÓÇ°ÃæµÄÄÇÒ»¾ä»°£¬¶¼¿ÉÒÔ˳Àû±àÒë°¡¡¡ |
4Â¥2011-07-17 23:32:05
snoopyzhao
ÖÁ×ðľ³æ (Ö°Òµ×÷¼Ò)
- ³ÌÐòÇ¿Ìû: 16
- Ó¦Öú: 157 (¸ßÖÐÉú)
- ¹ó±ö: 0.02
- ½ð±Ò: 18844.7
- ºì»¨: 29
- Ìû×Ó: 3803
- ÔÚÏß: 1422.4Сʱ
- ³æºÅ: 183750
- ×¢²á: 2006-02-13
- רҵ: ÎÛȾÉú̬»¯Ñ§
¡¾´ð°¸¡¿Ó¦Öú»ØÌû
¡ï ¡ï
jjdg(½ð±Ò+2): ¸Ðл²ÎÓë 2011-07-18 14:51:45
lijunjiexm(½ð±Ò+28): Thanks.Ï£ÍûÒԺ󻹿ÉÒÔÏòÄúÌÖ½ÌFortranµÄÎÊÌâ¡£ 2011-07-19 08:15:56
jjdg(½ð±Ò+2): ¸Ðл²ÎÓë 2011-07-18 14:51:45
lijunjiexm(½ð±Ò+28): Thanks.Ï£ÍûÒԺ󻹿ÉÒÔÏòÄúÌÖ½ÌFortranµÄÎÊÌâ¡£ 2011-07-19 08:15:56
|
Ŷ£¬¿´µ½ÎÊÌâÁË£¬ÄãÔÚµÚ 19 ÐÐÉùÃûÁË NUMEL£¬µ«ÄãÔÚµÚ 14 ÐÐÓÃÁË NUMEL ËùÒÔ£¬Äãɾµô 19 ÐÐÖÐµÄ NUMEL£¬È»ºó°´ÎÒµÚÒ»ÌûÖеķ½·¨ÔÙÊÔһϡ¡ |
5Â¥2011-07-17 23:35:19














»Ø¸´´ËÂ¥