FEAP 子程序及功能PPT格式课件下载.ppt
- 文档编号:13885083
- 上传时间:2022-10-14
- 格式:PPT
- 页数:21
- 大小:47KB
FEAP 子程序及功能PPT格式课件下载.ppt
《FEAP 子程序及功能PPT格式课件下载.ppt》由会员分享,可在线阅读,更多相关《FEAP 子程序及功能PPT格式课件下载.ppt(21页珍藏版)》请在冰豆网上搜索。
FEAP子程序及功能,MASTERMINIFEM-主程序,设置一维数组大小,调用子程序PCONTRSUBROUTINEPCONTR-控制数据输入与求解计算流程SUBROUTINESETMEM-检查一维数组大小是否够用SUBROUTINEZERO-将数组置零,SUBROUTINEPMESH-数据输入SUBROUTINEGENVEC-网格自动生成(通过线性插值)SUBROUTINEPROFIL-计算总体刚度矩阵元素的定位(去掉约束,带宽,主对角线元素位置)SUBROUTINEPMACR-引导与控制求解计算,SUBROUTINEELMLIB-单元库(计算单元刚度矩阵)SUBROUTINEADDSTF-集合总体刚度矩阵SUBROUTINEPLOAD-形成荷载列阵FUNCTIONPROP-计算当前时间的荷载乘子SUBROUTINEPFORM-计算单元刚度矩阵和和荷载列阵,并组集成总体矩阵,SUBROUTINEACTCOL-对称线性方程组三角分解SUBROUTINEUACTCL-非对称线性方程组三角分解SUBROUTINEPEIGS-计算特征值(动力问题)SUBROUTINEEUPDAT-更新速度和加速度(动力问题)SUBROUTINEPRTDIS-输入节点位移,SUBROUTINEPRTREA-输出节点反力SUBROUTINEPROMUL-矩阵乘C=C+A*BFUNCTIONFUNCTIONDOT-两个矢量相乘SUBROUTINENORM-将矢量单位化LOGICALFUNCTIONPCOMP-比较字符变量,单元库,SUBROUTINEELMT01-平面线弹性单元库SUBROUTINESHAPE-计算形函数及其导数SUBROUTINEPGAUSS-给出GAUSS点坐标及积分权函数SUBROUTINESHAP2-增加四节点等参单元插值点(增加中点),FEAP*QUADRANTOFACIRCUIARDISK(EXAMPLEPROBLEM)19111224COOR110.0.5050.610.2.1004.58282.1110.4.1403.04.1504.03.1600.5.1700.754.94341801.54.76971902.254.4651,ELEM111276151671211191111217161BOUN111-1500165-1016010FORC1600.-5.MATE110.30.0211END,线性求解,MACRTANGFORMSOLVSTREDISPSTOP,非线性求解,LOOP10TANGFORMSOLVDISPNEXTSTREDISPSTOP,TANGLOOP10FORMSOLVDISPNEXTSTREDISPSTOP,刚度保持不变,DTPROPLOOP10TIMEFORMSOLVDISPNEXTEND,增量荷载法,MASTERMINIFFMC-SETPROGRAMCAPACITYMAXMUSTAGREEC-WITHDIMENSIONOFMCOMMONM(2000)COMMON/PSIZE/MAXMAX=2000CALLPCONTRSTOPEND,LOGICALPCOMPCOMMON/CDATA/O,HEAD(20),NUMNP,NUMEL,1NUMMAT,NEN,NEQ,IPRCOMMON/LABEL/PDIS(6),A(6),BC
(2),DI(6),1CD(3),TE(3),FD(3)COMMONM
(1)DIMENSIONTITL
(2),WD(3)DATAWE/4HFEAP,4HMACR,4HSTOP/C-READACARDANDCOMPAREFIRST4COLUMNSC-WITHMACRDLIST1READ(5,1000)TITL,IF(PCOMP(TITL
(1),WD
(1)GOTO100IF(PCOMP(TITL
(1),WD
(2)GOTO200IF(PCOMP(TITL
(1),WD(3)RETURNGOTO1C.READANDPRINTCONTROLINFORMATIONDO1011=1,20HEAD
(1)=TITL
(1)READ(5,1001)NUMNP,NUMEL,NUMMAT,NDM,NDF,NEN,1NADWRITE(6,2000)HEAD,NUMNP,NUMEL,NUMMAT,NDM,1NDF,NEN,NAD,C.SETPOINTERSFORALLOCATIONOFCATAARRAYSPDIS
(2)=A(NDM)NEN1=NEN+1NST=NEN*NDF+NADN0=1+NST*2*IPRN1=N0+NEN*NDM*IPRN2=N1+NEN*IPRN3=N2+NSTN4=N3+NST*IPRN5=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.CHECKTHATSUFFICIENTMEMORYEXISTSCALLSETMEM(N13)CALLPZERO(M,N12),C-CALLMESHINPUTSUBHOUTINETOREADANDPRINTALLC-MESHDATAIII=0CALLPMESH(M(N2),M(N5),M(N6),M(N7),M(N8),M(N9),1M(N10),M(N11),NDF,NDM,NEN1,III)C-ESTATLISHPROFILEOFRFSULTINGEQUATIONSFORC-STIFFNESS,MASS,ETCCALLPROFIL(M(N12),M(N7),M(N9),NDF,NEN1,NAD)C-SETPOINTERSFORSOLUTION,ARRAYS,CHECKFORC-SUFFICIENTMEMORYN13=N12+NEQN14=N13+NEQ*IPR,NE=N14+NUMNP+NDF+IPRCALLSETNEM(NE)GOTO1C-CALLMACROSOLUTIONMODULEFORESTABLISHINGC-SOLUTIONALGOTITHMCALLPMACR(M,M(N0),M(N1),M(N2),M(N3),M(N4),M(N5),1M(N6),M(N7),M(N8),M(N9),M(N10),M(N11),M(N12),1M(N13),M(N14),M(NE),NDF,NDM,MENT,NST,NE)GOTO1,CALLPMACR(M,M(N0),M(N1),M(N2),M(N3),M(N4),M(N5),1M(N6),M(N7),M(N8),M(N9),M(N10),M(N11),M(N12),1M(N13),M(N14),M(NE),NDF,NDM,MENT,NST,NE)PMACR(UL,XL,TL,LD,P,S,IE,1D,ID,X,IX,F,T,JDIAG,1B,DR,CT,NDF,NDM,MENT,NST,NE),
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- FEAP 子程序及功能 子程序 功能