我郵箱[email protected]" E8 @9 z: u3 _& c
施加彎矩扭矩的方法其實有很多種方法,在這里介紹其中的5種,,并進行比較:. M( A( m( J; N
6 G! U j0 \) p3 [9 T
1.將矩轉換成一對的力偶,,直接施加在對應的節(jié)點上面,。: d5 @, |: a/ S4 T" \& Q
& l( K! S4 V9 L+ l) P5 A) d( Y9 }. u 2.在構件中心部位建立一個節(jié)點,定義為mass21單元,,然后跟其他受力節(jié)點耦合,,形成剛性區(qū)域,就是用cerig命令,。然后直接加轉矩到主節(jié)點,,即中心節(jié)點上面。7 Y4 x+ r% P- | Y
) }0 b" n" @/ \7 ?, S
3.使用mpc184單元,。是在構件中心部位建立一個節(jié)點,跟其他受力節(jié)點分別形成多根剛性梁,,從而形成剛性面。最后也是直接加載荷到中心節(jié)點上面,,通過剛性梁來傳遞載荷,。
1 A, h! @( Q8 E4 g6 D" l2 L
' C# ~2 d9 e; F8 ]! C/ r( N! X+ I 4.通過rbe3命令。該方法與方法2很接近,。
# Y# [$ @. W; ~ x/ k" H$ t% n
0 L! G" x- i' i8 t: d; D: B 5.基于表面邊界法:主要通過定義一個接觸表面和一個目標節(jié)點接觸來實現(xiàn),,彎矩荷載可以通過在目標節(jié)點上用“F”命令施加。
G' E2 h; |4 p7 A) ~, o: F* [& [* c: O7 e7 I$ p
對于方法1,,通過轉換為集中力或均布力,,比如施加扭矩,把端面節(jié)點改成柱坐標,,然后等效為施加環(huán)向的節(jié)點力;而施加彎矩,,可以將力矩轉化為端面的剪切均布力;但這種方法比較容易出現(xiàn)應力集中現(xiàn)象;, w% f2 x7 b" w! V3 D2 E0 |; P" c) t
4 R/ g" j& I1 w$ s' T! ?, E# Y
方法2,定義局部剛性區(qū)域,,施加過程venture講的很詳細,,這里就不在贅述。根據(jù)他的例子,,我在下面給出了一段命令流,。該方法有個不足,它在端面額外的增加了一定的剛度,,只能適用于小變形分析,。
& w( M9 Z$ S2 `7 O8 A" X2 b
0 ]; W6 T& W4 O* Y 方法3 ,相對方法2來說,,采用剛性梁單元,,適用范圍更廣一些,對于大應變分析也能很好的適用,。但在小應變分析下,,方法2和方法3沒有什么區(qū)別。
* d9 O( d, a$ i
$ ^+ {4 h) A. W9 F5 E 方法4,,定義一個主節(jié)點,,施加了分布力面,,應該說跟實際比較接近一點,,但端面的結果好像不是很理想,,結果有點偏大,,在遠離端面處的位置跟實際很符合,。
2 _6 O" F* Z8 ?, L9 e- g5 ]+ u6 N$ e( [4 R
方法5,它具體的受力形式有如下兩種:
9 I. [! \1 ^6 Y# V/ M- ^# @4 f9 _" w* o% t0 y! p p" w( j! V
剛性表面邊界(Rigid surface constraint)-認為接觸面是剛性的,,沒有變形,,和通過節(jié)點耦合命令CERIG比較相似;/ q' ]9 l$ f8 _6 X& X4 k+ f* H9 I
; h6 g3 G4 s/ j8 W1 S 分布力邊界(Force-distributed constraint)-允許接觸面的變形,,和邊界定義命令RBE3相似,。
1 U8 R6 x: |; A: D7 e4 _) | }! B g8 ~3 t' N" v
使用這種方法,需要用KEYOPT(2) = 2打開接觸單元的MPC(多點接觸邊界)算法,,
/ q2 Z! K+ s% X. y& }
/ b! |2 N. t/ ~+ [: ^9 b 下面針對venture給出的例題,,用不同的方法來實現(xiàn)的命令流。
1 u1 X E" U, y- F& ?. h, H! T5 X) o& b- \: K+ x
方法1不介紹了,,方法2:& r, K# U4 J6 t9 ^5 D
[( O, Q( @) g, [0 k /PREP77 [+ A {& M* H/ ^; \4 s0 e/ `( u
) i3 T% l! C# I0 ^/ M% Z ET,1,95
6 @ F- b) }+ n9 G" L0 I. Y9 j1 Z$ E; e( @5 K) a5 d! |% q
ET,2,21$ ?, q4 k3 |1 o" H& H0 k. F
" R% a0 l U/ l0 a4 X
KEYOPT,2,3,0
3 A; t' F$ t) J1 {2 l
% [% A/ f P6 |8 R R,1,1E-6
6 w( a3 W* J# Z' F3 x
6 j1 b; Z3 Z) q4 P$ W9 M) P MP,EX,1,2.01e5
; W, O+ C% i9 D" f1 \) e
O) Z) f; @9 E MP,PRXY,1,0.3' ^& j$ D; X& W0 f5 ?" [. Y0 `
G6 \1 e h, q) G CYLIND,15,10,0,200,0,360,, r- n7 y/ Q- ^5 v1 H- K
- d2 J0 N/ f( z! q wpro,,90,0 w/ m( P' D- s! w$ \
; k! Q& ~* E) |- j8 e- R+ {% m+ w5 I
vsbw,all
3 b; o0 ?5 ~) C) Z- B& A9 z. S o# p h% g8 M6 h. q1 d
wpro,,,90
# S5 z, }/ S9 k& \. C/ m, c2 i( w z+ y
vsbw,all% J" E& q. ?3 Z Z
' c5 F5 l+ g, t" e' Q+ E1 y
WPCSYS,-1,0
9 I7 I0 U7 f1 s0 V/ ~6 W
! d- x+ U# [6 D' O9 ~( X% } K,17 , , ,210
0 y7 K- _3 ?# h/ ?' H/ s
% V2 [7 u! F" p* Z, U! E5 k5 C lsel,s,,,13,16,17 f6 Y- c% z% b; _' j7 J, O
* \% H1 |- p9 ]2 B# g lesize,all, , ,8, , , , ,1$ z! g: u0 X; ?% d3 v
+ P7 \! K# v# I. ` lsel,s,,,22
7 @/ j# g% Z% W+ y% l" o& N b
0 f( E( f+ q+ E% y9 L lesize,all, , ,4, , , , ,1
4 j; i9 v; |7 V0 X) M4 ]. i1 q6 Q) n( _/ ~$ Q% L0 P: c
lsel,s, , ,17,20,19 ]; X: _. Z6 K2 j* b
7 r, n& J( N6 u4 @" K
lsel,a, , ,26,27,1
0 k% ~: O. u0 r6 ?% c, [; D+ Z9 R8 r: R j& Z# V1 e" e+ q0 t
lsel,a, , ,30,31,1 U% E H, J, Z7 v0 W
8 O5 ^4 m7 e! Q4 Q/ Y' ~
lesize,all, , ,20,0.4, , , ,16 b2 e- l" W* c3 k$ q
: V6 v0 Q- O: X. B$ D0 r7 I
alls, A$ M9 l* q O2 P# [ J
8 s. z# n( d& h$ J' \! h vmesh,all1 }# m4 I0 s% w: |& X2 P
; b6 ` Y7 H4 r' X r( C5 W* n
!!!!!下面一段開始各個方法有所不同,,由于前面的建模一樣,后面的例子就不再給出
- ~& v) ~4 f8 Y: d$ p- i/ r+ [$ `% M
ksel,s,,,17* P9 v: z2 E+ ?4 {2 }. S2 S" B8 ^
# @! \4 e: c D& I+ n' t
type,2! r) @& h2 E* i0 Y; E- t+ k+ c5 n* [$ z
3 Q2 P/ v# d0 c( g
real,1
7 e/ H7 n( w! [5 \. F5 w( f4 e1 X; g. K( L, T8 M+ j# l# w
kmesh,all$ n: V7 x3 m! W& n- C: y
% _5 {8 e: h+ p4 a8 W
allsel/ R9 X& L3 m' W& M5 A6 }4 D+ f
# M1 {6 g6 ?8 a
nsel,s,loc,z,200,2101 ^3 a, b" z [4 s
) U/ M5 \/ s7 N/ z6 D npolt1 K+ Q1 g0 G% ]( \" {
1 o' A6 `4 z! T2 T8 w- l CERIG,node(0,0,210),ALL,ALL, , , ,
e d* q# x# |8 O" S; k+ Z- p! o# k) k) y& U& q
!!!!!CERIG命令定義局部剛性區(qū)域7 c# r+ y( S4 I2 l
) G$ P, G2 H: `2 ?6 Z- m allsel
2 O" @' A5 ]' v; ?- c" h
. H% a" J) a5 i) O( |) S /SOLU
, e E4 q6 T( }' Q v0 v
0 S# T! _7 o( U* @' @: X f,node(0,0,210),mz,10e5$ \" J" y4 X: z3 d5 C
0 u' @% T4 w( {' \ d FINISH
3 a5 N; t5 K1 `: {6 V& D1 p
1 x8 B& `* U# g% y" V4 ]- i9 M !!!!!以下一段邊界條件的施加各種方法一樣,,后面例子也不再贅述
7 ^$ z& O% n8 d% |8 R3 y
- C" J, W* q# a+ V: Q; s7 e& B /SOL) h ]3 S1 }, H5 Q2 B7 E+ [8 D
3 ]0 b8 d8 \" d! c+ b nsel,s,loc,z,0
& P! n& l C6 S2 N' J/ v ~$ \
3 l* [# b: z: Z b, U d,all,all+ A7 h3 N$ Z& |' }1 G O
- i/ |* V/ Z4 Y+ c7 u' A% }. V+ d3 a allsel
- h+ M: H( v! |
8 e. W$ z1 z, m& K% M2 V solve5 {- Z4 \/ n8 i$ ~* y5 q
1 |; P) d6 n& ?. A! W 方法3:使用MPC184單元定義剛性梁3 G( H- F! F8 w
# o$ G# Q# c9 T" L3 C1 g ……6 l) i2 u) |. O: z; c6 i* w7 J
. {% v! ]8 c9 U et,2,184$ @% a$ K3 d/ O& u
4 D# p2 w' i7 q
keyopt,2,1,1
, S) R& y. ~- U- R6 Y: k9 X9 D3 X: \; U" l" q5 e& \
nsel,s,loc,z,200
* w% l/ L U8 w5 m2 k
+ ?# {, _0 h8 [* s. O: L n,15000 ,0,0,210( q$ n/ b8 J0 o5 z- H/ ^
- ]3 \1 a9 ]6 `: Z) i. E- i- | type,2
: I7 _( J" K$ e3 V+ q# I- h, r, ^* k; K* x. {' v* d% ^" t! X
*get,nnum,node,0,count X& l3 ~7 B6 g f* O
' e. P9 w3 \# k6 H6 \! ]
*get,ND,node,0,num,min0 _) q8 _% t# R! I0 Q/ ~
2 \5 m4 Q6 x$ U5 w2 \7 |( L! A
*do,i,2,nnum( j, `: Y* E" I7 n4 X1 v0 B+ v9 Q
6 i. |) g* O8 v6 F( l0 P9 m" ]% S
!!!!節(jié)點個數(shù)是nnum,,只需要生成nnum個mpc單元
- m# x( u; }' C* Y& v/ E* u6 S
3 R& c7 F! M2 c* Z/ A3 A, y E, 15000,ND8 Z! e, D3 U3 `9 I
7 O* D& D: D5 c: }9 y8 f" ~ ND=NDNEXT(ND)# u4 z0 j$ y5 I( E6 H9 e2 G7 T8 X
+ q4 C. `0 m/ z. N8 Y *enddo
+ @2 U) S2 Z! D7 @
1 d/ i8 r# U- s6 q allsel
5 w5 u; Y: V7 R( M2 T( R
' {7 ?0 f9 }4 B /SOLU0 r$ ~/ r! |& C1 }
, N( V; N; @. z
f,node(0,0,210),mz,10e5
# s- p( \; s. w0 v1 m' U2 U% K/ {
FINISH: v$ \4 f! ~/ J% B8 ^% f- p
6 u1 L) i( b$ u2 k
……3 Y. f, b5 A$ y7 z& Z5 v% H1 M4 o
7 |7 p! R0 ^! B4 L7 J0 W$ g8 k& K 方法4:rbe3命令$ B4 s! r; a) K- f0 [
( [' l' s$ \4 b; c2 w ……9 _# [7 Q9 ^' o$ m+ M& g. U
! N; v& M s) v: U5 \* I0 n
ET,2,21
' P6 Y3 e) ^8 A6 n+ I, h) M' V
5 `3 m7 I4 f; N) B7 w% |7 h KEYOPT,2,3,0
" _$ u( y j% I2 d% D! Z) E0 @0 {, S. C( d/ `
R,1,1E-69 d' i& A, A4 K/ |. H
# o5 u0 X5 M- I" h( j) L% ]6 j
K,17 , , ,210; s$ P3 e, a; \
& t0 W+ ~' D8 o$ `" P
ksel,s,,,17
, g8 t* v9 b0 v1 |' Z
% u9 C$ y3 C* q7 a$ A$ _/ h$ i type,2
; ?( Z- O% D, H! m) u) ?& x7 e! B+ ^$ K
real,1
2 B) Y' {7 u! U% E( Z- L2 K2 t
, ]4 \* }- k' r" P7 ~ kmesh,all
$ h2 g2 H5 u8 p9 ?) C4 ^
- F+ _$ `# k5 f; c# Z& A allsel, I& E! S" ?" V; |8 ]
" J( l ^5 y" j' V; [ z6 ~- A% T nsel,s,loc,z,2008 b2 M2 [4 ^- T& I
$ e0 V/ \6 Z7 U* u8 K" |
*get,nnum,node,0,count
: s o+ }$ [# _4 z& w) v9 N3 p) ?( A/ a3 W
*get,ND,node,0,num,min
: f) V$ B/ {& A. q
. K( ^4 P6 ^4 q8 o; w3 G5 C *dim,sla,array,nnum
) s) f$ j7 |$ ]+ _( O/ ]; }3 t) J2 _% _% N! o- z
*dim,sla2,array,nnum
( j% n- p( z9 P" x1 l# i; `5 Y6 T6 ]- y F7 e
*do,i,1,nnum. ~! V" `: f5 ]# f
3 \* Z& f; q A: Z7 T- A sla(i)=ND i2 d m; Z0 R Q. ~4 |2 K! Z
% m$ F& l8 g8 `2 t0 q sla2(i)=ND
+ x8 T: t4 s! w9 y8 v
: X* x. ^5 g+ C0 {1 s# B( n ND=NDNEXT(ND)" Q) h. b) C4 Y. F
& y: M. W$ ~ d& m! K; n& o& [ *enddo, H. i9 l2 m: }5 h9 `- Y
& {4 D! y( W+ G6 y2 M+ x5 Q
allsel
9 ?- a- `! k5 o% H r* [! ]
/ [; j) o" ^3 S; w0 z rbe3,node(0,0,210),all, sla,sla2
" W! @% @! a# q R5 L- l2 x, A! M8 `4 m
/SOLU' ?7 S9 ?3 |* ?& u! {1 K" I
* O- o- H* _+ d0 g: O- K, _ f,node(0,0,210),mz,10e5
; v q% |( d* J0 D. G; ?; n
" {5 d" ~" }/ [: R8 t2 l' G FINISH
4 x# q1 F# c0 V- G
3 Z7 x, Y, m4 b* X ……! e! T& N C9 o4 Q- Q2 @# J# ~
0 f& R* k' u% c1 K9 T: m8 _ n
方法5:定義剛性接觸面, v* c% u1 i3 p3 ^3 E
: m/ m" T$ U! ?1 u$ n7 \0 G( L
……
8 a; M6 ~9 U) B( F/ \" q: i* ~) u9 h& B+ u a9 V
n,15000 ,0,0,200
/ ~" W2 Z9 o# X+ I% \! _& O
) c/ Q; ]7 l5 c+ u' C6 l( B MAT,1
+ s0 q5 w+ t' N" n4 A, m. ?; v/ g* G/ v" [: ~4 T8 R8 Q
R,3
8 g, F, w4 K. U- Q/ @4 @( Y& t+ L4 Z! H( w Z
REAL,3
: z, Z: C8 t* K% P$ s7 v1 k3 V6 X
, a9 c) X1 c3 S' K, a } ET,2,170$ V& s3 r* z8 w `0 ?& _7 d8 }
4 y R# I' ^0 A8 U( @9 M
ET,3,175" }3 X) K1 ~) J- t P1 u
\- Q/ o( s$ v6 i4 L3 Q KEYOPT,3,12,5
6 j j7 O, I8 A! w
& c0 |9 L1 G; `8 u \ ? KEYOPT,3,4,10 X, P% o* S& m. q2 Q, X* U
) [; V; V0 I/ f J* C; o4 l" [8 n KEYOPT,3,2,2
/ k6 u% N, p! Z" I, j
& }3 \' C! \+ z) `6 y KEYOPT,2,2,0
1 e! t8 ]6 Z5 q6 W; q( e' a' T" O( S E5 k4 A! G* o5 R
KEYOPT,2,4,111111
9 ?5 i7 H& v! j& U8 ]5 }" P3 F
1 H0 e& C# G; d9 ]+ V$ C7 V TYPE,2+ O1 ?+ X0 g5 ]+ Y0 U
$ C0 ]7 h8 C4 [, E3 T9 w
! Create a pilot node) S6 U+ }/ ]/ E$ A
5 \, n" C/ i# J1 Y6 ^( g
TSHAP,PILO
7 g0 l+ L* P% w& H
, B" r- T9 ^& _% r- p3 V+ ~0 E6 j E,15000
" l3 p! s" L" W+ L! M. G* I/ T' U) [" ^
! Generate the contact surface, c7 R) F: d4 `/ [( K0 y8 m0 {
, R5 e. n& V6 b) F5 N ASEL,S,,,14( S) H. r9 J N, k
8 e# g* Y; p5 R1 D$ H- x( q6 y ASEL,A,,,19
/ ]& |6 q& _ i- I3 J6 V; i, @4 {( ]7 Q* u2 t' f# v: Z( I
ASEL,A,,,24( ]& `- {$ [& s" `3 [4 h
* V0 s# O/ w$ p( T/ U( y' c ASEL,A,,,28
- ]1 {( @! A% \5 C. x9 j
) R7 S9 g6 p7 S K% _ CM,_CONTACT,AREA: |0 d6 B r$ d5 ]/ w
+ ?( F+ D$ e# y6 r1 v
TYPE,3
5 L% @4 K* \3 I, K% _9 q5 h$ u; v( W8 L) g8 S
NSLA,S,1
6 ]0 P. Y, R; ~! _8 s7 a7 _( f+ a9 G2 Z" n
ESLN,S,0
7 V) h1 k$ y2 I1 s
0 h1 \# Y- F& n ESURF
0 b& d f9 F* f9 Q1 h
; }6 l- A! r' o8 V: A0 F. D0 q+ X ALLSEL S2 t; A, K) G' \3 F
7 b* r7 y3 b9 L3 z
allsel
. o- ] h9 \5 _# T" l$ p8 c4 \9 m7 N+ {* W$ I u; i4 W4 m
/SOLU
- @9 c$ U# @0 ?: t2 i' k3 d
! g7 A$ w7 o/ m# m f,node(0,0,210),mz,10e5
2 O3 J6 S, Q. C, D) s
1 X% K" b8 L+ S1 `: K2 O FINISH
. }) Q2 y' W# q; t4 ?$ p4 Z/ x5 K: v& m+ H
…...
: J! k' z( X0 \( D& r7 L3 j |