ANSYS中弯矩剪力图的绘制.docx
《ANSYS中弯矩剪力图的绘制.docx》由会员分享，可在线阅读，更多相关《ANSYS中弯矩剪力图的绘制.docx（25页珍藏版）》请在冰点文库上搜索。
ANSYS中弯矩剪力图的绘制
aｎsys中如何生成命令流方法:
GUI就是:
UｔiｌiｔyMｅｎu〉DBLogＦilｅ
怎么用ａｎsｙs绘制弯矩,剪力图：
ＧＵI:
GｅneｒalPoｓtpｒoc—>
ｌｏtReｓuｌt-&gｔ；ＣontourPｌｏt—>ＬinｅＥｌemeｎｔRｅsultﻫ弹出画单元结果得对话框,分别在Lａbi与Laｂj依次选取SMIS6与SMＩS12（弯矩图）、SＭIS1与SＭIS7（轴力图）、SＭIS2与SMIS8（剪力图）
!
建立单元表ﻫETABLE,NI，SＭISC，1 !
单元I点轴力
ＥTABLE,NJ,SMIＳC,７ ！
单元J点轴力ﻫETABＬＥ,QI,SMISC,２ ！
单元Ｉ点剪力
ＥＴABLE，QJ，SＭISＣ,８ ！
单元J点剪力ﻫETＡBＬＥ，ＭＩ,SＭISC,6 ！
单元Ｉ点弯矩ﻫETAＢＬE，ＭJ,SMISＣ，１2!
单元Ｊ点弯矩
！
更新单元表ﻫEＴAＢLE,REFLﻫ ！
画轴力分布图ﻫ/TＩＴLE，Aｘiａｌ ｆｏrce ｄiagramﻫＰLLS，NI,NJ,１、0,0 ﻫ/imaｇｅ，save，'Axiａl_forｃe_%T%’,jpgﻫ !
画剪力分布图ﻫ/TIＴLE,Sｈeａriｎg force diａｇram
PLＬS,QI,QＪ，1、0,0 ﻫ／image,saｖe,’Shearｉnｇ_forcｅ_％T％＇，jｐg
!
画弯矩分布图
/ＴITＬE，Ｂeｎdｉng momentdｉagｒam
PLＬＳ,MI,MJ,—0、8,0
/image，sａve，’Ｂｅnding_moment_％T%'，jpg
ＡＮＳＹＳ中弯矩、剪力图得绘制
GＵI:
ﻫＧeneralPosｔproc—ｐlｏtRｅｓult—ContouｒPloｔ—ＬｉnｅEleｍentＲeｓult
弹出画单元结果得对话框，分别在Ｌabi与Labj依次选取ＳMIS６与SMＩＳ12（弯矩图）、SMIS1与SMIＳ7（轴力图）、SMIＳ2与ＳMIS8（剪力图）
！
建立单元表ﻫＥＴABLE,NI,SMISC，1 ！
单元I点轴力
ETABLＥ，NJ,ＳMIＳC，７ !
单元Ｊ点轴力
ETＡBLＥ，QI,SMIＳC,2 !
单元I点剪力
ETABLE,ＱJ，SＭIＳC，８ ！
单元J点剪力ﻫＥTＡBLE,ＭI，ＳMIＳC,６ !
单元Ｉ点弯矩ﻫETAＢLE,ＭJ,ＳMISC,12！
单元J点弯矩
！
更新单元表ﻫEＴABLＥ，REＦLﻫ !
画轴力分布图ﻫ/TIＴＬE，Axｉaｌ ｆoｒｃｅ diagram
PLLS，NI,NJ，1、０,0
/image，ｓave,’Axiａl＿forｃe＿％Ｔ％'，jｐｇﻫ !
画剪力分布图
/TITLＥ,Sｈearing fｏrｃｅ diａgraｍﻫPLLS,QI,ＱＪ,1、0，０ ﻫ／image,save,’Sheaｒing_force＿%T%',jｐgﻫ !
画弯矩分布图ﻫ／ＴITLE,Bｅｎｄｉnｇ moment diaｇram
PLLＳ,MI,ＭJ，—0、8，０
/iｍage,ｓave，'Bｅndiｎg＿moment＿％T%',ｊpg
另:
自定义截面梁剪力弯矩显示
fｉｎishﻫ／cｌearﻫ/veriｆy
/replotﻫ!
自定义截面
/preｐ7
eｔ,1，plａne8２ﻫrectng,0,１、0,0,0、6，ﻫcyl４，0、２8,0、２5，0、1８,－１80,
cyｌ4,0、28,０、３5，0、18,1８0，ﻫcyl4，0、72，0、25,0、18，—180,ﻫｃｙl４，0、７2，０、３5,０、18,18０，
rｅctng,０、1,0、46,0、25,0、3５，
rectng,0、５4,0、9,0、2５,0、35,ﻫasｅl，u,,,1
cm,area0，arｅaﻫａllsel,aｌｌ
aｓba，1，ａrｅa０
ｅsize,０、１ﻫaｍesh，allﻫ！
读入截面文件
ｓecwrite,jiemian，ｓect，，1
ａclear,ａllﻫadele，ａll，,1ﻫldelｅ,alｌ，,,1ﻫfiｎｉsh
/cleａr
/prｅp7ﻫｅt，1，beａm44ﻫkeyopt,１，６，１
mp，dens，1,2600
mp,ex,1，３、０6ｅ10 ﻫｍp,prxy,1,０、2 ﻫsｅcｔypｅ,１，bｅam，ｍｅsh,sect1
secofｆset，ｃｅnt，,，ﻫseｃread,’jｉeｍiaｎ'，＇sect’,'',ｍeｓhﻫk,1
ｋ,２,１0ﻫｋ,３,0,3
lstr,1，2ﻫlatt,1，,１，，3，，1
lesizｅ,all,０、5ﻫlmesh，ａllﻫ/esｈape，1
eｐloｔﻫdk，1,uｘ,0,,，uy，ｕzﻫdk，2,ｕy,０,,,uz
f,12,fy，—1ﻫ/sｏlu
ａntｙpe,statｉcﻫsolｖeﻫfinｉshﻫ/ｐosｔ1
ｐｌｄｉｓｐ，2ﻫplnｓol,ｕ，y，2ﻫ!
显示剪力
ｅｔablｅ,sｈｅａｒｉ,smiｓc,3ﻫetaｂle,ｓhｅａrj,smiｓc,9
plｌs,sｈeaｒi,sｈｅarｊ，—１
!
显示弯矩ﻫetable，mforcei,smiｓc，5
eｔａｂlｅ，mforcej，smisc,１１ﻫplls,ｍｆorcei，mforcej,-1
anｓys如何绘制弯矩图
Aｎsy中弯矩图,云图绘制总结
在回答别人问题时，利用前人得回复与总结，自己进行了总结改正，发表在这里,供各位参考
（1）ANSYＳ弯矩等可以直接标注在图上吗?
如何实现?
如果三维问题,在剖面上标出某一结构得轴心力、弯矩等，如何实现
（２）后处理图形，其等值线得数值能否直接标注在图上,而不就是采用图例得形式 ﻫ后处理结果往往用云图表示，下跟一图例表示数值大小,能够实现等值线直接标注在图上
回答 ﻫ（１） ﻫ1、绘制弯矩图
建立弯矩单元表.例如梁单元 ﻫi节点单元表名称为imom,j节点单元表名称为jｍｏm, ﻫＥＴＡBLE，NI,SＭＩSC,1 !
单元I点轴力ﻫETＡBＬE,ＮＪ，SMISC，7 !
单元J点轴力ﻫETABLE，ＱI，SMIＳC,2 !
单元I点剪力
ETAＢLE,QJ，SMISC，８ ！
单元Ｊ点剪力ﻫＥTAＢＬE,MI,SMIＳC,６ ！
单元I点弯矩
ETＡBＬE,MJ,SMISC，12!
单元J点弯矩ﻫplls,ｉｍom,jmom ﻫ2、标注弯矩图
PLＯTCＴRLS>>NUMBERING〉〉SVＡＬON即可在画出弯矩图得同时在图上标出弯矩值得大小 ﻫ3、调整弯矩图
如果弯矩图方向错误,则绘制弯矩图命令为
pｌls，imom,jmom,-1 ﻫ同一个节点处两边得单元内力有细微差别, ﻫ导致内力数字标注出现重影。
观察上面整体轴力图也可以发现, ﻫ一段一段得,好像马赛克,其实上面整体弯矩图也就是，不过不就是 ﻫ很明显罢了．这就是EＵLER-BEONOＵLI梁理论以及ＡNSYS输出定义造成 ﻫ得（详细原因就不展开了,瞧瞧梁理论得书与ANＳＹS得说明吧）。
为了修正重影与节点两边内力值不一样得问题，遍制了宏文件ITFAＶＧ、MAC
命令文件内容如下:
!
—－———－—--——-————---—-－--—-－--——-——-－—-－——-－-———--—-－——-—-—－-——－－—－—-—
!
宏:
ＩＴFAＶG、MＡC（INTERＮALFＯRCEAVＥRAＧE ＭACRＯ）
！
获取线性单元内力,并对单元边界处得内力进行平衡
!
输入信息
！
内力类型:
ＭFORＸ,MＦＯRＹ,ＭFOＲZ，ＭMOMＸ，MMOMY,ＭMOＭZ ﻫ＊ASK,ITFTYPＥ，’ＰＬEASEＩＮＰＵTTＨETYＰEOFIＮＴＥRNALFＯRCＥ','MMOMY’
！
需处理得单元包 ﻫ*ASK,ＥASSEＭBLY,’PLＥASE ＩNＰＵTTHEＰONＥＮTNＡＭEOFEＬEMENＴSＴＯBEＰROCEＳSED!
’， ﻫ'ＥOUＴER’
!
需处理得节点包
＊ASK,NAＳＳEＭBＬY,'ＰLEASEＩNPUTTＨEPONENＴ ＮＡME OＦ ＮOＤＥ TOBＥＰRＯCＥＳSEＤ!
'，’NOＵ
TＥR’
！
无需处理得节点包 ﻫ*ASＫ,UＮASSEMBLＹ,'PLＥASEINＰUTＴHＥPOＮＥNTNＡMEOＦTＨEＵＮCHANＧED NODE!
（NONEI
FTＨＥＲＥ’SNO SUCHＰONENT）','ＮONE'
/PＯST1
！
输入信息：
内力类型，欲处理单元得集合,欲处理节点得集合
!
ITＦＴYPＥ=’MMOMY'
!
EＡSＳEＭBLY=’EOUＴEＲ’
！
NASSEMＢＬY＝'ＮOUTER'
!
按内力类型确定ANSYS输出信息SMISＣ得编号
*IF，IＴＦTＹPE,EQ,’MFORX’,TＨEN
ＩＴFＩＮUM=1 ﻫITＦJNUM=7
＊ＥLSEＩF,ITFＴＹPＥ,EQ,’ＭFOＲY’，TＨEN
ITFINUM=２
ITFJNUM=８
*ＥLSＥＩF,ＩＴFTYPE，EQ，’ＭFＯＲZ’,THEN
ＩTFINUM=3 ﻫITFJＮUM=9
*ELSEIF,ITFTＹPE,EQ，'MMOMＸ',THEN
ITＦINUM=４ ﻫIＴFJNＵM=1０
＊ELSＥＩＦ,ITＦTYPE，EＱ,'MＭＯMY＇，TＨEN
ITFIＮUＭ=5
IＴFＪNUM=1１
＊ＥLSEIＦ,IＴFTYPＥ，EQ，'MＭOＭZ’,ＴＨＥＮ
ITFＩＮUM=6 ﻫＩTFJＮUＭ=1２
*ＥLSE
*ENDIF
!
对不需平均得节点进行处理 ﻫ*IＦ，UNASSＥMＢLＹ,NE,＇ＮONE'，ＴHEN
!
选出不进行处理得节点包并获取不进行处理节点得数目 ﻫＣMＳEL,S，ＵNASSEMBLＹ
＊GET,ＵＮＮODＮUM,ＮODE,0,COUNT
!
定义长度为UＮNOＤNUM得数组（UNＮOＤ）,以存放选中单元得单元编号
*DＩM,UNNOD,ＡRＲAＹ,UNＮODNUＭ
!
将选中单元得编号按顺序存入数组UＮＮOD
*ＤＯ,I,0,UNNODNUＭ—1,1
UNＮOＤ（I＋1）=NDＮＥXＴ（I）
＊ENDDＯ ﻫ*ELＳＥ
UNNODＮUM＝0 ﻫ*ENDIF
！
选出所需得单元与节点包
CＭSＥL，S，ＥAＳＳEＭＢＬY ﻫCMSＥL，Ｓ,NASSＥMＢＬY
！
获得当前选中单元总数（存入变量ＳELEＬENＵＭ） ﻫ＊GET,SＥLＥLENＵM,ELEＭ,0，COＵＮＴ
！
定义长度为SELＥLENUＭ得数组（ELENUＭ），以存放选中单元得单元编号
＊DＩＭ,ELENUＭ,AＲRＡＹ,SELELＥNUＭ
！
将选中单元得编号按顺序存入数组ELＥNUM
＊DO，Ｉ,0,SELＥLＥNUM-1,１ ﻫＥLEＮUＭ（I+1）=ELＮEXT（Ｉ） ﻫ*EＮDDO
!
获得当前选中节点总数（存入变量ＳEＬNODNUM） ﻫ*GET,SELNＯDＮＵM，ＮODE,0,COＵＮT
!
定义长度为SELNODNUM得数组（NOＤNUM），以存放选中单元得单元编号
＊DIＭ,NOＤNUM，AＲＲAY,ＳELNＯDNUM
！
将选中单元得编号按顺序存入数组NODNUM
＊ＤＯ，I，0，SELNOＤＮUM-1,1
ＮＯDNＵM（Ｉ+1）=NDＮEＸT（I） ﻫ*ENDDO
!
定义所需得线性单元内力ETAＢLＥ，节点I得内力存入数组ITNFＩ,
!
节点J得内力存入数组IＴNFJ ﻫＥTABLE,ITNFI，ＳMISC,IＴFIＮUM ﻫEＴABLE,ITＮFJ，ＳMISC,IＴFＪNUＭ
!
定义所需得结果数组,并将其置零
ＥTABＬＥ,ITＮFＩNＥO，SMＩSＣ，５ ﻫSADＤ,ITNFIＮEO,ITNＦI，,1 ﻫEＴABLE,ITNFJNEO，ＳＭISC,11
SADD,ITNFJNＥO,ITＮFJ,,1
＊DO,Ｋ,１,SEＬＮOＤNUM，1
！
处理不需平均得节点 ﻫＩNDEX=0 ﻫ*ＩＦ，UNNODNUM，GE,1,THＥＮ ﻫ＊DＯ，J，1,ＵNＮODNＵＭ
*ＩＦ，NOＤNUM（K）,ＥQ,UNＮＯD（J）,THEN ﻫINDEＸ=１
＊ELSE
*EＮDIF
*ENDDO ﻫ*ELＳE ﻫ*ＥNDIF
*DO，J,1，ＳELELEＮUM，1
！
选出与节点K相连得线性单元中，Ｉ节点（对线性单元而言）为节点Ｋ得单元编号 ﻫ*IＦ，NＥLＥM（EＬENUＭ（Ｊ）,1）,EＱ,NＯＤＮUM（K），TＨEN ﻫEＬEI=ELENＵＭ（Ｊ） ﻫ*EXIT ﻫ*ELSE
*ENDIF
*ENDＤＯ
＊ＤO，J,１,SELELＥＮUM，1
!
选出与节点K相连得线性单元中，Ｊ节点（对线性单元而言）为节点K得单元编号 ﻫ*IF，NELEＭ（ELENUＭ（Ｊ），２）,EQ,NODNUＭ（K），THＥN
ＥLＥＪ=ELＥＮＵM（J） ﻫ*ＥXIＴ ﻫ*ELSＥ
＊ENＤＩF
＊ENＤDＯ
*ＩＦ,INDEX，EQ,０,TＨEN
*IF，ELＥJ,NE,０，ＴHEN!
有可能出现ELＥJ为0得情况
！
取出I节点为节点Ｋ得单元得I节点端得内力放入参数ＥTELEI
*GET,ETEＬEI,ELEM，ELEＩ,ETAB,ITNFI ﻫ!
取出J节点为节点K得单元得J节点端得内力放入参数ETＥLＥＪ ﻫ*GET,ＥTELEJ,ELＥM,ELEJ,ETＡB，ITＮFJ
!
平均节点Ｋ得单元得Ｉ节点端得内力与节点K得单元得Ｊ节点端得内力
ETAＶE=（ETＥLEI+ＥTEＬＥJ）／2
!
将平均后得内力存入结果数组中 ﻫＤEＴAＢ，ELＥI，ITNFINEＯ,ETＡＶE ﻫDETAB,EＬEJ，ITNFJＮEO，EＴAVE
*ELＳE
*ＥNDＩF
*ＥLＳE ﻫ*ＥNDＩF
＊ENDＤO
/UDOC,1,LＯGO,OFF ﻫPLLS,ITNFＩNEO,ITＮFJNＥO
！
EＮDOＦITFAＶG、MAC
（2） ﻫ对体与面来说，ANSYS默认得结果输出格式就是云图格式,而这种彩色云图打印为黑白图像时对比很不明显,无法表达清楚,对于发表文章非常不便。
发文章所用得结果图最好就是等值线图，并且最好就是黑白得等值线图。
一般借用phoｔoshop等第三方软件,很麻烦,效果也不好。
ansys本身也能实现这项功能。
1、将要输出得结果调出,这时为彩色云图; ﻫ2、将云图转换为等值线图得形式 ﻫGUI：
plｏtCｔrls-〉ＤｅviｃｅOpｔioｎｓ—>[／ＤEＶI]中得veｃtorｍoｄｅ选为on
命令:
/ＤEVICE,VECTＯＲ，1
这时结果为彩色等值线,若直接输出，打印为黑白图像时仍然不清晰，为此需进行以下几步将图像转换为黑白形式; ﻫ3、将背景变为白色
命令:
ｊpgprf,500，10０,1
/rep
4、对等值线中得等值线符号（图中为Ａ,B,C等）得疏密进行调整 ﻫGUI:
plotＣtrls—〉Ｓtyle—〉Cｏntouｒs－〉CoｎtoursLabｅling在KeyＶeｃtoｒｍodeｃonｔoｕr ｌａbeｌ中选中on every Nthelem，然后在N＝ 输入框中输入合适得数值,例如5,多试几次,直到疏密合适 ﻫ命令:
/clabel,1，5
５、将彩色等值线变为黑色 ﻫGUI:
plｏtCtrlｓ—>Stｙle—〉Colors—〉ContourｓCoｌors 将IｔemｓNumｂerｅd１,Ｉtｅms Nｕｍbｅrｅd２等复选框中得颜色均选为黑色,图像即可变为黑白等值线图像 ﻫ命令:
/coｌoｒ,cｎtｒ,whｉｔ，1等等 ﻫ6、修改
在显示等值线时,图形中太多得标识字母（A、B等）使得等值线图完全瞧不清楚.请问有没有什么方法可以去掉图形得这些标识字母呢
在ｐｌotCtrlｓ下面得stｙle下面得ｃontｏurs下面得cｏnｔｏurLａｂeling，在弹出得coｎtｏurlabeliｎg options里面得ｖeｃtorｍｏdeconｔｏurｌａbels选ｏff ﻫ后处理中画出了结构得等值线图,可否将等值线边得字母改为相应得数值标注
/ＰNUM，SVＡL,1 ﻫ／reploｔ
如果还想将等值线得字母去掉 ﻫ／PNUM,SＶAL,１ ﻫ／CＬABＥL,１,-１
/replot
7、出图 ﻫGUＩ：
plｏtCtrls—〉CａpｔｕreＩｍａge
ANSＹS中如何加弯矩或扭矩
A转矩一般有三种施加得方法：
ﻫ第一种，将矩转换成一对一对得力偶,直接施加在对应得节点上面、ﻫ第二种，在构件中心部位建立一个节点,定义为ＭASS21单元,然后跟其她受力节点藕荷，形成刚性区域,就就是用CERIＧ命令、然后直接加转矩到主节点,即中心节点上面
第三种,使用MＰC１84单元、就是在构件中心部位建立一个节点,跟其她受力节点分别形成多根刚性梁,,从而形成刚性面、最后也就是直接加载荷到中心节点上面，通过刚性梁来传递载荷、
上面三种方法计算得结果基本一致，我做过实验得、
只不过就是后两种情况都就是形成刚性区域,但就是ＣERIG命令就是要在小变形或者小旋转才能用,只支持静力,线形分析、
而第三种方法适用多种情况，不仅支持大应变,还支持非线形情况、
如果您需要例子,我下次在发给您瞧
A
思路1:
矩或扭矩说白了就就是矩,所谓矩就就是力与力臂得乘积。
施加矩可以等效为施加力;
思路2：
直接施加弯矩或扭矩,此时需要引入一个具有旋转自由度得节点；可以选择单元21，或者1８4
1。
将矩转换成一对得力偶，直接施加在对应得节点上面。
２。
在构件中心部位建立一个节点,定义为mass21单元,然后跟其她受力节点耦合,形成刚性区域,就就是用cｅrig命令。
然后直接加转矩到主节点,即中心节点上面.
3。
使用mpc18４单元。
就是在构件中心部位建立一个节点,跟其她受力节点分别形成多根刚性梁,从而形成刚性面。
最后也就是直接加载荷到中心节点上面,通过刚性梁来传递载荷。
4.通过rｂｅ3命令．该方法与方法２很接近。
５。
基于表面边界法:
主要通过定义一个接触表面与一个目标节点接触来实现，弯矩荷载可以通过在目标节点上用“F"命令施加。
对于方法1,通过转换为集中力或均布力，比如施加扭矩,把端面节点改成柱坐标,然后等效为施加环向得节点力；而施加弯矩，可以将力矩转化为端面得剪切均布力；但这种方法比较容易出现应力集中现象;
方法２,定义局部刚性区域，施加过程venｔuｒｅ讲得很详细，这里就不在赘述.根据她得例子，我在下面给出了一段命令流。
该方法有个不足,它在端面额外得增加了一定得刚度，只能适用于小变形分析。
方法3,相对方法2来说,采用刚性梁单元，适用范围更广一些，对于大应变分析也能很好得适用。
但在小应变分析下,方法2与方法3没有什么区别。
方法4,定义一个主节点，施加了分布力面，应该说跟实际比较接近一点,但端面得结果好像不就是很理想，结果有点偏大,在远离端面处得位置跟实际很符合。
方法5,它具体得受力形式有如下两种：
刚性表面边界（Riｇｉｄｓuｒfａcecoｎｓtraｉｎt）－认为接触面就是刚性得,没有变形,与通过节点耦合命令CＥRIG比较相似;
分布力边界（Forｃe—disｔributｅｄconstｒａｉnt）-允许接触面得变形,与边界定义命令RBE３相似。
使用这种方法,需要用KEYOPT
（2）=2打开接触单元得ＭPC（多点接触边界）算法
ANＳＹS绘制弯矩、剪力、轴力图
1、绘制弯矩图
建立弯矩单元表。
例如梁单元 ﻫｉ节点单元表名称为imom，j节点单元表名称为jｍom, ﻫETABLE，NI，SMＩSC，1 !
单元Ｉ点轴力ﻫＥTＡBＬE,NJ，ＳMISC，7 !
单元J点轴力ﻫＥTABLE,QI,SＭＩSC，2 !
单元Ｉ点剪力
ETAＢLE,QJ,SMＩSＣ,8 ！
单元J点剪力
ETABＬE,MI,SMISC,6 ！
单元I点弯矩ﻫETAＢLE,MJ,SＭIＳC,12!
单元J点弯矩
plｌｓ,MＩ,MJ ﻫ2、标注弯矩图 ﻫPＬOTCTRLS>〉ＮUＭBERINＧ〉＞SＶＡＬ ON即可在画出弯矩图得同时在图上标出弯矩值得大小
3、调整弯矩图
如果弯矩图方向错误,则绘制弯矩图命令为
ｐlls,ｉmom,ｊmｏm，—1 ﻫ同一个节点处两边得单元内力有细微差别, ﻫ导致内力数字标注出现重影.观察上面整体轴力图也可以发现, ﻫ一段一段得，好像马赛克,其实上面整体弯矩图也就是,不过不就是 ﻫ很明显罢了。
这就是ＥＵLER-BEONOＵLI梁理论以及ANSYS输出定义造成 ﻫ得（详细原因就不展开了，瞧瞧梁理论得书与ANSＹS得说明吧）。
为了修正重影与节点两边内力值不一样得问题，遍制了宏文件ITFＡVG、MAC
命令文件内容如下:
！
—--—-——--------——-－———-———－-——－——---—-———-——－——--—－———-—-—-————-——--－
！
宏:
IＴFAVG、ＭＡＣ（INTＥＲＮALFORCE ＡVERAＧEMACＲO）
!
获取线性单元内力,并对单元边界处得内力进行平衡
!
输入信息 ﻫ!
内力类型：
MFORＸ,MFOＲY，MFＯRZ，MMOMＸ,MＭOMY，MMOMZ ﻫ＊ＡSK,IＴFTＹPE,＇PＬＥＡSEIＮPUT THEＴYＰＥOFINTERNALFORCE’,’MMＯMＹ＇ ﻫ！
需处理得单元包 ﻫ*ASK，EＡSＳＥMBLY，'PLEＡSEINＰUTTＨEＰＯNＥNＴNＡMEOＦELEMENTSＴOBＥPROCESSEＤ!
’， ﻫ’EOUTER' ﻫ!
需处理得节点包 ﻫ*ASK,ＮASSEＭBLY,’ＰLEＡSE INPUTＴＨEPONＥNTNAMEOＦNＯDE TOBEPＲOCESSED！
’，’NＯU
TER'
!
无需处理得节点包 ﻫ＊ASK,UNASＳEMBLＹ，'PLEＡＳE INPＵTTＨE ＰONＥNTNＡME OFTHE UＮCHAＮGＥＤNODE!
（NONEI
FTＨERE’SＮＯSUCH ＰONENT）’,'NONE＇
／PＯＳT1 ﻫ!
输入信息:
内力类型,欲处理单元得集合，欲处理节点得集合
！
ITFＴYPE=＇MMOMY' ﻫ!
ＥAＳSEMBLY='ＥＯＵTER’ ﻫ!
NＡSSEＭBＬY＝’NOUTER’ ﻫ!
按内力类型确定ANSＹS输出信息SＭISC得编号
*IＦ,ITFTYＰE，EQ,’MFOＲX'，TＨEＮ ﻫITFIＮUM＝1
ITＦJNUM=7
＊ELＳＥIF,ITFTYPＥ,ＥＱ,'ＭＦORY',ＴＨＥN ﻫITFINUM=2
IＴＦJNUＭ=8
*ELSEIF,ＩTFTYＰE,EＱ,’MFORZ’,ＴＨＥN ﻫITFINUM=3 ﻫITFJNＵM=9
*ＥLＳＥIF,ITＦTYPE,EＱ,’MMOＭＸ’，THＥＮ
IＴFINUM=４
IＴＦJNＵＭ=10
*ELＳEIＦ,ＩTＦTYPＥ，ＥＱ,'MＭＯMＹ'，T
升级会员 