24小时热门版块排行榜    

查看: 1577  |  回复: 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文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 01:04 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.O.55.1.O.5.4,科目全,可+急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:52 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:33 by 4qO8zhfjr2AX
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-29 21:33 by 4qO8zhfjr2AX
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 3/150 2026-08-29 21:32 by 4qO8zhfjr2AX
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见