| 查看: 1581 | 回复: 12 | ||
| 【悬赏金币】回答本帖问题，作者阳关不锈将赠送您 15 个金币 | ||
[求助]
【求助】使用fortran语言编写umat子程序想实现循环积分，但一直报错，求教
|
||
|
if ((elt_m .ge. 1.0).and.(sigma_N.ge.0)) then !---------matrix tensile cracking--------- Gnc=Gnct ftype=1.0 T0_eff=T_eff/sqrt(elt_m) D0_eff=D_eff/sqrt(elt_m) Gc=Gnc+(Gsc-Gnc)*(Gs/(Gn+Gs))**eta Df_eff=D0_eff Df_eff=statev(6) E_mn=0.5*D0_eff*T0_eff do while(E_mn.lt.Gc) Df_eff=Df_eff+0.001 statev(6)=Df_eff E_mn=E_mn+0.001*exp(T0_eff*D0_eff/(Gc-0.5*T0_eff*D0_eff)*(1.0-Df_eff/D0_eff))*T0_eff end do Df_eff=statev(6)+0.001 if(Df_eff.le.D0_eff) then dm=dgdf2 ftype=1.5 endif 这是包含循环部分的程序，请帮我看看do while部分有没有错误，谢谢了！ |
回复此楼» 猜你喜欢
南通大学学科教学（物理）接受调剂
已经有0人回复
-
香港科技大学全奖招收博士/博士后，光子、纳米材料、液晶显示方向
已经有15人回复
-
物理学I论文润色/翻译怎么收费?
已经有152人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
超导三元凸包图
已经有0人回复
-
最近遇到这样一个问题
已经有3人回复
2楼2017-08-16 10:49:49
BeverlyNWPU
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 860.6
- 帖子: 76
- 在线: 40.7小时
- 虫号: 3707887
- 注册: 2015-03-03
- 性别: MM
- 专业: 流体力学
3楼2017-08-17 17:33:01
lijf517
木虫 (小有名气)
- 应助: 39 (小学生)
- 金币: 3391.9
- 红花: 14
- 帖子: 293
- 在线: 390.9小时
- 虫号: 2249203
- 注册: 2013-01-18
- 性别: GG
- 专业: 力学
4楼2017-08-17 17:55:38
5楼2017-08-17 19:45:48
6楼2017-08-17 19:48:26
lijf517
木虫 (小有名气)
- 应助: 39 (小学生)
- 金币: 3391.9
- 红花: 14
- 帖子: 293
- 在线: 390.9小时
- 虫号: 2249203
- 注册: 2013-01-18
- 性别: GG
- 专业: 力学
|
编译器报的是个语法错误，与具体算法没关系。从贴出代码段看，第一行的if一句应该在最后一行加个endif语句来闭合。do while语句没有明显错误。但要注意该循环语句段中有一行特别长，如果编译器选项没设置好的话，这一代码行可能只有前面一百多个字符被识别，后面的会自动被编译器丢弃，也可能会造成语法错误。最后，你的循环语句需要优化下。那个长语句行中，T0_eff和D0_eff的值在循环中一直是固定值，所以请把与这两项相关的乘法除法都移到循环外面计算，不要再循环内一遍又一遍计算已知的值，这是无用功，会降低计算效率 发自小木虫Android客户端 |
7楼2017-08-18 10:01:50
8楼2017-08-18 14:20:37
 |
9楼2017-08-19 12:30:47
 |
10楼2017-08-19 16:49:23
1(金币+5)
