24小时热门版块排行榜    

查看: 1365  |  回复: 6

952339077

新虫 (初入文坛)

[求助] ANSYS里面apdl命令流求助!!!!!!!!!!急 已有1人参与

/POST1

Nsel,s,loc,x,-10,10
Nsel,r,loc,y,-9,9
Nsel,r,loc,z,163,177
Esln,s,1

*DIM,SSN,array,10,12

*DO,I,1,10

   SET,LAST
   SET,1,I,1,

*Get,maxn,elem,0,num,max

        ETABLE,ETOX,EPTO,X
        ETABLE,ETOY,EPTO,Y
        ETABLE,ETOZ,EPTO,Z
        ETABLE,ETOXY,EPTO,XY
        ETABLE,ETOYZ,EPTO,YZ
        ETABLE,ETOXZ,EPTO,XZ

        ETABLE,ETSX,S,X
        ETABLE,ETSY,S,Y
        ETABLE,ETSZ,S,Z
        ETABLE,ETSXY,S,XY
        ETABLE,ETSYZ,S,YZ
        ETABLE,ETSXZ,S,XZ

        *GET,SSN(I,1),ELEM,maxn,ETAB,ETOX
        *GET,SSN(I,2),ELEM,maxn,ETAB,ETOXY
        *GET,SSN(I,3),ELEM,maxn,ETAB,ETOXZ
        *GET,SSN(I,4),ELEM,maxn,ETAB,ETOY
        *GET,SSN(I,5),ELEM,maxn,ETAB,ETOYZ
        *GET,SSN(I,6),ELEM,maxn,ETAB,ETOZ      

        *GET,SSN(I,7),ELEM,maxn,ETAB,ETSX
        *GET,SSN(I,8),ELEM,maxn,ETAB,ETSXY
        *GET,SSN(I,9),ELEM,maxn,ETAB,ETSXZ
        *GET,SSN(I,10),ELEM,maxn,ETAB,ETSY
        *GET,SSN(I,11),ELEM,maxn,ETAB,ETSYZ
        *GET,SSN(I,12),ELEM,maxn,ETAB,ETSZ

*ENDDO

*MWRITE,SSN,E:\desktop\ansysoutput\TEST3D60R,txt,,JIK,12,10
(12E15.7)
请问有大神知道我这个怎么算出来结果文件都是0?
我用的四面体实体网格,载荷步是10步,而且,坐标区域内的是存在我需要的单元的,可是怎么结果都是0啊
回复此楼

» 猜你喜欢

疲劳可靠性设计
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggggsahala

金虫 (著名写手)

APDL忘光了
我是来打酱油的
2楼2016-10-19 10:22:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
952339077: 金币+10, 有帮助 2016-10-19 20:18:45
检查一下以下命令流:
  Nsel,s,loc,x,-10,10
  Nsel,r,loc,y,-9,9
  Nsel,r,loc,z,163,177
  Esln,s,1
是否确实选择到了单元,如歌它没有选择到单元,所有的结果就都是零。
3楼2016-10-19 17:03:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

952339077

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by shgao20 at 2016-10-19 17:03:20
检查一下以下命令流:
  Nsel,s,loc,x,-10,10
  Nsel,r,loc,y,-9,9
  Nsel,r,loc,z,163,177
  Esln,s,1
是否确实选择到了单元,如歌它没有选择到单元,所有的结果就都是零。

选到了,非常确定选到了,就想问问,这个程序是不是选择最大单元编号的每个子步的12个应力应变值?
疲劳可靠性设计
4楼2016-10-19 20:17:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
952339077: 金币+30, 有帮助, 回答的不错,但是问题还没有解决 2016-12-02 19:41:49
这个程序是在以下命令流:
    Nsel,s,loc,x,-10,10
    Nsel,r,loc,y,-9,9
    Nsel,r,loc,z,163,177
    Esln,s,1
  所选择的单元中,对具有最大编号的单元,将载荷步 1、子步 1~10, 的 12 个单元应力、应变存放到数组 SSN(10,12) 中。
  如果你是想把 10 个载荷步的最后一个子步的结果存放到数组中,你应该将命令:
    SET,1,I,1,
  改为:
    SET,I,LAST,1,
5楼2016-10-22 10:44:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

952339077

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by shgao20 at 2016-10-22 10:44:26
这个程序是在以下命令流:
    Nsel,s,loc,x,-10,10
    Nsel,r,loc,y,-9,9
    Nsel,r,loc,z,163,177
    Esln,s,1
  所选择的单元中,对具有最大编号的单元,将载荷步 1、子步 1~10, 的 12 个单元应力、应变 ...

能解释一下如何加载1-6个载荷步,每个载荷步下20个子步,总共120个子步,这种如何编写呢?上面的确实是加载了一个载荷步下的10个子步
疲劳可靠性设计
6楼2016-12-02 19:44:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

你的意思是说,你的结果文件中共有 6 个载荷步,各 20 个子步的结果吗?如是,可以试试如下命令流:

  /POST1

Nsel,s,loc,x,-10,10
Nsel,r,loc,y,-9,9
Nsel,r,loc,z,163,177
Esln,s,1
*Get,maxn,elem,0,num,max

*DIM,SSN,array,6,20,12
*cfopen,E:\desktop\ansysoutput\TEST3D60R,txt, ,

*DO,I,1,6
*do,j,1,20

   SET,i,j    !  读取第 i 个载荷步,第 j 个子步的结果

        ETABLE,ETOX,EPTO,X
        ETABLE,ETOY,EPTO,Y
        ETABLE,ETOZ,EPTO,Z
        ETABLE,ETOXY,EPTO,XY
        ETABLE,ETOYZ,EPTO,YZ
        ETABLE,ETOXZ,EPTO,XZ

        ETABLE,ETSX,S,X
        ETABLE,ETSY,S,Y
        ETABLE,ETSZ,S,Z
        ETABLE,ETSXY,S,XY
        ETABLE,ETSYZ,S,YZ
        ETABLE,ETSXZ,S,XZ

        *GET,SSN(I,J, 1),ELEM,maxn,ETAB,ETOX
        *GET,SSN(I,J, 2),ELEM,maxn,ETAB,ETOXY
        *GET,SSN(I,J, 3),ELEM,maxn,ETAB,ETOXZ
        *GET,SSN(I,J, 4),ELEM,maxn,ETAB,ETOY
        *GET,SSN(I,J, 5),ELEM,maxn,ETAB,ETOYZ
        *GET,SSN(I,J, 6),ELEM,maxn,ETAB,ETOZ      

        *GET,SSN(I,J, 7),ELEM,maxn,ETAB,ETSX
        *GET,SSN(I,J, 8),ELEM,maxn,ETAB,ETSXY
        *GET,SSN(I,J, 9),ELEM,maxn,ETAB,ETSXZ
        *GET,SSN(I,J, 10),ELEM,maxn,ETAB,ETSY
        *GET,SSN(I,J, 11),ELEM,maxn,ETAB,ETSYZ
        *GET,SSN(I,J, 12),ELEM,maxn,ETAB,ETSZ

*VWRITE, i, j, (ssn(i,j,k),k=1,6)
(f2.0,f4.0,2x,6E15.7)

*VWRITE, (ssn(i,j,k),k=7,126)
(8x,6E15.7/)

*ENDDO
*EDNDO

*cfclose

如果不使用上面的 *cfopen,*vwrite 和 *cfclose 命令,也可以在循环之后使用如下命令:

*MWRITE,SSN,E:\desktop\ansysoutput\TEST3D60R,txt,,KIJ,12, 6, 20
(12E15.7)

因为我无法验证,只能请你自己验证并修改完善了。
7楼2016-12-03 08:59:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 952339077 的主题更新
信息提示
请填处理意见