|
发表于 2013-7-25 15:09:12
|
显示全部楼层
我邮箱1962940284@qq.com- Y: D( |6 A) W3 O# p
施加弯矩扭矩的方法其实有很多种方法,在这里介绍其中的5种,并进行比较:
) n z. w0 h9 `! C, Q' ]0 O" J' p, \- l* Q
1.将矩转换成一对的力偶,直接施加在对应的节点上面。
8 B* |. ~5 P% j* Z/ N/ S- S+ b0 S- Z7 k/ N
2.在构件中心部位建立一个节点,定义为mass21单元,然后跟其他受力节点耦合,形成刚性区域,就是用cerig命令。然后直接加转矩到主节点,即中心节点上面。2 Q0 s& u( d/ v! V
% e6 X$ N! @9 N2 f
3.使用mpc184单元。是在构件中心部位建立一个节点,跟其他受力节点分别形成多根刚性梁,从而形成刚性面。最后也是直接加载荷到中心节点上面,通过刚性梁来传递载荷。
( z, |5 u( G( T( t k7 p! H8 m; f( O; ~8 R$ w/ Z7 e
4.通过rbe3命令。该方法与方法2很接近。* r* j* D- ^; w8 Q' n4 d
. o f1 E+ k @ 5.基于表面边界法:主要通过定义一个接触表面和一个目标节点接触来实现,弯矩荷载可以通过在目标节点上用“F”命令施加。8 L7 e8 G$ B, W: Y3 g' `
7 i$ L, G" c9 w8 c
对于方法1,通过转换为集中力或均布力,比如施加扭矩,把端面节点改成柱坐标,然后等效为施加环向的节点力;而施加弯矩,可以将力矩转化为端面的剪切均布力;但这种方法比较容易出现应力集中现象;% q. \% u" Y" |: B
3 I7 G# Z# E" C) X 方法2,定义局部刚性区域,施加过程venture讲的很详细,这里就不在赘述。根据他的例子,我在下面给出了一段命令流。该方法有个不足,它在端面额外的增加了一定的刚度,只能适用于小变形分析。
. I/ i2 k0 ~% N( b- l j: W$ a# M y: _& M" J$ L' m
方法3 ,相对方法2来说,采用刚性梁单元,适用范围更广一些,对于大应变分析也能很好的适用。但在小应变分析下,方法2和方法3没有什么区别。
% _# R$ B( N3 x7 U2 v: s. S/ F3 `* G- H' h- g# g- d5 D% W
方法4,定义一个主节点,施加了分布力面,应该说跟实际比较接近一点,但端面的结果好像不是很理想,结果有点偏大,在远离端面处的位置跟实际很符合。' y5 y, z: p; c7 V
- l* y' L+ D( ]. M+ ]0 a 方法5,它具体的受力形式有如下两种:" @* A7 E* n3 b5 o
8 ], G) P/ j, c
刚性表面边界(Rigid surface constraint)-认为接触面是刚性的,没有变形,和通过节点耦合命令CERIG比较相似;
* ?3 \' ~8 \ y2 k' \
. R4 G; H% x1 q- ~0 _# F 分布力边界(Force-distributed constraint)-允许接触面的变形,和边界定义命令RBE3相似。/ o6 }9 J. s t* d1 ~
2 x) _ |4 v$ V% i) l
使用这种方法,需要用KEYOPT(2) = 2打开接触单元的MPC(多点接触边界)算法,: M i3 i; f' N$ b* o1 ^3 j) O( ]4 n
# M7 P( C+ O7 O! U
下面针对venture给出的例题,用不同的方法来实现的命令流。; u: y- u$ i: K# O, c' N6 F
2 b8 n- U C: v; ^. c# { 方法1不介绍了,方法2:2 P0 [2 W( d1 ?$ y
+ L( J& m3 l# J
/PREP7! D8 ~% F# e+ p
7 z6 a% C% M a" [
ET,1,953 c0 s) \6 I, [7 E5 Y3 R: z. ^ I
3 | o! d. A1 e4 o3 a; r
ET,2,21
J) ?' b" l0 k9 x- Y( s
" f) C+ @$ c; z4 @- W5 h KEYOPT,2,3,0
0 `& Z/ V; q; o. t" } \) G7 [$ o! E: m5 f2 W
R,1,1E-6; M5 ~& K, N5 u# e
3 ~; M. Y$ f+ O8 o& L2 y MP,EX,1,2.01e5
7 l% d( i* C6 f" ?2 l6 U" I0 |5 g3 T, e. A
MP,PRXY,1,0.3/ [- S! J( c5 {
2 d. t# p) }' e: N+ Q CYLIND,15,10,0,200,0,360,
( A( Z3 G- m r8 [9 a! ^& W. h
8 r! t) @. [5 l3 l/ y- Q/ D wpro,,90,
/ V5 } f/ G5 S' ]! P j3 W7 R
/ U& `, _1 x; P" i- ~+ {* q8 a4 T vsbw,all
5 \/ }& M# I+ Q+ E# L* Z# A0 L c4 C" j6 h4 S0 Q9 K7 g
wpro,,,90
9 b, ^" P) ~) w# C
! Y' @/ M, g* Q/ E. O; v vsbw,all) _6 W1 e5 C0 ~7 a
( d3 M c+ E5 z7 h2 y- q* L n) L WPCSYS,-1,0
~+ g. o2 ~- }0 Z' a7 S# z; B# y7 @/ g4 M- @
K,17 , , ,210
1 [- z+ I; U" }# f! t* I+ ^- D! x! H" i# Q* e
lsel,s,,,13,16,1
1 F; r8 u" |0 }- K' ?0 |+ @
. K; l5 G: J' Y, i) g, D lesize,all, , ,8, , , , ,1$ O, ]- }3 _! {" E; ?; W4 Y+ t
) q4 E. m3 u: C* }& t6 Z4 t8 r lsel,s,,,22
% U B3 Q* u$ B5 O- @6 k
6 f. m# x3 E1 j! O, q2 z8 I lesize,all, , ,4, , , , ,1
9 n9 E6 e9 Q% V3 m/ }
3 u5 T, a5 V- m+ E9 [/ f lsel,s, , ,17,20,1
" e) }- Z8 Q8 c% ?$ r+ _8 C3 s+ W% p- Y, X' X3 K7 ^( z
lsel,a, , ,26,27,1+ Q, E- E u/ e; e
: R) w# Z9 G/ U# S. X lsel,a, , ,30,31,1
" p" q) A3 }6 K/ c& V6 }$ R4 u2 W" O% N$ d
lesize,all, , ,20,0.4, , , ,1" L, [- N; V( [+ Z( B& f) N5 C
- f) v" |" D, g2 m# }3 U8 _% \
alls
& @* C" s8 X" @8 N
3 T% F- k5 W# T6 e u5 E. m" j+ y: ` vmesh,all
- f2 Q& P+ w* J7 E/ G1 W
% g) c4 T, x) X' b( K# { !!!!!下面一段开始各个方法有所不同,由于前面的建模一样,后面的例子就不再给出
7 |5 x D# l) b% R4 z6 t* s3 ^) b7 z3 P+ L3 Q( i3 t/ y
ksel,s,,,17
) |, p2 R- ~$ s. E) |4 A$ M# M6 A5 p
type,29 e0 t0 H* y' {+ m
/ L, a$ E" |6 K, z, T1 W0 R real,1
7 n J3 l; {! B7 ?+ a0 X! V |5 ~+ b9 K( j6 |! z+ X; I: d2 U
kmesh,all" H& c$ n' ?# p% ?% m
* _- ~5 R7 ?+ x3 m3 ^9 R# A( b allsel
4 g2 a0 T4 d( f9 J6 t7 I- u0 P, I6 P9 i& ] d5 _
nsel,s,loc,z,200,210
c" C) Y. |' P- }! D5 o: y6 g" t
npolt0 M# \$ g V% f, a2 V
4 J/ p4 m# m( S% w
CERIG,node(0,0,210),ALL,ALL, , , ,
0 U8 ~/ U( f- R, {; c3 ?8 ~
6 \+ i/ }8 d& @8 b0 J0 j !!!!!CERIG命令定义局部刚性区域* V) `+ [& L1 w+ t* Y, T5 F. ], h
8 N# h* _- z4 [
allsel$ X D+ K( W: J- T
1 D/ e6 {( ?/ m7 {/ M& A
/SOLU
0 ^4 A2 ?" }6 n# p+ q# ~, H8 ]3 f3 R" B- q' E9 \
f,node(0,0,210),mz,10e5' a8 O: d# E8 h: A# k, D3 k* ~
' ]* j# R. s T9 w* P) @; ^& Y
FINISH
) Y, l7 d+ p) j: }& T9 m) f* g8 q
, |$ F+ h' G( X6 Z9 F/ E !!!!!以下一段边界条件的施加各种方法一样,后面例子也不再赘述
/ }7 ^- o; H, C5 ]( [, g3 k* L4 V D
/SOL% D% H3 c* E- W9 z, | g l- \
! m6 O( |1 t( g) z y nsel,s,loc,z,09 `; h: T3 ]8 Z+ w* w& K
, x5 {/ S- `; Q1 g0 U& E+ g$ b
d,all,all
! h# N) E! g8 E. \. [; {( W7 }) F6 E3 u7 h
allsel
: @! I/ H- p6 g: d) a( _
3 `- I7 ]% e' s( k' }& L) D solve3 `$ w# N+ M! e' P% D, n
9 ?" Q7 N0 \5 H 方法3:使用MPC184单元定义刚性梁
/ [3 V4 |. G- \4 ^$ ?# m4 B' j8 ]3 @# ^4 L& }. P0 i. @
……. D8 N# N; W8 L. A7 Q2 Y
7 m5 S2 z, ~8 l3 |. ^+ g
et,2,184' v" K" X1 X, {8 m( }+ w
% l! u+ b6 S% b$ m" m% F keyopt,2,1,1
! w3 ^0 h- V& _3 O& f8 l
5 F$ M' H4 d0 s) ~: t/ j. H nsel,s,loc,z,200' Q) c& Q; L) A6 h e/ `
4 c1 T* }% B. M, B
n,15000 ,0,0,2104 ~" o" l. n7 H
4 x% Q$ e1 h6 O/ m type,2
* ?# B. b. }8 s4 Y
1 g! r$ \* d$ v0 k2 T4 [1 ~0 s [ *get,nnum,node,0,count
' y5 a( Q/ v H; p8 a+ F# p. t9 h" y, Q
*get,ND,node,0,num,min& n# g" \/ P; z+ t/ w) t; G
7 k* L: A: G# F0 ]7 G1 L( P
*do,i,2,nnum0 `# B/ d( c7 N0 ]1 s3 K2 N9 N
/ G2 X0 B3 I9 @) g0 a !!!!节点个数是nnum,只需要生成nnum个mpc单元
. I; ?; o& [! H1 U$ M A
7 ^4 J/ Y# v A5 Z3 Q. {$ ]1 v E, 15000,ND3 o. I& X A9 M( H# f7 G% _
( P, d a, E& }0 a# K5 y
ND=NDNEXT(ND)
% s8 R$ c/ d$ x! v' _5 V. N+ ^5 q7 U3 A7 _4 D+ s
*enddo
, x. B! y c# B" R% I- r) }7 w8 S2 F2 i
allsel
5 f( a6 {5 j1 W2 C
0 @ D: V( F3 O$ Y/ Z /SOLU
" Y1 R8 G8 b, u8 Y# D3 K0 D3 @- i" J3 q! ?3 E4 `1 M; t& b
f,node(0,0,210),mz,10e55 v' z2 I# o- G1 V& j: `& E
' P8 q$ ?7 v% O. L4 Q" g FINISH
! Z( D3 o# w2 g0 o
w, s5 `8 g9 }6 }' w' @/ K ……- z& a2 q* i. n6 ]: b K9 L
( t, v, P* M9 a" |& p& B 方法4:rbe3命令
: _4 r) h& F5 U5 S' ?# D9 G7 Q. E" U q1 n5 s" ~, l+ e
……
- I6 P1 h# `$ u) G
5 F) z E+ |4 E/ { ET,2,21
5 Z# H; {+ i* x% D$ l/ g; V6 _3 z2 D+ M) U- S* I
KEYOPT,2,3,0! V3 J$ r/ {# d
4 \+ c$ T1 z. {8 B0 @
R,1,1E-6
: W, ~( J/ \) \' [: F8 F# o, L& C* q$ e
K,17 , , ,2109 m# n$ k: h0 r5 w6 h3 e0 F& G
- t- Z1 G# C+ r6 t, U3 [ ksel,s,,,17
* ~% q+ ]5 ^1 y8 T+ L+ | ^( G, N! ?* O& ~. F5 k" v% o" ]
type,2
7 l' z. t3 l6 D: e7 q# M7 _! @- Q, ^2 V r0 h
real,1
$ t! U8 B; z5 W/ m @+ u; T/ B+ q( g$ F; Z. G0 [& V+ z! ~" x/ g
kmesh,all
- A- {5 ?% A1 Y2 j% F$ C" _
+ v: Z9 }+ v( t allsel3 x! m# Z" p5 k- t8 r, p
6 O5 ]2 n, _) ` Y1 M3 |* f
nsel,s,loc,z,200
- E7 N4 S$ t$ x( c& b d1 \! L2 d+ U- B' [' {9 S4 Q0 Z
*get,nnum,node,0,count) t7 M9 Y& A5 N* O( ?
5 L4 e4 k0 d" r0 ` *get,ND,node,0,num,min, M7 o) S; \. r
3 d! n/ Q0 Q, F4 }) w *dim,sla,array,nnum5 x1 H$ g+ q# i( k; n. a0 V* ^
9 u( D$ T9 o& v6 d *dim,sla2,array,nnum
' B. o' f0 T+ J2 p! o+ r6 u
! r+ M& a) x) o2 y+ w, P8 ]% D *do,i,1,nnum! s+ y+ O; t- q3 B- {4 @4 U( _- @
6 [- J& ?% z6 A0 {9 J/ X
sla(i)=ND
: @6 z; S: e2 A: q+ ?8 P7 J% f# |; w8 Y0 j$ a! E6 c/ w. D: O
sla2(i)=ND
, ^8 ~9 J$ s# f; t8 H% y+ n3 `1 V5 \" E# |9 C+ w V7 q
ND=NDNEXT(ND)
! z) V( h" E6 n: b1 C% Q R1 ?; B% ]$ k/ `6 E! K9 q( @% {# k
*enddo
$ O/ U7 K! _5 A+ [8 t+ K% p1 J j) l) {. A8 I
allsel9 i8 p& X. r" p* y
3 M' }/ V/ `+ q( A
rbe3,node(0,0,210),all, sla,sla2
. W. [$ K- }* t- Y6 O- Z9 y
7 v* F9 ?- q0 x2 p /SOLU+ C- ?& Z: k g% q! k) J, C: P
X1 a0 ^6 p* c/ `+ F& D0 X
f,node(0,0,210),mz,10e5
3 M( K$ W1 s( }6 V/ [2 A
5 \% n7 z) B _! M3 C& `1 M1 q FINISH* T( ?: g8 n+ i* C9 Z j
1 {! v* e2 s' x6 @% ^3 q1 P
……
7 E- X, r6 \2 ^# T
+ b+ L# Q; J% G9 G3 i0 H 方法5:定义刚性接触面4 k8 s- y. f5 s( a. b# S) s
1 [0 w) @6 C" N6 W. @* j' c ……: d6 G5 E# V& K3 _$ r0 i0 g; O* Z
+ i8 j9 D0 k9 G
n,15000 ,0,0,200
. n1 m( s0 `5 Z! E( M1 b+ K* e; m$ r' s" @' }
MAT,1- b$ U2 \4 U( u ^4 H( a. ]
( u! P3 X( f! k; `% O R,3
& _- w8 G# N% O4 M, p+ g6 z5 c: I) z; z/ J% f9 [0 i
REAL,3' I5 m* V$ M0 g0 O
0 s$ ^4 }* @: I: ^ ET,2,170
5 X2 W% i3 ]/ C) ?+ Z* W/ v5 J& S$ d: [
ET,3,175
7 s& S4 L; j. L0 Y! K
* W" h% J* _' h KEYOPT,3,12,5* P: G5 l& j' b
. L+ w9 X& A* D6 m* Y# Z
KEYOPT,3,4,1/ r3 W" o4 Y y! v# [ o$ Z; ?
^( c g3 Z4 x2 e
KEYOPT,3,2,28 ~% }5 h' Y/ ~2 `9 v5 A- y( p: A
* b* a' z' c6 E R% {
KEYOPT,2,2,0
5 p$ D' S3 a f& u) k0 |7 s
& I# k5 {1 `! I3 q, L' X KEYOPT,2,4,111111
/ C+ B- n1 K& x
& F; H1 `7 D2 V0 S: h4 n3 A0 ^9 A TYPE,2
) [# p j2 h( V# Q7 ~* n% g/ t# K
! Create a pilot node! n+ A# e# @% J- d$ ~
n' X' D5 H X TSHAP,PILO
( H7 H- a+ L# F2 Y9 Z0 r E$ N
+ c9 j0 | c% T9 j0 q' H' Z E,150001 M" V2 _. a2 @5 h
1 d( O3 m/ H2 \( J4 y" B
! Generate the contact surface+ E# r5 A) X6 @# G3 w2 L5 P
+ ^" V0 I* W" l8 D* [0 D5 i: z
ASEL,S,,,14
( {) L1 E7 { v3 x: U2 |# Y/ o4 b/ c; t) F. }
ASEL,A,,,19
8 G9 R U. p5 }- R" H: B5 {. w D) h7 G" t
ASEL,A,,,24
N4 L& | b F6 B l
+ a; M B6 l% V9 _2 l6 `- F$ {$ A ASEL,A,,,28
* n1 U; a$ d% D/ u0 U9 n& u# |$ \: j6 A1 ` o/ D
CM,_CONTACT,AREA8 g$ L, Y/ q& ?. g5 R* i
9 L& Q7 b5 z! S8 U* J$ F
TYPE,3 x1 d9 J5 r9 M, X& h) [
4 i0 l. E9 h+ N% x! T( {$ r7 a
NSLA,S,1+ X" K2 I$ {+ Q, w6 I
- @0 u/ j5 o( V7 K ESLN,S,0* f: J/ e, F% x; L4 @3 W0 i# [3 E
1 d5 S5 l: }/ l9 `, m$ o
ESURF
0 ^! M& R, I8 W/ B$ ^
7 C& c) y5 f6 J5 B2 w ALLSEL0 z3 o& ?9 Q2 k2 N1 E) R
$ y+ u, b& ~1 R" ?
allsel3 ?# a. G6 t) Q: K) ~
& T+ X6 z p& r8 j$ o /SOLU2 w/ J' O# ^/ n" z6 B. G1 g
& f5 b/ N1 [7 M- \6 d
f,node(0,0,210),mz,10e57 u6 r/ g) Z8 ^3 |9 X
' i6 `- v3 r2 A1 u# i6 O FINISH
2 }1 Z5 C2 q' r6 F W4 J
6 w3 P# @+ J! w. O( x& j. }! o( I3 U …...$ t0 @- n. O* S$ t7 u
|
评分
-
查看全部评分
|