应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【求助】ANSYS 中命令流编程（fortran）
121/2返回列表12下一页
查看: 4930  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yaojinnan

金虫 (正式写手)

[交流] 【求助】ANSYS 中命令流编程（fortran） 已有11人参与

大家好，我是做焊接模拟的，用的是ANSYS软件。现在想做一段程序，其中想做个循环交替加载。
比如：载荷作用时间总共为10s，有两个载荷。
          其中，0-1s 给选定物体加载第一个载荷，
                    1-3s给选定物体加载第二个载荷
                    3-4s给选定物体加载第一个载荷
                    5-7s给选定物体加载第二个载荷
                    7-8s给选定物体加载第一个载荷
                    8-10s给选定物体加载第二个载荷
其实简单的可以这么理解。
但是我的载荷是一个由ANSYS函数编辑器定义的一个表载荷，
并且这两个载荷都是移动的。所以不清楚应该怎么用FORTRAN来编程。希望大家能给予帮助。定有回报。
下面就是我的命令流程序：
finish
/clear,start                                    !清除数据库以开始新的分
/filname,reverse-shuang tuo qiu,0
/title,hanhanahanhan
/units,si

/prep7

ET,1,SOLID70
ET,2,SOLID90
!et,2,surf152
!keyopt,2,4,0
!keyopt,2,5,1
!keyopt,2,8,3

MPTEMP,1,20,100,200,300,400,500,
MPTEMP,7,600,700,800,900,1000,1500
MPDATA,KXX,1,1,152,159,164,165,168,196,
MPDATA,KXX,1,7,209,162,162,162,162,162
MPDATA,DENS,1,1,2680,2662,2646,2620,2600,2580,
MPDATA,DENS,1,7,2436,2384,2384,2384,2384,2384
MPDATA,C,1,1,956,963,1047,1130,1224,1308,
MPDATA,C,1,7,1412,1084,1084,1084,1084,1084
MPDATA,HF,1,1,8.22,11.0,13.7,23.2,33.4,46.8,
MPDATA,HF,1,7,58.0,68.5,68.5,68.5,68.5,68.5

lx=0.08
ly=0.2
lz=0.005

r=0.002
D=0.002！焊缝宽度
ah=D/2！双椭球热源参数定义
chf=4*ah
chb=6*ah
bh=0.003
rf=2*chf/(chf+chb)
rb=2*chb/(chf+chb)
K=0.5 功率有效系数
I=130
U=30
q0=I*U
q=q0*K  !有效功率
pi=acos(-1)
step=0.002
v=0.003
!count=ly/v
count=20
!***********************************
!定义单元类型
!***********************************
!建立几何模型
!***********************************
wpstyle,,,,,,,,0
block,0,0.006,0,ly,0,-lz
block,0,0.012,0,ly,0,-lz
block,0,0.024,0,ly,0,-lz
block,0,lx,0,ly,0,-lz
vovlap,all
/view,1,1,1,1

!定义焊缝区单元尺寸0.002
lesize,7,0.001
lesize,2,0.001
lesize,12,0.001
lesize,11,0.001
lesize,6,0.001
lesize,3,0.001
lesize,8,0.001
lesize,1,0.001
lesize,4,0.001
lesize,9,0.001
lesize,5,0.001
lesize,10,0.001
!定义过渡区的单元尺寸
lesize,15,r
lesize,18,r
!定义远离焊缝区的单元尺寸
lesize,34,1.25*0.001
lesize,35,1.25*0.001
lesize,47,1.25*0.001
lesize,46,1.25*0.001
lesize,39,4*r
lesize,42,4*r
lesize,57,4*r
lesize,60,4*r
lesize,27,4*r
lesize,30,4*r
lesize,58,4*r
lesize,59,4*r
!***********************************
!划分网格，生生有限元模型
!***********************************
type,1
mshkey,1
mshape,0
vmesh,1
vmesh,7
type,2
mshkey,0
mshape,1
vmesh,5
vmesh,6
numcmp,area
save
!***********************************
!循环加载、求解
!***********************************
!施加载荷
/solu
*do,i,0,count
disy=i*step
t=(disy+step)/v
antype,4
trnopt,full
tunif,293        
outres,basic,last      
time,t
nropt,full,,on        
autots,1
nsubst,1,,,1     
kbc,1
lnsrch,1
tsres,erase
!TINTP,,,,1,,,,,
timint,1,therm
!***********************************
!焊接过程开始
!***********************************
nsel,s,ext,
nsel,u,loc,x,0       !------改
sf,all,conv,200,293    !------改
allsel,all
!***********************************
!***********************************
!双椭球热源的加载
!***********************************
!前半球节点的选取
!***********************************
LOCAL,11,2,0,disy,0.000,,,,chf/ah,bh/ah,
NSEL,S,LOC,X,0,ah   
NSEL,R,LOC,Y,0,90   
NSEL,R,LOC,Z,-90,0
CM,FRONT,NODE   
LOCAL,13,0,0,disy,0.000,,,,1,1,
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*DEL,'TFRONT'   
*SET,_FNCNAME,'TFRONT'  
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),q   
*SET,_FNC_C2(1),rf  
*SET,_FNC_C3(1),ah  
*SET,_FNC_C4(1),chf
*SET,_FNC_C5(1),bh  
*SET,_FNCCSYS,13
*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS%
%_FNCNAME%(0,0,1)= 0.0, -999
%_FNCNAME%(2,0,1)= 0.0
%_FNCNAME%(3,0,1)= %_FNC_C1(1)%
%_FNCNAME%(4,0,1)= %_FNC_C2(1)%
%_FNCNAME%(5,0,1)= %_FNC_C3(1)%
%_FNCNAME%(6,0,1)= %_FNC_C4(1)%
%_FNCNAME%(7,0,1)= %_FNC_C5(1)%
%_FNCNAME%(0,1,1)= 1.0, -1, 16, 3, 0, 0, 0
%_FNCNAME%(0,2,1)= 0.0, -2, 0, 6, 0, 0, -1
%_FNCNAME%(0,3,1)=   0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,4,1)= 0.0, -1, 0, 1, -3, 3, 17
%_FNCNAME%(0,5,1)= 0.0, -2, 0, 1, -1, 3, 18
%_FNCNAME%(0,6,1)= 0.0, -1, 0, 0, 0, 0, 0
%_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, 0, 0, -1
%_FNCNAME%(0,8,1)= 0.0, -4, 0, 1, -1, 2, -3
%_FNCNAME%(0,9,1)= 0.0, -1, 0, 3, 0, 0, -4
%_FNCNAME%(0,10,1)= 0.0, -3, 0, 1, -4, 3, -1
%_FNCNAME%(0,11,1)= 0.0, -1, 0, 2, 0, 0, 2
%_FNCNAME%(0,12,1)= 0.0, -4, 0, 1, 2, 17, -1
%_FNCNAME%(0,13,1)= 0.0, -1, 0, 2, 0, 0, 19
%_FNCNAME%(0,14,1)= 0.0, -5, 0, 1, 19, 17, -1
%_FNCNAME%(0,15,1)= 0.0, -1, 0, 1, -4, 4, -5
%_FNCNAME%(0,16,1)= 0.0, -4, 0, 2, 0, 0, 3
%_FNCNAME%(0,17,1)= 0.0, -5, 0, 1, 3, 17, -4
%_FNCNAME%(0,18,1)= 0.0, -4, 0, 2, 0, 0, 20
%_FNCNAME%(0,19,1)= 0.0, -6, 0, 1, 20, 17, -4
%_FNCNAME%(0,20,1)= 0.0, -4, 0, 1, -5, 4, -6
%_FNCNAME%(0,21,1)= 0.0, -5, 0, 1, -1, 1, -4
%_FNCNAME%(0,22,1)= 0.0, -1, 0, 2, 0, 0, 4
%_FNCNAME%(0,23,1)= 0.0, -4, 0, 1, 4, 17, -1
%_FNCNAME%(0,24,1)= 0.0, -1, 0, 2, 0, 0, 21
%_FNCNAME%(0,25,1)= 0.0, -6, 0, 1, 21, 17, -1
%_FNCNAME%(0,26,1)= 0.0, -1, 0, 1, -4, 4, -6
%_FNCNAME%(0,27,1)= 0.0, -4, 0, 1, -5, 1, -1
%_FNCNAME%(0,28,1)= 0.0, -1, 0, 1, -3, 3, -4
%_FNCNAME%(0,29,1)= 0.0, -1, 7, 1, -1, 0, 0
%_FNCNAME%(0,30,1)= 0.0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,31,1)= 0.0, -1, 16, 3.14159265358979310, 0, 0, 0
%_FNCNAME%(0,32,1)= 0.0, -2, 0, 3.14159265358979310, 0, 0, -1
%_FNCNAME%(0,33,1)= 0.0, -4, 0, 1, -2, 3, -1
%_FNCNAME%(0,34,1)= 0.0, -1, 0, 1, -4, 3, 20
%_FNCNAME%(0,35,1)= 0.0, -2, 0, 1, -1, 3, 19
%_FNCNAME%(0,36,1)= 0.0, -1, 0, 1, -2, 3, 21
%_FNCNAME%(0,37,1)= 0.0, -2, 0, 1, -3, 4, -1
%_FNCNAME%(0,38,1)= 0.0, 99, 0, 1, -2, 0, 0
BF,FRONT,HGEN,%TFRONT%
CMDELE,FRONT
ALLSEL,ALL

!!!后半球节点的选取
LOCAL,12,2,0,disy,0,,,,chb/ah,bh/ah,
NSEL,S,LOC,X,0,ah   
NSEL,R,LOC,Y,-90,0  
nsel,r,loc,y,-90,0  
CM,BEHIND,NODE  
!!!
LOCAL,14,0,0,disy,0,,,,1,1,
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*DEL,'TBEHIND'   
*SET,_FNCNAME,'TBEHIND'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),q   
*SET,_FNC_C2(1),rb  
*SET,_FNC_C3(1),ah  
*SET,_FNC_C4(1),chb
*SET,_FNC_C5(1),bh  
*SET,_FNCCSYS,14
*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS%
%_FNCNAME%(0,0,1)= 0.0, -999
%_FNCNAME%(2,0,1)= 0.0
%_FNCNAME%(3,0,1)= %_FNC_C1(1)%
%_FNCNAME%(4,0,1)= %_FNC_C2(1)%
%_FNCNAME%(5,0,1)= %_FNC_C3(1)%
%_FNCNAME%(6,0,1)= %_FNC_C4(1)%
%_FNCNAME%(7,0,1)= %_FNC_C5(1)%
%_FNCNAME%(0,1,1)= 1.0, -1, 16, 3, 0, 0, 0
%_FNCNAME%(0,2,1)= 0.0, -2, 0, 6, 0, 0, -1
%_FNCNAME%(0,3,1)=   0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,4,1)= 0.0, -1, 0, 1, -3, 3, 17
%_FNCNAME%(0,5,1)= 0.0, -2, 0, 1, -1, 3, 18
%_FNCNAME%(0,6,1)= 0.0, -1, 0, 0, 0, 0, 0
%_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, 0, 0, -1
%_FNCNAME%(0,8,1)= 0.0, -4, 0, 1, -1, 2, -3
%_FNCNAME%(0,9,1)= 0.0, -1, 0, 3, 0, 0, -4
%_FNCNAME%(0,10,1)= 0.0, -3, 0, 1, -4, 3, -1
%_FNCNAME%(0,11,1)= 0.0, -1, 0, 2, 0, 0, 2
%_FNCNAME%(0,12,1)= 0.0, -4, 0, 1, 2, 17, -1
%_FNCNAME%(0,13,1)= 0.0, -1, 0, 2, 0, 0, 19
%_FNCNAME%(0,14,1)= 0.0, -5, 0, 1, 19, 17, -1
%_FNCNAME%(0,15,1)= 0.0, -1, 0, 1, -4, 4, -5
%_FNCNAME%(0,16,1)= 0.0, -4, 0, 2, 0, 0, 3
%_FNCNAME%(0,17,1)= 0.0, -5, 0, 1, 3, 17, -4
%_FNCNAME%(0,18,1)= 0.0, -4, 0, 2, 0, 0, 20
%_FNCNAME%(0,19,1)= 0.0, -6, 0, 1, 20, 17, -4
%_FNCNAME%(0,20,1)= 0.0, -4, 0, 1, -5, 4, -6
%_FNCNAME%(0,21,1)= 0.0, -5, 0, 1, -1, 1, -4
%_FNCNAME%(0,22,1)= 0.0, -1, 0, 2, 0, 0, 4
%_FNCNAME%(0,23,1)= 0.0, -4, 0, 1, 4, 17, -1
%_FNCNAME%(0,24,1)= 0.0, -1, 0, 2, 0, 0, 21
%_FNCNAME%(0,25,1)= 0.0, -6, 0, 1, 21, 17, -1
%_FNCNAME%(0,26,1)= 0.0, -1, 0, 1, -4, 4, -6
%_FNCNAME%(0,27,1)= 0.0, -4, 0, 1, -5, 1, -1
%_FNCNAME%(0,28,1)= 0.0, -1, 0, 1, -3, 3, -4
%_FNCNAME%(0,29,1)= 0.0, -1, 7, 1, -1, 0, 0
%_FNCNAME%(0,30,1)= 0.0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,31,1)= 0.0, -1, 16, 3.14159265358979310, 0, 0, 0
%_FNCNAME%(0,32,1)= 0.0, -2, 0, 3.14159265358979310, 0, 0, -1
%_FNCNAME%(0,33,1)= 0.0, -4, 0, 1, -2, 3, -1
%_FNCNAME%(0,34,1)= 0.0, -1, 0, 1, -4, 3, 20
%_FNCNAME%(0,35,1)= 0.0, -2, 0, 1, -1, 3, 19
%_FNCNAME%(0,36,1)= 0.0, -1, 0, 1, -2, 3, 21
%_FNCNAME%(0,37,1)= 0.0, -2, 0, 1, -3, 4, -1
%_FNCNAME%(0,38,1)= 0.0, 99, 0, 1, -2, 0, 0
BF,BEHIND,HGEN,%TBEHIND%
CMDELE,BEHIND
ALLSEL,ALL
SOLVE
BFDELE,all,HGEN
*ENDDO
SAVE

其中从双椭球热源加载到最后就是加载的过程，这是一个载荷，还需要另外一个载荷，这里没有加载，就像上面我所说的那样循环交替那样应该怎么加载，希望大家给予帮助。。应该用到怎么样的语句呢？？
回复此楼

» 收录本帖的淘帖专辑推荐

硕博精华

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
只有有了优秀的对手，你才会变得更加优秀
1楼 2010-09-09 20:40:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

PZtech

银虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
jjdg(金币+1): 感谢参与 2011-10-22 00:34:34
把命令先写在文本文档中,然后拷贝到命令行中.
赞一下(3人) 回复此楼
6楼2011-10-21 21:06:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Liv_Guo

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主有没有用函数表示热源的，不用table表格
赞一下(1人) 回复此楼
10楼2018-04-11 16:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

graduate7303

至尊木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
学习学习了。。。
回复此楼
2楼2010-10-18 14:06:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dgh338

新虫 (初入文坛)
这么多的程序，看不懂
赞一下 回复此楼
3楼2010-10-19 15:45:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by yaojinnan at 2010-09-09 20:40:33:
大家好，我是做焊接模拟的，用的是ANSYS软件。现在想做一段程序，其中想做个循环交替加载。
比如：载荷作用时间总共为10s，有两个载荷。
          其中，0-1s 给选定物体加载第一个载荷，
                   ...

ANSYS怎么打开命令流编程？
赞一下(1人) 回复此楼
好好学习，天天向上。
4楼2010-11-12 19:48:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leixiaokuan

新虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖
jjdg: 先google一下吧 2011-10-19 06:42:39
我是做焊接的，想学Ansys模拟焊接过程，有这方面的教程吗？若有请传我一份，邮箱：leixiaokuan_0512@1363.com
赞一下(2人) 回复此楼
结交广大焊接工程师
5楼2011-10-18 13:41:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

16443

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
使用多载荷步实现
赞一下(1人) 回复此楼
~~SCI---StupidChineseIndex
7楼2011-11-18 11:17:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

包子小妹

新虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你好，我想问一下设置求解选项怎么弄的？分为表格和载荷步，但是我不知道这两个的区别，包括哪些命令该用，比如trnopt，trses，这些，命令什么时候会用到呢？
赞一下 回复此楼
8楼2014-07-02 14:41:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

热七月流火

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主，有没有解决，我也遇到这种问题，求分享啊
赞一下 回复此楼
9楼2017-11-06 11:39:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yaojinnan 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-16 01:21 by Xttdmn
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 3/150 2026-03-15 17:32 by 小物理化学
[考研] 0703化学调剂 290分有科研经历，论文在投 +3 腻腻gk 2026-03-14 3/150 2026-03-15 17:28 by 小物理化学
[基金申请] 国自科面上基金字体 +4 iwuli 2026-03-12 5/250 2026-03-15 17:07 by 风云无泪
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 293求调剂 +5 上班不着吉 2026-03-09 5/250 2026-03-14 02:37 by JourneyLucky
[考研] 化学工程321分求调剂（南京工业，浙江工业） +3 大米饭！ 2026-03-09 4/200 2026-03-14 02:34 by JourneyLucky
[考研] 332分材料工程调剂 +3 莓好时光海苔 2026-03-09 3/150 2026-03-14 02:03 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
[考研] 0857 资源与环境 285分 +6 未名考生 2026-03-09 6/300 2026-03-11 21:08 by 30660438
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[基金申请] 提交后的基金本子，已让学校撤回了，可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
关闭
请填处理意见
关闭