理解Femap的后处理

assassin_007

我们费了那么多的时间终于做出了一个模型，可是结果怎么样呢？
8 ]' R; t8 h7 H2 c7 Q, d. Y. _2 R我们如何评判结果呢？
0 {4 m7 I3 J8 i) `5 p- L* ]: x在Femap的结果中又有那么多的选项？又是节点，又是单元，还有最大最小，平均值，还有什么Corner？
# z- j5 E4 j& p) v9 \; j  P9 r+ Z我该怎么来选择呢？
这里我们会以一个静态分析为例，很容易的，这样大家就可以理解后处理的作用，以及如何评价你的有限元模型。

assassin_007

现在我们看下面的这个例子：
! I7 |, L- \. `$ H% R+ P
这里面有四个六面体单元，单元的边长为1.什么单位？没有单位的！Femap是无量纲计算，忘了吗？
约束：单元1（最右面）下面的四个节点为固定。载荷：单元1上表面8个点上，四个角点为15，四个中间点为10，合力为100
! P- j+ G& d) ?! @* t/ s计算结果的考察（节点17，上表面正中的节点）：

assassin_007

打开.f06文件，查看节点17的应力数据

单元		应力
1 / l2 ^/ n7 w: B  g& V' f# x	Center * e2 b! R. j: C" U" g	41.72927
	G17 9 S# `+ V2 c- n7 A* z. p" [- ]	122.8854 / C4 n( w  R& ]1 u  m/ k
	Max（G3） % f( P7 e8 ]% i3 a- b) g5 E	318.3802
	Min（G1） 5 D% ]* N" q& d1 w	36.97896
2	Center ) w" R$ A! d$ c; M& q	1.731009 4 G+ Z! e7 K( L6 R
	G17	84.71432 . A' X* a. y! ^' F
	Max（G6）	113.5346 9 b4 s  n5 \' b( ]) \- v
	Min（G16） 5 U5 }8 d" v* c, ^- g* s( j5 f/ N	14.36555 6 [& ^7 ?5 [) Z/ ^2 H: P0 B# j; t
3 $ l" x6 n4 E$ {' o8 W5 j	Center 1 _5 H2 ?+ l  _8 v: W) E	1.731009 & R0 x* V- |8 G: v) G- }# o
	G17 8 [: K2 U  c  D	84.71432
	Max（G6）	159.1384 9 Q$ b, x7 C# p6 E0 x* \
	Min（G10）	14.36555 5 E- H" k6 M9 B
4	Center	1.527067 6 j8 a' L5 ^$ e  W2 m1 V( s
	G17 & U- R6 \8 _( q; Y+ ]+ [5 k	66.49532 - ?0 W& @2 y( F1 V$ k
	Max（G17）	66.49532 9 V. l/ }; o: |* w0 C# r+ _5 u
	Min（G15）	11.21177

assassin_007

怎么了，找不到.f06文件？很简单呀，在你的Femap创建的文件保存的地方就可以了。
: _% I1 @; |7 S% `, z找到后用文本文件打开，里面就是计算的结果了，找到这一行：
5 E6 g( \6 L+ Z  kS T R E S S E S   I N   H E X A H E D R O N   S O L I D   E L E M E N T S   ( H E X A )
下面就是计算的应力结果了。什么，看不到这个？没有关系，有一点你肯定是忘记了。
在Femap里面创建分析集的时候，最后一步，对话框的右下角是什么？
' S& P0 L! R- g/ n3 S1 `结果目标嘛！
默认的是仅仅用于后处理，那么我们要选择打印和后处理，也就是第四项。
; Y- N% A' Z6 ^; T1 l再运行一遍，肯定有了

assassin_007

我们把相关的部分看一下，节点3的应力值最大，为381.3802

单元1的应力值最大，为41.72927
不是吗？
好了，我们现在来看看到底怎么回事。

assassin_007

鼠标点右键，选择Post Data，对话框上面再选择Countour Options，看见没有，选择Nodal，Average，Use Corner Data。
好了吗？
Ok，再次Ok
看一下屏幕啊，
最小值：11.21
最大值：381.4
# Y1 R% H9 f9 B7 n+ e. C7 R节点17上面的值呢？
找到选择工具条，第一个按钮，下面选择Element，第二个按钮，下面选择Show ToolTips，
% a, Y" d$ ]* V6 P% o0 m6 G然后呢，把鼠标放到第17个节点上吧，显示的值就是89.70234。
那么，这个值代表什么意思呢？为什么上边让大家从f06文件里把那些数据抄出来呢？
$ L6 ?! b2 j( w! z7 d, o1 @( ^) m算一下吧，就是节点17上面的四个值的平均值。
9 H) e( o' C+ `: k* U4 l明白了吗？
3 M  g% e7 m7 t+ l

assassin_007

下面，我们就不那么罗嗦了，大家挨个试吧，得到如下结果：

  o/ g0 t6 \! zNodal， Average
4 X! Z, \" \0 e( y/ P6 A, YMin：1.527
Max：41.73
/ b( G2 a" e( e$ {( ?4 oG17：11.67959
& ]' ^9 d. g7 T" I  c. z& n* WG17的值为四个单元中心的应力的平均值

assassin_007

5.        Nodal，Max，Corner Data
Max：381.4
  A' X8 e3 a$ aMin：11.21
G17：122.8854
; O/ Q% s( z$ u  T8 V" OG17的值为四个单元上G17的最大值
% ?* W$ m5 `, }  j2 L' \+ O4 ^3 j6.        Nodal，Max
+ |* }) p3 [( p& G9 ]; RMax：41.73
Min：1.527
1 E' R+ R+ l3 x6 ?- A  |3 h: |& g9 AG17：41.7297
: M5 Y9 c6 \5 P! [% m! nG17的值为四个单元的应力值中最大的应力
7.        Nodal，Min，Corner Data
% k% `. T% G0 z8 u) y/ k* R8 EMin：11.21
6 b5 l7 B3 z3 `& ZMax：381.4
G17：66.49532
G17的值为四个单元上G17的应力值中的最小值
8.        Nodal，Min,
2 k% D' D3 y+ ?* G4 ?8 M: e9 iMax：41.73
. N% v) }$ h  L( D7 Y7 m0 k; kMin：1.527
G17：1.527
8 E) h' n: o% V( k. _8 [1 `" _" R3 V- qG17的值为四个单元的应力值中最小的值

assassin_007

那么，单元的呢？通常情况下，在单元连续，所谓连续主要是单元类型相同，材质相同，选择节点的值和单元的值没有什么大的区别，你会得到和上面一样的结果。然而，如果单元不连续，换句话说，就是单元的类型不同，比如，采用的单元是六面体和四面体混合，或实体单元与板壳单元或线性单元混合，就需要注意了，又或者单元类型一样，但是材料不一样，也需要注意了。特别是各种情况混合的时候，所以呀，一般还是看单元的值吧，要不要Corner Data自己看了！

assassin_007

从这个意义上讲，结果的判定需要注意了。理论上一个节点的值只应该有一个，可是，由于有限元计算采用的是积分的近似计算，导致了同一个节点上会出现不同的值，那么最大值和平均值越接近越好了，对不对？如果差别很大，那么，需要调整一下网格了？差别大的地方密一点吧！也牵扯到了取舍的问题。所以，有限元结果的判定呢，还是需要有一定的工程经验才可以，不是点两下鼠标就OK了。不是吗？

assassin_007

希望大家喜欢这个结果！

assassin_007

十分抱歉，最近一直在忙Solid Edge的事情，有些耽误了，还请大家原谅！

heyantian

很好很通俗，接着来，顶你了

狼跃冲顶

嗯，看来这个软件不错
有不少人已经开始用它了

jimfly

没记错的话，Femap是第一支拥有Windows化下拉menu的 CAE软件，
当时就支援那时期的多支CAE软件，如Ansys、NASTRAN(MSC & NE)、Abaqus
. l- g% a) z5 X" P. H9 Z) x(可能当时是在windows下直接发展，MSC还曾经以此出PC版NASTRAN)，
; E0 x9 P' ^+ y4 b9 z& Y3 l- _哪个时期，许多知名软件都是以UNIX起家，
如MSC/Patran、Ansys、Ideas、Abaqus及HyperMesh，都是以X-Windows权充，
6 Y4 @3 ]  b1 Q& U, J9 s可惜至今 Femap虽然在画面操作上有进步(v9之后有CAD的影子)，
但是在mesh功力上好像进步不多，而其他产品也陆续推出正统PC版，
) N4 ^* O  e& p+ m导致 FEMAP换了几手老板，User好像也渐渐流失了
& }% `- j6 H' Y# T, v
感谢分享

heyantian

很通俗，这样容易掌握些。

hyping0218

楼主好久没来了

saaya

是呀，楼主哪去了

banlizl

楼主讲得很好啊,希望多发这种贴啊!~

hellowww

谢谢，很好，希望继续

xxdong

楼主忙着抱孩子那
7 Y# r5 C) \3 `) C  ]2 s哈哈
内部消息，嘘嘘 titter titter

bosieng

是啊,给楼主发信也不见回,好久不见他了.