|
发表于 2013-7-25 15:09:12
|
显示全部楼层
我邮箱1962940284@qq.com
# ~" k- x* D J# H! r/ w/ W施加弯矩扭矩的方法其实有很多种方法,在这里介绍其中的5种,并进行比较:
! }. K7 }$ _' Y$ a
% ?$ l! X( u/ T" D1 [$ e 1.将矩转换成一对的力偶,直接施加在对应的节点上面。
# I8 [8 K: I: I3 x7 T, x5 _" ^; k" ^" k( }6 m
2.在构件中心部位建立一个节点,定义为mass21单元,然后跟其他受力节点耦合,形成刚性区域,就是用cerig命令。然后直接加转矩到主节点,即中心节点上面。
* T; k( N: P# E" ?- W" t `1 z# N& N i& R0 V* @
3.使用mpc184单元。是在构件中心部位建立一个节点,跟其他受力节点分别形成多根刚性梁,从而形成刚性面。最后也是直接加载荷到中心节点上面,通过刚性梁来传递载荷。, ?4 J% N* R9 i* L* \* m
8 ^0 U/ |% Y6 X- B" V; d* r2 j2 F 4.通过rbe3命令。该方法与方法2很接近。# T6 M1 i5 o9 Y" t+ j( c# j
* J0 ^+ M5 c* J. j! | 5.基于表面边界法:主要通过定义一个接触表面和一个目标节点接触来实现,弯矩荷载可以通过在目标节点上用“F”命令施加。6 l2 Z3 u5 O8 u, O) Y
0 C' e" x1 A+ _7 Y5 m6 P 对于方法1,通过转换为集中力或均布力,比如施加扭矩,把端面节点改成柱坐标,然后等效为施加环向的节点力;而施加弯矩,可以将力矩转化为端面的剪切均布力;但这种方法比较容易出现应力集中现象;6 E9 a# i' | d! m4 ^, ~$ M% B
! q$ ?/ n6 O7 u2 m, c+ p% ^
方法2,定义局部刚性区域,施加过程venture讲的很详细,这里就不在赘述。根据他的例子,我在下面给出了一段命令流。该方法有个不足,它在端面额外的增加了一定的刚度,只能适用于小变形分析。! Y" N# H9 Z! u& g: N
/ Q- E% w$ h& ]6 M, U | 方法3 ,相对方法2来说,采用刚性梁单元,适用范围更广一些,对于大应变分析也能很好的适用。但在小应变分析下,方法2和方法3没有什么区别。
0 W' m( f8 |+ ?! N
: F+ {% i, A' B- t. @) h 方法4,定义一个主节点,施加了分布力面,应该说跟实际比较接近一点,但端面的结果好像不是很理想,结果有点偏大,在远离端面处的位置跟实际很符合。5 r+ u- N9 D# Q" ?- y: O
8 y; b$ S. \+ A3 Q9 w1 m- h 方法5,它具体的受力形式有如下两种:8 z) G0 E# _$ q! P- O" d' Q
9 A" X) T8 P3 T
刚性表面边界(Rigid surface constraint)-认为接触面是刚性的,没有变形,和通过节点耦合命令CERIG比较相似;6 _; b L5 U c6 A* \" i0 [" u* Q0 |
) h/ p" M7 w' s- I+ O1 y2 g8 B6 _; A 分布力边界(Force-distributed constraint)-允许接触面的变形,和边界定义命令RBE3相似。
5 a- E0 _$ I# j" \" ~7 F* ]7 n9 i
7 b9 A) g+ @% d2 O) X 使用这种方法,需要用KEYOPT(2) = 2打开接触单元的MPC(多点接触边界)算法,
; f& T; E8 J# @8 e$ C8 c. \: `
2 g Z- V, j4 u$ b 下面针对venture给出的例题,用不同的方法来实现的命令流。% c9 p- {6 C Y0 D' t+ T9 ]
) ^( B# F" d1 m# j2 W% R0 u 方法1不介绍了,方法2:
+ _4 b/ x, y7 [5 _3 g, u1 r6 V5 q q, Q$ s+ K4 W& ?
/PREP7
5 [5 p) [8 L0 Y, u! n& r! C6 `) {; o% n" j* c* M6 H% u
ET,1,95
3 [, _& _$ n4 w$ s, _- r/ r
5 Q8 }- |! ~+ j2 w% U5 E3 u ET,2,21# F4 j9 h8 ] j6 p
% s! W% y/ `* w5 _
KEYOPT,2,3,0
. j" S2 g# D) W% ]9 ^& w% b( R
/ I* N! `; c- @( y R,1,1E-6$ O! T0 q2 m. |* v4 \) j
+ I7 E+ j# S$ [1 w7 T1 @+ o% g MP,EX,1,2.01e5; ?& {3 Z% V& _( \, X: ^+ |
- v* D" w; S( Q. j MP,PRXY,1,0.3
1 l0 p" O; e a2 n1 O, e
( l8 A p$ x0 [ CYLIND,15,10,0,200,0,360,. Z! p5 W$ F& N$ T) L, m
2 r, z M9 j, u" Y3 s4 H0 t7 x) _ wpro,,90,% F, D" X. E6 Y. O! g
" a' g' C* n" D* J* \: i
vsbw,all1 X3 x3 U0 ~& k$ p4 y
# H+ I. z3 Y6 o/ x2 C: o6 o- F- L$ e wpro,,,90
' x( e( s2 N! Z8 N+ T; ?
5 J1 R2 W' F- j- K vsbw,all# n7 I0 ]5 p$ I- z, ]+ i. M5 e* n
' E% {7 U8 @2 c) j3 Z9 Q+ z
WPCSYS,-1,0
! Y( X( @+ Z5 B4 c2 }- y" u- P8 `3 J _# [9 K& s: `7 \
K,17 , , ,210
& G) W6 ]' z+ j6 x& f
C$ {$ |8 D" H- P lsel,s,,,13,16,1; l. l2 F L, b7 R) K6 b8 @
% I; B: j8 ]2 V) ]% C" D
lesize,all, , ,8, , , , ,1* S u+ l- b! b- \* ?: L O
6 E. {4 _& x# `- y G- F7 t; L( N1 h
lsel,s,,,22. z1 w/ H8 g8 g: P) G
. f9 v% M. a$ x- S I" C
lesize,all, , ,4, , , , ,1/ ~+ f( ~3 |8 ]" ^1 I" P7 q
) v0 R1 Y. o6 t6 l1 x, t% }; Y lsel,s, , ,17,20,1
# b. M. K8 z% x& i6 W3 A7 x' f; Z2 X$ p2 X7 y3 b# d
lsel,a, , ,26,27,1& T0 C+ h8 Y! p. T5 S- ?9 D# Q
% v# L& T4 p/ T
lsel,a, , ,30,31,1
- A' V8 U6 E$ I( u+ V4 o9 O$ E: X5 F" g0 }, L
lesize,all, , ,20,0.4, , , ,1' f& W2 k4 j, b; s4 U4 h
8 o2 p. D! Q% B7 P. o2 Y9 c
alls
* t* _* W! H( I% O- F% _* N+ k5 Q( I( S
vmesh,all
/ K$ W+ P1 ?- `, d- F K
# `) g* f, X. ]* y5 K7 G# u8 S5 t/ _ !!!!!下面一段开始各个方法有所不同,由于前面的建模一样,后面的例子就不再给出
+ \$ N3 _! Q4 ?. N2 ~
( P1 J; X: A" `% U7 c/ K( \ ksel,s,,,17
8 c+ A, f; u& z; F5 i# d2 F# O
type,2
, T" U) C) R+ a* w" m# g( S% _
2 }0 I/ g: j# W9 b @9 i. \ real,1
6 n& G! g2 z `4 `5 @" \- U
, i5 t' Y0 n, M2 w$ _/ R kmesh,all4 N* V' r' |6 v+ X$ ~( d6 p9 ~
$ y2 s0 x( U# r% J m8 Q; W allsel
1 A, T0 I, e/ h0 w6 N: E1 A6 a `
' i4 M% R) p. _9 e3 d: n* t$ g nsel,s,loc,z,200,2100 l* a7 `- H* j9 Y- n4 Q1 b
) v& U) W" ~& k npolt
# Z: W/ p* p0 R3 _ c/ J( K. B7 F
CERIG,node(0,0,210),ALL,ALL, , , ,4 ~$ p& G/ }# m2 P) v( c8 Z1 F+ a, X
8 N. G$ k8 Z# |/ b4 [& U' i7 ^
!!!!!CERIG命令定义局部刚性区域" `: ? J' i. K4 ^
# N# y) i& v- Q& ~ V+ { allsel) b! \* p( `, H: w, k
5 m2 T3 w7 @( X) U: _
/SOLU" h }. M9 g( m4 L3 x" _
& W. u5 Z: ~# o% I }( y f,node(0,0,210),mz,10e5
& n& `9 j f7 ~5 y7 M/ ]
. @: N8 u' Z' r0 R) k r FINISH
3 q5 ^( [) k* ]
# ?0 o( `8 B/ y% r !!!!!以下一段边界条件的施加各种方法一样,后面例子也不再赘述- k ]& X" S* T T2 B9 {; _6 C
) {8 _$ ?1 c# R; K; v! o* m# D
/SOL+ H. z9 X7 J" `
/ z' w5 B9 A, B5 K5 G
nsel,s,loc,z,0
' ]" T: n q+ D- B
3 O- V8 o6 `+ H7 }. W* x- b' S d,all,all
( w x- O0 Q& y3 c
5 e9 D7 p9 k6 @! i8 O allsel# M4 L/ E/ r4 H9 }) r( O+ ?, p
, u# q- v! a$ K7 c! _" S* @
solve
5 B, r4 q5 V, m- \4 {
. ~* A1 M1 Q& y4 b( i 方法3:使用MPC184单元定义刚性梁. h1 U9 \6 \# o/ o; C
7 n5 Z- Z" ~* u9 M5 R7 k3 D6 a
……
9 f) j6 s) S9 |) J4 G. p! c/ c9 R9 e# P8 O, L1 L
et,2,184
# v) [0 H- b; {
, k. C% Q) q! \4 L' S+ { R9 ^ keyopt,2,1,17 J: g s3 \; E
3 s( [$ w" e) O0 ]1 V nsel,s,loc,z,200
) ?+ w& |6 S8 \" E5 o" b
5 @6 A0 Z3 J. c: r! z n,15000 ,0,0,210
! G5 A4 {6 l( G* x; z! F8 z6 p( C% A
type,2$ m" V; q0 g( Q' B! Q
1 u( j3 p% C. p& m+ m) N# h. f *get,nnum,node,0,count
. G, D/ ^: j, c
7 Z: K5 F4 o2 `, R) }4 n- X *get,ND,node,0,num,min+ T6 [% [( Y E+ C
' L) Q; x# k5 E, C% | D5 S
*do,i,2,nnum
* h! U- z3 l' v/ g
: ], Q2 B9 k+ ~+ P !!!!节点个数是nnum,只需要生成nnum个mpc单元
" n. f2 o/ u- B7 o5 W/ ?9 [* T8 Q/ [3 e
E, 15000,ND
, V: a- j- _" M1 e$ \# y& N0 q* K m$ s6 G& h7 x2 `+ L2 |
ND=NDNEXT(ND)
6 Y+ |. b- I4 {7 A8 }' L$ x- Y1 a3 z9 P* s5 h
*enddo) p6 Y- _; K+ w$ b' Q) I
- y! [9 a8 X+ [( ^5 Q1 B, n' t( ` allsel
* D, }% c$ i! W) E' S- S9 v4 j1 O8 J; f
/SOLU. E: E& F5 ^' d: h
1 l0 K" x# s0 t# c# F f,node(0,0,210),mz,10e5. q; O* v, [: X: O9 d/ k& O
3 h1 l6 N& H% T3 \: C
FINISH
+ _& ?$ {, d, k" e! a& Q' h2 C; u% D4 z% {/ a- L4 S
……$ C% I: L8 N0 F" T8 m
! H A' p- g/ b; S. Y
方法4:rbe3命令
% S) N6 |3 A: i' Y5 _8 Y8 N! o8 V" b1 J2 T5 }7 u4 u+ {
……
/ I0 y' t' j& q& x; {( ]: a! [4 `0 H* E8 g
ET,2,21* k# s% s* h) k. S$ i o0 x
" C2 d# C, \; G3 _+ Z
KEYOPT,2,3,0
+ j$ f+ g1 i9 O; Q" n* w8 c
) R' Z, c5 \. g& h1 z; i& | R,1,1E-6
; J) ~5 u( e3 Y8 D8 z+ E/ n/ z( Z; Q6 X2 ?) s l7 f; d' n3 Y
K,17 , , ,210
$ S2 k. l$ `4 M$ q/ U1 w; x! F& f- Y+ c2 @, M6 P; e
ksel,s,,,17- ?- T6 e' p/ r7 ~- V1 z; M
3 j7 _3 q3 n3 d/ l) m8 G" ~% F& ]! U
type,22 m: S7 D( P+ F2 I* I, n! ~$ U- n
2 ~' L; P/ `. _7 r; u e' H: n real,1, Q& i [& |' U' z
! B& O4 d" s ^( Q0 z kmesh,all; K. J+ g; v0 D
+ l5 S1 z6 S3 ^/ s( z( G allsel- f) f: g+ D) {1 k; ^' D. N
& I! E! C% k: J5 G! T, J6 C- r4 X
nsel,s,loc,z,200
& P& V5 F$ U' S1 N2 T; g/ J" m
2 S5 J4 Z% z1 G- B *get,nnum,node,0,count
C8 _+ [6 c/ z1 E( V, `4 y( d$ w/ ^" I X8 ^, Q) z1 p) Q3 r- ]
*get,ND,node,0,num,min s# o3 Q8 I4 H+ P$ y0 ?6 X
: Z+ Y. x; X/ }+ U; L5 y( z- e *dim,sla,array,nnum5 d" L" u$ Y* J1 v$ @5 W/ L
+ b5 x! ^8 F+ C! b' W" ^" n" { *dim,sla2,array,nnum B. d" X* |" j& M
+ z5 ^! c- F; F$ ]5 a: m. f: y *do,i,1,nnum
/ l# q- h3 q* d n
# p6 `2 X- x1 | sla(i)=ND) i: f& f1 J9 I: d! N0 }% q
" c( v, l8 h! P/ V+ n0 N& }
sla2(i)=ND
3 H- S' W. ^5 q, ]' J
" u, j6 y$ |# b. V! r ND=NDNEXT(ND)
. @: x. T' g' }" J) o3 h3 K+ B# x3 d
*enddo- |8 f& d3 |' n* o/ `
2 e* V) X( X" `
allsel( C. ]9 D4 _2 e F+ K% x. r! l
; Z5 d) H4 P4 l# b* S2 r: x4 I! u rbe3,node(0,0,210),all, sla,sla2
( T% m$ j. ~! i+ q7 B) Q2 o; s& W
5 @2 U( T& ]% D; ?) J @3 N/ S /SOLU
1 d. I" a; F$ V" T9 }
* M( z3 ?! @# Z5 t3 | f,node(0,0,210),mz,10e5! b3 h( x* e4 b, r, K
# H% P/ I& n& P
FINISH$ Z s8 B8 n" n# @
# ~" f. G# y; F+ k. ]7 K2 Y* H. B
……
1 s. b+ C a$ x* u
) W, }' v2 _* o! @ 方法5:定义刚性接触面4 ~6 m2 Z! d) t" I/ h7 [$ c
% p. z$ g/ k3 K1 J- S7 n$ J2 k' q* b
……
! A+ U+ c0 y# o8 R$ r1 Q- [
5 @# l( k4 A/ A n,15000 ,0,0,200
, |% {; f- y7 R- F' X5 u q' a0 z5 B V2 p6 h
MAT,1 P y2 S5 D0 b, v$ }# N* q/ W
( Q/ R% r/ a& k5 b8 m% T R,35 Z! A, J) O' p4 n7 `0 f
+ C# |; x t3 f7 U
REAL,3
4 I @. _$ ]3 s
% V1 |$ e8 ] m& l- q8 O+ L+ d ET,2,1709 s% p* R0 l$ S& j6 F
/ {/ _9 d: ]) ^" u5 Z' Z ET,3,1758 q4 A. [0 H! H* Z+ {8 ~) a
. G# j: H) m! x8 m- H8 F6 e
KEYOPT,3,12,5
2 r: L5 w; n4 G' `& A% }1 N! t8 L$ j% A. ]/ g: b; R
KEYOPT,3,4,1
, k0 w: j& y' S$ M( B
5 T/ ~8 E# U( j KEYOPT,3,2,2; m5 a7 S) {+ p9 m5 ~+ @
. t# z$ K, C7 `
KEYOPT,2,2,0, J+ h6 M% }' D" z8 L+ [+ u
; p( a7 x9 X# a' F9 z6 X! E6 J( x# c KEYOPT,2,4,111111( x( P9 O. [* r# q$ `) }" i B
# Z( e+ h( ^0 G2 d TYPE,2
& e& H# [/ p/ s; \" G3 J& M% R3 q5 q% J# S: w, @: g4 i
! Create a pilot node
5 q( P# _4 C6 P. A$ l# W: R: `6 S2 a1 S: h' I/ c5 X* Z
TSHAP,PILO
8 g) a) _, k7 {- G
6 K* S( B- Y4 q1 s( E5 x! O. X E,15000
c1 t! |9 j% M8 z7 x3 K8 P) x$ a: k. O. a6 f) ?5 T! }
! Generate the contact surface
- v5 P0 \/ s$ l! ?- o) ]5 T) b7 n
# Z T8 N* x+ J ASEL,S,,,141 e* ?+ M& t1 Q: D" F) E
) a) n. [- H. d0 X& \ ASEL,A,,,19& @' b j% G4 t; \5 i1 y% E
! S0 q7 y/ Z" t* a+ H ASEL,A,,,245 M( I& Y" P: x$ ~! q9 }
/ h, `) V4 e/ M! ~! k" h* G) ]( D ASEL,A,,,28
) [0 \8 H2 d% H9 b, q( ?
* Y. |2 G1 h; w2 u CM,_CONTACT,AREA
. V2 |$ m( N q8 I% w' j3 E! d# p* c* G
TYPE,3
% j! g5 Q6 x) S( i& s. B, O1 {& B% X. C- x
NSLA,S,1
' {# f$ Q# k, m. }/ {6 j X' Y
1 { t" j5 g( l1 B& R8 t- b {, L ESLN,S,0
. U; @) q1 O4 n4 w- Y6 `6 A0 E3 v* Q
+ x) V6 j, l7 U7 x4 z- V8 B! ~8 v ESURF" g) @7 _1 d- t/ Q9 U! Z- J. H
3 D' B# T- S: a6 Y! r1 |! { ALLSEL# U Z1 y* A) F, A# T" f9 K
$ m6 M& I- c0 B. S allsel
2 q7 f. z% ?" S4 f% v) Q, s! Y6 L( Z) ]: B9 `
/SOLU
3 I7 _2 B( ]$ }9 e+ ~& |/ D* G
" h3 V1 A# g# ^2 m f,node(0,0,210),mz,10e5
( m$ [' b$ K8 g. t% Q8 E2 R
# k7 ~7 o/ i2 J! H+ O0 N! h* H2 v FINISH, e' C3 b, U% v0 ~5 i
. ?+ Q* ?% t. j( m …...: Z! B( [6 a7 E" c3 q% q: ?
|
评分
-
查看全部评分
|