声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1233|回复: 2

[Fluent应用] 关于fidap的subroutine的编写 请教

[复制链接]
发表于 2007-7-20 10:11 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
以下是本人写的一个subroutine 是想计算一个依存于表面流体流速而变化的htransfer值 但是下面这个subr应该有很多的错误  请高手指教 谢谢 流速和htra的关系式是现成的htra=VEL**0.8 还有就是怎么把这个SUBR代回到主程序里呢 比如主程序的
写的是htra(set=1,subroutine) 好像也不对吧  很急
SUBROUTINE USRGAP (NELT,NE,NG,COEF,VARI,DVARI,NDFCD,LDOFU,SHP,DSDX,XYZL,PROP,TIME,NPTS,NDP,MNDP,IERR,ZIGN,ZIGT,GAP,TEMP2)
include "IMPLCT.COM"
include "PARUSR.COM"
      DIMENSION COEF(NPTS), GAP(NPTS), TEMP2(NPTS)
      DIMENSION SHP(NPTS,MNDP), DSDX(NPTS,NDFCD,MNDP), XYZL(NPTS,NDFCD)
      DIMENSION PROP(*), VARI(NPTS,*), DVARI(NPTS,NDFCD,*),LDOFU(*)
      DO 100 I=1, NPTS
UX=VARI(I,LDOFU(KDU))
UY=VARI(I,LDOFU(KDV))
UZ=VARI(I,LDOFU(KDW))
VEL=SQRT(UX*UX+UY*UY+UZ*UZ)
IF(VEL.GT.1)
htra=VEL**0.8
COEF(N)=100
IF(VEL.LE.1)
COEF(N)=10
100 CONTINUE
RETURN
END
回复
分享到:

使用道具 举报

 楼主| 发表于 2007-7-21 08:47 | 显示全部楼层

回复 #1 frank1981 的帖子

flunt公司的人给了我些建议 说是htra=VEL**0.8这个语句不能在subr里面用 而且主程序里stra(set=1,subr,reft=,,,,)的reft要做设定  我现在就是不知道htra=VEL**0.8这个语句到底怎么表达出来了  谁能给个建议阿。。。
发表于 2007-7-21 11:57 | 显示全部楼层
估计这里用fidap的人不多
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 07:10 , Processed in 0.067828 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表