找回密码
 注册
Simdroid-非首页
查看: 1414|回复: 25

[软件学习与使用] 疲劳分析中的几个概念及其软件中的用法

[复制链接]
发表于 2009-4-17 21:11:35 | 显示全部楼层 |阅读模式 来自 重庆沙坪坝区
本帖最后由 apprent 于 2009-4-23 12:04 编辑 9 t3 y* ^: R& m3 b! T" I

$ G  D# ?' @( \# Y* j2 z+ g3 Z题目有点大,其实我是来问问题的,而且是很基础的问题。请各位不吝指教。我自己也会继续翻资料学习。
: B( N8 _/ K. h4 A$ F; i: q' m7 ~
~~~~~~~~~~~~~~~~~~~~~~~
; d. t7 Y5 q! l# c讨论背景:  Z# p& {9 F6 T
使用的软件为MSC.Fatigue 2005% {% s' R2 h! D/ a+ Y  t
做S-N分析。; u7 V5 `( }7 }+ f1 P2 {. I3 j
& k( E$ a3 `  u1 a: |5 N) ^2 U
借用浙大的一个PPT中的图(稍有改动):/ v4 Y) y8 l7 X; W" r0 ~, `2 f7 A5 G
+ z2 u& q5 F4 d2 g. D6 ]; T
5 I1 h0 y9 _& K1 z
参考讨论帖:
2 z; w3 y7 L) m& C' t' uhttp://forum.simwe.com/viewthread.php?tid=769124(初学MSC.fatigue问题总结); T; x3 n, J: B0 R
http://forum.simwe.com/viewthread.php?tid=146161(对存活率、M.s-n与C.s-n做了讨论)7 @/ g; t' O3 [- A1 h% n% j( ?0 V8 h0 k
http://forum.simwe.com/viewthread.php?tid=819125(存活率,疲劳分析的意义): y9 a. s2 h7 G
0 m7 L# g# Z* I" j7 @# |5 g
~~~~~~~~~~~~~~~~~~~~~~~) P) y0 h4 x1 d& b. V
问题:9 f  k' a7 o- a0 g7 j

2 [* B3 {7 r& [( X. [, @4 e$ |( V1. 平均应力修正(Mean Stress Correction)+ M! c2 G$ @: C. \. S
4、各应力修正方法中,Goodman简单且偏于安全,在工程上应用较多,Gerber偏于危险,None就是无平均应力修正了,它适用于平均应力为零很小的情况。Max.Abs.Principle是用第一强度理论,von mises是用第四强度理论。4 T% y4 M6 N' G9 A1 {
Beans110 发表于 2007-3-26 22:17

7 y8 n$ U9 i9 O) M5 }7 j3 U那么,当我在MSC.fatigue中加载的载荷谱(力/力矩)的平均值不为零的时候,只要在这里选一个应力修正模型就可以了吗?
% j' }' N7 x+ x是否MSC.fatigue在计算时会先把这个载荷谱进行雨流计数,得到其Range、Mean和Cycle,然后套用所选的应力修正模型?( Q" g. ~2 U& E2 T; Q2 _0 S
, k# e. z6 |% `, S8 f* Z

8 R+ O7 M1 ~8 k" H1 G3 E, T! |4 a/ a4 n- p5 i7 e" ~- j
( B; j* v* ^* w( D0 v" Q$ U# i

2 L1 o7 R/ |7 s3 j' ]9 k( a* P, f& F! L

: F) Z% i+ P3 O# e$ z2. 存活率(Certainty of Survival)2 d6 }: O4 F- \8 J3 f1 ^; H
疲劳试验用了十根试样,其寿命值因材料、试验等原因导致寿命值存在分布,即得到十个寿命值。用此数据预测你的零件寿命时,你用那一个呢?50%即为平均值,即十个试样中,有五个低于此寿命值,有五个高于此寿命值!且在Msc/fatigue中是起作用的!
% u: x* l. O% A! i2 _" h+ vflycomet 发表于 2003-11-14 12:25
这个值的选取不同,得出的寿命还是不同的
! D- y6 u$ e1 C' s我对我们的一个正在进行中的试验,先预测其寿命 ' n1 U  @3 U; j3 p9 M
选96%时,寿命为63万次,选50%为180万次
; ], m9 u" G% i8 f* a- M$ ^, PJerry 发表于 2003-11-14 14:25
从安全性考虑,如果一个如航天构件,生存率取值要高;如果对一个农用运输车而言,生存率则可取低些!总之,根据你的需要进行选取!8 P0 l( Z' w  |7 u
flycomet 发表于 2003-11-14 19:01
7 S% A# P9 ]4 I% f
……- e$ H% q% I9 B* G* o7 H
不同的存活率得到的寿命不同。此时你最好是求出1%-99%的一条分布曲线来。
8 T" X, g3 Y% F* L" tkkkyyy 发表于 2003-11-14 21:22
不能说疲劳没必要,是因为要做好疲劳很难,难在很多数据很难得到,因为要做很多次实验比如S-N曲线,很多都是根据经验得到的,或者直接通过极限抗拉强度,屈服应力估计的,' d4 t. g7 Z1 R# L/ i/ }1 l6 }1 O# ~
实验室里面,用特定的样品形状,施加特定的载荷来做出来的,而且是通过概率统计经过处理,
一般的软件给出的S-N曲线在是存活率50%时候的,但在工程上不允许那么低,起码90%以上,这就要对S-N曲线进行校正
  P" t' W. A  @9 @" A4 r+ d6 uhuanghui01 发表于 2008-3-4 16:09

2 S! {9 p* G; c, _  S# d那么,我这里的计算有要求,存活率使用97.7%,我就在这里直接把存活率设为97.7%就可以了吗?4 s1 b7 y* k2 G1 t3 n, A
存活率修正的公式是什么呢?(既然大家的S-N大都是用UTS和E估计出来的)。前述Goodman、Gerber的公式教科书上是都有的。
3 a* o+ @& y" r+ a8 S- E! x; Y
1 J! L$ H; C3 X4 O9 v
# W3 l; B$ B0 v, w
5 y% I& U! T$ W6 V
3. 取节点或单元的应力
- E, [+ x+ o: u( j
: e8 x+ x6 h( H

, F. x  P. k% Q# e' e  g. n& ]" p" @5 b4 F; ^3 u
这个选项的含义我大致是明白的。主要就是涉及到从积分点应力向节点应力插值的情况。) ~: g/ e2 h; v- d' @7 M8 u( m# p/ A
对于有4个积分点的quad4单元,这里选Node的话,那么使用的是从积分点向4个节点线性外插的应力值。
9 E% X7 u5 y9 }; b' W选Element的话,使用的单元的应力,即4个积分点的平均值。' _8 Q; U- z. t0 x1 s) l
我在一些计算中发现,选择这两种方法计算出来的疲劳寿命相差是比较大的……
" O: F5 J$ R4 l# K* ~2 J* L猜想是不是FEM模型在局部还不够精确,应力梯度比较大(网格不够密,应力分布不平滑),造成两种应力计算方式有比较大的差异,因此造成疲劳结果差异较大。  j. A! R* n" H& V' C
7 A4 A7 h; R% e. o( r/ o

/ j6 H( i4 T0 R2 z4 |% ]1 x4. 加工工艺% K' F% J5 P, W/ n; C
+ n7 k+ u# q( l4 L: P$ X

& N5 K4 V1 q! x+ N& w$ Z" \4 q- y1 _/ s, j: U
7 f4 Q$ R+ L$ _; @: w) N# ^7 x& d8 K
表面处理、热处理工艺对疲劳寿命的影响很大,软件中选择不同的工艺计算出的结果差异明显,那么这里的各选项应该只是在某个地方考虑了一个系数吧?哪里可以查到这个系数(模型、经验公式、理论)呢?
) T! }" j# d8 m  B1 p8 m/ t
$ u, Y8 I- W; |9 N  ]( n* C! E
4 K( s  [9 G% c. c, L! @: {=========================================
! \  f7 M  }5 Yps:给最近发的几个相关帖子做个链接
, N6 l3 m8 q( k3 r, U, ~% p) n/ Y5 H/ Chttp://forum.simwe.com/thread-875979-1-1.html: o' \4 h$ L* Z/ q
http://forum.simwe.com/thread-876186-1-1.html

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×

评分

1

查看全部评分

 楼主| 发表于 2009-4-17 21:20:19 | 显示全部楼层 来自 重庆沙坪坝区
Simdroid开发平台
本帖最后由 apprent 于 2009-4-18 17:51 编辑
9 K* x' k/ K* D8 Z8 k1 y; c
" X% Y$ D2 X) j8 A1 @) M占一楼以备添加。# r1 H  r3 t' q0 P
~~~~~~~~~~~~~~~
/ U5 N4 r) ^8 E  J- m
6 M3 F  n1 m- u$ w感觉理论基础还是差得多,而软件又都过于黑箱了……& C' _* }4 r8 T& e0 q
估计如果自己用数学软件编一套计算程序后才会非常明白。, B0 [" s! b7 A% l6 ?% l' v, I
7 v# J$ p/ ^' t5 F
~~~~~~~~~~~~~~~1 v3 P4 r8 \3 z! _8 z
' m7 h. n: l- d$ P5 j) u1 u
我像是来这里做笔记的,呵呵- U/ U( @: a+ x. n) @
希望还是能对软件新手们有点帮助。请众位高人不要鄙弃。
回复 不支持

使用道具 举报

发表于 2009-4-18 10:00:08 | 显示全部楼层 来自 山西太原
和楼主同感?有的时候知其然不知其所以然。
回复 不支持

使用道具 举报

 楼主| 发表于 2009-4-18 17:07:27 | 显示全部楼层 来自 重庆沙坪坝区
本帖最后由 apprent 于 2009-4-18 17:55 编辑 / Y$ j4 t- w; Q: @4 _

/ h; k  k- S' a; x+ E6 r) l. p& t/ k9 d对于——$ c2 ^2 q- f6 ~& t5 ^6 U. S

' i4 `! A: |7 H6 u3. 取节点或单元的应力
$ M& i9 @; D, t; P# z% g8 E
- H* P+ ?. E: W6 t1 T2 `4 G9 i
3 t1 N- m7 I7 Z/ T, [( R" c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! Y! A4 R6 i- N! z" s来自《MSC.Fatigue User’s Guide 》的解释:3 Z& @& P. Q) ~1 ~$ \

" t8 m; n! a& |( V  [7 w4 fResults Location:
; G6 k; t* Q! l# f8 vThis parameter tells MSC.Fatigue whether to expect Nodal stress/strain results or Elemental centroid stress/strain results.
/ J: r2 k) q4 F; v! t4 Y) I% n- ~This dictates whether the user is setting up a global multi-node or global multi-element fatigue analysis.
. m+ r4 i5 U1 ^1 S$ uSubsequent parameters, results file types, and results displays are dependent on whether nodal or elemental data is being considered.
- o0 x& e$ S, |: dIf nodal data is being considered, the resulting fatigue lives are reported at the nodes. & `0 v5 R: y8 w9 O4 I% ?% N' ]" _
Conversely, if elemental data is being considered, the fatigue lives are reported at the element centroids.
# U# G6 O0 ~5 X: |/ {- H- nFatigue cracks invariably occur at free surfaces, and hence when a crack initiation method is used, node points results are usually required. ( t2 h' l4 H9 ^( Y+ ~3 C0 j
The exception is when a shell model is used, element centroid results may be extrapolated to the top or bottom surface.
+ a, G9 ?2 ~/ l9 C! }) R% `This is useful when there may be some doubt as to the accuracy of the node point results due to extrapolation and/or nodal averaging practices. 6 \+ H' W4 T. @7 @
The spot weld analyzer uses forces and moments from both nodes and elements.
( D3 d& o- w& c. a. E. J" D& lThe SEAM-weld analyzer takes stresses from both the top and bottom surface and needs both nodes and elements for this.' ?- I' c* O3 K3 s7 O  d3 o

; V  f2 f7 a7 e/ N  k翻译:
# S, b- x0 U2 c  L8 v: B这个参数告诉MSC.Fatigue去提取节点的应力应变结果还是单元中心的应力应变结果。
! W" {# ^" `, Z6 Y' K) @这决定了用户将要进行全局多节点还是全局多单元的疲劳分析。# h+ y. f2 R8 ^5 z  Q
其后的设置参数——结果文件类型和结果显示,将依赖于上述选择——使用节点还是单元的数据。
# h- q7 H1 Y) i  ^5 f如果使用节点数据,则疲劳结果将被呈现在节点上4 v8 a5 _( K: U
相反地,如果使用单元数据,则疲劳结果将被呈现在单元中心上& M# L. v6 Q: ?: e
疲劳断裂总是发生在自由表面,因此当进行“裂纹扩展”分析时,通常要求使用节点的数据。
/ ]* `2 J5 B2 x2 A; n! I6 V5 E例外的情况是,当使用壳单元模型时,单元中心的结果可能要被外插到顶层和底层表面。
; P# Q2 H( ?8 `9 z此时这个选项便发挥作用了,因为外插和/或节点平均操作会给节点结果的精确性带来疑问。( S( u. l0 g; h
点焊分析要使用节点和单元两者的力合力矩数据。2 l+ I6 R2 b* u7 Y
缝焊分析要使用节点和单元两者的从顶面到底面的应力数据。7 v  ~( @7 p- ^

6 I6 D1 P# i" J1 o# e) q3 jNodal Averaging:
! F) O' T1 Y. K8 ^* zDepending on how the finite element results are defined, nodal averaging of the stresses or strains may take place. ) y. Y' n! O( n! W& j4 _$ h
If grid point stresses exist and are selected, no averaging will occur.
) B4 g  `  a+ JHowever if the stresses or strains selected for the fatigue analysis are elemental based, such as results at integration points or elemental nodal values such that each element has a different value at the shared grid points, then nodal averaging will occur.
) s( j, d" j% |+ YThis averaging is done on a global basis such that every contributing element surrounding a particular node will be used in the averaging. : v; D9 F3 O" W. n- Q& G& {  |
The exception to this is if Group is selected in which case only elements in the Current Group will be used in the averaging. * c$ y! m# O3 z, w4 ]4 u# f
For the SEAM-weld analyzer only the current group can be used and no choice will be given.+ U0 V: B9 ]* P
7 }# Q1 n- J, h* H0 L0 o" X
翻译:
- T' D5 R  C+ S( a* H+ o根据有限元结果的定义方式(使用节点还是单元数据),应力或应变的节点平均操作可能被采用。: E. F0 A2 N% {. f2 c4 q
如果节点的应力存在并且被选择使用,那么节点平均就不会发生(不需要了)。3 e7 {, o; |1 `& j
然而如果基于单元的(比如积分点上值的或各单元在公共节点上有不同值的节点值)应力或应变被选中来进行疲劳分析,那么就需要进行节点平均操作。' k* D' l5 ]$ R# k4 X: Q3 C% o
节点平均操作是全局性的,也就是说,所有环绕着某一节点并对其有贡献的单元(在该节点的值)都将参与平均操作。
/ _1 H# m% a& a5 H1 O% _如果Group被选中,那么只有在当前组中的单元参与平均操作。: V: i& X. E" z. _9 f
对于缝焊分析,只有当前组选项,因此这个选项变为不可改动。
7 ~4 U# Y- H& d8 M: V
# d; R3 C8 S% D3 o0 Z: Z2 k
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! K7 S7 C* ?4 ^3 o/ m" U2 h# P8 ^" P5 F! @- G
我的模型中使用的单元的结果,节点平均选择了Group,因为我给模型用了“表面覆壳”,表面加了一层0.0001mm厚的壳单元,设为一个单独的组,疲劳计算时只针对这一个组的壳单元。& `* {. ?) C* J& [4 o8 N
因为:0 h% G! _* F2 i+ m/ r! d; p3 h7 r
1. 疲劳破坏/最低寿命一般都是发生在结构表面;
# I0 z/ r3 Z6 E8 F5 U2. 使用壳单元提取出的应力都是结构表层的,而用实体单元的话,表面节点的应力是从表层下的积分点外插上来的,相对而言,使用壳单元提取的应力更精确;
* U# m8 E& g+ k  {9 P( |" z& J3. 显然,只计算壳单元一个组,能减少计算量;8 j% Z- M. G. w* v, y: Y' S, i
4. MSC.Fatigue手册中还给出一个理由,由于理论知识的不足,我不是非常明白:The shell element stresses are normally recovered in element coordinate systems which have the advantage of being in the plane of the surface. This makes multiaxial assessment easier.(壳单元的应力坐标系有一个优点:在结构表面的平面上。这个特点使多轴评估更容易。)
0 n! W# G4 o! G: Z2 O-----------------------------------------3 y, J! h+ D, h1 G. I' G% l
) Z: |* H% ^" ~% y5 A+ r. L$ ]% ]
我想知道:6 r5 i* ]. W: Z- U+ e4 w
上面用颜色标出来的 global multi-node fatigue analysis or global multi-element fatigue analysis 是什么意思?% n" E9 C/ I" Z5 Q* I
, g7 R& p4 q6 w- U& b7 ?
resulting fatigue lives are reported at the nodes fatigue lives are reported at the element centroids 又有什么区别?(这个知道了,见下图——)
& ~9 ]" c6 e9 @0 p9 ]5 v3 s6 L) u. i4 ~1 y  u* a% J  k& m
刚刚做了两个例子试了试,发现基于我的FEM应力结果,选单元和选节点,疲劳计算结果差别不大(如上图),云图分布趋势一致,而且用单元计算的,云图更加圆滑。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
回复 不支持

使用道具 举报

发表于 2009-4-28 15:20:38 | 显示全部楼层 来自 北京丰台
学习过程中,也有同样的疑问,尤其是问题1和2.
" T: C% f5 S. G" P$ t向楼主学习。
回复 不支持

使用道具 举报

发表于 2009-5-5 19:14:27 | 显示全部楼层 来自 辽宁大连
1. 平均应力修正(Mean Stress Correction)2 G+ I8 Z0 C+ [4 M
那么,当我在MSC.fatigue中加载的载荷谱(力/力矩)的平均值不为零的时候,只要在这里选一个应力修正模型就可以了吗?
$ v: m1 m; Y: ?6 f0 k是否MSC.fatigue在计算时会先把这个载荷谱进行雨流计数,得到其Range、Mean和Cycle,然后套用所选的应力修正模型?
/ s9 A9 v0 N3 T. T2 p6 H* H3 C' c# r. e
时间关系,先只讨论这个问题。
2 M  i( H! b7 y& s2 w    我一般只适用疲劳软件里的雨流处理,然后自己分析,因此对一直做到底的软件内容不太了解。但我是这么认为的:; W% A6 B$ |; N1 Z9 f# N! J
    如果选择应力修正模型,就必须给出平均应力,不给出平均应力的话,还是不进行修正。所以,既然软件已经雨流出了应力范围、和平均应力的结果,只要用户选择了修正模型,软件就应该自动把雨流出来的平均应力带入到模型里面去的,而不用用户再行专门带入。
回复 不支持

使用道具 举报

发表于 2009-5-5 19:16:57 | 显示全部楼层 来自 辽宁大连
第二个问题将在“关于GL Wind中对S-N曲线的规定”的帖子讨论,不过今晚没时间了
回复 不支持

使用道具 举报

发表于 2009-5-6 15:42:22 | 显示全部楼层 来自 广东中山
楼上的是做风电的吗,对风机结构疲劳计算应该很牛吧
回复 不支持

使用道具 举报

发表于 2009-5-9 13:57:44 | 显示全部楼层 来自 湖南长沙
学习了,谢谢讨论!
回复 不支持

使用道具 举报

发表于 2009-7-16 19:03:22 | 显示全部楼层 来自 重庆沙坪坝区
学习了,谢谢先!
回复 不支持

使用道具 举报

发表于 2009-8-28 16:54:27 | 显示全部楼层 来自 湖南长沙
学习软件中,谢谢讨论!
回复 不支持

使用道具 举报

发表于 2011-1-3 21:12:36 | 显示全部楼层 来自 辽宁大连
一般情况下,是取节点力的
回复 不支持

使用道具 举报

发表于 2011-2-15 09:58:40 | 显示全部楼层 来自 辽宁沈阳
正准备做疲劳相关的分析,谢谢楼主的分享,关于上面谈到的有一点疑问想请教一下:6 O$ E4 I" X4 f3 d$ Z
“我的模型中使用的单元的结果,节点平均选择了Group,因为我给模型用了“表面覆壳”,表面加了一层0.0001mm厚的壳单元,设为一个单独的组,疲劳计算时只针对这一个组的壳单元。% ”5 u  b1 j# B; N) ]
  这个怎么实现的?
回复 不支持

使用道具 举报

发表于 2011-2-22 21:28:41 | 显示全部楼层 来自 陕西咸阳
做风电疲劳的人不少。
回复 不支持

使用道具 举报

发表于 2011-4-17 15:01:37 | 显示全部楼层 来自 河北秦皇岛
好帖子,先mark了再仔细研究
回复 不支持

使用道具 举报

发表于 2011-4-18 18:55:37 | 显示全部楼层 来自 陕西西安
要好好研究一下
回复 不支持

使用道具 举报

发表于 2011-4-18 18:59:46 | 显示全部楼层 来自 陕西西安
要好好研究一下
回复 不支持

使用道具 举报

发表于 2011-4-27 11:36:31 | 显示全部楼层 来自 河北秦皇岛
标记一下,慢慢学习
回复 不支持

使用道具 举报

发表于 2011-4-29 15:09:15 | 显示全部楼层 来自 台湾
学习了,谢谢!!
回复 不支持

使用道具 举报

发表于 2011-5-1 07:27:11 | 显示全部楼层 来自 江苏泰州
谢谢啊。这么好的资料。。
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-28 11:49 , Processed in 0.074782 second(s), 17 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表