理解Femap的后处理

assassin_007

我们费了那么多的时间终于做出了一个模型，可是结果怎么样呢？
' _8 i. z% A7 W  a, \4 X" P- ~我们如何评判结果呢？
% g; I0 @' t6 ^, X: E, B, |在Femap的结果中又有那么多的选项？又是节点，又是单元，还有最大最小，平均值，还有什么Corner？
我该怎么来选择呢？
这里我们会以一个静态分析为例，很容易的，这样大家就可以理解后处理的作用，以及如何评价你的有限元模型。

assassin_007

现在我们看下面的这个例子：

这里面有四个六面体单元，单元的边长为1.什么单位？没有单位的！Femap是无量纲计算，忘了吗？
. X; J4 I3 x6 S- K: F约束：单元1（最右面）下面的四个节点为固定。载荷：单元1上表面8个点上，四个角点为15，四个中间点为10，合力为100
计算结果的考察（节点17，上表面正中的节点）：

assassin_007

打开.f06文件，查看节点17的应力数据
& @% j3 L% G* B6 W  t/ L

单元 8 T2 V0 h6 D" }		应力
1	Center " _3 y$ U. }" r* s( _0 y	41.72927
	G17	122.8854
	Max（G3） : Q. L  _0 F0 b2 O	318.3802 4 \$ A( F2 b  t
	Min（G1） & q9 [1 Z5 M/ @( p7 f( y6 Z	36.97896 2 i8 ?9 C3 U! Q0 L) ~3 w9 u
2 " y: ^  d. l5 G, g; K# ?* O	Center	1.731009
	G17	84.71432 % n( H6 y( c' ?0 R, n, E! C& Q& G
	Max（G6）	113.5346
	Min（G16）	14.36555
3	Center ( F9 T- q+ P! R! G8 n	1.731009
	G17 / ~" A& M5 }& z9 k3 S: c	84.71432 : A6 X' e" s; R+ D& c
	Max（G6）	159.1384
	Min（G10） / R' ^5 O1 O: [* E! X8 y	14.36555
4 4 g7 H5 u1 n  C# X	Center	1.527067
	G17	66.49532 ; v4 ^, @  I: D2 S7 {4 `' }5 k$ ^
	Max（G17）	66.49532
	Min（G15） ) I& ^' h0 Q3 T	11.21177 . q- q% Y& c7 K  K

assassin_007

怎么了，找不到.f06文件？很简单呀，在你的Femap创建的文件保存的地方就可以了。
找到后用文本文件打开，里面就是计算的结果了，找到这一行：
S 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 )
% O, O- [6 S3 u- l$ v下面就是计算的应力结果了。什么，看不到这个？没有关系，有一点你肯定是忘记了。
/ @" b' a4 U7 u在Femap里面创建分析集的时候，最后一步，对话框的右下角是什么？
+ s' P/ ?8 b! n) n9 O2 Q2 b结果目标嘛！
默认的是仅仅用于后处理，那么我们要选择打印和后处理，也就是第四项。
再运行一遍，肯定有了

assassin_007

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

5 i# \& ]8 K6 z4 \9 `2 R. v单元1的应力值最大，为41.72927
不是吗？
好了，我们现在来看看到底怎么回事。

assassin_007

鼠标点右键，选择Post Data，对话框上面再选择Countour Options，看见没有，选择Nodal，Average，Use Corner Data。
" ?) O$ r6 ~, d7 \3 s; B8 Y好了吗？
1 p$ t& @# X* V: B* I  sOk，再次Ok
8 X; y5 ?) @& l& d' x8 n9 r( w看一下屏幕啊，
最小值：11.21
9 K9 A- ~, }5 c最大值：381.4
节点17上面的值呢？
/ L* q8 J! r1 _找到选择工具条，第一个按钮，下面选择Element，第二个按钮，下面选择Show ToolTips，
然后呢，把鼠标放到第17个节点上吧，显示的值就是89.70234。
3 e$ `* L( h8 R. n2 M; a( Q! k5 b$ e那么，这个值代表什么意思呢？为什么上边让大家从f06文件里把那些数据抄出来呢？
算一下吧，就是节点17上面的四个值的平均值。
& L* ^6 f: y7 e" c! l  i明白了吗？

assassin_007

下面，我们就不那么罗嗦了，大家挨个试吧，得到如下结果：
5 K) k7 [! o5 ?/ L/ V
Nodal， Average
Min：1.527
Max：41.73
G17：11.67959
$ X* O5 O7 z# ^1 iG17的值为四个单元中心的应力的平均值

assassin_007

5.        Nodal，Max，Corner Data
, s  U& M/ E% o9 @  @Max：381.4
Min：11.21
: b' `- [/ S/ d  M+ lG17：122.8854
G17的值为四个单元上G17的最大值
$ X6 J1 a6 T2 H2 p- }- o6.        Nodal，Max
- Y# Y$ E+ O8 L$ J+ g/ bMax：41.73
' w& y9 O" n& J+ `6 {9 C" BMin：1.527
- }( I2 S% s6 P" R  O/ |  a) xG17：41.7297
G17的值为四个单元的应力值中最大的应力
7.        Nodal，Min，Corner Data
Min：11.21
4 x- p7 r2 [7 eMax：381.4
G17：66.49532
G17的值为四个单元上G17的应力值中的最小值
8.        Nodal，Min,
% m6 r- k0 H5 N9 H% H' k1 FMax：41.73
Min：1.527
G17：1.527
G17的值为四个单元的应力值中最小的值

assassin_007

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

assassin_007

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

assassin_007

希望大家喜欢这个结果！

assassin_007

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

heyantian

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

狼跃冲顶

嗯，看来这个软件不错
2 v+ n* [8 o7 w/ V有不少人已经开始用它了

jimfly

没记错的话，Femap是第一支拥有Windows化下拉menu的 CAE软件，
& d/ g) m) ^* b4 P# m当时就支援那时期的多支CAE软件，如Ansys、NASTRAN(MSC & NE)、Abaqus
1 s  d, g( U0 ~4 w(可能当时是在windows下直接发展，MSC还曾经以此出PC版NASTRAN)，
哪个时期，许多知名软件都是以UNIX起家，
9 b0 ]* s. z4 `7 ^7 |( l3 o$ s5 _如MSC/Patran、Ansys、Ideas、Abaqus及HyperMesh，都是以X-Windows权充，
可惜至今 Femap虽然在画面操作上有进步(v9之后有CAD的影子)，
但是在mesh功力上好像进步不多，而其他产品也陆续推出正统PC版，
$ Y  v% o$ O, {' ]导致 FEMAP换了几手老板，User好像也渐渐流失了
# x; f# s( y) X7 E2 e  {6 P
. b; L/ G2 K/ v4 j8 r感谢分享

heyantian

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

hyping0218

楼主好久没来了

saaya

是呀，楼主哪去了

banlizl

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

hellowww

谢谢，很好，希望继续

xxdong

楼主忙着抱孩子那
/ C! E3 ]0 t/ o9 S( A1 G. P3 m8 h- F哈哈
6 w$ H* v& K. C  O) @内部消息，嘘嘘 titter titter

bosieng

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