网站大量收购独家精品文档,联系QQ:2885784924

工程电磁场数值计算FDM程序(GMC).ppt

  1. 1、本文档共1页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
* READ IN THE PROBLEM FROM FILE GMCIN.DAT READ (NINPT,50) VTOP 50 FORMAT (F6.2) READ (NINPT,55) NODHS, NODVS 55 FORMAT (2I3) * GIVE THE INITIAL VOLTAGE VALUE OF ALL NODES DO 1 I=1,NODHS DO 1 J=1,NODVS VOLTG (I,J) = 0.0 1 CONTINUE DO 2 I=2,NODHS-1 VOLTG (I,NODVS) = VTOP 2 CONTINUE * SUCCESSIVE OVERRELAXATION METHOD(SOR) NTIM=0 100 MUN=0 NTIM=NTIM+1 DO 10 I=2,NODHS-1 DO 10 J=2,NODVS-1 VNOD=VOLTG(I,J) VNOD1=0.25?(VOLTG(I+1,J)+VOLTG(I,J+1)+VOLTG(I-1,J)+ VOLTG(I,J-1)) VNOD=VNOD+ALFA? (VNOD1-VNOD) DEV=ABS(VNOD1-VNOD) IF (DEV.GE.DF) MUN=MUN+1 VOLTG(I,J)=VNOD 10 CONTINUE ( 7行 ) 程序说明行 ( 8-11行 ) 读入: ? VTOP — 顶盖电位值(10.0) ? NODHS, NODVS — 水平、垂直方向的节点数(5, 5) ( 12行 ) 程序说明行 ( 13-16行 ) 循环赋值语句 φ (i,,j) 赋值 : ? 各内节点的零初始值 ? 侧壁与底面的给定边值 ( 17-19行 ) 循环赋值语句 顶盖给定边值 φ (i,,5) ( 20行 ) 程序说明行 ( 21行 ) 迭代次数(NTIM)赋零值 ( 22行 ) 不满足收敛指标的离散点数(MUN) ( 23行 ) 迭代次数 + 1 (24-32行 ) DO 10 循环体 一次迭代计算过程

文档评论(0)

132****9295 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档