24小时热门版块排行榜

返回列表

zhuqw

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 1217.4
帖子: 77
在线: 74.6小时
虫号: 2408038
注册: 2013-04-08
性别: GG
专业: 半导体光电子器件

[求助] 模态分析求解出错问题已有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

回复此楼

» 猜你喜欢

1楼 2015-12-30 15:09:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

专家经验: +809
仿真EPI: 1
应助: 1311 (讲师)
金币: 13381.4
红花: 267
帖子: 3531
在线: 473.4小时
虫号: 3644138
注册: 2015-01-14
专业: 动力学与控制
管辖: 仿真模拟

【答案】应助回帖

感谢参与，应助指数 +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 (幼儿园)
金币: 1217.4
帖子: 77
在线: 74.6小时
虫号: 2408038
注册: 2013-04-08
性别: GG
专业: 半导体光电子器件

引用回帖:

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