CosmosWorks算法验证文件 9 (Buckling Analysis)

tigerdak

.
& y5 s8 F2 M% _5 O; s提示：屈曲分析（特征值法）。
6 d6 N5 @& N5 W
Title      Buckling of a Bar with Hinged Ends (Line Elements)
: z$ a" ]# ?3 Q8 d4 F
& P7 i  A* f. p! D" b# p9 sOverview

Reference:	S. Timoshenko, Strength of Material, Part II, Elementary Theory and Problems, 3rd Edition, D. Van Nostrand Co., Inc., New York, NY, 1956, pg. 148, article 29.
Analysis Type(s):	Buckling Analysis + k" g( w4 I3 O, D+ J3 dStatic   c. ?5 }8 {: d' S8 D% c, U9 f

# H& [( i- v4 \3 X# wTest Case

( s6 g6 S2 V* M+ v/ `2 a, [4 {Determine the critical buckling load of an axially loaded long slender bar of length L with hinged ends. The bar has a cross-sectional height h, and area A.
/ s0 n" }) P+ w$ Z9 t# q
+ o  s6 S$ f9 n8 R; C/ \Figure 127.1  Buckling Bar Problem Sketch
2 W. W. D7 ]( b5 m" y+ l: l
0 a8 Q( ~- L; u. R5 L! p$ b. X4 j

Material Properties E = 30E6 psi

Geometric Properties l = 200 in A = 0.25 in2 h = 0.5 in

Loading F = 1 lb

! ]' ~5 c; c/ i, VAnalysis Assumptions and Modeling NotesOnly the upper half of the bar is modeled because of symmetry. The boundary conditions become free-fixed for the half symmetry model. A total of 10 master degrees of freedom in the X-direction are selected to characterize the buckling mode. The moment of inertia of the bar is calculated as I = Ah2/12 = 0.0052083 in4 .

Results Comparison

	Target	ANSYS	Ratio
Fcr, lb	38.553	38.553 [1]	1.000

* L$ E( w: o4 o1 l8 D2 V) Q

• Fcr = Load Factor (1st mode).

  f) j/ n+ ~) m4 r2 l' i
[ 本帖最后由 tigerdak 于 2007-11-8 18:44 编辑 ]

tigerdak

屈曲分析方法分两种：特征值法和非线性分析法。特征值法偏向于纯理论计算，不考虑零件的缺陷和环境影响，所以非线性分析方法更准确。以下方法在杆的中间点加上微小力来模拟零件缺陷或环境影响。一般会比理论计算值低。
+ F$ Y; P. ?3 o( K' P
7 u/ Z& O5 r! y) Q1. 在杆的中间点加上1e-9 lb 初始扰动的非线性屈曲分析
! N) E9 z/ c, V2 ~" r% j8 k  ]


5 R1 {9 K) f  _+ J# {
, k0 \8 e" ^1 _$ s2. 在杆的中间点加上1e-3 lb 初始扰动的非线性屈曲分析


$ \, L8 n$ e; j3 ~% T6 h& v
[ 本帖最后由 tigerdak 于 2007-11-8 12:44 编辑 ]

tigerdak

对比CosmosWorks的计算结果

1. 特征值方法，即扭曲研究

Load＝38.553 lb
' q9 v  e: a- e2 f
1 U9 G5 Q& U' s5 K6 l# p6 d: B
- L; B3 u6 n1 l6 M

2. 加 FX＝0.001 lb 初始扰动的非线性分析，从图中可看到在0.732996秒发生扭曲，此时Load＝36.65 lb
. D% {8 l1 a. {& L

, f0 E6 ?% a' o8 {7 P1 b- l  v