| 查看: 4281 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13381.4
- 红花: 267
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
2楼2015-12-31 07:49:50
zhuqw
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1217.4
- 帖子: 77
- 在线: 74.6小时
- 虫号: 2408038
- 注册: 2013-04-08
- 性别: GG
- 专业: 半导体光电子器件
|
回流求解的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
|
本帖内容被屏蔽 |
4楼2017-08-27 22:35:45
5楼2018-07-29 16:50:39
30