24小时热门版块排行榜    

查看: 1542  |  回复: 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-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 10:39 by 3OOjAIS77qg2
[基金申请] fileCode有新解读? +8 Tide man 2026-08-08 16/800 2026-08-08 09:59 by Tide man
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 08:27 by 6vVgjDL4CnGu
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 KXLV3nuBVBY7 2026-08-07 4/200 2026-08-08 08:21 by 6vVgjDL4CnGu
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 关于代码变化问题,想知道的进来 +9 且听虎啸 2026-08-07 12/600 2026-08-07 21:53 by 理想三旬2025
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
请填处理意见