|
|
我郵箱1962940284@qq.com' c0 V; K4 A. q5 r$ \6 w
施加彎矩扭矩的方法其實有很多種方法，在這里介紹其中的5種，并進行比較：/ j2 C6 N8 [6 p5 |7 w. Z6 q
1 j% c: r, h4 @5 @0 `0 [" D
1.將矩轉換成一對的力偶，直接施加在對應的節點上面。
% X4 d2 a' w s( S4 H# v' d. Q4 y% r
2.在構件中心部位建立一個節點，定義為mass21單元，然后跟其他受力節點耦合，形成剛性區域，就是用cerig命令。然后直接加轉矩到主節點，即中心節點上面。
2 b! W' E8 s; B6 c* K& P/ l+ g; [$ w6 k: d1 `$ c8 h! ~% U3 f
3.使用mpc184單元。是在構件中心部位建立一個節點,跟其他受力節點分別形成多根剛性梁，從而形成剛性面。最后也是直接加載荷到中心節點上面，通過剛性梁來傳遞載荷。
! k2 Q6 F6 w' N
+ m* g/ u( G% d7 h! X' ? 4.通過rbe3命令。該方法與方法2很接近。0 c# d* T9 O3 a; ~. ^, j8 G. e; o
% @# J* m$ p# v8 d 5.基于表面邊界法：主要通過定義一個接觸表面和一個目標節點接觸來實現，彎矩荷載可以通過在目標節點上用“F”命令施加。& Y0 c5 i/ c2 p5 m4 \$ m( ]6 I
0 q. U$ |" b! j# p# G0 X0 l, F' @
對于方法1，通過轉換為集中力或均布力，比如施加扭矩，把端面節點改成柱坐標，然后等效為施加環向的節點力;而施加彎矩，可以將力矩轉化為端面的剪切均布力;但這種方法比較容易出現應力集中現象;6 ^1 S& I0 T: }
2 |6 o3 Q/ q6 C8 h3 L2 L* M1 w9 Q 方法2，定義局部剛性區域，施加過程venture講的很詳細，這里就不在贅述。根據他的例子，我在下面給出了一段命令流。該方法有個不足，它在端面額外的增加了一定的剛度，只能適用于小變形分析。( h6 X8 M2 O/ `2 i/ i, O# E
2 b" H( p, u V+ }
方法3 ，相對方法2來說，采用剛性梁單元，適用范圍更廣一些，對于大應變分析也能很好的適用。但在小應變分析下，方法2和方法3沒有什么區別。% e# I4 x+ W7 H" N
1 j" k- ~8 Y9 f6 w5 O, C, e9 x$ v
方法4，定義一個主節點，施加了分布力面，應該說跟實際比較接近一點，但端面的結果好像不是很理想，結果有點偏大，在遠離端面處的位置跟實際很符合。
- A$ B3 A: v8 Z0 M% d0 T L( t8 M8 s
方法5，它具體的受力形式有如下兩種：
3 h' B( d5 y$ _$ N) s; a
4 P3 y7 w) G8 ~1 s. _ 剛性表面邊界(Rigid surface constraint)-認為接觸面是剛性的，沒有變形，和通過節點耦合命令CERIG比較相似;1 w# o" g& ?, `) b& B8 ^
+ W! A! d, X5 L
分布力邊界(Force-distributed constraint)-允許接觸面的變形，和邊界定義命令RBE3相似。
' W+ i* M2 v9 n5 ?7 e
5 {' o9 Q" W- h, c' K8 }2 l" F# a 使用這種方法，需要用KEYOPT(2) = 2打開接觸單元的MPC(多點接觸邊界)算法，
4 ~. ], N" ?1 V, n- Y6 ^2 c, d
L7 Q4 ~( b$ Z- b; p7 c! T 下面針對venture給出的例題，用不同的方法來實現的命令流。
7 y7 D# m4 t) H5 @8 \ j. v& d. b& E( [5 A
方法1不介紹了，方法2：
0 o" _" p) M7 u; W. o. V* ]
9 z/ R3 R4 ~" b& E /PREP7
( G; O i1 Y( g: t! D6 A1 d- W, g( `( q3 K4 H7 G
ET,1,95" u4 C/ l; p) [
$ t8 {. \: g% Y4 u ET,2,21' e2 M" n0 F$ x5 W4 _
0 l' z0 Z/ w, E Q, e7 S0 i+ e! g$ x KEYOPT,2,3,0
" o, U) d" R, C0 `9 @) e, ]3 e4 _. o, u7 v9 I, l
R,1,1E-65 L) n F$ S! d* _$ v! l* W
. ?- S- b; x/ _) q) t/ v MP,EX,1,2.01e5
4 L1 }7 R- Q! K7 R/ {# i1 R# X# n3 q. D* J A* M9 ~$ ?
MP,PRXY,1,0.3
4 \ ?5 e8 S- C0 P0 a( [: N+ j8 r! A, v' d2 V X) _1 ]
CYLIND,15,10,0,200,0,360,( s& P x, H9 r" h$ p6 O' P, T
- n7 J; b$ z# ], M2 g! u+ h$ S d wpro,,90,4 d. C/ E$ e: Z" r- Q- c' O+ F
8 ~% c; p1 B6 ?* c I vsbw,all
; L0 i: |6 Z% U8 P0 p1 X
. d- i7 q+ [9 l& T3 I wpro,,,90" Z8 }( j9 H6 X+ e
) U" q5 m. x9 k! {+ N
vsbw,all
/ R! k5 w+ t. P5 j2 C9 M f9 F) m+ g$ Y0 Z; I+ |
WPCSYS,-1,0
. G2 w/ `% S( l- ~7 `% U [
k3 w2 M5 T) n2 ^ K,17 , , ,2107 X; G' ~# k. v1 b( T
, N" k( @& Y8 j' a
lsel,s,,,13,16,1! S1 @4 O$ q4 a* @- k* B
1 N9 g/ a; x# {- d9 S lesize,all, , ,8, , , , ,1
+ G8 n) ~! l4 b. ~# k; c" L5 u
3 C) M" V% w0 X" ~0 v0 d lsel,s,,,225 |: k- r Q. }* x8 g+ ~
/ t7 q4 G2 ]! e$ \/ E. o; s lesize,all, , ,4, , , , ,11 ?/ v# \0 H* [, W
7 R8 u" g) ?! s0 r1 n8 d; |
lsel,s, , ,17,20,1
) q: @& U9 i) y6 T
7 f* d9 M/ N1 I* Q G5 c7 W lsel,a, , ,26,27,1/ y6 a2 h/ k( o' ]2 C& b- z. l1 }
/ n" L2 D' ~( {5 N2 p lsel,a, , ,30,31,1' l2 V: V4 Q7 ~ S
9 Z. }, h- J6 C4 l$ B: K lesize,all, , ,20,0.4, , , ,1
* b* y- O, |% u6 j0 |" I# h, g7 l- i2 \. n& ~5 e
alls9 X9 @4 W! P) W. `6 t
, d; i% S: t3 q' G
vmesh,all3 D' k0 q- ? C8 U) a
$ I5 A2 Z6 `/ `' ~
!!!!!下面一段開始各個方法有所不同，由于前面的建模一樣，后面的例子就不再給出
& W/ o/ a% @$ |2 r3 g1 L
& }+ }- E4 h' c- S! F ksel,s,,,17) @/ v2 t, j$ y! t; T
- |3 o! K E% A; Z* ?' Y4 F# c
type,2- u, n d; M# l$ s; h
0 a, G# ?8 a# _3 H. j0 x. @: F
real,1
. a# i3 V( G2 o& ^, t3 L1 k8 u! ^" T5 [1 |$ B- Z( \& V! Y
kmesh,all4 H, S# E& n% c7 ~
6 Q% `) x) j, d3 |; ]& W/ ~. g+ @* j allsel) n* O* O# Q \; G
5 \, _, E0 U% I1 V/ o
nsel,s,loc,z,200,210
3 m6 O$ y* S# x' E2 Y# a4 ?" {2 y+ Y" p% J& c& i
npolt
& l& E3 b+ \& t* P% {; ]( ? G) b
) {6 Z" ~5 ?( u N CERIG,node(0,0,210),ALL,ALL, , , ,
8 f$ ]* I2 j' h# X z: d3 I0 T+ M# R
7 {# b7 R8 \6 U( Y !!!!!CERIG命令定義局部剛性區域
5 q& ^- V; S4 Y( y2 W! p' s! T9 j+ ^. a/ K/ z! H2 ~/ l- {- y( i
allsel5 P; l. x* g* c
- `; ~, N _! ]! z0 r /SOLU
3 {6 ]0 j/ u8 w
, D0 u8 f' O0 w& ^ f,node(0,0,210),mz,10e5
* U- D4 g) b& d }: e1 s; B# e+ z3 Y' f, d: h7 w' ?
FINISH
6 X5 u6 f& ~5 e0 k2 M: _1 f1 h1 e7 G2 K; g5 x
!!!!!以下一段邊界條件的施加各種方法一樣，后面例子也不再贅述
+ P* h2 X/ r) w! f1 z# k' Z. r8 U
/SOL6 ~6 H" e/ d9 x; Z
9 _7 \0 W2 i; j( a& ] V; ~. i" y
nsel,s,loc,z,0
5 Z$ c( R7 K+ s' z( i0 D) f
, f% w5 H+ L" z y9 Y9 j d,all,all
. r8 V% P# s0 ]% a! q
8 s0 S. ], g: I: B1 y- d+ F allsel
. Y; l' f7 [9 P* R8 L# g5 z
8 i+ f% [4 e3 l" w2 g; I& Q solve
- q, v2 C: e2 p `0 b$ }. s Y; M7 {
方法3：使用MPC184單元定義剛性梁
! u0 i4 V7 o+ D9 x: C8 E: Y( G
. i R2 r3 O/ i$ w ……3 w1 \8 z: {& \
9 S+ R# Z8 N1 {3 T4 R* O* D( R) r et,2,184* m z: L5 X) [: e7 ^4 O5 `! P
- Y- y Y, C) A9 U keyopt,2,1,1
! {. g! k3 j+ O# `' X& i8 Q/ y1 o; h% L. [! c" ~% z( L
nsel,s,loc,z,200* B" A5 Z$ M, h, T, y8 ~
1 f: K0 Q. A: ], J% k1 k n,15000 ,0,0,210
1 N( G* w: a( x4 p5 x5 K
. `1 k" U0 ^$ `9 |0 v, k type,2
/ i; R/ E( ~, P( X& G% N) S. N( ~( Z
*get,nnum,node,0,count! f3 \! _: B) ] `5 T q
8 ]( J6 S0 F' |
*get,ND,node,0,num,min
_! |+ O% o8 O# o% b, Z3 d7 W; q( e( P3 [! d6 i
*do,i,2,nnum2 ?4 }2 H" f5 x5 s" A* w/ b
& W! A' ~1 I5 r: m4 I6 B7 Z8 s !!!!節點個數是nnum，只需要生成nnum個mpc單元: |& t; y# K6 h
% Q2 S: d$ r H) p- N
E, 15000,ND+ R5 N1 `, u" k0 E
Q8 d2 q4 o0 h8 n1 R. _
ND=NDNEXT(ND)
+ s7 Y" \- C% {' V6 {$ F
( g1 A- C" F1 ?5 i *enddo
, v6 ?" h4 D- R" f: T/ Q" x, ?2 S; Q/ ?, [8 P( Q5 w3 a, \
allsel
+ a, i# ]) \5 x% ]
0 j, [& |6 {/ W9 N: s /SOLU3 H6 ~1 U, ?' n) ^3 Q! F
0 ~9 z/ \8 W) {% G* m) P' [ f,node(0,0,210),mz,10e5
- X& A: r; P3 t; Z
6 t' ^5 {7 ]" B- P7 }% X% L" ]: V FINISH) D- }; S" d/ v# R+ _
4 X6 k8 v4 c5 t8 T9 R; \1 E
……' K/ ^& B. p3 l. x4 @
1 d& z& _! S3 k; w3 H 方法4：rbe3命令8 Y5 r# x( E3 @% W) A! p& C
; H$ j3 |$ c) \6 l8 i5 F( V
……
3 C- C" |+ S- J# p A+ Q$ a! j
. {4 |- m2 O0 F {4 ~ ET,2,21
. q1 F* J# b/ E0 R9 u
/ f# u1 t6 f& B8 i. w8 | KEYOPT,2,3,0) B' X# F" P" x
* c+ o# u+ N) Z) Z: s R,1,1E-67 r$ e& K" m1 m' r0 V' I" m
/ ]. S/ u4 E7 W: ?3 R
K,17 , , ,210& \+ q* V- R0 \. P, ^7 H8 W. d3 A
, ]9 k4 t: {1 C
ksel,s,,,17' |6 F5 i. f$ ]
* k( h, p1 ~ u type,2
$ `' Q9 `1 Y E5 I; m- |. x' L2 ?
/ @; q& d, n. V! Q% v real,1
: B9 m7 C' b3 C% s5 ^
5 N6 k9 d1 A/ a& O! @: A+ X kmesh,all
! a- a: M. r8 N% @' r5 |* d0 I# R6 w Y' e
allsel
1 i+ P& ]8 x8 o2 N8 e$ u; V# D. B' m6 m- E0 h
nsel,s,loc,z,200/ I$ R5 N1 \7 d: _% O6 R6 y% U
W2 u* i1 c# F, d+ _; R" d
*get,nnum,node,0,count
0 Z# A7 Y) X6 W/ q# P m0 w. ]' P ]% U( w) o% p* T& u
*get,ND,node,0,num,min: {1 t% ]8 n6 s" J7 R
/ x* Z1 O, B3 Z* _, J% f$ r
*dim,sla,array,nnum
2 D2 d% p- _) p) y& D
# ?, ]9 Y2 y% i: D2 A" k4 A *dim,sla2,array,nnum
. U) d3 f( {/ @ l$ S/ t- U$ X* v+ }1 I+ L& Z5 x' @
*do,i,1,nnum
$ Q3 c4 Q3 ]* f8 {0 E, N7 A
6 k( k A8 { z! X8 K sla(i)=ND
+ N1 X; o! f; ]) k; m7 C
& J/ s* ?/ D0 e! Y1 u& ?2 I sla2(i)=ND
' K9 M# c4 s r5 e$ d6 J" P) X
; {! G0 ]( W% v8 d( \ ND=NDNEXT(ND); f- S# `3 r; u+ N
+ u6 B7 ]$ j5 f9 g& l *enddo* I- g4 x) i9 p; w; ~' O) U3 j* n
$ L& I& e$ c9 z# X
allsel
& D7 n! h0 D# ?* B$ w8 |2 @' ]+ x8 m9 y
rbe3,node(0,0,210),all, sla,sla25 B2 u/ g2 o" K7 w- g
0 Z7 l6 \6 [' @: W$ N# l /SOLU3 U) @: x- A T2 \# @0 f
. V1 S1 B$ x6 t r( [: w t f,node(0,0,210),mz,10e5
* v9 m: J; G6 m: E7 u
& U# l, g0 S* _6 Y# H FINISH ~7 g5 Q) u' z# y9 x) S% o
- E+ ^; D+ s+ q0 \' \: q7 N1 _* l
……- w+ W( ~# Q! C8 q/ n* Y
2 Y N# f- j0 h4 q+ {; H) c$ Q
方法5：定義剛性接觸面
1 L# m5 O& ]; J+ [+ Q4 K q$ L D& x+ k9 ]) C5 `
……
1 G) G9 D0 {0 r
( S8 v5 H4 S! G) ]2 P n,15000 ,0,0,2001 H( m( c! N6 m5 Y7 \! N! p: k
; v( ^4 @' e7 t
MAT,1
: s9 }/ B( Y' k
+ m' Y6 m4 @! A' S/ `7 H R,3* R) X- `: D. ^- d7 I L
, L: i- u6 f* i REAL,37 _, ~: G5 f$ S. t
2 c x9 M& c/ S6 J$ i ET,2,170
- n0 }& {1 C. |8 j% ]8 {5 I7 l" \' Q1 O1 h
ET,3,175
* W! c& Z. n* [3 B# p0 f8 I4 l. @* T1 z3 {% ~. F" a/ v& p7 |* H3 P" G
KEYOPT,3,12,5
- h2 P( Z0 n$ i R) D5 |; B
+ d: p0 A" M4 y. N# t$ m KEYOPT,3,4,1
: g+ q) H- N7 a3 {' y0 q, f, p
1 S, c# n# Z, O KEYOPT,3,2,2) g* ^5 H. N; H9 e
4 R U* ]1 M5 D0 \
KEYOPT,2,2,0- B! T+ H! i5 @5 J4 y1 a6 Y* o' N
$ [2 j" y" m* j KEYOPT,2,4,111111 M3 y+ D2 P. I: k6 a* q& v* R
! y/ I: c: t) X5 [% x# E TYPE,23 f4 n9 J# G. u
- A8 m; ~( n' d/ w1 X/ u7 I
! Create a pilot node1 g: X* o7 _2 f/ `& {) `
$ Q( i5 d! i" V% E; T; |
TSHAP,PILO
* y! } C: r) U8 D" A5 H! ]) h! P j
E,15000
- t* i8 d: }/ N3 }; K) b; ~6 U
' A8 }% B% L G ! Generate the contact surface
4 J% m$ ^6 d* R* f/ i2 J8 e9 B3 _$ a% k' k. k* O$ v P
ASEL,S,,,14
5 c! ~8 D5 P1 _+ O1 R0 g
2 a; }6 V! R7 G8 C! E ASEL,A,,,19
2 Z+ p/ d5 W/ k R7 J Z, |$ t5 X& ?- y* c" B1 z/ [, i
ASEL,A,,,249 z* m5 A! A# x; g, ~ E
1 ]8 y% T0 o. ~: Q& R; m0 X, [- Q ASEL,A,,,28" [: o/ V/ z/ Q a
/ e( g6 x, n2 i
CM,_CONTACT,AREA
! I1 l# ?! t: Y# |8 ?% S7 n# K" E$ A. k2 s2 S0 Y: q
TYPE,3
, q! C: i: B6 h C
9 g# I% D1 X* P* u/ q NSLA,S,1- c6 S" ?% ]8 |/ [
8 ~4 @8 l9 b6 W. {1 P' O ESLN,S,0
/ v: D' S' b- P1 y: `$ J7 \6 ^, l9 i4 d+ w- a1 `2 H
ESURF
; |8 n9 x6 S5 c
: h5 W2 \: G* \% F ALLSEL9 o( ]( V+ S: R7 V
2 L/ q( u2 R2 l1 v+ d5 y allsel
: f5 {( G( E% j- X% ?3 [7 C- M% v, V# \9 O
/SOLU7 Y4 Q+ d+ z$ X0 k7 X" |
* U, s z8 D1 w1 k
f,node(0,0,210),mz,10e55 G! u4 ~6 A& C: i+ ^) u: E; R% }8 Y8 W3 e
" P9 Q+ X. }8 O: D1 G, C X% d; k: A FINISH
" B: \ ^+ P2 {9 V6 ]9 l9 p
* b0 N( Z4 V! u' c# I" h7 q4 F …...4 A' t$ i$ X. ]/ F, I
|
評分
-
查看全部評分
|
發表于 2013-7-25 14:59:05
QQ好友和群
QQ空間
