|
发表于 2013-7-25 15:09:12
|
显示全部楼层
我邮箱1962940284@qq.com5 d% l. k5 v9 m
施加弯矩扭矩的方法其实有很多种方法,在这里介绍其中的5种,并进行比较:
# x; e/ B3 d( y9 Q' H" `! ?# }/ a. c, Y2 I
1.将矩转换成一对的力偶,直接施加在对应的节点上面。
) W8 K$ A. O4 I) ^) {
m ^2 ~& C U0 z* {0 k 2.在构件中心部位建立一个节点,定义为mass21单元,然后跟其他受力节点耦合,形成刚性区域,就是用cerig命令。然后直接加转矩到主节点,即中心节点上面。+ P& Y: h. I; t4 X+ j
7 P( T$ ^( i/ P/ z 3.使用mpc184单元。是在构件中心部位建立一个节点,跟其他受力节点分别形成多根刚性梁,从而形成刚性面。最后也是直接加载荷到中心节点上面,通过刚性梁来传递载荷。" V: b6 `. F5 [6 P2 }9 H
2 {/ z! n6 @$ V R: i6 F, S7 z) H. O! ^ 4.通过rbe3命令。该方法与方法2很接近。
) H9 _1 f4 l5 ]# Y, a
3 y" k. Y7 t9 r0 X" h. s T; T& d; D 5.基于表面边界法:主要通过定义一个接触表面和一个目标节点接触来实现,弯矩荷载可以通过在目标节点上用“F”命令施加。
0 f* d- e6 ~' K$ @0 m
" n) V) x; @, u4 J# r }' Z, m 对于方法1,通过转换为集中力或均布力,比如施加扭矩,把端面节点改成柱坐标,然后等效为施加环向的节点力;而施加弯矩,可以将力矩转化为端面的剪切均布力;但这种方法比较容易出现应力集中现象;$ \& V6 V1 i6 ?8 a
- }% z1 N1 _' R# D: z 方法2,定义局部刚性区域,施加过程venture讲的很详细,这里就不在赘述。根据他的例子,我在下面给出了一段命令流。该方法有个不足,它在端面额外的增加了一定的刚度,只能适用于小变形分析。# ?% X) H& x5 ]0 P ~2 |/ |' t& r8 v7 T2 _
. g* b* ?% X4 s) u; y0 B* r0 `+ B3 r
方法3 ,相对方法2来说,采用刚性梁单元,适用范围更广一些,对于大应变分析也能很好的适用。但在小应变分析下,方法2和方法3没有什么区别。' w+ G$ n9 G% x; ]( q+ l
' |+ U8 ^, k# m1 k 方法4,定义一个主节点,施加了分布力面,应该说跟实际比较接近一点,但端面的结果好像不是很理想,结果有点偏大,在远离端面处的位置跟实际很符合。
7 \* ^+ b8 |& t/ \3 W j7 S
& O- { S. P2 H, O6 b& X% G 方法5,它具体的受力形式有如下两种:
! S- a6 l$ ^4 D O {) S, X2 S9 z
刚性表面边界(Rigid surface constraint)-认为接触面是刚性的,没有变形,和通过节点耦合命令CERIG比较相似;" t6 r! S. U' R- R
, ]' [ b1 F! q8 \9 q1 f2 v 分布力边界(Force-distributed constraint)-允许接触面的变形,和边界定义命令RBE3相似。5 y" s; z$ F# Y* F' [
8 V, M% ^8 L/ b3 u! f* E' i
使用这种方法,需要用KEYOPT(2) = 2打开接触单元的MPC(多点接触边界)算法,; @0 ^) N0 w$ Q: E/ O) _: a( ]9 c
! l: S: \1 P4 G# }8 h8 y9 u4 n: P 下面针对venture给出的例题,用不同的方法来实现的命令流。
, m8 O, C" g# O; R. m6 R8 v! U$ Q4 h, k8 c! A
方法1不介绍了,方法2:( o! I8 o& j' [4 ~& ~
* ?, ~# R/ N$ O* {1 U! h' ^1 Y
/PREP7
& t5 S6 u4 l- U2 [& H# k7 G6 _# D0 Y9 d, D2 D7 M: p' w3 v
ET,1,95. n" q& x0 k) A& P+ P, n' {9 f: C
! T/ D% O4 A! q0 [+ n/ D
ET,2,211 d# ^# ^& J, j8 O
7 d' E, ]0 [# L. T9 G KEYOPT,2,3,0+ O) G( }: J1 V& Q+ a% ~
! `+ Q0 o; B2 S7 h' a9 }) B' ~ R,1,1E-6
8 Q7 h5 M/ }# ?) S$ d. v. E. [& j& d$ D' s1 E) s- Z* y5 j
MP,EX,1,2.01e5* y4 t# F# c+ [, B$ C" A4 J
! h2 m( f G+ W. _ MP,PRXY,1,0.3( `+ m B/ M1 R7 u! c3 V
; m. w4 g- @( L8 i
CYLIND,15,10,0,200,0,360, W# Q; b! t3 G2 Y$ @
' R4 a, x; C8 | E5 E) D m0 I2 ~0 @+ @ wpro,,90,
$ t; L8 o* s3 e1 o( M
, e! h3 L* Y: n- _5 V* s& c1 H vsbw,all
* |% ~6 T% U9 o! _ p8 J6 x4 c- o9 F# c' I$ _" r" V4 d) R) |, w
wpro,,,90- K0 ]& N: g$ u8 I! ~
+ n) [* d0 M8 _6 [( f' U. F; }( S3 b vsbw,all8 ?* _: { }$ E" [) O# T/ c0 K' t
/ t' ?: W3 v, ^) }5 u9 U* ]9 W, X
WPCSYS,-1,09 g% H: ^' r$ t9 @- K6 H6 A1 |$ E
" F4 n" d' P! G( w$ q# R; J K,17 , , ,2100 E2 C/ K0 _- q
& {* V- M) }. b% f0 o$ n
lsel,s,,,13,16,11 k5 X0 |& g$ O9 e
0 J1 m" N- T) d% b* s! \ lesize,all, , ,8, , , , ,1
4 _* H8 }( _# R! Q. ?
; O9 P0 D( s/ A+ G: z, s$ c% a lsel,s,,,22( S$ K; G$ ^! C
# m4 f5 r. e, O: m
lesize,all, , ,4, , , , ,19 r1 Q p0 k( i# k& [3 N# O) ?
J. \( J3 y$ K3 m% W7 R: j' ?) f0 G
lsel,s, , ,17,20,1
. ~: b( k5 z( f+ I0 O4 @4 b' ]6 {5 t' h0 @# V
lsel,a, , ,26,27,1* S9 t& j% u4 r2 m
) C' |' |8 b7 a0 E$ K
lsel,a, , ,30,31,1, Z$ K7 X& I+ M0 G2 `
. c8 F0 i* E& d: C0 ^
lesize,all, , ,20,0.4, , , ,1
+ @; o6 C9 M. t- H, c0 T% |4 U6 I, p1 t( ^
alls( e, u. v E, c- H" I; I
0 ~- c9 X3 u! n2 V% d1 {5 k vmesh,all, q2 \1 a' n6 U/ @1 A/ E! I! |7 K3 t
% w* t' b2 m; b0 p/ B' j/ @) j$ Y !!!!!下面一段开始各个方法有所不同,由于前面的建模一样,后面的例子就不再给出
: V$ s( V3 ~- x1 X+ o, Z' w5 U1 g, ~1 @5 g
ksel,s,,,17
4 I4 ^4 N9 a k$ _. P4 J4 B# `- A5 W2 L% S( R
type,2
5 a, i$ u; r5 T W2 P
# y( m% r7 ^! Q2 k1 d3 z2 L real,1/ Z/ G, u' S1 u0 w
/ f" ^6 ?1 X5 p7 R: ^) N% K kmesh,all u" G' {; s: E4 x2 O5 U' Q. A: L
: n6 p% H$ L9 @0 m# Y, L$ I allsel
# t. J! c7 G# f* |% |% {
+ b3 t+ h e, c, v nsel,s,loc,z,200,210
/ d( ?- Q! u: ?% z4 E( @- c( G4 `* A) M% p$ ]
npolt' }* G% K s4 _3 G0 k
* b, O, a( }" O8 s1 G# Z
CERIG,node(0,0,210),ALL,ALL, , , ,' I# z1 w1 X0 g' m8 H5 v8 ]$ |: X6 N5 V
+ O t3 Y. ?* D !!!!!CERIG命令定义局部刚性区域
3 @( N9 w& p; R- T$ N: `6 }+ B% C+ L4 m
allsel
; z/ V$ n, A4 m* [" j6 A0 g) X
% S( K/ U/ f* F" | /SOLU
% x2 n3 {0 E9 [+ \# M* ~. i# c. k( \5 o! X: v' Z
f,node(0,0,210),mz,10e5
. e4 ]3 z+ J2 X% C
+ y! \) C+ |/ J# S' f FINISH% C7 P' j3 K* m! P7 l, R" @! X
/ p( [6 j# u% H$ Z) F% \3 c
!!!!!以下一段边界条件的施加各种方法一样,后面例子也不再赘述
4 ?0 y, ?8 ~0 U$ B0 R. y: R) k4 s
# d# B0 ^2 x# I$ O& h X /SOL4 M9 ~9 g. a4 S, M: s6 o3 F: E
1 `- j0 ^. d) A$ ~, q nsel,s,loc,z,0
' a7 j% d# Y: i( n, y# D& W4 w- M) m4 y
9 i8 a9 |* w" w4 o K; a d,all,all
/ F" ]1 k+ ?0 f+ y
7 ^) p- I/ O- | allsel
/ I( @& \+ E2 X m& G1 n/ g+ P& f% p% R K" v; e: \
solve
5 Q7 x+ j) O/ B0 p/ E s, h9 U o+ I5 g Z; _! p- ?9 W
方法3:使用MPC184单元定义刚性梁- p! F5 k" S1 O( g5 \! a
& z8 v& ?( o9 T V: C- g
……, y9 w: k' p+ T1 E1 v3 u8 T. x3 F }
- X2 Y: b+ n$ q; u1 p et,2,184" G/ \% i0 A1 t' P+ _8 a0 m& v) |
8 ]1 d+ a( p: _3 X! H
keyopt,2,1,1
! c6 v( f+ H7 F; c5 `% u! V& }
: H- I% r0 P& l1 \# ]" a( T nsel,s,loc,z,2001 O9 ^# e$ @, C& N+ [& d4 k) X; b$ c
6 p2 Q& i" S. J
n,15000 ,0,0,210
3 |$ R+ O# k- S8 f* y- n8 q% i! r+ x& ^# H' X3 m+ F, U% b
type,2
! Q" O: c4 j' k; B
$ t( V/ M1 h! V! D/ X1 `7 ]* I *get,nnum,node,0,count
0 Q7 E2 }% b: f% f5 r, r% x2 y% \# t* g* R; i" m1 W( J3 f& u( B
*get,ND,node,0,num,min
, h* u; `9 j: a& G5 q
1 @0 X @5 y" w9 Y8 b *do,i,2,nnum
) k3 p4 q; z9 w; w: B/ T' \
8 T( p' Y! }! q( W+ |( D x !!!!节点个数是nnum,只需要生成nnum个mpc单元& P, K7 i7 H0 D4 g: T4 r
- O; S: q+ j2 _% a! {! j E, 15000,ND
! g- u- ?3 N P1 R! ^7 Y. v
0 M& u- J6 @' d) J0 ]- E5 g ND=NDNEXT(ND)
i' c" `: Z2 B* R9 U. l% ~0 V* ^: m" W- o) ?
*enddo
0 m' r% L* b" o* Z1 M( o. J$ o# q* s3 U" s( `
allsel" A6 d+ L- D: H- b. {4 S
4 y9 O, X4 L/ n /SOLU
+ l1 F! [# x7 F0 N
7 j: j) O0 ~0 c# B f,node(0,0,210),mz,10e59 N e* k6 G* y$ b* D
- v5 [5 M0 @2 s FINISH$ i# q6 l& o# S
( G: G4 |; U, [3 P# i& U
……3 A/ H" o) g' x5 x! \, D j& u
1 Z6 X1 m/ K4 n* b: d
方法4:rbe3命令
0 c' i: e! G: T! Q& J6 r- _/ H$ Y6 p+ [* |: s) m1 j
……
+ C4 D) Z# [% L4 k
1 s2 Z# X! h0 ?$ r7 p ET,2,21- |- i m' } G; s6 L
! O, U8 G8 C$ n j) v1 q KEYOPT,2,3,0& ~: h: l9 s6 r/ c
, N/ b D D; M. o6 Y8 v7 k' N8 t
R,1,1E-6! E2 ]6 x) t4 k4 a, @- q
; M: I7 K! a* ]0 A C
K,17 , , ,210
8 B) w( x( \ X) a# J
P( C1 J$ i( m ksel,s,,,17$ F. i1 }; @# V W/ Q5 @
. W& z* @" O- C8 c" t
type,2
+ }3 }# R+ V; P& v; }+ h2 k
9 x0 w: \% s' s8 F# Q real,14 y& t- s: C ?# k* T$ z4 w/ K/ X
: t* v' F6 ?9 P( r5 v9 S
kmesh,all$ \8 {8 K7 n$ n: ~" ?9 T% I
( W% I; F# e4 m* p0 v1 ]
allsel& m9 }. a( \3 `- J. L2 \: A
$ O3 I- D& `/ }% ? Z nsel,s,loc,z,200
* @2 B* |, \, \; d7 y# ?# i
- N: z' T& P: \+ E *get,nnum,node,0,count
9 m, C1 B5 Y' D+ A5 q& I3 V# X) {- A
*get,ND,node,0,num,min
) P2 i; c9 s, F0 O* D+ z6 `/ T
1 m' j3 T$ h) m: r- T: a *dim,sla,array,nnum
# X- O: v8 W5 m" L- q0 j) o% ]$ r4 @' j
*dim,sla2,array,nnum
9 F$ O1 d* }' h& p0 f/ [; n0 e9 V, G" n- A
*do,i,1,nnum
8 a! N% Y6 D2 K- d
; E) @% `3 o+ P' v% E9 y, S0 O sla(i)=ND, Q/ j6 o9 G8 f5 U
9 Q* w: R) X1 j$ E$ x4 w7 \: Q
sla2(i)=ND
1 z( P( y# v7 f, H6 ~7 Q0 p; ]% j& c
ND=NDNEXT(ND)
$ q1 u1 c5 S. ?4 P/ ^3 \. w3 N. w3 B, P: e& L' n
*enddo
0 K- W3 j' J- S% h/ x1 |1 p
2 E& I, u6 C3 Q7 e2 W allsel8 p* @ c" b7 W T" ]
9 p/ u- c4 D I7 e ^
rbe3,node(0,0,210),all, sla,sla29 }( l7 l$ `" u8 O0 Y4 z* x
/ }3 O: P5 E, ?- q+ V1 {7 i
/SOLU
) b$ v c$ C" W0 u; F' a! F2 h/ O% C8 |" k; Q; U8 c. L4 F. H- H
f,node(0,0,210),mz,10e59 [' x( {' ]/ M7 j
0 K d' F/ j4 [7 F2 W+ G FINISH8 u% i2 x2 g* X0 t5 \
# k9 B8 S9 F" f5 _9 |" a& x ……
" o' i/ B: q* c" K& d1 Q4 w
* w! h, Z1 i' j7 r7 J% ]; p 方法5:定义刚性接触面( y( @# B0 C6 I) e# t
3 w" d f; F8 ^$ P7 l* }+ I8 E$ a4 P
……
S. T4 X, ?1 [7 d; [1 J
& d) x/ B! _9 o" E/ G' R n,15000 ,0,0,200- a v3 `* J- _ k/ e
; U( d- s- {2 |- Z MAT,1: |2 [% T) o/ {# C
; p% }: m6 x% n& E
R,3" v7 G) H+ {* R0 W9 `
4 Y4 @: G r/ [; `% \1 n# } REAL,36 d$ m9 X) w# R& h/ ]5 c ]* y
" g: _4 l8 k& G5 j4 O3 O
ET,2,170
3 ~) `, ~1 _* c- c; p* q8 _3 x4 y4 y7 b5 h
ET,3,175
) Q" K, S) ^4 `. Z
) S2 f8 C/ l8 x. Q KEYOPT,3,12,5
& A: I6 U- q' A" J! p2 M- l8 U7 ~- R+ W- M( n2 E+ _
KEYOPT,3,4,18 v, q% e5 d" E0 d, F" I
8 s7 K, ^+ ]8 }2 m7 d, A3 |+ s KEYOPT,3,2,2
! A4 _2 z6 \1 z% g+ G# w7 ]- R9 N8 B+ [, Z
KEYOPT,2,2,00 Q. O- }2 r8 M0 t9 b* }/ k
; W( n3 j9 J. m3 Y/ q* \" a KEYOPT,2,4,111111
1 m& [, C B7 Y7 X- |6 H: Q: W0 g
" Z+ X/ p9 E9 y3 N3 N) ^) q. o9 Z: @ TYPE,2- R' U, \3 s; [2 f; t0 j
: f9 p/ t) i2 j% ]- [" {* u; W ! Create a pilot node3 x ^ j: z9 k9 `' L) _. h
/ y8 \6 M: g. b2 O1 g% Y TSHAP,PILO9 w1 r3 z0 l5 ~ x, q% O; h
) w, e5 n) ^+ L/ S/ ~7 ~
E,15000
- X" Y! M. R) ^$ U4 f
4 ^& q" c& {" V/ c$ q& g, @ ! Generate the contact surface3 R, r* M# z3 y& g: H1 W! s- _
: `* ]; A( I. V% ~0 b
ASEL,S,,,14
6 O) ~ n- Z1 L" X
% W4 s( s. T& I6 [. F4 R ASEL,A,,,19* P" @# C4 c+ z7 n( s) F
& t4 n& g; z, N( d+ Q$ ]
ASEL,A,,,24" S0 U: |8 U) `$ A' c* L
1 o) A( z0 @' U
ASEL,A,,,28
- i9 T b- ]- A3 E# m3 a# B3 E( [$ @
CM,_CONTACT,AREA0 q. v0 `* a8 o+ w$ g$ _& C
1 O& i" }- o" X. m4 A) Q2 ]. ?( e) r TYPE,3
2 C& g5 T3 j; O, q7 b0 s4 M9 C- g# R
NSLA,S,1
" D1 }' a- z7 L9 r/ k
$ o+ v4 P+ ~+ a# r% W* {# ~' U ESLN,S,01 ]% t6 i+ v8 c8 y3 f
% N d" W2 Z* K, r1 R7 y ESURF
6 S2 r9 Z% u+ o6 g# ~' \' Q9 i
& x3 o9 ^/ a# h5 y4 H X6 n ALLSEL
+ p7 {" f& ~" s% C6 P' P) a/ S
: b" C/ ^1 M7 L& o0 g allsel' z% z* k. J; Z2 h" F1 X; r) R
/ i& I& ^. H) R- v /SOLU
( K- Z: `) t) i. Y7 _- I/ ^6 ?/ N5 H- M; @1 b" f* c, v" K7 F
f,node(0,0,210),mz,10e5, _6 F# `9 `' H6 _( I
0 i3 b- V; z3 u e
FINISH
8 t% o$ I, n1 q. _6 ~6 i! b2 k3 i* @( k+ H
…...
+ f4 m8 }+ @5 e, G$ c4 D, }) Q |
评分
-
查看全部评分
|