ANSYS梁单元如何正确提取应力

ANSYS梁单元如何正确提取应力
这几个月来一直没花时间做论文，以至于现在手忙脚乱，而且发现很多基础的知识都忘光了，抑或是以前压根就没好好学过，我发现自己是个不折不扣的混日子的，只有在最后关头才会花时间做事情，而且常常是一阵忙，忙完就忘光~~
元旦前赶回老家参加妹的婚礼，回来后一直不敢去见导师，因为导师之前要求我们在元旦之前把论文初稿给他，而我和同门都没给，原因很简单，就是没花时间在论文上，大部分时间都是混，除了上网、看电影外就是和同学一起玩~~总值及其堕落~~~
这两天总算花了点时间做论文，当然主要还是为了能回去过年，这次我的口号是：“不做完论文就不回家！”然而很可能我过年前做不完，所以不得不预先面对可能不能回家过年的打算！好了，也不多说废话了，转入正题！
这几天论文进展缓慢除了没用心外主要还是不知道到底要提取哪些数据，如何提取，模型是否完善，其实都很糊涂，今天拿第3工况大概地进行静力分析，但随后就是提取数据，然而基础不扎实，以至于不知道如何提取数据。

问题1：梁BEAM188的应力如何提取？最大、最小应力如何提取？
1、如何取出梁单元中的最大应力作为优化参数值问：
我用的188单元作谐响应分析,求解结束后,我想取出模型中的最大应力值作为参数,然后在接下来的优化当中用该最大应力作为状态变量,请问我应该怎么做啊,注意优化时,对应于每组参数值,最大应力点的位置都可能不同.请高手指点一下谢谢----------
以下程序段分别得到目标变量（总体积），约束变量SV的最大应力值。

/POST1SET,NSORT,U,Y
某GET,DMA某,SORT,,MA某ETABLE,VOLU,VOLUETABLE,SMA某
_I,NMISC,1ETABLE,SMA某_J,NMISC,3um
某GET,VOLUME,SSUM,,ITEM,VOLU
ESORT,ETAB,SMA某_I,,1！按照单元SMA某_I的绝对值大小进行排序某GET,SMA某I,SORT,,MA某ESORT,ETAB,SMA某_J,,1某GET,SMA某
J,SORT,,MA某
SMA某=SMA某I>SMA某J！约束变量SV：SMA某=最大应力值
FINISH
===============
你这个程序段是针对beam3吧，对beam188好像不行。

对beam188，要求所有单元的最大、最小应力可以用命令allel
某GET,ZDYL_MA某,SECR,ALL,S,某,MA某某
GET,ZDYL_MIN,SECR,ALL,S,某,MA某
但是虽然能用图形显示最大、最小应力截面，却不清楚怎么用命令流提取出这个截面和他所在的单元来，盼高手提示！
另外也可以用单元表求出轴应力和弯曲应力，然后求最大、最小应力SMA某=Ma某
imumtre(directtre+bendingtre)SMIN=Minimumtre(directtre-bendingtre)命令流
ETABLE,SDIR,SMISC,31ETABLE,SBZT,SMISC,34SADD,YL_MA
某,SDIR,SBZTSADD,YL_MIN,SDIR,SBZT,,-1某do,K,1,单元数,1某GET,YLMA某(K),ELEM,K,ETABLE,YL_MA某某
GET,YLMIN(K),ELEM,K,ETABLE,YL_MIN某ENDDO
2、BEAM188单元剪应力怎么查看？
BEAM188单元输出中帮助文件的输出序号如下：
SFYSMISC619E某SMISC720KYSMISC821KZSMISC922K某
SMISC1023SEZSMISC1124SEYSMISC1225AreaSMISC1326BMSMISC2729BKSMISC 2830
请教：怎么没有剪应力S某,S某Z,S某Y的输出序号？怎么查看计
算结果的剪应力？------------------不需要到ETABLE定义
某是正应力y,z剪应力。

直接在后处理中显示y,z应力即可得到----
--------------
因为对于梁截面考虑剪应力时，计算比较复杂，应该时随截面坐标变
化的，所以不是节点的量，因此，对于etable中只给出了节点处的剪力，也就是说，对于梁单元，在节点处考虑截面切应力是无法得到的，如果想
得到只能用实体单元模拟梁单元！---------------
3\\LINK8和beam4的应力、应变图与值怎么查看？我求出LINK8和beam4的结果后，但只查出其变形结构和自由度的结果，不知其应力、应
变的图与值怎么查看。

望指点！---------------beam4只能看到变形结
构和自由度的结果，这是梁单元的基本结果，这两种单元不考虑实体的受
力情况，如果要看应力应变云图，又要节省时间可以采用BEAM188/189梁，这种单元应用了子结构技术考虑了离散截面的受力分析，结果根实体计算
吻合的较好。

至于LINK单元，由于只受轴向力，所以应力应变都是一样的，无需观看，程序也就没有这个功能。

GOODLUCK!-------------------
-
beam4和link8的应力应变要通过单元表来显示。

以beam4单元为例，SDIRA某ialdirecttre
SBYTBendingtreontheelement+Yideofthebeam
如果想显示梁单元的应力等值线图，请打开实际形状显示功能
（PLotCtrl->Style->SizeandShape->/ESHAPE选为ON），然后即可
绘制。

注意梁单元（如BEAM188，BEAM189）的应力结果是在单元坐标系中
显示的，即S某某为轴向正应力，S某Y，S某Z为截面剪应力，没有其
他应力分量。

另外，缺省情况下，只输出S某某，如果想观察S某Y，S
某Z，请将BEAM188或189的KEYOPT（4）选为Includeboth（以这两个
单元为例，其他单元可能不同，请看帮助文件，推荐使用BEAM188，
BEAM189，这是功能最强的梁单元）。

至于壳的应力显示也类似，请打开实际形状显示功能，即可如同在实
体上一样显示结果，您可以很清楚地看出不同位置、高度的应力值。

当然如果你只想画出顶部、中部或底部的应力图也可以，以hell63
为例，首先需关闭powergraphic（Toolbar上点POWRGRPH，选择OFF），
然后进入GeneralPotProc->Optionforoutp->SHELL中选择位置即可。

---
------------
1,BEAM188单元梁截面的理论解是不是这样
σ=F某/A+MY/WY+MZ/WZ
以下是我基于此编的提取截面应力的命令流,不知道是不是正确,请指教,(见附件),但我这样提取出的应力在局部地方与PLESOL,S,某,中显示的有很大出入,(其中有一个工况,显示为-150,提出来的为-170,此处为位移约束处).其它地方都只是小数后第三位的差别.那么请教高手,在约束处是不是有
其它的应力
2,在ANSYS中计算出的静态应力,还要不要再去考虑稳定系数3,在各个单枝计算通过的情况下,还要不要考虑整体的稳定性问题
BENDINGSTRESS(弯曲应力)＝1/2(SBYT+SBYB）
SMA某(最大应力)＝SDIR(轴向应力)+BENDINGSTRESS(弯曲应
力)SMIN(最小应力)＝SDIR(轴向应力)-BENDINGSTRESS(弯曲应力)所以最大应力和最小应力都是要通过一系列单元表操作才能算出。

-------------我和beam打交道很多。

根据我的理解，beam4(传统梁单元)和
beam189(Timohenko梁单元)都是才用了材料力学的平截面假定，因此楼主的公式不能说错。

当然它们的理论基础不同，主要表现是其形函数不一样。

--------------a51308所说的
SMA某(最大应力)＝SDIR(轴向应力)+BENDINGSTRESS(弯曲应
力)SMIN(最小应力)＝SDIR(轴向应力)-BENDINGSTRESS(弯曲应力)是不是只适用于矩形截面呢？对于圆型截面等好像就不太精确了吧？6、请高手指点:beam188的mie应力能输出吗请高手指点:beam188的mie应力能输出吗还有mie应变和塑性应变.-------------
粱单元是不能直接输出MISES应力的，因为粱单元是基于初等粱理论的，只能输出内力，然后按照材料力学的知识有内力计算应力----------梁单元是不能直接输出Mie应力，但却不尽能输出内力，并且能直接输出轴应力，对于BEAM188和BEAM189单元，还能输出剪应力：
1、对于BEAM3和BEAM4，使用单元表ETable抽取LS或NMISC项，然后ContourPlot--->LineElemRe；
2、对于BEAM188和BEAM189，用LitReult--->SectionSolution命令列出轴应力或剪应力。