24小时热门版块排行榜    

查看: 4243  |  回复: 4

zhuqw

金虫 (小有名气)

[求助] 模态分析求解出错问题 已有1人参与

目前在做半导体激光器热应力耦合分析,涉及到模态分析,遇到些问题,卡了半个月了不知道怎么解决。我做的程序分析顺序是这样的:
1.先进行回流时候的应力分析
2.在回流的基础上,也就是以回流后的变形和应力情况作为初始条件,加在模型上,进行工作情况下的热应力分析。

    回流和工作分开的热应力我都能做出来,但是现在加在一起就出现了些问题,程序运行到求解工作时的应力时出现了以下错误提示,然后程序就自动关闭了:

*** ERROR ***                           CP =     196.218   TIME= 14:41:33
The size of the displacement vector in the .ESAV file is 27914.  It     
does not match with the size of the displacement vector in the current  
analysis which is equal to 83742.  Please send the data leading to      
this operation to your technical support provider, as this will allow   
ANSYS, Inc to improve the program.                                      

*** ERROR ***                           CP =     196.748   TIME= 14:43:51
File file2.mode does not exist.                                         

*** WARNING ***                         CP =     197.513   TIME= 14:43:54
The current solution may have been produced using different model or   
boundary condition data than is currently stored.  POST1 results may   
be erroneous unless you perform a new solution using the stored data.   

*** WARNING ***                         CP =     199.478   TIME= 14:44:06
The requested EPPL data is not available.  The PLNSOL command is        
ignored.                                                               

    我以为是因为回流应力求解的数据存在file.ESAV文件下了,所以在求解工作温度的时候,在求解前重新命名一个文件,file2,就是给命令流做了以下修改:

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算工作温度和应力!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/FILNAME,file2,1
/PREP7
ETCHG,STT
ALLSEL
FINISH

    可是问题没有解决,依然出现上面的错误提示框。

    接着我改变了求解的方式,结果就是程序直接崩溃。
*** WARNING ***                         CP =     209.587   TIME= 20:35:04
The modal analysis has non-zero displacements defined, which will cause
a load vector to be generated.  The load vector may be used by the      
mode superposition analyses.  This feature is not documented.           

*** ERROR ***                           CP =     209.681   TIME= 20:35:04
An unexpected error ( SIG$SEGV ) has occurred...  ANSYS internal data   
has been corrupted.  ANSYS is unable to recover and will terminate.     
Previously saved files are unaffected.  Please send the data leading   
to this operation to your technical support provider, as this will      
allow ANSYS, Inc to improve the program.                             



我最终的求解代码是:
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!求解回流!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/SOLU
ANTYPE,STATIC
EQSLV,JCG
TOFFST,273.15
BFV,ALL,TEMP,25
ALLSEL,ALL

TREF,156

NSEL,S,LOC,Y,0
D,ALL,UY,0
NSEL,S,LOC,Z,-12445
D,ALL,UZ,0
ALLSEL
SOLVE
FINISH


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!提取回流结果!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/POST1
ALLSEL,ALL

SET,1
*GET,MXNID,NODE,,NUM,MAX
MAX_NODE=MXNID
*DO,I,1,MAX_NODE,1
  *IF,I,LE,MAX_NODE,THEN
    NODE_SELECTED=NSEL(I)
    NODE_SELECTED=1
  *ENDIF
*ENDDO  
      
*VGET,NL,NODE,,NLIST                 !Return the list of selected nodes
*VGET,EL,ELEM,,ELIST
*VGET,SX,NODE,,S,X
*VGET,SY,NODE,,S,Y
*VGET,SZ,NODE,,S,Z
*VGET,NUX,NODE,,U,X
*VGET,NUY,NODE,,U,Y
*VGET,NUZ,NODE,,U,Z


!Generate Ist File
*cfopen,node-based-ist.dat
*VWRITE
('INISTATE,SET,NODE,1')
*VWRITE,NL(1),SX(1),SY(1),SZ(1),NUX(1),NUY(1),NUZ(1)
('INIS,DEFI,',F6.0,' , , , , ',E,' , ',E,' , ',E,' , ',E,' ,',E,' ')
*CFCLOS


/PREP7
UPGEOM, , , ,file,RST       !Adds displacements from a previous analysis and updates the geometry of the finite element model to the deformed configuration
CDWRITE,COMB,file,CDB
FINISH

/SOLU
!DDELE,ALL,ALL
!D,ALL,ALL
NALL
ANTYPE,STATIC


!*Read in Node Based IST Data
/inp,node-based-ist.dat
ESEL,S,ELEM,,1,9,1
/out
INISTATE,LIST            !读出初始应力状态
ALLSEL,ALL
PSTRESS,ON        !Thermal prestress
EMATWRITE,YES

SOLVE
FINISH

!*Perform prestress modal analsis
!*计算工作温度

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算工作温度!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/FILNAME,file2,1
/PREP7
ETCHG,STT
ALLSEL
FINISH

/SOLU
ANTYPE,STATIC
DDELE,ALL,ALL
EQSLV,SPARSE
BCSOPTION,,INCOERE !运行核内计算
TOFFST,273.15
NSEL,S,LOC,Y,0
D,ALL,TEMP,25
ASEL,ALL

VSEL,S,LOC,Y,8255,8405

*SET,TOTAL_OUTPUT_POWER,20
SINGLE_AVERAGE_POWER=TOTAL_OUTPUT_POWER/19
SINGLE_VOLUME=(150E-6)*(10E-8)*(2E-3)
TEMP_HEAT_GENERATION=SINGLE_AVERAGE_POWER/SINGLE_VOLUME
HEAT_GENERATION=TEMP_HEAT_GENERATION*(10E-6)
BFV,ALL,HGEN,HEAT_GENERATION
ALLSEL,ALL
SOLVE
FINISH

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!计算工作应力!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/PREP7
ETCHG,TTS
TREF,25
ALLSEL
FINISH

/SOLU
DDELE,ALL,ALL
LDREAD,TEMP,,,,,,RTH
NSEL,S,LOC,Y,0
D,ALL,UY,0
NSEL,S,LOC,Z,-12445
D,ALL,UZ,0
NALL
ANTYPE,MODAL
PSTRESS,ON
EQSLV,SPARSE,1.0E-3
MODOPT,LANB,6,
!PSOLVE,CGSOL
PSOLVE,EIGLANB
ALLSEL
SOLVE
FINISH
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
如果是找不到 file2 的有关文件,在回流求解完成,你将工作文件名改为 file2 后,可以将 file.mode 和 file.esav 分别拷贝为 file2.mode 和 file2.esav,前提是两个 file.... 文件存在;如果两个 file.... 文件本身就不存在,则是你的求解过程设计有问题。
2楼2015-12-31 07:49:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuqw

金虫 (小有名气)

引用回帖:
2楼: Originally posted by shgao20 at 2015-12-31 07:49:50
如果是找不到 file2 的有关文件,在回流求解完成,你将工作文件名改为 file2 后,可以将 file.mode 和 file.esav 分别拷贝为 file2.mode 和 file2.esav,前提是两个 file.... 文件存在;如果两个 file.... 文件本身 ...

回流求解的file....文件只有file.esav,因为那个是稳态求解,没有mode文件。
然后我现在把中间提取回流数据部分删除了,直接在回流部分开启预应力效果,在工作应力部分也开预应力,最后依然是算到工作应力部分程序崩溃。没有提示说找不到文件啥的

*** WARNING ***                         CP =     133.740   TIME= 10:17:39
The modal analysis has non-zero displacements defined, which will cause
a load vector to be generated.  The load vector may be used by the      
mode superposition analyses.  This feature is not documented.           

*** ERROR ***                           CP =     133.849   TIME= 10:17:39
An unexpected error ( SIG$SEGV ) has occurred...  ANSYS internal data   
has been corrupted.  ANSYS is unable to recover and will terminate.     
Previously saved files are unaffected.  Please send the data leading   
to this operation to your technical support provider, as this will      
allow ANSYS, Inc to improve the program.
3楼2015-12-31 10:26:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuan699

禁虫 (初入文坛)

本帖内容被屏蔽

4楼2017-08-27 22:35:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kelianbaba

新虫 (小有名气)

请问楼主如何解决的?遇到了同样的问题
5楼2018-07-29 16:50:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhuqw 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] fileCode有新解读? +8 Tide man 2026-08-08 16/800 2026-08-08 09:59 by Tide man
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 08:27 by 6vVgjDL4CnGu
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 KXLV3nuBVBY7 2026-08-07 4/200 2026-08-08 08:21 by 6vVgjDL4CnGu
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见