24小时热门版块排行榜    

查看: 4285  |  回复: 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的回帖

zhuan699

禁虫 (初入文坛)

本帖内容被屏蔽

4楼2017-08-27 22:35:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

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的回帖

kelianbaba

新虫 (小有名气)

请问楼主如何解决的?遇到了同样的问题
5楼2018-07-29 16:50:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 11:32 by l0VvVHGBGRLv
[考博] 找导师 +4 yuanjiabao 2026-08-29 5/250 2026-08-30 09:52 by onesupeng
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 07:05 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 06:57 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见