彈塑性力學(xué)講義-FEAP 子程序及功能

1、FEAP 子程序及功能MASTER MINIFEM -主程序，設(shè)置一維數(shù)組大小，調(diào)用子程序PCONTRSUBROUTINE PCONTR - 控制數(shù)據(jù)輸入與求解計(jì)算流程SUBROUTINE SETMEM - 檢查一維數(shù)組大小是否夠用SUBROUTINE ZERO - 將數(shù)組置零SUBROUTINE PMESH - 數(shù)據(jù)輸入SUBROUTINE GENVEC -網(wǎng)格自動(dòng)生成（通過線性插值）SUBROUTINE PROFIL - 計(jì)算總體剛度矩陣元素的定位 （去掉約束，帶寬，主對(duì)角線元素位置）SUBROUTINE PMACR - 引導(dǎo)與控制求解計(jì)算SUBROUTINE ELMLIB - 單元庫（計(jì)

2、算單元?jiǎng)偠染仃嚕㏒UBROUTINE ADDSTF - 集合總體剛度矩陣SUBROUTINE PLOAD - 形成荷載列陣FUNCTION PROP - 計(jì)算當(dāng)前時(shí)間的荷載乘子SUBROUTINE PFORM -計(jì)算單元?jiǎng)偠染仃嚭秃秃奢d列陣，并組集成總體矩陣SUBROUTINE ACTCOL - 對(duì)稱線性方程組三角分解SUBROUTINE UACTCL - 非對(duì)稱線性方程組三角分解SUBROUTINE PEIGS - 計(jì)算特征值（動(dòng)力問題）SUBROUTINE EUPDAT - 更新速度和加速度（動(dòng)力問題）SUBROUTINE PRTDIS - 輸入節(jié)點(diǎn)位移SUBROUTINE PRTREA

3、-輸出節(jié)點(diǎn)反力SUBROUTINE PROMUL - 矩陣乘 C = C + A*B FUNCTION FUNCTION DOT - 兩個(gè)矢量相乘SUBROUTINE NORM - 將矢量單位化LOGICAL FUNCTION PCOMP -比較字符變量單元庫SUBROUTINE ELMT01 -平面線彈性單元庫SUBROUTINE SHAPE - 計(jì)算形函數(shù)及其導(dǎo)數(shù)SUBROUTINE PGAUSS - 給出GAUSS點(diǎn)坐標(biāo)及積分權(quán)函數(shù)SUBROUTINE SHAP2 - 增加四節(jié)點(diǎn)等參單元插值點(diǎn)（增加中點(diǎn)）FEAP*QUADRANT OF A CIRCUIAR DISK （EXAMPLE

4、PROBLEM）19 11 1 2 2 4COOR11 0. 0.50 50.61 0.2.10 0 4.58282.11 1 0. 4.14 0 3.04.15 0 4.03.16 0 0.5.17 0 0.754.943418 0 1.54.769719 0 2.254.4651ELEM11 1 2 76151 6 7 12 11191 11 12 17 161BOUN11 1 -150 0 165 -1 0160 1 0FORC160 0. -5.MATE11 0.3 0.0 2 11END線性求解MACRTANGFORMSOLVSTREDISPSTOP非線性求解LOOP 10TANGF

5、ORMSOLVDISPNEXTSTRE DISPSTOPTANG LOOP 10FORMSOLVDISPNEXTSTRE DISPSTOP剛度保持不變DTPROPLOOP 10TIME FORMSOLVDISPNEXTEND增量荷載法 MASTER MINIFFM C-SET PROGRAM CAPACITY MAX MUST AGREEC-WITH DIMENSION OF M COMMON M（2000）COMMON/ PSIZE/ MAX MAX=2000CALL PCONTRSTOPEND LOGICAL PCOMPCOMMON/CDATA/ O,HEAD(20), NUMNP, NU

6、MEL, 1 NUMMAT, NEN, NEQ, IPRCOMMON/ LABEL/ PDIS(6), A(6), BC(2), DI(6), 1 CD(3), TE(3), FD(3)COMMON M(1)DIMENSION TITL (2), WD(3)DATA WE/ 4HFEAP, 4HMACR, 4HSTOP/C-READ A CARD AND COMPARE FIRST 4 COLUMNSC-WITH MACRD LIST 1 READ (5,1000) TITLIF (PCOMP (TITL(1),WD(1) GO TO 100IF (PCOMP (TITL(1),WD(2) G

7、O TO 200IF (PCOMP (TITL(1),WD(3) RETURNGO TO 1C.READ AND PRINT CONTROL INFORMATION DO 101 1=1，20 HEAD (1)=TITL(1) READ(5, 1001) NUMNP, NUMEL, NUMMAT, NDM, NDF, NEN, 1 NAD WRITE(6, 2000) HEAD, NUMNP, NUMEL, NUMMAT, NDM, 1 NDF, NEN, NAD C. SET POINTERS FOR ALLOCATION OF CATA ARRAYSPDIS (2)= A(NDM)NEN1

8、= NEN+1NST = NEN*NDF+NADN0 = 1+NST*2*IPRN1 = N0+NEN*NDM*IPRN2 = N1+NEN*IPRN3 = N2+NSTN4 = N3+NST*IPR N5 = N4+NST*NST*IPRN6=N5+NUMMAT N7=N6+10+NUMMAT*IPRN8=N7+NDF*NUMNPN9=N8+NDM*NUMNP*IPRN10=N9+NEN1*NUMELN11=N10+NDF*NUMNP*IPRN12=N11+NUMNP*IPRN13=N12+NDF*NUMNPC. CHECK THAT SUFFICIENT MEMORY EXISTSCALL

9、 SETMEM (N13) CALL PZERO (M, N12) C- CALL MESH INPUT SUBHOUTINE TO READ AND PRINT ALL C- MESH DATA III=0 CALL PMESH (M(N2), M(N5),M(N6), M(N7) , M(N8) , M(N9) , 1 M(N10) , M(N11), NDF, NDM, NEN1, III)C- ESTATLISH PROFILE OF RFSULTING EQUATIONS FOR C- STIFFNESS, MASS, ETC CALL PROFIL (M(N12), M(N7),

10、M(N9), NDF, NEN1, NAD) C- SET POINTERS FOR SOLUTION, ARRAYS, CHECK FOR C-SUFFICIENT MEMORY N13=N12+NEQ N14=N13+NEQ*IPRNE=N14+NUMNP+NDF+IPRCALL SETNEM(NE) GO TO 1C-CALL MACRO SOLUTION MODULE FOR ESTABLISHING C-SOLUTION ALGOTITHM CALL PMACR (M, M(N0), M(N1), M(N2), M(N3), M(N4), M(N5), 1 M(N6), M(N7), M(N8), M(N9), M(N10), M(N11), M(N12), 1 M(N13), M(N14), M(NE),NDF, NDM, MENT, NST, NE)GO TO 1 CALL PMACR (M, M(N0), M(N1), M(N2), M(N3), M(N4), M(N5), 1 M(N6), M(N7), M(N8), M(N9), M(N10