关于 combin14 的使用问题
关于combin14 的使用问题,为此做了一些探讨,现将探讨结果于大家分享,希望能有所启
发。
如何正确使用 Combin14 单元的一些探讨
由于不少网友的帖子提到如何正确使用 Combin14 单元的问题,下面通过一些算例的比较,对该问题做了一些探讨。希望起到抛砖引玉的作用。
例 1
45º 直线,长度 100,划分两个 beam3 单元,
两端 (下端为节点 3,上端为节点 4) 用 combin14 连接到重合节点 (下端为节点 1,上端为节点 2),
. ]7 Q f8 k, X- y' K9 F( @! L, B
节点 1 – 3: x,y 方向两个 Combin14 单元;
3 ` Y2 I( m2 b# m
节点 2 – 4:仅 x 方向的 Combin14 单元
对非梁节点 (节点 1,2) 约束三个位移自由度以消除刚体运动。
梁截面 A = 1, I = 1, H = 1
Combin14单元 K = 1
梁中间节点 (节点 6) 受垂直于直线的横向力 1.0×sqrt(2)' @. h6 v; X8 F( t7 G1 B (x 和 y两个方向各 1.0)
各节点坐标如下:
NODE9 r0 l6 j) A$ C2 f* U# s
X
Y$ h7 P* {8 E2 F
4 S" _) O
5 t; M. |5 A
Z
1, m* N( B1 s7 R g+ Y& T
0.00009 {& o+ ~7 ~" x0 m0 q
0.0000
0.0000
1 b7 n. q) `3 Y) \
2+ F& a$ X/ c; V' [6 Q( ?6 D
100.00
8 y/ |1 _ d; {+ e' G5 c# I0 r3 e
i# U( _! q2 z) ]9 @3 _; y, ~
100.00
0.0000
3. s$ L& E9 R f' N/ k
0.0000
0.0000
0.0000
4
100.00
100.00# x: |6 R" H7 [# B% J( |5 ^ 0.0000
5. x) L$ Q7 a( s6 k J
25.000" J* T5 `$ c4 V1 O
25.000
0.0000
6
50.0004 g Q# ~) v- D
50.000' I' E. ]+ k o3 m
0.0000
6 n3 [ Z. L: R6 ~; l6 H+ x
7
75.000
75.000. n, b t* c% [, |" a4 c/ K 0.0000
计算结果:
位移结果:
NODE: W/ }' G9 F" i6 ?7 ]$ R4 A4 U UX
UY
UZ
USUM
1 l) n" O5 f+ X: [/ [$ I1 m# U; L 0.0000
0.0000' u3 w. T) Y4 `& H
0.00003 n3 N: }6 t3 ?# f* q, |- C 0.0000- c0 l% o( x% y# V
2
0.0000
0.0000" T+ R$ D8 K" M/ s+ h8 B& V
0.0000
0.0000
39 |6 W7 O, @2 l; L; m
0.51548E-12 -1.0000
0.0000
1.0000
9 q* Z3 P# o7 |7 Q7 q9 t, Y) L& ]
4- g6 Q x- B( |% @4 E9 ?: [- \
1.00005 R) d' F! P; i1 u. h
-2.0007! H9 n( t. i( f( R' ~8 \. Z
0.0000
2.23674 ~7 b; c' k# y1 L: s
: ^* R' S8 v% d; _% g
5
0.44291
-1.4431. i3 T4 w: L+ ^* N- B8 [/ n! I" @# z 0.00007 J+ q$ E( R6 S9 b% P% ?
1.5095& ]7 u% [ O# c- J- j7 d
, \/ R9 L/ E! M+ [/ q
6' q9 `6 d3 `# m: |
0.78060
-1.78097 E. l! i3 V( t; S
0.0000
1.9445
7
0.94291
-1.94346 E8 L! v, W- e% w
0.0000
2.1601
1 @" m* F* y# c% Q: l/ a6 \
MAXIMUM ABSOLUTE VALUES
9 {6 ?1 Y) D4 A' `
NODE
46 z0 g. Q4 L* l, z( F& N
4
02 i. d$ Y5 o( r
4
; d/ O0 }8 H8 U
VALUE
1.0000
-2.0007/ \" s0 k% [( s1 t$ j/ R) C7 r) n 0.00002 K; N$ x9 C6 W: s
2.2367
反力结果:
NODE( B" |8 C8 m, }" X% g9 }# x
FX' E* R5 t9 E8 R4 r* h1 Z/ d
FY- Y* b6 \9 ]1 _; k2 l! z
FZ
1 F5 q" l' X: V1 c9 ]( D
1 -0.51548E-12
1.0000* w" o0 J' P: M' X8 \' B
0.0000
2% Z; [# _+ E6 g" b( j
-1.0000
0.0000
0.0000
! l# i. }! G% S0 S/ z( V- i# l9 r( E TOTAL VALUES
6 i1 r
7 C! y9 y. L5 {
VALUE2 u# j& F. d& F/ ?- o/ c; _, a
-1.0000
1.00002 x$ P/ `7 q; d6 z+ s+ o
0.0000
弹簧单元内力:
ELEM% T4 W: H% U5 `6 X+ ^* K* i5 b SMIS1' b$ y& \- ?) {" e# c* {
5
0.51548E-12+ R$ t! O: A2 I
= (UX1-UX3)*1≈ 0
% C( W- ~& E x/ l. W: t
4 G1 s8 Y. {( v+ f6 a* ~
67 e0 t+ H/ w7 H! B
1.0000
/ g' T4 Q3 u6 J
. c0 V; ]9 w) g4 k* F' Q$ |
= (UY1-UY3)*1
7! k# V% L% C: K$ c$ F
-1.0000
= (UX2-UX4)*1
可以看到,合反力与所加的载荷相同。单元 5 (连接节点 1、3 的 x 方向弹簧) 的内力为零,即该弹簧单元基本不起作用。
例 2
结构同上,但改用两个 combin14 单元,即三个 combin14 单元中任意去掉一个,结果如下:
(1)删除单元 5 :连接节点 1、3 的 x 方向弹簧
位移结果
NODE
UX& V3 l* S9 T) [7 C
UY * k6 }# E1 Q9 d9 A& f3 Z
UZ
USUM
4 d2 U0 y- a% h, @! \
5 A. q
1
0.0000
0.0000( V) Z$ i0 {( s4 R" {7 O
0.0000
0.0000' n2 s1 [# ~6 r- F
2 V& t+ p0 l# P4 h! C/ X+ c8 R7 J
0.0000
0.0000
0.0000
0.0000
6 l4 u' @' H/ W8 D
3$ S5 A1 w m% S" p
0.88890
-1.00002 [9 e) R3 Z& l+ L4 _2 I5 W
0.0000
1.3380% }6 `/ ]9 O+ H R% G: Y
4
1.0000
-1.1118
0.0000* _+ ~9 k0 n& Z! c6 t
1.4953 7 u% Y. n5 s6 u- u y2 [
5) U1 Y; Y, d+ ?% p& K
1.1096; j( `( P' x$ q8 _ a6 |( i -1.2209# i* C) ^& q; n; N2 _
0.0000
1.6497
( m$ V4 L7 Y; b
6' x* F# E" C6 ]& i
1.2250
-1.3365
0.0000
1.8130
7: n; G5 k% h) S" x+ Q P0 k0 M
1.1651
-1.2767
0.0000
1.7285
MAXIMUM ABSOLUTE VALUES 6 {$ P* \; R- j0 E }% D
NODE7 O! G8 f8 {5 b# w- J, v/ L! J+ K 6$ k9 g9 L/ k3 v4 j+ h: \, J {
6
6
/ v+ I" L3 g+ k$ V0 K& n
1 k, A8 _ v& ^5 l" l G
VALUE
1.22500 z$ g& i( @$ J& A$ d1 E
-1.3365) }* Q5 j% ~. P3 w+ n2 v, L
0.0000
1.8130
反力结果
NODE5 m6 z2 p |0 [- \
FX
FY
FZ7 ]& Z2 D0 G! ~" g
1) o( y: x% k: v8 e& i$ x
0.0000
1.0000
0.0000' w- @" i( h+ B0 z
2
-1.0000/ C+ L7 p& O$ v
0.0000$ J2 S; l+ T8 i& F
0.00001 v: t1 a' N$ f, e0 s- X: P
7 \3 h* z6 ?7 p, k$ x' J* h
TOTAL VALUES
2 A
3 ^8 [3 W: L3 |$ a: Z: d
4 H/ V
VALUE
-1.0000" H, |& X. o# L2 {* f
1.0000# v# L; h# {& v3 C2 b' V* Z, }1 n- X 0.0000
弹簧单元内力:
ELEM$ U4 U! k4 M# ]/ |6 v6 A! h g( o SMIS10 K6 C% u: t# A& a X7 M
6, d/ F' _( H6 |. ~3 J, V
1.0000% ~* \# X6 f( P
= (UY1-UY3)*1
7
-1.0000
= (UX2-UX4)*1
(2)删除单元 6 :连接节点 2、4 的 x 方向弹簧
(3)删除单元 7 :连接节点 1、3 的 y 方向弹簧
都不能求解。原因是约束不足。
例 3
结构同上,但改用两个 combin14 单元,单元方向的定义改用 KEYOPT(2) = 0,然后用 KEYOPT(3) 定义为 2D 单元。结果也不能求解。原因是单元长度为零报错。可见,对于 KEYOPT(2) = 0 的情况,combin14 单元的节点不能重合。
例 4
同例 3,但将节点 1 和 2 从节点 3 和 4 重合的位置沿 x,y 各移开 1。
不再出现单元长度为零的错误,但不能求解,原因是约束不足。
例 5
在例 4 的基础上增加约束节点 6
的 ROTz,仍然约束不足。
原因是只有沿同一直线方向的两个弹簧,不能约束转动和沿该直线横向的位移。
例 6
在例 4 的基础上增加两个约束:节点 3 的 Uy 和 ROTz,可以求解。
位移结果:
NODE
UX$ {8 Z; J" U) g2 w, B
UY7 s) V; u% `4 W6 Z. Z
UZ
USUM4 D W: n0 z; a5 F7 N
$ ~ Y: p5 k* ~; b) H6 Z
1/ B0 G [* _; C1 X4 C- ^( |$ g" U
0.0000
0.00007 D9 O, G& k5 J3 o5 [
0.0000
0.00008 t& H/ X- k q2 y1 _
5 V, J, d5 F( e2 ^
22 Q* K, J J# m# h
0.0000# o3 v5 ?$ M% L) E( E
0.0000' q. X) c; `4 K; H, }, ?5 n. x 0.0000
0.0000
3# A; J5 o) m5 s' f/ Z
1.00038 B! `& t [0 z$ m
0.0000
0.00006 C4 Y3 N- ^2 E, Z
1.0003, e- e4 `: i9 h) T. q" ]8 s; Z: b
4
2.4030
-1.4033. e j! Q$ B& t4 z' o
0.0000- b) Z2 U- T* }, x% y. C: u! c 2.7827; a7 [& v( T C3 M! U; D
: W. X9 v0 z. z3 j
5
1.1756' V1 p+ _; X, F
-0.17546+ m- J$ i" U/ B0 Y* _
0.0000# e: b3 ]. X0 y1 k
1.1886* B6 E& T5 `! a8 D l
0 h0 x- t; R& E9 K" Y7 v
6
1.56146 q* n ?% D8 |
-0.56136
0.0000
1.6592
: P5 N3 b& s8 J
7
1.9822/ Y2 Y, n: n2 v6 h- {( M' F
-0.98235+ M6 W5 P( A: N2 G
0.0000
2.2122
; H% t3 {) I; Y7 x8 o$ y: ^ MAXIMUM ABSOLUTE VALUES
NODE, \1 X" A. l* P
4
4
09 G8 }3 r$ ^: t! b& g
4
' x# L1 p, c$ Y
VALUE& Q+ J1 s H1 v5 V% X
2.4030. ?+ L9 b* a5 L, X2 J$ Z
-1.4033
0.00008 Z3 `2 C( \# i" C6 Z8 x* P' M
2.78272 H6 _) ~8 ]/ c
反力结果如下:
: H; k% ]6 s2 v. `& W+ ^" v
NODE
FX6 g( W# w" k# U% B! T6 p
FY
FZ
MZ: r& `1 c! W: c0 z, Q1 F
: Z- L+ y8 T8 |# A4 C# a
1 -0.50017. {& M5 M
2 ]* s) E6 A
-0.500175 q* K& L( ]$ U7 T, E
0.0000- G" X+ y) ~ }% J5 ]5 K, V
2 -0.49983: ^- i+ Y* E5 ~6 k& |2 |. r8 N -0.49983
0.0000" Q4 m* ~0 H' n
* B5 r; I n, {0 P+ P, v
3
8 y9 o# ~! B0 \* v: ]& ^
, {* N, j7 w$ J e; `! T, m" e
2.00004 p1 ]9 p2 K$ E. h2 n
100.00' h/ K$ G+ z: X- J: r5 Q
TOTAL VALUES
VALUE + T: [% @3 o8 b
-1.0000# B" Y9 G, o) g( ~4 X) H) Q6 Q: V
1.0000
/ q0 t5 H5 k: N, W6 ~ t/ ]! {; p7 d
0.00003 a! w- a$ T- ~" T5 T
3 t1 Q
4 A, `7 ~( w8 T
100.00
弹簧单元内力:
ELEM' I7 q1 E* K: l
SMIS14 e+ D: u2 H! A4 ?) I4 t/ v. f
% P6 R2 h7 b/ k. G; D; V: C! J( c
5. z" t6 o4 P- ~. T- x
0 V) Z* L5 P. j( c( v1 S
0.70711- k* S) A& }- h# i8 \- J2 x
= ((UX3-UX1)*cos(45) + (UY3-UY1)*sin(45))*1
A- @9 V/ d' W5 |$ D
6- S4 v- `( g: j4 \
-0.70711# D3 `) J- Y7 t# S! v
= ((UX2-UX4)*cos(45) + (UY2-UY4)*sin(45))*1. ]4 _$ I& K8 n0 g- |' I3 ~
注意:由于约束节点 3,所以反力中出现了外力相对节点的弯矩 100.0。这说明,combin14 单元只具有沿其轴向 (节点 I 到 J 的方向) 的刚度
例 7
在例 4 的基础上增加两个约束:节点 6 的 Uy 和 ROTz,可以求解。
位移结果:
NODE. T; H) w1 U" {/ @2 M$ T
UX
UY; ?7 S7 F" ^2 T9 g' x
UZ/ l: K# V5 t# F7 W/ d
USUM3 f! R5 B' o3 W3 ?1 o
1
0.0000. z- e* X1 J7 ?# G6 u
0.0000% B8 v. M, a; n1 {9 I+ C6 d
0.0000% C+ m% R5 _" }' w
0.00004 A8 }9 _/ S2 [% Q
2
0.0000
0.0000! H5 ^' ^0 M, t6 Q& |1 j
0.0000' s. G M3 ]# n5 y: w6 ~3 s+ b 0.0000
+ |$ p+ f! e6 G6 l9 p. q
3
1.0002
-0.16836E-03
0.0000* r. y9 a7 Q' a; S1 m
1.0002# v8 v5 j' p% i5 _; c
( F5 j: v; Q& H2 A0 h7 s1 ?
4
1.0002' s4 x7 A( Z) M M% s
-0.16836E-03
0.00001 G$ [+ m; X/ X' o7 r Q! g! p
1.0002- R! A$ x: p# {% \- d0 t& r
% S- l1 R' Z, w2 [; |
51 f* L$ O$ y6 [1 S( f
1.0003
-0.84179E-045 C$ t4 w; E5 G- q9 C
0.0000* a+ G; |9 C1 j% O, B7 L% ]9 z) w
1.0003
! i6 _5 p% U- l: g. [- y
6
1.0003+ O/ ]8 B0 @( k$ S$ w6 N! N/ Z 0.0000
0.0000" g% o$ o4 g( R: r7 R
1.00030 `% B4 L. {3 M% K' |2 @% X
$ r U. ~" H V- V' o
7; B: R2 X' w* `9 `$ {) ~/ d4 A8 ]
1.0003
-0.84179E-04$ d2 ]4 F% R) A2 C
0.0000
1.0003
( U4 W" j' W- M6 w( p- E
MAXIMUM ABSOLUTE VALUES
/ @1 n1 v7 J/ [; U! L( Y5 [
NODE0 X% g4 ~& [ N0 Z) L
69 d9 w1 Z% Y4 b9 o( l
32 H+ r; r/ K4 g$ ]2 _& V
0# `/ z" ?- `# [/ H% h) }
6
' `5 |. L* e6 P! ]$ r8 \* i
VALUE; P# N8 F4 \( j6 {# f# D+ d) V% Q 1.0003
-0.16836E-034 v/ R6 u- B9 r' b- `* r
0.0000
1.0003
反力结果:
NODE
FX+ G' b! n( A% t7 t4 Z0 O
FY( M% @) W2 R! U2 l
FZ4 ]/ w& H$ U3 }) ~4 h1 u
MZ6 l( I8 n9 I) A& P
. L7 G/ p0 o1 M, q! C
1 -0.50000
-0.50000
0.0000
2 -0.50000
-0.500005 S- _6 i/ a4 b1 D7 ~: K9 v 0.0000
% L s; b( P E
6" z) ]+ e( W" Z" X/ L. F
2.00001 W. d- K: Y! f
* ?+ S9 E& `( ^: |
$ h% ] j( y# T. t: r$ x( q2 v
0.15827E-10
TOTAL VALUES
VALUE" z! R. Z* i2 Y. H' j3 I! u }: x; l
-1.0000
1.0000$ |& i" ~& y5 F# F0 U
0.0000
' b2 g' v, F0 R3 s
! f9 x$ W% |; k
-0.15827E-10
弹簧单元内力:
ELEM- G! ]3 y# q: O: M% i5 ^& l
SMIS1
5
0.70734
= ((UX3-UX1)*cos(45) + (UY3-UY1)*sin(45))*1
: I$ W" F$ Q- W9 b% o- h
6( _1 F/ D1 I% q% @3 k
-0.70687
= ((UX2-UX4)*cos(45) + (UY2-UY4)*sin(45))*14 n8 |8 p8 o- j/ c" l
注意:由于改为约束节点 6,所以反力中外力相对节点的合弯矩基本是零。
同样,combin14 单元只具有沿其轴向 (节点 I 到 J 的方向) 的刚度
例 8
同例 3,但将节点 1 和 2 从节点 3 和 4 重合的位置沿 x,y 分别移开 1 (节点 1 沿 x 方向,节点 2 沿 y 方向),可以求解。
位移结果:
. B! u$ w0 Q# @0 [' H
NODE8 Z8 x( @4 ^2 b! u j
UX$ z( h7 P2 w7 i7 y
UY4 l K% _4 f( [3 w* G$ t' C0 u; T
UZ- F$ `9 D' C3 C Z
USUM
) m. l9 {! u1 @6 r$ n( a' @% m2 Q# @
1
0.0000/ V; R; r& C) b) q
0.0000 P7 V/ B1 U8 b9 L- G0 _8 E3 X 0.0000' Y0 F0 t' ^: R$ v1 T
0.0000
2
0.0000
0.0000
0.0000$ @3 S0 x$ U3 _
0.0000
6 T) g; y& s! N* z8 z" r0 i5 a
37 S/ ]5 ?, _1 i" Q5 V* D1 E
1.0000: V: l4 r" R* C' F0 j0 D
-0.532396 F2 i5 T$ O9 o
0.0000
1.1329
4
1.4683
-1.0000; j$ _$ x i4 b" Q3 [
0.0000
1.7765
5
1.3100
-0.84220/ F9 @2 p; w0 L- b
0.0000
1.5574
64 M# Z) X1 k; h( r. q1 j# ~$ \, _ 1.5147, U, N: ^$ M1 T5 N7 N. ~
-1.0468
0.0000
1.8413
3 P6 k- G( b7 V! G
4 R: b F
7
1.54414 C; m, F8 J1 ~$ ^+ Y: j4 I
-1.0760# z7 ]5 ?4 J( W% o: ^8 J
0.0000# n" N+ w0 E6 z
1.8820, d/ H E. ^: l g2 A4 e
MAXIMUM ABSOLUTE VALUES NODE
7& {9 a8 t5 i8 a; q
7
7
4 y8 i1 @9 y3 `) o/ S
4 i8 n* I* O0 [* t
VALUE
1.5441
-1.0760
0.0000
1.88206 h+ u9 @6 I/ |/ L
反力结果:
NODE
FX7 N: {. T$ a) k R$ I
FY# \# }% B+ a+ F# H t
FZ2 s# k9 Z; i# g4 E) Y
1' H5 [4 |# F7 B% R/ m
-1.0000* T% x- G7 o7 }& j
0.0000) k4 t& D- u" p; Z4 ?9 N0 i
0.0000
, O" R) _, A: i6 Q. ?( u! O
2% w% M1 d3 A2 `3 q, {# r% l
0.0000
1.00007 Y% D$ c6 W6 E7 a2 ?2 t! u% z! k 0.0000
TOTAL VALUES
9 N! ?1 ?4 z0 C2 V+ C
VALUE
-1.0000
1.0000
0.0000, u& j% ^+ A) ?, N/ M5 o
弹簧单元内力:
ELEM
SMIS11 u; s' [2 _( Z0 Q* s& B- b) A
5
* m" c5 g* V. R) w/ o- {' s
1.0000 ; A" E5 K: j% d% P" D4 _* E2 A$ _
= (UX3-UX1)*1
8 \/ g& [% H/ V6 c
6: D3 A, n W. w$ e( T% w
-1.0000
= (UY2-UY4)*1
例 9
同例 3,但将节点 1 和 2 从节点 3 和 4 重合的位置沿 y,x 分别移开 1 (节点 1 沿 y 方向,节点 2 沿 x 方向),可以求解。
位移结果:
5 ~, y% c3 j/ w5 o, X
NODE
UX
UY
UZ0 O5 ^' @3 v0 K. V+ o
USUM3 }' r! S& a% F1 \; t- Q( b
13 I) H& p/ f: \
0.0000
0.0000
0.0000* l! m V4 i) M l) m! b
0.0000( w1 V, [! {) i, f9 {" c
2
0.0000- c. Z, l' ?% V: r5 r/ f
0.0000
0.0000
0.0000
- J% M8 y7 C4 o/ N! g3 t' x8 M- }
3/ N+ X; }' N" e
0.88890
-1.0000
0.0000
1.33803 @2 ^' L& a& q9 C A E0 k6 I) Y
0 M6 \5 y) P% p# ?! n. O( [* N1 E
4- z; k# c! ^& w+ K1 D# n' A
1.0000; U! C) Y% T* C. O( |
-1.1118
0.0000# y9 `# N0 C9 I" c) A; U/ \4 J, S- ~
1.49537 Q; J$ k" e3 Z) T
' P& K1 k0 J; C4 |" m: ?" M
5+ C3 D9 z) k( }: _+ m0 g
1.10969 d& D: W9 o/ v5 l
-1.22099 n0 D R$ v$ r5 l8 w8 N
0.00007 v; S3 e: w {5 o3 c
1.6497
6& n: j' J4 O/ E% ^! ^- |% W) w& a! Z$ k2 c 1.22509 C/ \4 U. x, p! j
-1.33653 z4 X8 S& J3 r1 V7 @! ~
0.0000
1.81303 {5 B! G( P! F& ~
) [1 P- u* ?: f0 x
7! W1 n9 v4 d$ [) _
1.1651
-1.2767
0.0000
1.7285) N6 C9 h* V$ ^9 L8 \* D% c7 D MAXIMUM ABSOLUTE VALUES
NODE, V; P7 N5 {4 J& o
6
63 \) M# \' G6 X4 d( C
6
VALUE+ \+ [8 |$ N2 l- p2 B
1.22508 F. Y l2 j. l' m
-1.3365
0.0000 l# T5 v3 J+ A! I& p8 l# K
1.8130
反力结果:
& C2 L0 B4 D6 _* q
NODE
FX
FY& }- }* J* u7 `
FZ
' k* {( `, G; T8 o
( g+ `0 n) S- y
! W4 l# K$ h; ^$ I s y
1
/ o" O% @9 P, Z0 w
0.0000
1.00007 z8 @6 z0 o$ F9 G& n
0.0000
21 }3 \3 ~) H7 u# r1 u- W
* _1 ^0 H* p. k W2 z* b- B1 x o) R- e# n
-1.0000- p% E S8 O8 e# `2 O; w! {
0.0000
0.00004 n; I# q( o, L+ e( X
TOTAL VALUES
VALUE, q+ y5 `9 R% U: i( L, j' o0 G& n' f
-1.00000 S* b! ^# W2 y _8 L' M5 i; u" H
1.00006 E# j2 }) M/ e2 y3 ~' o" d, H: v; Q3 q
0.0000
弹簧单元内力:
ELEM
SMIS1, Z$ m5 E% V$ C Y
S w9 ]8 w9 t+ o: Z# w# H
5/ \% M; |+ M) O# v
1.0000$ {1 o% f$ w" r7 o+ e
) Q. Q! O- _+ b( H1 O
= (UY1-UY3)*1
7 l6 x$ j! p: m) {2 D* g1 x
6* E" j! f/ }0 E4 e* W/ ^) o) ?4 v
: J7 l3 t7 g4 |% e9 b
) J, E0 Q) w$ i5 z- }8 n! [# d4 {( B" Y
-1.0000* X% w5 `6 R H0 ~* s
1 r: e- s9 ^
2 m# T# ]
3 O
= (UX2-UX4)*1
经验算,上述各例的合反力与所施加的外载荷的合力是一致的。弹簧内力也与相关节点的相对位移与弹簧刚度之积一致。
下面再看看,如果 KEYOPT(2)>0,但节点不重合的情况。
例 10
基本同例 1,但将节点 1 和 2 沿 45º 方向移动,x、y 方向各移动 1 (节点 1 移动 -1,节点 2 移动 1)。
位移结果:
2 _! r' ~" C( r7 @6 U9 l; c
ANSYS软件介绍与实例讲解
一简述ANSYS软件的发展史。 1970年,Doctor John Swanson博士洞察到计算机模拟工程应该商品化,于是创立了ANSYS公司,总部位于美国宾夕法尼亚州的匹兹堡。30年来,ANSYS 公司致力于设计分析软件的开发,不断吸取新的计算方法和技术,领导着世界有限元技术的发展,并为全球工业广泛接受,其50000多用户遍及世界。 ANSYS软件的第一个版本仅提供了热分析及线性结构分析功能,像当时的大多数程序一样,它只是一个批处理程序,且只能在大型计算机上运行。 20世纪70年代初。ANSYS软件中融入了新的技术以及用户的要求,从而使程序发生了很大的变化,非线性、子结构以及更多的单元类型被加入到子程序。70年代末交互方式的加入是该软件最为显著的变化,它大大的简化了模型生成和结果评价。在进行分析之前,可用交互式图形来验证模型的几何形状、材料及边界条件;在分析完成以后,计算结果的图形显示,立即可用于分析检验。 今天软件的功能更加强大,使用更加便利。ANSYS提供的虚拟样机设计法,使用户减少了昂贵费时的物理样机,在一个连续的、相互协作的工程设计中,分析用于整个产品的开发过程。ANSYS分析模拟工具易于使用、支持多种工作平台、并在异种异构平台上数据百分百兼容、提供了多种耦合的分析功能。 ANSYS公司对软件的质量非常重视,新版的必须通过7000道标准考题。业界典范的质保体系,自动化规范化的质量测试使ANSYS公司于1995年5月在设计分析软件中第一个通过了ISO9001的质量体系认证。 ANSYS公司于1996年2月在北京开设了第一个驻华办事机构,短短几年的时间里发展到北京、上海、成都等多个办事处。ANSYS软件与中国压力容器标准化技术委员会合作,在1996年开发了符合中国JB4732-95国家标准的中国压力容器版。作为ANSYS集团用户的铁路机车车辆总公司,在其机车提速的研制中,ANSYS软件已经开始发挥作用。 二节点﹑单元﹑单元类型的基本概念。 节点:几何模型通过划分网格,转化为有限元模型,节点构成了网格的分布和形状,是构成有限元模型的基本元素。 单元:有限元模型的组成元素,主要有点、线、面、体。 单元类型:根据实体模型划分网格时所要确定的单元的形状,是单元属性的一部分,单元类型决定了单元的自由度,包括线单元(梁、杆、弹簧单元)、壳单元(用于薄板或曲面模型)、二维实体单元、三维实体单元、线性单元、二次单元和P–单元。 三用ANSYS软件进行分析的一般过程。 1建立有限元模型 (1)指定工作文件名和工作标题。 该项工作并不是必须要求做的,但是做对多个工程问题进行分析时推荐使用工作文件名和工作标题。
纷享销客用户使用手册 管理员第二版
纷享销客用户使用手册-管理员) 第二版 (
1 1
目录 创建企 业 (3) 邀请同 事 (4) 纷享管理后 台 (6) 部门与员 工 (7) 权限管 理 (11) 帐号安 全 (14) 应用管 理 (15) 企业设 置 (16)
外勤签到设置 (19) 考勤签到设置 (21) 报数系统后台管理 (23) 版本说明: 1 纷享销客用户使用手册版本:V4.7.0 日期:2015.11.19 版本更新: 文档版更新日更新概 新人V2.0 2015/11/19 王创建文档 学琳
纷享销客用户使用手册 创建企业企业第一个注册用户即为企业后台管理员,企业创建成功后,登录前往管理后台进行管理操作。 进入纷享销客官方网页,点击立即注册,进入注册流程,通过手机注册进入企业创建界面,填写企业名称,填写当前注册用户姓名和密码,点击创建。企业创建成功后,系统生成企业帐号方便登录。 纷享销客用户使用手册 邀请同事点击中部黄色指引系统指引会引导管理员前往邀请同事,管理员登录进入企业,条中的[邀请同事],弹出窗口提供给管理员两种邀请方式,一是通过手机号码或邮箱发送邀请,二是直接生成邀请链接复制黏贴给同事。.纷享销客用户使用手册
如果管理员没有通过引导方式邀请,也可在移动端设置-邀请同事,通过微信、短信、邮箱的方式调取本机的应用快速发出邀请。管理员可在此界面进行用户注册同事接受邀请并进行注册,提交信息给管理员,审核,审核通过的同事就正式加入当前企业,等待权限等其他设置。.纷享销客用户使用手册 纷享管理后台和应用,管理管理员登录进入企业后,网页导航栏包括功能模块:工作、CRM员还会拥有管理后台入口,点击[管理]进入管理后台,需要再次输入密码,确保后台管理的安全。企业信管理员可通过管理后台进行部门与员工维护、权限设置、账号安全管理、纷享提供用户分类管理员设定,公告管理、息维护、删除信息管理、日志统计等;不同权限的管理员分门别类的对各自模块继续维护和管理。.纷享销客用户使用手册
有限元ansys静力分析的一个小例子
有限元 学院:机电学院 专业: 姓名: 学号:
一、问题描述 如图所示的平面,板厚为0.01m,左端固定,右端作用50kg的均布载荷,对其进行静力分析。弹性模量为210GPa,泊松比为0.25. 二、分析步骤 1.启动ansys,进入ansys界面。 2.定义工作文件名 进入ANSYS/Multiphsics的的程序界面后,单击Utility Menu菜单下File中Change Jobname的按钮,会弹出Change Jobname对话框,输入gangban为工作文件名,点击ok。 3.定义分析标题 选择菜单File-Change Title在弹出的对话框中,输入Plane Model作为分析标题,单击ok。 4.重新显示 选择菜单Plot-Replot单击该按钮后,所命令的分析标题工作文件名出现在ANSYS 中。 5.选择分析类型 在弹出的对话框中,选择分析类型,由于此例属于结构分析,选择菜单Main Menu:Preferences,故选择Structural这一项,单击ok。 6.定义单元类型 选择菜单Main Menu-Preprocessor-Element Type-Add/Edit/Delete单击弹出对话框中的Add按钮,弹出单元库对话框,在材料的单元库中选Plane82单元。即在左侧的窗口中选取Solid单元,在右侧选择8节点的82单元。然后单击ok。 7.选择分析类型 定义完单元类型后,Element Type对话框中的Option按钮被激活,单击后弹出一个对话框,在Elenment behavior中选择Plane strs w/ thk,在Extra Element output 中,选择Nodal stress,单击close,关闭单元类型对话框。 8.定义实常数 选择菜单Main Menu-Preprocessor-Real Constants Add/Edit/Delete执行该命令后,在弹出Real Constants对话框中单击Add按钮,确认单元无误后,单击ok,弹出Real Constants Set Number 1,for Plane 82对话框,在thickness后面输入板的厚度0.01单击ok,单击close。 9.定义力学参数 选择菜单Main Menu-Preprocessor-Material Props-Material Model 在弹出的对
ANSYS中简支梁的模拟计算
1 E c ; / E c lE s _2卜+僅 12 (5-30) 通过大型有限元软件ANSYS 对简支梁进行模拟计算 下面以钢筋混凝土简支梁的 ANSYS ①程序数值模拟的应用实例,对ANSYS ⑧程序的应用方法及 模拟效果进行验证,梁的尺寸、配筋及荷载如图5-9所示。钢筋采用H 级钢,混凝土强度等级为 C30。 2.1单元类型 i )混凝土单元:采用 ANSYS ①程序单元库中 SOLID65单元。 (ii )纵向钢筋:PIPE20 (iii )横向箍筋:PIPE20 2.2材料性质 i )、混凝土材料 [16~ 19] 混凝土立方体抗压强 度 f cu ( N / mm 2) 弹性模量E c 2 (N/mm ) 泊松 比 V 单轴抗压强度f c ' 2 (N/mm ) 单轴抗拉强度f r (N/mm ) 裂缝间剪力 传递 系数P t 张开 闭合 30 24000 0.20 25.0 3.1125 0.35 0.75 ?单轴受压应力-应变曲线(二-;曲线) 在ANSYS @程序分析中,需要给出混凝土单轴受压下的应力应变曲线。在本算例中,混凝土单 轴受压下的应力应变采用 Sargin 和Saenz 模型[17,18]: ①22①22 E20 ①22 RCBEAM-03 图5-9 2①82①82①8 2①82①82①8 ① 8@75@75@75 2①22①22①22 150 150 150 150 RCBEAM-01 150150150 RCBEAM-02 (b )、梁断面图 梁尺寸、配筋及荷载示意图 f ①24 ①24 ①22 150 150 ■4- ------------- P P 125 1200 600 (a )、梁的几何尺寸及荷载示意图 600
ANSYS塑性变形模拟例子
/一个周边简支的圆盘,其中心受到一个冲杆的周期作用(假定冲杆是刚性的),需要进行圆盘在冲杆的周期作用下的塑性分析。本实例的模型简图如图19.1所示,材料特性如下所示,塑性时的应力-应变关系如表19.1,载荷历史如表19.2所示。 弹性模量:EX=70000,泊松比:NUXY=0.325 /PREP7 /TITLE,Circular Plate Loaded by a Circular Punch - Kinematic Hardening !* 下面定义建模分析时需要的参数 EXX=70000 RPL=65 RPU=5 H=6.5 STS1=55 STN1=STS1/EXX STS2=112 STN2=0.00575 STS3=172 STN3= 0.02925 STS4=241 STN4= 0.1 NEX=15 NET=2 NEX1=nint(0.8*NET) NEX2=NEX-NEX1 !* ET,1,42,,,1 !定义单元PLANE42,设置为轴对称 !* MP,EX,1,EXX !定义材料属性 MP,NUXY,1,0.325 !* TB,KINH,1,1,4, !定义多线性随动强化准则 TBPT,,STN1,STS1 TBPT,,STN2,STS2 TBPT,,STN3,STS3 TBPT,,STN4,STS4 !* 创建节点 N,1,RPL,,,,,, N,2,0,,,,,, N,3,,H/2,,,,, !* 创建关键点 K,1,,-(H/2),, K,2,RPU,-(H/2),, K,3,RPL,-(H/2),, KGEN,2,ALL, , , ,H, ,3,0 !复制并平移关键点
七个ansys经典入门实例
“有限元分析及应用”课程有限元分析软件ANSYS6.xed 上机指南 清华大学机械工程系 2002年9月
说明 本《有限元分析软件ANSYS6.1ed:上机指南》由清华大学机械工程系石伟老师组织编写,由助教博士生孔劲执笔, 于2002年9月完成,基本操作指南中的所有算例都在相应的软件系统中进行了实际调试和通过。 本上机指南的版权归清华大学机械工程系所有,未经同意,任何单位和个人不得翻印。
目录 Project1 简支梁的变形分析 (1) Project2 坝体的有限元建模与受力分析 (3) Project3 受内压作用的球体的应力与变形分析 (5) Project4 受热载荷作用的厚壁圆筒的有限元建模与温度场求解 (7) Project5 超静定桁架的有限元求解 (9) Project6 超静定梁的有限元求解 (11) Project7 平板的有限元建模与变形分析 (13)
Project1 梁的有限元建模与变形分析 计算分析模型如图1-1 所示, 习题文件名: beam。 NOTE:要求选择不同形状的截面分别进行计算。 梁承受均布载荷:1.0e5 Pa 图1-1梁的计算分析模型 梁截面分别采用以下三种截面(单位:m): 矩形截面:圆截面:工字形截面: B=0.1, H=0.15 R=0.1 w1=0.1,w2=0.1,w3=0.2, t1=0.0114,t2=0.0114,t3=0.007 1.1进入ANSYS 程序→ANSYSED 6.1 →Interactive →change the working directory into yours →input Initial jobname: beam→Run 1.2设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK 1.3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete… →Add… →select Beam 2 node 188 →OK (back to Element Types window)→Close (the Element Type window) 1.4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural→Linear→Elastic→Isotropic→input EX:2.1e11, PRXY:0.3→OK 1.5定义截面 ANSYS Main Menu: Preprocessor →Sections →Beam →Common Sectns→分别定义矩形截面、圆截面和工字形截面:矩形截面:ID=1,B=0.1,H=0.15 →Apply →圆截面:ID=2,R=0.1 →Apply →工字形截面:ID=3,w1=0.1,w2=0.1,w3=0.2,t1=0.0114,t2=0.0114,t3=0.007→OK
ansys二次开发及实例
ansys二次开发教程+实例 第3章ANSYS基于VC++6.0的二次开发与相互作用分析在ANSYS中的实现 3.1 概述 ANSYS是一套功能十分强大的有限元分析软件,能实现多场及多场耦合分析;是实现前后处理、求解及多场分析统一数据库的 一体化大型FEA软件;支持异种、异构平台的网络浮动,在异种、异构平台上用户界面统一、数据文件全部兼容,强大的并行计算功能 支持分布式并行及共享内存式并行。该软件具有如下特点: (1) 完备的前处理功能 ANSYS不仅提供了强大的实体建模及网格划分工具,可以方便地构造数学模型,而且还专门设有用户所熟悉的一些大型通用有 限元软件的数据接口(如MSC/NSSTRAN,ALGOR,ABAQUS等),并允许从这些程序中读取有限元模型数据,甚至材料特性和边 界条件,完成ANSYS中的初步建模工作。此外,ANSYS还具有近200种单元类型,这些丰富的单元特性能使用户方便而准确地构建出 反映实际结构的仿真计算模型。 (2) 强大的求解器 ANSYS提供了对各种物理场量的分析,是目前唯一能融结构、热、电磁、流体、声学等为一体的有限元软件。除了常规的线性、 非线性结构静力、动力分析外,还可以解决高度非线性结构的动力分析、结构非线性及非线性屈曲分析。提供的多种求解器分别适用于 不同的问题及不同的硬件配置。 (3) 方便的后处理器 ANSYS的后处理分为通用后处理模块(POST1)和时间历程后处理模块(POST26)两部分。后处理结果可能包括位移、温度、应力、应变、速度以及热流等,输出形式可以有图形显示和数据列表两种。 (4) 多种实用的二次开发工具 ANSYS除了具有较为完善的分析功能外,同时还为用户进行二次开发提供了多种实用工具。如宏(Marco)、参数设计语言(APDL)、用户界面设计语言(UIDL)及用户编程特性(UPFs),其中APDL(ANSYS Parametric Design Language)是一种非常类似于Fortran77的参数化设计解释性语言,其核心内容为宏、参数、循环命令和条件语句,可以通过建立参数化模型来自动完成一些通用性强的任务;UIDL(User Interf ace Design Language)是ANSYS为用户提供专门进行程序界面设计的语言,允许用户改变ANSYS的图形用户界面(GUI)中的一些组项,提供了一种允许用户灵活使用、按个人喜好来组织设计ANSYS图形用户界面的强有力工具;UPFs(User Programmable Features)提供了一套Fortran77函数和例程以扩展或修改程序的功能,该项技术充分显示了ANSYS的开放体系,用户 不仅可以采用它将ANSYS程序剪裁成符合自己所需的任何组织形式(如可以定义一种新的材料,一个新的单元或者给出一种新的屈服 准则),而且还可以编写自己的优化算法,通过将整个ANSYS作为一个子程序调用的方式实现。 鉴于上述特点,近几年来,ANSYS软件在国内外工程建设和科学研究中得到了广泛的应用。但这些应用大多局限于直接运用ANSYS软件进行实际工程分析,对利用ANSYS提供的二次开发工具进行有限元软件设计却很少涉及。本文首次利用ANSYS软件的二次开发功能,以VC++6.0为工具,运用APDL语言,对ANSYS进行二次开发,编制框筒结构-桩筏基础-土相互作用体系与地震反应分析程序。 3.2 程序设计目标 针对某一实际工程问题,ANSYS所提供的APDL语言可对ANSYS软件进行封装。APDL语言即ANSYS软件提供的参数化设计 语言,它的全称是ANSYS Parametric Design Language。使用APD L语言可以更加有效地进行分析计算,可以轻松地进行自动化工作(循环、分支、宏等结构),而且,它是一种高效的参数化建模手段。使用APDL语言进行封装的系统可以只要求操作人员输入前处理 参数,然后自动运行ANSYS进行求解。但完全用APDL编写的宏还存在弱点。比如用APDL语言较难控制程序的进程,虽然它提供了 循环语句和条件判断语句,但总的来说还是难以用来编写结构清晰的程序。它虽然提供了参数的界面输入,但功能还不是太强,交互性 不够流畅。针对这种情况,本文用VC++6.0开发框筒结构-桩筏基础-土相互作用有限元分析程序(简称LW S程序)。
纷享销客使用规定3.0
关于纷享逍客的使用规定3.0 一、日志(每日工作计划) 1、总部:总部人员每日下班前应填写工作日志,填写内容为当日的主要工作内容及完成情况、第二天的工作计划内容,发送时间不得迟于当天24:00,逾期发送乐捐人民币20元/次。 2、营销部:各项目经理应于当天将工作总结及次日工作计划发送至“企信项目组”,逾期发送乐捐人民币20元/次。 二、周计划 1、填写人员:各项目经理 2、填写时间:每周日下班前填写完毕,发送时间不得迟于当天24:00,逾期发送乐捐人民币50元/次。 3、填写格式:填写本周主要工作内容及完成情况及下周的工作计划内容,工作内容要求客观详实,条理清楚。周表报统计时间为周一至周日,原“销售周报表”及“销控表”以附件形式添加并发送。 4、抄送范围:余总、林总、营销总监、项目策划、行政点评人:营销总监 三、月计划 1、填写人员:各项目经理 2、填写时间:每月初填写完毕,发送时间不得迟于5日24:00,逾期发送乐捐人民币50元/次。 3、填写格式:填写本月主要工作内容及完成情况及下月的工作计划内容,工作内容要求客观详实,条理清楚。月表报统计时间为1-31日,原“销售月报表”、“绩效考核”及“销控表”以附件形式添加并发送。 4、抄送范围:余总、林总、营销总监、项目策划、行政点评人:营销总监 四、案场会议记录 1、各项目经理应组织销售人员召开早会/晚会,并将早会视频发送至微信“阳光项目管理群”及晚会会议记录拍照发送至“纷享逍客企信的营销经理群”;经理休息或开会需交待其他人员主持早会,未完成视为经理监督不到位。 2、发送时间:早会视频在9:30前发送,全体上班成员均应参加,特殊情况未参加且未报备说明原因的个人及经理各乐捐50元/次,晚会记录应在当天 24:00前发送,逾期发送经理乐捐50元/次。 五、请假审批 1、所有人员请假需按照格式填写请假的起止时间及请假事由。 2、总部人员请假:抄送范围为总部,第一审批人为其直属上级,超过天数者第一审批人审批后应再次提交上一级审批人。 3、销售人员请假:抄送范围为“行政部”+"所在项目",请假时间1天内审批人为“项目经理”,超过1天者,由项目经理审批后需再次提交“营销总监”审批; 4、项目经理请假:抄送范围为“行政部”+"所在项目",审批人为“营销总监”,并由营销总监视具体情况决定是否需要提交上级审批。 5、请假必须提前1天填写《请假条》,获得批准后方可休假。未经批准视为无效请假,擅自休假的视为旷工。
ansys工程实例(4经典例子)解析
输气管道受力分析(ANSYS建模) 任务和要求: 按照输气管道的尺寸及载荷情况,要求在ANSYS中建模,完成整个静力学分析过程。求出管壁的静力场分布。要求完成问题分析、求解步骤、程序代码、结果描述和总结五部分。所给的参数如下: 材料参数:弹性模量E=200Gpa; 泊松比0.26;外径R?=0.6m;内径R?=0.4m;壁厚t=0.2m。输气管体内表面的最大冲击载荷P为1Mpa。 四.问题求解 (一).问题分析 由于管道沿长度方向的尺寸远大于管道的直径,在计算过程中忽略管道的端面效应,认为在其长度方向无应变产生,即可将该问题简化为平面应变问题,选取管道横截面建立几何模型进行求解。 (二).求解步骤 定义工作文件名 选择Utility Menu→File→Chang Jobname 出现Change Jobname对话框,在[/FILNAM] Enter new jobname 输入栏中输入工作名LEILIN10074723,并将New log and eror file 设置为YES,单击[OK]按钮关闭对话框 定义单元类型 1)选择Main Meun→Preprocessor→Element Type→Add/Edit/Delte命令,出现Element Type 对话框,单击[Add]按钮,出现Library of Element types对话框。 2)在Library of Element types复选框选择Strctural、Solid、 Quad 8node 82,在Element type reference number输入栏中出入1,单击[OK]按钮关闭该对话框。 3. 定义材料性能参数 1)单击Main Meun→Preprocessor→Material Props→Material models出现Define Material Behavion 对话框。选择依次选择Structural、Linear、Elastic、Isotropic选项,出现Linear Isotropic Material Properties For Material Number 1对话框。 2)在EX输入2e11,在Prxy输入栏中输入0.26,单击OK按钮关闭该对话框。 3)在Define Material Model Behavion 对话框中选择Material→Exit命令关闭该对话框。 4.生成几何模型、划分网格 1)选择Main Meun→Preprocessor→Modeling→Create→Areas→Circle→Partail→Annulus出现Part Annulus Circ Area对话框,在WP X文本框中输入0,在WP Y文本框中输入0,在Rad1文本框中输入0.4,在Theate-1文本框中输入0,在Rad2文本框中输入0.6,在Theate-2文本框中输入90,单击OK按钮关闭该对话框。 2)选择Utility Menu→Plotctrls→Style→Colors→Reverse Video,设置显示颜色。 3)选择Utility Menu→Plot→Areas,显示所有面。 4) 选择Main Menu→Preprocessor→Modeling→Reflect→Areas,出现Reflect Areas拾取菜
ANSYS应用实例:钢筋混凝土简支梁数值模拟
(ii )纵向钢筋:PIPE20 (iii )横向箍筋:PIPE20 2.2 材料性质 (i )、混凝土材料 表5-4 混凝土材料的输入参数一览表[16~19] ·单轴受压应力-应变曲线(εσ-曲线) 在ANSYS ○R 程序分析中,需要给出混凝土单轴受压下的应力应变曲线。在本算例中,混凝土单轴受压下的应力应变采用Sargin 和Saenz 模型[17,18]: 2 21??? ? ??+???? ??-+= c c s c c E E E εεεεε σ (5-30)
式中取4' 4')108.0028.1(c c c f f -=ε;
断面图配筋图断面图配筋图断面图配筋图RCBEAM-01 RCBEAM-02 RCBEAM-03 图5-12 各梁FEM模型断面图 (a)单元网格图(b)钢筋单元划分图 图5-13 算例(一)的FEM模型图 2.4 模型求解 在ANSYS○R程序中,对于非线性分析,求解步的设置很关键,对计算是否收敛关系很大,对于混凝土非线性有限元分析,在计算时间容许的情况下,较多的求解子步(Substeps)或较小的荷载步和一个非常大的最大子步数更容易导致收敛[2]。在本算例中,设置了100个子步。最终本算例收敛成功,在CPU为P41.6G、内存为256MB的微机上计算,耗时约为8小时。 2.5 计算结果及分析 2.5.1 荷载—位移曲线 图5-14为ANSYS○R程序所得到的各梁的荷载-跨中挠度曲线,从图中可以看出: (i)、梁RCBEAM-01:曲线形状能基本反映钢筋混凝土适筋梁剪切破坏的受力特点,而且荷载-跨中挠度曲线与钢筋混凝土梁的弯剪破坏形态非常类似,即当跨中弯矩最大截面的纵筋屈服后,由于裂缝的开展,压区混凝土的面积逐渐减小,在荷载几乎不增加的情况下,压区混凝土所受的正应力和剪应力还在不断增加,当应力达到混凝土强度极限时,剪切破坏发生,荷载突然降低。
混凝土参数表
混凝土结构设计规范GB50010-2002第条。 C30混凝土受压和受拉时的弹性模量为:(10)4 N/mm2,即30KN/mm2. 000N/m2=3*1010pa=3*104Mpa=30GPa 2500Kg/m3 泊松比为 梁采用C40混凝土,弹性模量E=33GPa,密度γ=2500kg/m3,泊松比为 2500Kg/m3 )))))))))))))))))))))))))))))))))))))))))))))))))))))) 管片结构采用C50混凝土,弹性模量为35 GPa泊松比为,, 密度γ=2500kg/m3 混凝土强度等级为C25,重度r= kN/m3,弹性模量E= GPa,泊松比μ= 2500Kg/m3 $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 混凝土为C20,弹模E=26000Mpa,,,泊松比, 2500Kg/m3 1Gpa=1000Mpa 1Mpa=1000000pa 1Gpa=00pa ━━━━━━━━━━━━━━━━━━━━━━━━━ 名称弹性模量 E 切变模量 G 泊松比μ GPa GPa ───────────────────────── 镍铬钢 206 合金钢 206 碳钢 196-206 79 铸钢 172-202 球墨铸铁 140-154 73-76
灰铸铁 113-157 44 白口铸铁 113-157 44 冷拔纯铜 127 48 轧制磷青铜 113 41 轧制纯铜 108 39 轧制锰青铜 108 39 铸铝青铜 103 41 冷拔黄铜 89-97 34-36 轧制锌 82 31 硬铝合金 70 26 轧制铝 68 25-26 铅 17 7 玻璃 55 22 混凝土 14-23 纵纹木材 横纹木材 橡胶 电木 尼龙 可锻铸铁 152 拔制铝线 69 大理石 55 花岗石 48 石灰石 41 尼龙1010 夹布酚醛塑料 石棉酚醛塑料 高压聚乙烯 低压聚乙烯 聚丙烯
ANSYS 中使用函数加载的一个简单例子
ANSYS 中使用函数加载的一个简单例子 本文将通过一个具体实例说明在ANSYS 中如何使用函数加载,后续将通过该实例在分析过程中遇到的一个问题提出自己的一点看法。 实例的具体说明: 一个1/4 圆柱,内半径30 mm,外半径42 mm,长度100mm,如图1 所示: 所用材料为双线性弹塑性材料,其机械性能为: 弹性模量 E = 201000 Mpa;泊松比μ=0.3
屈服应力σ= 200 Mpa;切线模量Et = 2010使用单元类型solid185 (8 节点六面体单元)。 取整体单元边长4 mm,然后可以直接对该几何模型划分MAP 网格,划分网格结果如图2: 约束条件为: 轴向两个截面为对称边界条件;一个端面约束轴向位移Uz。 载荷条件为: 在外表面施加变化的压力载荷,载荷函数为: P (y) = 8e7 + 7E7 * (Y/42)
即: X = 0 ,Y = 42 (最高点) 时,P = 15E7; X = 42,Y = 0 (最低点)时,P = 8E7。 我们采用函数方式来施加这一压力载荷,首先定义函数: 在Solution 模块中,点击菜单路径: Solution > Define Loads > Apply > Functions > Define/Edit 将会弹出一个函数编辑器,可以在其中定义所需的函数。 在函数编辑器中,函数类型选择为Single equation,即单值函数;计算函数值时使用的插值坐标系( (x,y,z) interpreted in CSYS) 选择0,即总体直角坐标系,如图3 所示:
然后,在函数编辑器中间位置的“Result = “ 小窗口中输入要定义的函数表达式,如果表达式中有x, y, z, time 等变量(供定义函数时使用的“自变量”),可以用{X},{Y},{Z},{TIME} 等的形式输入;或者点击下面一个小窗口右边的小箭头,会出现一个下拉列表,列出可以选择的变量,然后从该列表中选择某个自变量,则该自变量会按照上述格式写入函数中,如图5 所示:
格构梁的ANSYS有限元模拟分析实例运用
龙源期刊网 https://www.360docs.net/doc/d018448025.html, 格构梁的ANSYS有限元模拟分析实例运用作者:张少剑刘真 来源:《城市建设理论研究》2013年第10期 摘要:本文通过一工程实例运用ansys模拟计算。针对格构梁的研究,合理地简化模型,取出1.5米宽的土体、梁和面层单元,两边加对称约束,从而达到模拟空间结构梁的目的。本文还模拟了基坑的开挖过程的时空效应,共分七步,土体在自重应力作用下的沉降为第一步,梁与面层的激活、力的施加和土层杀死共分六步。梁的最大受力状态并不发生在最后一步完成后,而是在第六工况。 关键词:格构梁有限元分析模拟分析 中图分类号:K826.16 文献标识码:A 文章编号: 1 土体、梁、锚索和混凝土面层共同作用 基坑支护的受力机理是土体的土压力作用在格构梁和混凝土面层上,混凝土面层的力传递到格构梁上,格构梁再把它受到的力传递到和它相连的锚索上,锚索则和被支护土体嵌固为一体,格构梁和混凝土面层除起到承受土压力外,格构梁还起到平均弯矩和变形的作用,喷射混凝土面层则有保护土体表面,防止土体表面非格构梁作用部位坍塌的作用。 2模型简化及技术处理 根据基坑开挖深度,根据实际的土体性质建立土体模型。格构梁的作用是承受弯矩的,可以选用Beam4梁单元,考虑到钢筋混凝土格构梁中有钢筋的作用,其弹性模量、泊松比等设置有所调整。在建模时,如果混凝土面层的长宽与厚度的比都大于5,所以在有限元分析中采用板壳单元可以全面地反映其变形特征和应力分布规律。混凝土面层用Shell63单元模拟,其参数的取值和梁单元相同。 由于格构梁的受力性状,锚索的模拟对格构梁的受力影响较小,本模型忽略考虑锚索的模拟。预应力锚索的作用简化为作用在纵横梁交点处的集中力。 对于格构梁和土体、混凝面层之间的接触,模型采用节点耦合,以实现共同变形和受力。 3.1ANSYS有限元模拟计算 3.1.1模型的参数 1.土体的参数见下表:
纷享销客用户使用手册
纷享销客用户使用手册 快速上手版-管理后台 外勤签到设置.................................................................................. 考勤签到设置.................................................................................. 报数系统后台管理 ............................................................................ 版本说明: 版本:V4.7.0
日期: 版本更新: 创建企业 管理员登录进入企业后,网页导航栏包括功能模块:工作、CRM和应用,管理员还会拥有管理后台入口,点击[管理]进入管理后台,需要再次输入密码,确保后台管理的安全。
管理员可通过管理后台进行部门与员工维护、权限设置、账号安全管理、企业信息维护、删除信息管理、公告管理、日志统计等;纷享提供用户分类管理员设定,不同权限的管理员分门别类的 对各自模块继续维护和管理。 [1]管理后台首页展示员工统计数据,待审核的邀约,待分配部门的员工,未激活的员工,今 日登录的员工数; [2]安全数据统计,显示设备数据,和异常帐号数据 [3 [4] [1] [2] [3]添加员工:点击进行员工新建,填写姓名,昵称(默认与姓名一致),手机,登录帐 号(默认与手机号一致),设置初始密码,一个员工可从属多个部门,填写
其职务,选择直属上级(审批审核),新建用户默认为办公版,如当前企业购买了营销版 帐号,则可进行帐号设置。 [4]导入:点击导入,按照弹窗中的操作步骤提 示进行成员的批量导入操作,下载模版,补 充完成后上传即可;还可勾选是否发送短信 通知已创建的用户,方便员工及时收到通知[5] 传说中的小黄条 [6]员工列表信息包括:姓名,帐号,部门,职务,手机,创建时间及账户版本类型 (批量)选择员工记录时,可进行如上图操作:禁止登录,允许登录,批量编辑归属部门,重置密码,停用帐号, 激活提醒,设置版本类型。
几个ansys经典实例(长见识)
平面问题斜支座的处理 如图5-7所示,为一个带斜支座的平面应力结构,其中位置2及3处为固定约束,位置4处为一个45o的斜支座,试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台,分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2,3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题,使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu:Preprocessor →Coupling/ Ceqn →Constraint Eqn :Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu:WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2,0,0,THXY:45 →OK ANSYS Main Menu:Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu:Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流; !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系,进行旋转,施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end
c50混凝土弹性模量表格
竭诚为您提供优质文档/双击可除c50混凝土弹性模量表格 篇一:c50混凝土配比计算书 混凝土配合比试验计算单 c50混凝土配合比计算书 一、设计依据 tb10425-94《铁路混凝土强度检验评定标准》 tb10415-20xx《铁路桥涵工程施工质量验收标准》jgj55-20xx《普通混凝土配合比设计规程》tb10005-20xx《铁路混凝土结构耐久性设计规范》tb10424-20xx《铁路混凝土工程施工质量验收标准》gb/t50080-20xx《普通混凝土拌合物性能试验方法标准》gb/t50081-20xx《普通混凝土力学性能试验方法标准》 gb/t50082-20xx《普通混凝土长期性能和耐久性能试验方法标准》设计图纸要求二、技术条件及参数限值设计使用年限:100年;设计强度等级:c50;要求坍落度:160~200mm;胶凝材料最小用量360kg/m3;最大水胶比限值:0.55; 耐久性指标:56d电通量<1000c;
三、原材料情况 1、水泥:徐州丰都物资贸易有限公司,p·o42.5(试验报告附后) 2、粉煤灰:中铁十五局集团物资有限公司,F 类Ⅱ级(试验报告附后) 3、砂子:(试验报告附后) 4、碎石:5~31.5mm连续级配碎石,5~10mm由石场生产;10~20mm由石场生产;16~31.5mm由石场生产;掺配比例5~10mm为30%;10~20mm为50%;10~31.5mm为20%(试验报告附后) 5、外加剂:山西桑穆斯建材化工有限公司,聚羧酸高性能减水剂(试验报告附后) 6、水:混凝土拌和用水(饮用水)(试验报告附后) 四、设计步骤 (1)确定配制强度 根据《普通混凝土配合比设计规程》jgj55—20xx、《铁路桥涵工程施工质量验收标准》tb10415-20xx,混凝土的配制强度采用下式确定: fcu,0fcu,k1.645501.6456.059.(9mpa) (2)按照《铁路混凝土结构耐久性设计设计规范》 tb10005-20xx规定,根据现场情况: 1、成型方式:混凝土采用罐车运输,混凝土泵送施工工艺。 2、环境作用等级:l2、l 3、h3、h 4、t1、t2。 3、粉煤灰掺量要求:水胶比≤0.55,粉煤灰掺量要求
一个经典的ansys热分析实例(流程序)
/PREP7 /TITLE,Steady-state thermal analysis of pipe junction /UNITS,BIN ! 英制单位;Use U. S. Customary system of units (inches) ! /SHOW, ! Specify graphics driver for interactive run ET,1,90 ! Define 20-node, 3-D thermal solid element MP,DENS,1,.285 ! Density = .285 lbf/in^3 MPTEMP,,70,200,300,400,500 ! Create temperature table MPDATA,KXX,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12 ! 指定与温度相对应的数据材料属性;导热系数;Define conductivity values MPDATA,C,1,,.113,.117,.119,.122,.125 ! Define specific heat values(比热) MPDATA,HF,2,,426/144,405/144,352/144,275/144,221/144 ! Define film coefficient;除144是单位问题,上面的除12也是单元问题 ! Define parameters for model generation RI1=1.3 ! Inside radius of cylindrical tank RO1=1.5 ! Outside radius Z1=2 ! Length RI2=.4 ! Inside radius of pipe RO2=.5 ! Outside pipe radius Z2=2 ! Pipe length CYLIND,RI1,RO1,,Z1,,90 ! 90 degree cylindrical volume for tank WPROTA,0,-90 ! 旋转当前工作的平面;从Y到Z旋转-90度;;Rotate working plane to pipe axis CYLIND,RI2,RO2,,Z2,-90 ! 角度选择在了第四象限;90 degree cylindrical volume for pipe WPSTYL,DEFA ! 重新安排工作平面的设置;另外WPSTYL,STAT to list the status of the working plane;;Return working plane to default setting BOPT,NUMB,OFF ! 关掉布尔操作的数字警告信息;Turn off Boolean numbering warning VOVLAP,1,2 ! 交迭体;Overlap the two cylinders /PNUM,VOLU,1 ! 体编号打开;Turn volume numbers on /VIEW,,-3,-1,1
ansys实例-正确地模拟过盈配合
过盈配合在机械产品的装配中使用的相当普遍。比如轴与轴承、轴与轴瓦、汽车的制动盘等,都是通过一定的过盈量来使两个装配部件紧密地连接起来。 下面讨论如何在ANSYS 中正确地模拟过盈配合。 过盈配合在有限元分析中是一种典型的非线性接触行为。在有限元分析中设定了接触,从本质上来讲就是对相互接触的两个部件施加了某种约束,不同的接触算法对于接触约束的处理方法有所不同。接触约束的理论算法的选择,在ANSYS 中是通过设置contact 单元的KEOPT(2) 选项来实现的。 在ANSYS 中目前主要有5 种接触约束算法: KEYOPT(2)=0 Augmented Lagrangian - 加强的拉格朗日算法,这是ANSYS 的缺省选择; KEYOPT(2)=1 Penalty function - 罚函算法; KEYOPT(2)=2 Multipoint constraint (MPC) - 多点约束算法; KEYOPT(2)=3 Lagrange multiplier on contact normal and penalty on tangent - 接触法向采用拉格朗日乘子,接触切向采用罚函数的综合算法。 KEYOPT(2)=4 Pure Lagrange multiplier on contact normal and tangent - 法向和切向均采用拉格朗日乘子算法。 各种不同的约束算法各有其优缺点,各有各自最适用的场合,具体情况需要具体对待。大部分情况下,默认选择KEYOPT(2)=0 就够用了。 过盈配合所致的接触分析的难点在于如何确定初始接触状态。初始
接触状态设置得不对,会导致错误的计算结果或者不准确的计算结果,下面举两个例子来说明。ANSYS仿真计算代做:模态分析,瞬态动力学,谐响应分析和谱分析、械结构的疲劳、损伤,CFD流体;结构的强度评估和优化;企鹅:690294845 例1.两个圆柱体在几何上是刚好接触,划分网格后有限元模型有间隙。如图1 所示。 这两个圆柱体,在几何上是刚好相切的,即处于几何上刚好接触的初始状态。划分网格后,由于在圆周上用小段直线代替了弧线,两个圆柱体之间产生了一定的间隙,两个圆柱体的有限元模型的初始状态不再是接触的。此时,如果接触参数设置不当,就会因为初始约束不足,圆柱体出现刚体位移,得到错误的结果。