|
|
我郵箱1962940284@qq.com! z2 x+ ^/ O' T5 p, T
施加彎矩扭矩的方法其實有很多種方法，在這里介紹其中的5種，并進行比較：
, ~9 \9 G d9 ]$ L
7 S% f" i: P# j 1.將矩轉換成一對的力偶，直接施加在對應的節點上面。5 u/ k6 b7 s! ]) A% _8 b
5 m/ o) l9 X& i. T# \) T! P2 N
2.在構件中心部位建立一個節點，定義為mass21單元，然后跟其他受力節點耦合，形成剛性區域，就是用cerig命令。然后直接加轉矩到主節點，即中心節點上面。
, k) @) B9 |& t# \- J
( N# i9 V) @: t: I6 S0 q 3.使用mpc184單元。是在構件中心部位建立一個節點,跟其他受力節點分別形成多根剛性梁，從而形成剛性面。最后也是直接加載荷到中心節點上面，通過剛性梁來傳遞載荷。
+ d& C7 q1 l4 E9 p& [1 q
: m1 O+ R# \3 }6 M 4.通過rbe3命令。該方法與方法2很接近。$ L! x" D7 m; B" y+ w
) t0 u5 d5 j9 e. P; `' M 5.基于表面邊界法：主要通過定義一個接觸表面和一個目標節點接觸來實現，彎矩荷載可以通過在目標節點上用“F”命令施加。3 H y+ g) B3 g' I1 p0 j
. e" ? x7 }2 {- ?% _$ I
對于方法1，通過轉換為集中力或均布力，比如施加扭矩，把端面節點改成柱坐標，然后等效為施加環向的節點力;而施加彎矩，可以將力矩轉化為端面的剪切均布力;但這種方法比較容易出現應力集中現象;
/ ?9 C+ A6 A: k4 q4 l2 a& I3 F5 `/ e5 b1 ^9 n d8 g: Q' `8 {
方法2，定義局部剛性區域，施加過程venture講的很詳細，這里就不在贅述。根據他的例子，我在下面給出了一段命令流。該方法有個不足，它在端面額外的增加了一定的剛度，只能適用于小變形分析。
4 `& F/ n/ b. ]& j( {2 P8 H
! d$ P# V, v7 ^$ `; j 方法3 ，相對方法2來說，采用剛性梁單元，適用范圍更廣一些，對于大應變分析也能很好的適用。但在小應變分析下，方法2和方法3沒有什么區別。
& }( z) {* l% F. u0 }4 O8 H7 h& T: e0 a: \/ h
方法4，定義一個主節點，施加了分布力面，應該說跟實際比較接近一點，但端面的結果好像不是很理想，結果有點偏大，在遠離端面處的位置跟實際很符合。
3 \% L F9 I4 Q3 l$ X9 R) ?* K6 ^/ ~2 d: Z
方法5，它具體的受力形式有如下兩種：4 W7 r$ x0 M2 w6 z! X* l E+ a+ \
0 J2 R, a9 `' Q0 l9 R; N
剛性表面邊界(Rigid surface constraint)-認為接觸面是剛性的，沒有變形，和通過節點耦合命令CERIG比較相似;1 B+ @& t- _4 E7 _
& m. ?0 \% @ C7 Y# @
分布力邊界(Force-distributed constraint)-允許接觸面的變形，和邊界定義命令RBE3相似。' [; k4 G$ B+ G0 A1 n
. t7 |' ~" v; c2 o9 p
使用這種方法，需要用KEYOPT(2) = 2打開接觸單元的MPC(多點接觸邊界)算法，! V; t$ s4 q8 Q4 H
/ k( Q- y, d7 K
下面針對venture給出的例題，用不同的方法來實現的命令流。
: r7 z- w# y, a) H+ r# w! g" C( k; E
8 ? U3 D# P. r! T& R, l& m 方法1不介紹了，方法2：2 q1 u: W- s& C* y
2 r8 O) P) P' Q# ]0 J0 M% M+ J" d /PREP78 h" g# E7 L! o3 z% q
% r! U/ \, q' \6 n! ]; a9 V: h ET,1,95
! K9 e/ A+ _4 C5 q, M3 ?& X3 F) f4 {+ L7 N1 }
ET,2,21" `9 g: F; d! v, W) o
; o9 }6 S2 Y( Q+ W. {: Y KEYOPT,2,3,00 ~6 U; r( x5 ~9 z$ F
. K* e( h p/ ?2 S2 i9 U
R,1,1E-6
# [& I2 d- r: v/ R7 b
! m" i' M0 U7 b% ^8 V; g MP,EX,1,2.01e5
$ f, | l9 l5 `+ }, Z ]* T
7 c/ u4 S( o5 m MP,PRXY,1,0.3! L0 Q% H7 ?3 K! e
2 a0 W% h7 ^, X. R# |- Y$ D1 X5 M CYLIND,15,10,0,200,0,360,
- J3 g' x, L4 G) A" o
1 r; p" [/ o5 v wpro,,90,6 h2 \! h5 ~! Z9 o2 D' F8 t
5 \) Y" u8 A1 q- `$ X
vsbw,all
+ Y1 s: a8 Y! D' w
9 f! {' R4 S' h, d k wpro,,,902 N6 m. Z* Y7 C8 e0 X
J7 D% P$ Q" V
vsbw,all
3 i, E) l' X; W' \
. O' ^3 N6 ^9 s& _) \2 ^, N5 e WPCSYS,-1,0
6 {7 a: ], V1 M( I7 y- s5 `% @' U
K,17 , , ,210
# C' b5 U8 Y5 E9 [7 I- x4 [, j, ]. t9 J: z$ s# Y- F+ Z6 \0 ?' q
lsel,s,,,13,16,1, g; e; S8 v' X! R, Q9 S0 J* |2 {
. a' w+ D0 U: w lesize,all, , ,8, , , , ,1
; q" g9 o8 x a7 I. I; e; G0 b* r ?) o! _% {
lsel,s,,,22
' s8 M) {" ` d1 ~2 t. H9 h$ t/ b0 r: s; s, D
lesize,all, , ,4, , , , ,1" ]9 h- a8 q" z, a% D- w9 O
; o4 H* f( _& d1 D lsel,s, , ,17,20,1; \ \2 E0 }$ s5 ~0 |
$ V1 B& ]6 ?+ Z) t
lsel,a, , ,26,27,1
9 w4 v: R, V& n+ R
" k# e; B# L0 i8 z' J0 O lsel,a, , ,30,31,1) @1 B' c! G6 z2 j# x" t' }: c
5 Z/ j5 O5 B% P: @* l" Q: ~. }
lesize,all, , ,20,0.4, , , ,1+ R w9 u" U0 [6 K$ z
% P l2 ~: |6 C( ` } alls
& D i( F$ q9 F( g* I7 T. T4 }
4 i3 T0 T% V! ]% a" x vmesh,all! o( E6 c4 b( |4 K# Z
8 A: H1 ^& P" X( r: }8 t' y
!!!!!下面一段開始各個方法有所不同，由于前面的建模一樣，后面的例子就不再給出
. U; u5 @. y2 M* J8 e3 z5 [ w) v p' v2 a: o1 q
ksel,s,,,17/ l% G, R) `2 P# m! T* S. E; K
( g2 f h' A2 `1 [) b type,29 W8 z7 E7 e0 v
( _5 }! X/ v' W- F
real,1
8 e! c3 S& d- l% k% o
7 H" p) T8 O: |, C5 M2 _ kmesh,all
( T7 i1 P( v, g- B1 g1 ?* e5 N9 B8 j! W5 Q4 _) L
allsel
, G8 k6 D) V# K" R1 }
2 Y, O5 n- {0 A3 z nsel,s,loc,z,200,2101 d* D) A% \; {1 u
; t9 n6 t# \& L+ J) m: y w npolt1 j3 o8 D1 d! m# I
5 q* ?, m) b t4 r1 I CERIG,node(0,0,210),ALL,ALL, , , ,: @0 Q2 Q$ @ p! N" l, B
/ H0 P/ C r' y, ~- L !!!!!CERIG命令定義局部剛性區域
, }8 V; V ]' D7 A7 k6 Q
; o) Z+ c5 C( W( c" m5 Q- ] allsel$ m: H& O# z6 P `. M
! n/ g+ G$ `/ ~
/SOLU
& |8 O. K S: l' W+ H9 _$ K7 [- b8 o% V
f,node(0,0,210),mz,10e5$ d& V% _- T: i* x- g
2 i5 v5 p$ k, i- x+ b0 ]
FINISH- O0 [; `9 d& T4 i7 i1 |
6 U. K8 r$ @8 X; h
!!!!!以下一段邊界條件的施加各種方法一樣，后面例子也不再贅述( f1 ~+ ? e2 \/ [8 t2 K
7 ]; Z# {! B& i) _4 ]) I9 b /SOL, K% w7 a' ?4 B) y
5 ?' a* Y" E# s: `, R$ I& }
nsel,s,loc,z,05 P% N$ B5 d5 h
6 w* s& Y, \/ C2 p" g
d,all,all
. R0 E; }+ t* _! t1 O+ S: U& t' m: }/ Q
allsel0 a8 g5 O$ t6 D$ V
. h `1 \& O& O solve
, F; w0 d U0 X! W
* s1 H+ Y0 @3 M. q" v8 y; Z& M 方法3：使用MPC184單元定義剛性梁
. s; }7 K2 C. V; p" R
! K: H+ K+ {1 e3 G3 I) `3 s9 p ……' [+ l& f* r9 o8 B: O
, @6 Y/ q# C0 J* b
et,2,184! ~4 N9 ~; I" U# J7 Y
; u$ M. ~/ @$ {6 Q) |8 L% D; q keyopt,2,1,1
6 z' ^# d& L4 b0 M! h9 l7 Z! y3 Z R
nsel,s,loc,z,200
) N$ o* Q! ?* \; s, C6 e
* F9 t) d( z3 L* ]2 l n,15000 ,0,0,210
1 B4 L" ~- t6 ~/ K" W
- s+ |9 Q! B/ r. D& o( W1 u type,2+ J! v% g! G' r7 r, ?( n
8 t3 m. Z) c1 }9 P
*get,nnum,node,0,count( G, X( X& l* `/ y1 G9 W
, j, g5 ?3 `1 O- C% m *get,ND,node,0,num,min
* s2 d6 i5 I2 s6 F8 r, E" U" c& g: H0 {/ g. L- t* C4 I
*do,i,2,nnum8 O) v- f4 M! [* {7 \
! v; ?# G6 C' B% t, K !!!!節點個數是nnum，只需要生成nnum個mpc單元0 m7 W3 `* k) X4 I
/ F0 }% g! E& W9 v7 R0 T
E, 15000,ND
: k: ?: P( y3 e$ e/ X/ |) B Z7 W# g: o# h( J% u$ j. l8 k( b0 J$ G
ND=NDNEXT(ND)
$ I$ b" c0 ]3 _8 Y, S" l. E* K% S$ a* Y E; E
*enddo
. v1 d+ O0 f2 C& t/ L1 p/ c% L1 @$ X9 c3 W0 ?3 M9 B3 Q
allsel
( H6 |. ^9 _0 U/ c7 y/ |! v
/ b; Q7 l9 Y7 [4 Y% `1 ]( W. A m! T /SOLU
, E! v b, m4 l
1 _& G8 Z0 K( w) p f,node(0,0,210),mz,10e5
( s' G4 [- s9 q. O
2 `4 u9 S) D0 L/ Z0 T+ ] q FINISH
; J8 J- y! |0 N/ V" J
: \( i; d. Z2 }9 S4 @$ ^' x ……
1 C S% ]( [2 j( _. z) P7 x. z* k8 {6 H+ V- L8 Z& I( q
方法4：rbe3命令
+ B+ }" {' \& v' \' D# e6 i# L+ K0 u/ E8 q$ W$ a9 u# X
……
& P4 Q+ S% z1 B' _% P$ G5 `
5 h8 ~2 l; }) W+ C. H+ p ET,2,21
/ N" u0 ~4 M; e% [4 F0 f- x/ q
+ \0 b! q4 ^/ W5 W3 P KEYOPT,2,3,0
, N7 R A$ c* O/ o' u+ P
2 t# M5 J2 F9 s) l5 O5 ~ R,1,1E-6
/ V Y! v8 q: H; F. f1 [- G5 {5 }1 I' J1 j- q
K,17 , , ,210. x% K7 y( T+ X w. T3 z
, Z6 n+ m1 W9 o% v" S
ksel,s,,,17
G- u* a3 C7 u
+ F, T* x! N1 i# `" A( Z type,2
! v/ e- P. j1 w' A% r o
# y4 G0 I' C; H+ K! ?/ b8 w4 y2 J real,1
) R* t2 l H* h1 ~) e" r9 u0 e' P9 v0 K0 S. y1 \7 d
kmesh,all
' R* f$ }6 X, S7 C1 e2 z) [5 L! x$ _. N
allsel4 S) p0 z0 ^9 B/ j% _7 W" U: u
" }( e. b V) K: g* r8 ] nsel,s,loc,z,200: Y( x, w( W1 i7 d, `% |
) A6 P1 f9 }: ]" b8 [& [& d) Z$ z
*get,nnum,node,0,count
# j. d( b9 W; {! u" [6 S
' C" R, J3 E( @- ` N T% O5 j! y *get,ND,node,0,num,min
$ Z% W; [* t1 a% h6 ?) L8 H, N2 ]- ~. w# {! B T/ o; j/ ~
*dim,sla,array,nnum
: r$ n6 W/ p0 H
! A1 A% S Y7 y *dim,sla2,array,nnum8 l5 P( U n" a
2 ^3 m2 {2 w) Y/ F" I1 ^
*do,i,1,nnum
; Y1 o3 [1 d5 v7 J. K
- P# X+ ?* i, n* c9 |% D sla(i)=ND( ]1 x1 j8 {: q4 s0 X! ^0 {5 e7 G
! J9 d" G. g+ u. S' s W
sla2(i)=ND
% ~- N" Q Y4 n3 z, J5 i
! `$ O1 t* G1 |3 K/ u) f/ ]) r, J6 T ND=NDNEXT(ND)
0 f9 o1 _9 j) B+ D- f3 F3 T" E% l0 m4 s [5 A2 _. A+ D, {
*enddo; U0 k+ E/ x3 m' ?! T! A
9 w/ [/ Y) V! T/ H: M" j! d, h2 M
allsel7 z# s! c$ e9 p+ ?) e) L `
) ? w+ |3 B) P5 l4 J
rbe3,node(0,0,210),all, sla,sla2
+ a% ~1 P9 w; O+ ?5 v1 R: r/ e$ x8 A
/SOLU# f7 S1 w% s6 V X# d8 }4 ?
5 I( @: X. P! J2 Y' m
f,node(0,0,210),mz,10e5: g4 a% F# k' \( u w3 E7 [% m
! h! Q9 k2 ~5 P' h+ N7 D" u
FINISH
0 G; C/ X9 V5 k* {* i, I
% S, A w+ N% w9 ~ ……
8 K3 J$ C1 N. s' ~) P4 L; L
! J0 b* F1 d P# y, m+ g/ @ 方法5：定義剛性接觸面5 ~" g4 h+ \. D
$ [$ T7 S. P! ]$ U3 U, y
……
. Y- a6 V0 t$ l: U
, {# y! Z: d9 w* O n,15000 ,0,0,200& u% S& [$ x' s
% w; |( |; l* y0 z' u4 \8 x2 T) T
MAT,1: Z9 j2 R' _' C D: A p5 y/ o
' }$ V) X. [ z+ j9 L$ \
R,39 v/ F) D3 k$ B/ h8 Y
, z* r% g k. [$ d REAL,3$ U; a% n& [( `/ P: [, b
1 l" \" w" a1 O" E, {& Q' ?
ET,2,1702 o8 @! c0 J r
. d0 X; h1 b( f6 o w* Z4 m ET,3,175
4 D9 h1 W# y/ B- y& n
) b& i7 c8 l1 V, Q+ b4 t KEYOPT,3,12,5
* a7 M/ `! Z! I9 Z! P# g4 B% X8 f* ], E' U
KEYOPT,3,4,1
6 o( N( E7 z; t | H/ s! H
: b$ P+ o- Q6 n' d5 N3 J# q KEYOPT,3,2,2$ ]- B9 ~* V; h
- ] Y8 j0 A0 P) `& _ b% G
KEYOPT,2,2,0% f, T- K0 S% z; W0 y
/ H |2 n( I0 j6 f8 u7 p; Z KEYOPT,2,4,111111
1 y( ? c: Q0 Y* r2 H
$ l& A+ e* b. l+ A TYPE,2
% P% k1 C) m) E0 h
" F! [( X- M' Y0 n, ~! V5 P ! Create a pilot node
7 R4 O0 }& O0 p0 e3 L: k
" p- L4 L* R# l, i; B3 ^ TSHAP,PILO# C6 K5 m+ A: W S1 l
" ?: U" h S( [
E,15000/ |: K$ J( H) f! f F9 C* j6 u- W0 C
. G* \7 u" |& L1 p: V! e ! Generate the contact surface
& ]) {+ f6 |& M7 a) F& b& Y" ]/ k8 C3 ? H! r
ASEL,S,,,14
/ d9 X1 X) X1 E2 O9 N6 |7 ]( M6 L' ~7 Q; }. ?( O2 ~) u9 N
ASEL,A,,,19
2 e5 _) Z2 G% u @) g
- v9 N! x0 _3 B i* j% Y ASEL,A,,,24
2 D" ?7 k) B& c' q: j- z9 y, V) L& r$ g; J
ASEL,A,,,28: V) T6 X) g: k
0 B9 A0 t( e3 P# e+ v5 I" F+ x1 s CM,_CONTACT,AREA }6 s# i- I7 J. t: n/ k
! I, K) I3 o) Z5 T% B9 T
TYPE,3
2 h: w6 G4 Y* e9 `3 E/ v0 P2 W* e) {, q+ D0 a1 d+ [/ s
NSLA,S,1
+ Q) n8 p) Q# J8 E2 ^$ d$ f5 `& U, h5 K5 |4 m0 w i0 ]% @. x
ESLN,S,01 l" R, D7 \2 l3 ^( t! H. p; L
0 l$ D+ v1 K+ P- U
ESURF$ q: k5 ^1 M$ T" [
" p% ^& t6 @# ]. D ALLSEL
& F5 A( B( _; _3 e
7 u. [5 `( j5 ]5 v. I9 r; ^ allsel
6 ]- f/ t. ?0 Q6 T* F/ W/ `& q4 B4 X Q8 T4 ?/ T
/SOLU
( H% Z! f9 e& l
/ T& a. ^* o) j4 z3 u0 ?- Y f,node(0,0,210),mz,10e5
2 ~* j0 R i- m/ d5 X" R
7 B& d! P4 H: V: Q3 b* s FINISH
7 n# ~' h8 m; a$ X3 Z/ t! v5 v: B- t+ F6 X1 D8 a
…...
2 c* K% V1 U+ `/ i |
評分
-
查看全部評分
|
發表于 2013-7-25 14:59:05
QQ好友和群
QQ空間
