24小时热门版块排行榜    

查看: 1579  |  回复: 3

zhuqw

金虫 (小有名气)

[求助] 预应力求解问题

问题是这样的,我现在要做激光器封装应力应变分析,先进行回流,在回流求的形变和应力的基础上在求解激光器工作时候的应力和应变。单独进行回流和单独进行工作的都做过,都没问,但是现在整合到一块就各种报错,现在有一个问题,程序跑不下去了,错误日志的提示是这样:
*** ERROR ***                           CP =     143.209   TIME= 21:08:24
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 =     143.677   TIME= 21:08:58
File file.mode does not exist.                                          


我的命令流是这样的,建模和网格就不放了,因为模型是最简化的用来测试求解代码能不能运行的,所以网格很粗陋,错误很多(已经2000多K了错误日志),不过我觉得应该不影响后面求解命令的执行。求解命令:

*** ERROR ***                           CP =     143.209   TIME= 21:08:24
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 =     143.677   TIME= 21:08:58
File file.mode does not exist.                                          

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!求解回流!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!*
/SOLU
ANTYPE,STATIC
EQSLV,JCG
TOFFST,273.15
NSEL,S,LOC,Y,0
D,ALL,TEMP,25
ALLSEL,ALL
SOLVE
FINISH


/PREP7
ETCHG,TTS
TREF,156
ALLSEL
FINISH

/SOLU
ANTYPE,STATIC
LDREAD,TEMP,,,,,,RTH
NSEL,S,LOC,Y,0
D,ALL,UY,0
NSEL,S,LOC,Z,-12445
D,ALL,UZ,0
ALLSEL

SOLVE
FINISH

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!求解工作!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/POST1
*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,SX,NODE,,S,X
*VGET,SY,NODE,,S,Y
*VGET,SZ,NODE,,S,Z
*VGET,UX,NODE,,U,X
*VGET,UY,NODE,,U,Y
*VGET,UZ,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),UX(1),UY(1),UZ(1)
('INIS,DEFI,',F6.0,' , , , , ',E,' , ',E,' , ',E,',0.0,0.0,0.0')
*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
!*计算工作温度

/PREP7
ETCHG,STT
ALLSEL
FINISH

/SOLU
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,X,0
DSYM,SYMM,X,,
NSEL,S,LOC,Y,0
D,ALL,UY,0
NSEL,S,LOC,Z,-12445
NALL
ANTYPE,MODAL
PSTRESS,ON
MODOPT,LANB,5,
!EMATWRITE,YES
PSOLVE,ELFORM
ALLSEL
SOLVE
FINISH


命令是参考着ANSYS帮助手册的node-based initial写的,很多还不明白为什么,比如NALL是什么的缩写?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

zhuqw

金虫 (小有名气)

新出来的两条错误信息,是接在那两个信息之后的

*** WARNING ***                         CP =     149.090   TIME= 09:43:57
SOLVE is not a recognized BEGIN command, abbreviation, or macro.  This  
command will be ignored.                                                

*** WARNING ***                         CP =     149.636   TIME= 09:44:02
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.
2楼2015-12-23 09:47:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuan699

禁虫 (初入文坛)

本帖内容被屏蔽

3楼2017-08-20 12:37:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuan699

禁虫 (初入文坛)

本帖内容被屏蔽

4楼2017-08-20 12:38:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhuqw 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:25 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:41 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 4/200 2026-08-29 22:27 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +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
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
请填处理意见