|
发表于 2013-7-25 15:09:12
|
显示全部楼层
我邮箱1962940284@qq.com
% A5 S) ?, ~# Z1 O; Z2 J" K. G施加弯矩扭矩的方法其实有很多种方法,在这里介绍其中的5种,并进行比较:- S' C1 S: N1 v5 W% N* t8 e
( v# V& @7 b5 q7 d% j 1.将矩转换成一对的力偶,直接施加在对应的节点上面。 L2 ?, V! B% u, U8 O% {
5 M2 P7 B* ~) M* |
2.在构件中心部位建立一个节点,定义为mass21单元,然后跟其他受力节点耦合,形成刚性区域,就是用cerig命令。然后直接加转矩到主节点,即中心节点上面。
+ j9 h j, q6 V5 X8 Q
4 c5 l3 Q9 v) P8 E* o' t" M6 E# h 3.使用mpc184单元。是在构件中心部位建立一个节点,跟其他受力节点分别形成多根刚性梁,从而形成刚性面。最后也是直接加载荷到中心节点上面,通过刚性梁来传递载荷。
8 {7 c2 S H: d0 A' ]4 D! j+ h; O* O& g- C, q |( w
4.通过rbe3命令。该方法与方法2很接近。 p! o9 {% s- @. k, S
# U9 W0 a9 N% N5 h! Y7 b
5.基于表面边界法:主要通过定义一个接触表面和一个目标节点接触来实现,弯矩荷载可以通过在目标节点上用“F”命令施加。) }. Q" b* X0 h
) n' L1 c' ?" Q C: O/ A2 I
对于方法1,通过转换为集中力或均布力,比如施加扭矩,把端面节点改成柱坐标,然后等效为施加环向的节点力;而施加弯矩,可以将力矩转化为端面的剪切均布力;但这种方法比较容易出现应力集中现象;
3 ]2 e2 D) B6 C9 O2 m" Q) ^3 P
# d7 f* S9 H+ Y2 |. a: H9 B 方法2,定义局部刚性区域,施加过程venture讲的很详细,这里就不在赘述。根据他的例子,我在下面给出了一段命令流。该方法有个不足,它在端面额外的增加了一定的刚度,只能适用于小变形分析。( [% ^ }# z1 ]* z
, r+ j9 P' {) E: o6 u7 N% t' a$ h# T 方法3 ,相对方法2来说,采用刚性梁单元,适用范围更广一些,对于大应变分析也能很好的适用。但在小应变分析下,方法2和方法3没有什么区别。
! K. V" _! `+ G: Z. X7 p! i: V/ j% Z9 Z& l' e/ K4 b
方法4,定义一个主节点,施加了分布力面,应该说跟实际比较接近一点,但端面的结果好像不是很理想,结果有点偏大,在远离端面处的位置跟实际很符合。 r! p8 E2 ]' b2 D; C
% [: V7 r) z' w: o3 B, w/ P' W
方法5,它具体的受力形式有如下两种:
: Z- K. m2 r# |$ P* s! @4 H: N M7 {' k# H5 q
刚性表面边界(Rigid surface constraint)-认为接触面是刚性的,没有变形,和通过节点耦合命令CERIG比较相似;! f- s- A, k% d+ P; y8 L
+ I# {% P0 y7 @$ W
分布力边界(Force-distributed constraint)-允许接触面的变形,和边界定义命令RBE3相似。& } Q. A; M7 R: o
: ~+ o6 Q( |1 [' k& @1 f% m( P 使用这种方法,需要用KEYOPT(2) = 2打开接触单元的MPC(多点接触边界)算法,. Z+ K. R0 I1 I. f/ j3 e# T
0 _ l! r2 G z4 G
下面针对venture给出的例题,用不同的方法来实现的命令流。
( W8 R# t8 {& g* D" `
3 L1 [; K! k/ M" ^) O, d7 Q 方法1不介绍了,方法2:* t5 l/ o7 C0 U4 p" y
( g2 i# o9 I% V9 M" D5 ?* ? /PREP7 s$ ?. \' N# M, u/ ^( i, z/ k
% D0 z( b7 `4 T* c
ET,1,95
5 p; d1 a$ B* d0 n5 i2 F7 f
8 u/ ]0 N9 T/ h2 k ET,2,21
$ \9 C* \+ m1 s+ D: S" J3 K9 R7 \4 m! E3 q
KEYOPT,2,3,0
- e: ^: p0 ]* O7 R' X0 N( A: d/ F `" c( A7 V3 {- x
R,1,1E-6& u( Z6 |- f" _& k
5 l( a: T- _) D! ]" O% D MP,EX,1,2.01e5, f$ ^ ^9 F0 R V: w3 C9 h
. G# F1 x, }8 S MP,PRXY,1,0.3
! X- S& L4 j$ y3 _. A! N
* W; {" W+ p/ f7 I$ { CYLIND,15,10,0,200,0,360,
( q0 W) D' V0 x6 J; M! c h* W F9 U0 w. t5 S
wpro,,90,/ J. b- i% r" H+ ?7 F; ?
C! p. L" W& K+ V& j4 T1 U
vsbw,all
' w8 G2 p5 J, K/ H7 F$ p) X2 Y
6 {9 r) [) C9 E wpro,,,903 X& A! B& m( F& c+ n- S) G
& c; L6 W# R( F$ ^$ w
vsbw,all5 i" M, q, h0 r% s
) J! N' u+ ^6 b: |/ g WPCSYS,-1,0& d8 ~$ L3 a2 ^7 L7 y2 X/ i
+ F8 b, X. k* k( j K,17 , , ,210
$ s7 ?7 H! R9 E) u$ Y4 b
1 y- E7 x6 g( L; Y2 X$ } lsel,s,,,13,16,1: i/ X6 v. H0 K- E$ f# X& [8 _9 l
6 ^. f- D" i' d, m$ z6 [* D
lesize,all, , ,8, , , , ,15 r6 r8 P8 S) P
9 H$ B9 D# g) a$ D$ B: `( w
lsel,s,,,222 D3 g1 O7 ?( \
9 f9 T2 P- U2 C6 t5 g% ^; _ lesize,all, , ,4, , , , ,1 ~! n- {7 M' s! j7 q0 |
. z6 k0 ]% _+ _# ~# f5 U, ` lsel,s, , ,17,20,16 i6 V6 D2 d' Q# r) v* G7 x
. E' U/ G5 A) N
lsel,a, , ,26,27,1
; V' l0 x# Z5 h& u: ]
h! \1 Z) V9 g2 w, I lsel,a, , ,30,31,1
4 e ?6 f/ H' z7 \, l2 b2 P1 |
; f! ]/ y/ A) x8 \ lesize,all, , ,20,0.4, , , ,1
4 M( s" K! p$ T% F, I. X$ b% d" S" S& W
alls
/ r s/ f+ B' O6 K
; {0 l9 E2 [7 R' W: w+ N% c( i vmesh,all7 h5 c4 e+ R. @- O/ Z5 x& {" A
8 ^' [; ^9 A+ d
!!!!!下面一段开始各个方法有所不同,由于前面的建模一样,后面的例子就不再给出: G" s5 t+ J- V- M9 }+ t
, _, K5 y5 y& T+ @6 b
ksel,s,,,17% Q- R1 U2 z1 b! m4 f
& _- a: k0 B( D, P6 T) I8 c- v
type,2% u6 c/ t: ^% f% J$ Q
Z. h$ A2 P/ w) }7 p
real,1) P9 G& g- t }! B1 E" @. d: z6 z
5 E! ]1 Z1 L$ p6 J& H9 l; M8 ~; G9 f4 E
kmesh,all
& d8 A+ G* n# ]4 H3 Y8 ~
+ t$ j% E. p x! j/ I) f& D& I. { allsel& C1 V# \ ], ]
. E* a3 d8 U& t* v8 O8 U
nsel,s,loc,z,200,210
9 O, p3 I" t# J3 e+ f1 o: k# c! K* |, V
9 n6 h4 {* G" b( h. j4 d npolt
& {7 t% L- l" ^; \, d8 e$ O
& z( Q8 y7 J$ T: E CERIG,node(0,0,210),ALL,ALL, , , ,% D- l# A& @9 p }
& Y: {! A6 M& R6 U2 B @5 S k !!!!!CERIG命令定义局部刚性区域 Y5 e3 x1 b* l' y- C& y1 a! A
3 P4 c: Y6 b; p Q; n& n allsel& {2 u7 y# x4 V; D9 z g# V
8 c; q! p% B6 Y2 F
/SOLU
$ K. W) ]. w |. U8 I
* E: U' S( p* }5 a' H f,node(0,0,210),mz,10e5) L- ~% D2 U; }5 A, \9 e
5 j* T9 \/ E Y1 C+ I FINISH) m9 I, Y( m% G9 j& G4 F4 h
0 \% | _2 y* J0 b5 r* E1 w !!!!!以下一段边界条件的施加各种方法一样,后面例子也不再赘述; ?8 N- V/ B' c7 J4 s* @) j9 m
8 s; p: C$ G2 e( a# R
/SOL
/ e3 U8 d) Q9 U& e6 f' ~$ H
5 n- A5 a- J4 c+ I6 W nsel,s,loc,z,0! h% _6 `3 }9 i
8 o! s1 M6 `# w9 Z( M1 Z, o4 Z
d,all,all
- o& W0 c! C6 V! ?7 L( q# q: p5 l( B# a; y! C
allsel4 X% [1 x$ _6 N
+ {5 u8 X7 F# E3 J solve. l/ ]' `5 n; x$ ~* Q( |
4 K3 Q, h% v! u8 L
方法3:使用MPC184单元定义刚性梁
. `: Z' R: C# ]2 s- r) y$ R+ V/ s- p% ~
……
3 a0 _/ o# n" t. J' z$ B# ] E% n
$ r3 ?/ }6 A% _ et,2,184( G8 I9 _9 {, f9 ?7 u8 S
) U- Z4 b, L9 f/ P( k- R* t
keyopt,2,1,1
3 C" |7 K- \" ]. h- [5 ~- r' I+ U( d* u' y# Z
nsel,s,loc,z,200
, g$ P X6 f; B! K# g7 g- s! t0 \$ y% Y+ i# V
n,15000 ,0,0,210
* G S+ h7 {8 _! R$ W- f
$ A U+ M. z' j- K' Y) ^ type,2
$ o7 j/ M$ g; e5 e% B
' h* h6 t: x/ I3 @ *get,nnum,node,0,count
0 [- N l+ x) @# a* f- ~" f" s8 ]3 d* G- c
*get,ND,node,0,num,min; K Y2 \4 i9 d( R" S7 {( t
, n Z( G) F6 c5 w
*do,i,2,nnum; R" c& J: j0 W5 {4 A1 l, l
( y4 [% g, ~( O" f& V$ H! R
!!!!节点个数是nnum,只需要生成nnum个mpc单元1 _- M- o5 O7 S5 `( d z
3 |4 E$ l3 g! T9 D& ^9 |5 T E, 15000,ND
i5 T. d% F( Y' I( i8 v: e$ r
5 B4 a0 o; Q% \5 n ND=NDNEXT(ND)
' {* X" [9 M) I( W3 l; x; J Y! U% s
*enddo
) ?" q# E8 H8 h# m: _$ A/ X; j! z5 r1 G) t: I1 t4 P% y7 T
allsel
6 d/ H2 y' h' M) R, i( ]
& a0 j/ V$ `+ y7 h8 F" W. A/ u% V /SOLU$ q' D3 B. g% `, w/ ~0 q8 z
6 A. b k% _* {% y f,node(0,0,210),mz,10e56 e) q6 a1 i# A/ T
2 r+ X5 \( y; |+ N
FINISH
- J% Y! S. O; e# |" b& _! t
+ A4 b8 M$ C& c% V9 N ……
- O4 s# Y/ k4 X$ f U H8 q- p/ d& Z& s- ^
方法4:rbe3命令0 i8 m: H% ?; ^6 y+ I: X9 `
% l$ f" r8 X& }$ m ……
2 [+ b9 \3 {8 f' B1 z3 |. Z& x7 N" s. A. k5 h
ET,2,21* z' m9 v& f0 D- t, i
/ Q& x0 T& h% C0 Y" g( l& i
KEYOPT,2,3,0
0 g8 ]$ G3 L( o) v, |. w
8 c3 s* w- h& [" e8 _/ \9 w R,1,1E-63 U- @9 M/ N6 T$ W; k
% l- m5 C, K+ A8 i K,17 , , ,210
$ Z* l6 |+ z: {% j
0 Q" F' S! X0 Q o' |$ }3 X# ^. | ksel,s,,,170 q6 G* P7 c9 C1 v1 S
' A) N& \2 W! \1 N P type,2
- J6 g1 T! C J; H0 N+ w9 E' R( b5 r( w- c" B% a+ I! l
real,10 T' B2 s/ I5 r1 r
' c" O& Y7 x) u
kmesh,all/ ]) r( L+ J! G# j2 L
/ c2 A; C# k8 ^- s0 Q allsel3 v, F5 v: W& V+ k+ J
& P$ f' v7 S4 j3 e+ ]' d
nsel,s,loc,z,2008 F: t7 c3 j$ ^0 M- o) g; f) C- N
5 L' k7 [5 L! C1 G4 r: o# n; H4 H; q
*get,nnum,node,0,count
+ @/ o) @/ h0 `4 O: \9 x! S2 K G
% U* v) ?& F/ ]0 X) @" A *get,ND,node,0,num,min, s$ o- y" G; f6 v( r5 U, x
( i. r; X m. T- ?6 Y8 Y7 g+ k *dim,sla,array,nnum
+ |& I/ T! K; }& D3 u( X4 E5 x9 c" b0 L7 l
*dim,sla2,array,nnum$ G# P' C4 D, d0 _
/ [4 B9 u3 O; s0 j
*do,i,1,nnum m, d0 u7 N) l9 f7 d3 ?7 G
5 j6 ^; N6 G; ^5 r/ I) h' l! F9 A& \% g
sla(i)=ND4 Z' f" H3 p5 f D
8 q4 J4 }- K. U& | sla2(i)=ND
) q, m5 N- h1 M- E7 j r3 M/ ]) O/ z0 o8 ?& y; S
ND=NDNEXT(ND)
, v& | h' L: U! l' F- o' k' ^8 G$ P( k
*enddo
5 \$ @. y+ j" |0 P: {$ ]
- ~* V9 w; Y$ Z p6 e( v allsel
6 z6 g7 V) B; m3 B! _6 y: U. w- @; ?0 X; U- N* {, l
rbe3,node(0,0,210),all, sla,sla2* w! o$ A8 a2 f1 b/ T
* w4 ~& n. }2 \6 B# _# o. J
/SOLU
6 C1 b2 K. @ a* K9 n+ _. i6 M2 s2 ^
% P% l$ m( N& `( \ f,node(0,0,210),mz,10e5! c8 }% N. k9 ]& \: | _
; a9 y) I, D, L2 ?; w FINISH
/ C' Q( z6 d) d/ [. s0 e" x6 `# X
% d: p: c3 G! l% \; k7 ^; r ……6 Y7 J# X0 H& T J% b
, L+ T) r1 P+ [* Z- X5 q8 N 方法5:定义刚性接触面
( H! W7 |2 B. W6 N3 g0 w$ J1 R, t' [5 R! Y. o
……6 W$ e6 N, W0 g; f% T# H' e
: V1 P _: {9 V9 ~8 m1 a n,15000 ,0,0,200, v! X [# [% Q4 E
* R+ O: m. ]+ C% C
MAT,1# B( g' W$ i: Q1 W2 ~, [( R; O8 G
1 ~' f% G m- Y8 i: a
R,3; V1 k' ?1 i; @* n6 y
( l8 j9 s. ~! u- K1 b+ j0 m" |
REAL,3. `( ~6 Y1 L. C" i
5 f$ A9 h0 }) S+ Z ET,2,170% [8 p+ i( s+ ]8 P# {
5 h, l$ x$ z/ m$ T+ B$ O2 G: A ET,3,175' c# o9 T1 w3 r9 ]& {: ]
$ Y8 g- s( m" w E8 R& \3 J
KEYOPT,3,12,5
$ \9 B1 L+ T( t6 I9 ` [7 u$ t
0 [8 K; W2 ?. q* X+ H+ |, C, w KEYOPT,3,4,1
7 U, u0 ?6 n4 ?, R% E( j$ J& n4 r2 l/ s$ j
KEYOPT,3,2,23 B2 N1 D" r( E, |! I
) E+ G; ?& k8 j' f2 C$ A
KEYOPT,2,2,0& R2 U' y7 r+ f! ^
- x+ ^& B0 C; s, k" n4 w# O KEYOPT,2,4,111111, l9 |' D* t `" c! `; v1 m/ H
" Z5 O! F# v1 G( @2 O" ?4 M TYPE,2
! u5 b$ `3 U6 o# d2 M9 M; I7 o6 m- w! `" D. L
! Create a pilot node
, u& v R T. n" ?. y7 @
6 B$ B2 b5 m" y; J* U9 P1 i6 _ TSHAP,PILO
o4 T7 y3 r6 i. i! C) c3 q4 ?5 d
" x; I0 W& s9 W( F$ x2 d E,150004 _/ D1 c7 [, b+ Z D* V) f9 _
' c( b* Y' [4 _1 e. ^7 U8 y* i ! Generate the contact surface
- U1 @- Y* d y) _* |4 \9 B; P, W/ @7 E5 u1 T; |! j
ASEL,S,,,14. h, g/ `( N4 D$ L* S
) J. i" x/ V, {, d8 _* t. T Z ASEL,A,,,19
, j8 ^# f: k7 `) b. j; e( G0 R
, k# E) a( J$ `) I. L ASEL,A,,,24: n- Q# |* |8 r$ T* g# Y
2 Y% s* x2 H$ r. [. i. w0 V: Q3 P ASEL,A,,,284 i( G0 @' ` x$ R q8 ]- w
0 R2 ^8 I3 |7 t) j4 a
CM,_CONTACT,AREA4 |! g; B' Z: t1 q! ?5 ]% `6 y# m. F1 o
% W, m& L, b. ~" _& p1 g, \
TYPE,3
/ m$ d; B/ \' ? h6 c8 P$ K" k$ X" W+ d' d
NSLA,S,1; }! C* F1 x& @$ G1 Y I" B
8 [9 I, J+ t7 N/ M
ESLN,S,0
! H. O: ?) b! h' q0 a
# A* k0 }4 p7 D3 n% ~1 t" f ESURF/ K: {+ z1 @6 W* X
% Y& X2 p0 V! c: V# T ALLSEL
m5 W+ W9 ^* j( p2 F4 t5 }
u0 G) r3 v# E7 V4 ~0 y" F9 y$ G allsel
) v$ b3 `- r$ `2 A$ r8 _5 J% p+ U
/SOLU
+ W& _% V" I' V2 k3 ?3 m$ g0 [" }% k! b9 ^$ W- \
f,node(0,0,210),mz,10e5$ M& H: v6 n* b$ F3 d2 I
- g) v( L9 [! u" r$ a$ h FINISH
8 C3 R1 n; _0 B. @' c5 v7 F& l. J3 ~' L( F" x; P+ S$ d
…...
r7 h+ o5 c1 o3 A |
评分
-
查看全部评分
|