{"id":167542,"date":"2022-03-21T10:48:27","date_gmt":"2022-03-21T10:48:27","guid":{"rendered":"\/forum\/forums\/topic\/thermal-steady-state-analysis-with-ansys-mechanical-apdl-and-matlab\/"},"modified":"2022-03-23T17:13:26","modified_gmt":"2022-03-23T17:13:26","slug":"thermal-steady-state-analysis-with-ansys-mechanical-apdl-and-matlab","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/thermal-steady-state-analysis-with-ansys-mechanical-apdl-and-matlab\/","title":{"rendered":"thermal steady state analysis with Ansys Mechanical APDL and Matlab"},"content":{"rendered":"
\n
\n

Hi everybody,<\/p>\n

i'm launching a Steady State Thermal analysis with an Ansys Mechanical APDL macro code from Matlab in order to automatize a decoupled CHT caculus, but i'm facing with some issues. The matlab string used to launch the simulation is the following: <\/p>\n

<\/p>\n<\/p>\n

"!"C:Program FilesANSYS directoryANSYS192.exe" -b -i input.mac -o output.txt"<\/p>\n

<\/p>\n<\/p>\n

where the file "input.mac" has been created in ANSYS Mechanical with the command "write input file" and contains either all the informations about geometry, mesh, boundary conditions, and the command "solve" in \/solu environment, or the manually written macro code for the results export (written in \/post1 environment after the "solve" command).<\/p>\n

The main issue is that commands written after "solve" are not considered, no matter what i write (i tried to write some random digits but it doesnt give any errors\/warnings) included the macro code for the results export.T he output file doesn't show any errors and at the end "solution converged after XX iterations" is written.<\/p>\n

For a further test i tried to load the macro code on Ansys Mechanical APDL with the graphical interface. After the loading, solution data are not available and a manually written "\/solu solve" is needed.<\/p>\n

I report the \/solu and the \/post1 environments of "input.mac" here:<\/p>\n

<\/p>\n<\/p>\n

\/com,****************************************************************************<\/p>\n

\/com,*************************  SOLUTION    ********************************<\/p>\n

\/com,****************************************************************************<\/p>\n

\/solu<\/p>\n

antype,0          ! static analysis<\/p>\n

!eqsl,sparse        ! using sparse solver which is default so no eqsl command needed<\/p>\n

cntr,print,1        ! print out contact info and also make no initial contact an error<\/p>\n

rescontrol,,none         ! Do not keep any restart files<\/p>\n

dmpoption,emat,no       ! Don't combine emat file for DANSYS<\/p>\n

dmpoption,esav,no       ! Don't combine esav file for DANSYS<\/p>\n

\/com,****************************************************<\/p>\n

\/com,******************* SOLVE FOR LS 1 OF 1 ****************<\/p>\n

esel,s,type,,18<\/p>\n

nsle<\/p>\n

sf,all,hflux,5000<\/p>\n

!OTHER BOUNDARY CONDITIONS ARE NOT REPORTED<\/p>\n

nsel,all<\/p>\n

esel,all<\/p>\n

autots,on<\/p>\n

nsub,1,10,1<\/p>\n

time,1.<\/p>\n

outres,erase<\/p>\n

outres,all,none<\/p>\n

outres,nsol,all<\/p>\n

outres,rsol,all<\/p>\n

outres,eangl,all<\/p>\n

outres,etmp,all<\/p>\n

outres,veng,all<\/p>\n

outres,fflux,all<\/p>\n

outres,cont,all<\/p>\n

\/nopr<\/p>\n

esel,s,type,,22<\/p>\n

esel,a,type,,23,49<\/p>\n

cm,_elmisc,elem<\/p>\n

esel,all<\/p>\n

\/gopr<\/p>\n

outres,misc,all,_elmisc<\/p>\n

cnvtol,heat,,.001,,1.e-006  ! Program Controlled. Set MINREF equivalent of 1e-6 Watts<\/p>\n

! *********** WB SOLVE COMMAND ***********<\/p>\n

! check interactive state<\/p>\n

*get,ANSINTER_,active,,int<\/p>\n

*if,ANSINTER_,ne,0,then<\/p>\n

\/eof<\/p>\n

*endif<\/p>\n

solve<\/p>\n

!PROBLEM: COMMANDS WRITTEN BELOW ARE SKIPPED<\/p>\n

\/com *************** Write FE CONNECTORS ********* <\/p>\n

CEWRITE,file,ce,,INTE<\/p>\n

\/com,****************************************************<\/p>\n

\/com,*************** FINISHED SOLVE FOR LS 1 *************<\/p>\n

*get,_wallasol,active,,time,wall<\/p>\n

<\/p>\n<\/p>\n

\/post1<\/p>\n

\/INPUT,'postpr','dat','C:Users
\nafanelliDocumentsITER1A',, 0  !MACRO FOR POSTPROCESS<\/p>\n

xmlo,ENCODING,ISO-8859-1<\/p>\n

xmlo,parm<\/p>\n

\/xml,parm,xml<\/p>\n

fini<\/p>\n

\/gopr<\/p>\n

*get,_walldone,active,,time,wall<\/p>\n

_preptime=(_wallbsol-_wallstrt)*3600<\/p>\n

_solvtime=(_wallasol-_wallbsol)*3600<\/p>\n

_posttime=(_walldone-_wallasol)*3600<\/p>\n

_totaltim=(_walldone-_wallstrt)*3600<\/p>\n

*get,_dlbratio,active,0,solu,dlbr<\/p>\n

\/com,— Number of total nodes = %_numnode%<\/p>\n

\/com,— Number of total elements = %_numelem%<\/p>\n

\/com,— Element load balance ratio = %_dlbratio%<\/p>\n

<\/p>\n<\/p>\n

\/wb,file,end        ! done with WB generated input<\/p>\n

<\/p>\n

<\/p>\n<\/p>\n

I think there is a problem with the solution command but i'm not able to figure out.<\/p>\n

Thank you so much<\/p>\n","protected":false},"template":"","class_list":["post-167542","topic","type-topic","status-closed","hentry","topic-tag-ansys-apdl","topic-tag-apdl-command","topic-tag-APDLthermal-1","topic-tag-matlab","topic-tag-mechanical-apdl"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_btv_view_count":["1604"],"_bbp_likes_count":["0"],"_bbp_subscription":["251712"],"_bbp_topic_status":["unanswered"],"_bbp_status":["publish"],"_bbp_topic_id":["167542"],"_bbp_forum_id":["27791"],"_bbp_engagement":["88","23687","79912","251712"],"_bbp_voice_count":["4"],"_bbp_reply_count":["3"],"_bbp_last_reply_id":["209064"],"_bbp_last_active_id":["209064"],"_bbp_last_active_time":["2022-03-23 17:13:26"]},"test":"isacco-rafanelliunifi-it"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/167542","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/167542\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=167542"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}