|
发表于 2013-7-25 15:09:12
|
显示全部楼层
我邮箱1962940284@qq.com% ]- [# h' e- H. T- d4 S
施加弯矩扭矩的方法其实有很多种方法,在这里介绍其中的5种,并进行比较:/ M$ |( ^, b8 f8 v; c
! f9 ^" }6 m* p 1.将矩转换成一对的力偶,直接施加在对应的节点上面。: M4 b7 h. W% n$ _
/ v; U7 n& C; p/ K# g/ A6 q D 2.在构件中心部位建立一个节点,定义为mass21单元,然后跟其他受力节点耦合,形成刚性区域,就是用cerig命令。然后直接加转矩到主节点,即中心节点上面。
4 T6 Y& G2 {# Y6 E2 p
" Y2 w8 t/ H4 F* d H) b1 C- a( w 3.使用mpc184单元。是在构件中心部位建立一个节点,跟其他受力节点分别形成多根刚性梁,从而形成刚性面。最后也是直接加载荷到中心节点上面,通过刚性梁来传递载荷。0 K z+ s- M0 P0 a
% h0 ~7 R/ l) @ 4.通过rbe3命令。该方法与方法2很接近。
1 Z3 F7 }5 @$ l5 D
9 P. V6 K% @+ m 5.基于表面边界法:主要通过定义一个接触表面和一个目标节点接触来实现,弯矩荷载可以通过在目标节点上用“F”命令施加。, s/ P |! A! A# E/ J, z+ s
% x- s% X$ c, u0 M2 k 对于方法1,通过转换为集中力或均布力,比如施加扭矩,把端面节点改成柱坐标,然后等效为施加环向的节点力;而施加弯矩,可以将力矩转化为端面的剪切均布力;但这种方法比较容易出现应力集中现象;
) D2 d1 q9 [9 s1 }4 ~4 r8 |+ W2 D0 j& P" b$ r2 x6 V1 Q4 r
方法2,定义局部刚性区域,施加过程venture讲的很详细,这里就不在赘述。根据他的例子,我在下面给出了一段命令流。该方法有个不足,它在端面额外的增加了一定的刚度,只能适用于小变形分析。0 E4 N5 {1 l: H |; f1 g9 T
) t9 s2 m) F" |# o: h/ }0 e1 ]
方法3 ,相对方法2来说,采用刚性梁单元,适用范围更广一些,对于大应变分析也能很好的适用。但在小应变分析下,方法2和方法3没有什么区别。7 X7 x1 M" T( j
' I/ G8 E7 V- o$ x
方法4,定义一个主节点,施加了分布力面,应该说跟实际比较接近一点,但端面的结果好像不是很理想,结果有点偏大,在远离端面处的位置跟实际很符合。
a7 f- U- M! T6 x& F
" I2 X. ~& O0 ^8 o: n2 P7 v 方法5,它具体的受力形式有如下两种:
+ M& n5 |4 N2 n9 d: w' P6 i( Q/ H" u: _5 h% _
刚性表面边界(Rigid surface constraint)-认为接触面是刚性的,没有变形,和通过节点耦合命令CERIG比较相似;
8 @$ E- t! Q8 j
5 \3 g! j k) J5 F# U/ a, I 分布力边界(Force-distributed constraint)-允许接触面的变形,和边界定义命令RBE3相似。
7 A5 W4 l& C) G$ i' `0 Z' S% D; t/ E) d% j
使用这种方法,需要用KEYOPT(2) = 2打开接触单元的MPC(多点接触边界)算法,- P6 j" o8 ]4 b' A3 B0 H- e
/ l" a* U( K i, B% _* H; A
下面针对venture给出的例题,用不同的方法来实现的命令流。) I, P3 C2 k6 [7 w. r
" q( e! c0 E( A. C. X
方法1不介绍了,方法2:
" M! `1 K' @* [9 W$ u8 H/ [4 P) V e
/PREP7
2 w& F, G" Q' N! R0 Y) Z X+ x/ _# P1 R
ET,1,95
' d' C% @# x( ^
: l2 c8 j* _! B3 q6 b% q ET,2,21
3 h% W/ K# X$ d! ?6 q% J
: P* K1 s; G: ^( W KEYOPT,2,3,0
: F7 J( _6 Q! l0 u8 x: D a9 T9 w2 L- r+ ], c2 E- E
R,1,1E-6
7 T# c, ]5 M' J, j S) J
1 L2 C4 P4 i9 m2 i, d) n MP,EX,1,2.01e5
/ {1 @8 P; c/ T( b" ^9 F7 V! |* m0 [5 ]& B# t& o {' _3 q* T
MP,PRXY,1,0.36 E" o9 I1 |# G: a
% S/ A' v2 H1 V0 h4 Z) d9 d9 m CYLIND,15,10,0,200,0,360,
" \: E* [- K% _ K' p! B3 K5 _, f6 Z3 r" @
wpro,,90,2 W$ B+ E" I5 ^' K3 M
5 Z A: R3 a a; @$ {9 Y
vsbw,all
q, q# ~% |1 O4 d
. p: }6 n/ n+ b! l6 E. K wpro,,,90
: w. h- u& B! _# e4 W
6 m& \, D/ s: s! z vsbw,all
$ t( B3 U3 Q% E8 j: z6 x! k& o ]# v8 H/ }1 o5 J% a
WPCSYS,-1,0
$ ~" M% z7 N+ x7 T) s9 M5 i! z% C
; b& z) A7 J2 b0 v K,17 , , ,2108 e( P1 l) C( P' ~5 h: Z
- V! |! G+ Z* A9 h lsel,s,,,13,16,1
' x; o' p! D' s1 ]) _6 U/ s4 @- S! Y' N! e
lesize,all, , ,8, , , , ,1/ ]& P9 m& \2 |
+ p i! e/ t; A) ^ q
lsel,s,,,22
" F% s R5 ]/ N/ O
* ?5 \# }0 o+ W lesize,all, , ,4, , , , ,12 ^! m! F# p# X. _% j
8 i$ x0 B: M( k+ x/ y& m, }7 l
lsel,s, , ,17,20,1) o, L4 e* p+ m3 h
! \9 G4 |. [4 j3 I& q+ t* L lsel,a, , ,26,27,1
% W) |" W0 N1 g+ ]7 L( o4 v; }1 j5 N& e# S
lsel,a, , ,30,31,1
, K) m u$ z1 L8 w, P$ t( i
7 o3 T2 q8 u, J+ a lesize,all, , ,20,0.4, , , ,1
1 n1 o: R( V1 G9 u- \( Y. Q
$ ]$ Y, M/ A) F% n* p2 q alls, N+ g- [/ M5 [
4 @2 B* A8 |/ W, ]/ y vmesh,all3 Q3 I: t3 [4 k- i" s
7 V6 T9 Q) }+ `& I !!!!!下面一段开始各个方法有所不同,由于前面的建模一样,后面的例子就不再给出& O: ^% L9 S& k0 V6 Z+ b' @ y$ |' w
6 B- u2 o7 i1 g7 f% a- J6 f7 }* o ksel,s,,,17' a/ r; k$ H: N; R6 U
) N5 ~4 d5 q# O a8 d' ^9 j P: g type,2
# k% s+ ?' P2 z* s: y" P% N z: l9 T& h, s: ?$ J, B
real,1: l' m4 N9 b/ X4 G9 s9 w7 q, w
/ S0 O' M8 V N
kmesh,all/ @! m ~4 X- G! ^% X. c* K
5 y O5 r2 h; g" Z% q& v4 o& g, M6 z allsel6 C- T% a' R+ J0 r, Y+ q
9 P U+ B* @& c \8 [: a
nsel,s,loc,z,200,210
% P( L( n; O/ ?! E* C2 b% Y) Y J$ h" ]+ [7 \
npolt
+ A! K) p' K9 f, f6 U y
5 r- t% G% }+ k, i CERIG,node(0,0,210),ALL,ALL, , , ,5 @$ K6 e( w/ _; Y
5 Y- Q5 H3 L8 |. E0 O6 m) E. @ !!!!!CERIG命令定义局部刚性区域- O( q9 _) B; w. E
6 B+ p. Q$ p9 Z% b9 i
allsel& [# b/ Z- U3 d g9 k: O0 O
& B7 k* _1 `/ J L
/SOLU8 X3 M. k, a, n$ |/ x# W0 t2 x
5 B' \* x' y5 C( K
f,node(0,0,210),mz,10e57 u0 c; C/ h4 l
5 O; @( m4 h* E+ M8 P9 c FINISH
# T V/ o! h; X% m
- Z: y4 H$ s* U: j !!!!!以下一段边界条件的施加各种方法一样,后面例子也不再赘述
) c- I3 l* v! _ R0 I0 c" p* `7 u( [& Z7 _+ m9 B( j2 \0 I
/SOL! j( q! ~7 _+ d) b" K7 L' J, O
1 w$ l2 t7 K' w9 f$ ^& n nsel,s,loc,z,0
6 K& v2 X6 z; `* |7 U1 e- A1 p1 b: A2 R' F0 _
d,all,all
) g- H" s# C5 z' T/ v, T
. I+ E; e8 l( e6 T' e. w allsel
c5 H' L0 G# E" Y/ Z+ l" U3 z' W: T0 s: z v' @
solve; s- W& h3 }% `% z! s
8 j+ Q0 e3 c2 ]# i3 V, p
方法3:使用MPC184单元定义刚性梁- ^- U P8 J+ o
1 q/ W: t0 C$ ~% G
……
, l) _( B4 Z; P" t4 Z9 d% ~5 q: Y. F+ c# ?) ?9 ^
et,2,184
) f* q' Q5 P9 ~# Q2 z* g* {
; O4 v2 k3 z9 t8 d: S keyopt,2,1,1
, l. \$ j5 w/ O( f
5 k9 z# x8 ~! ] `' x; K nsel,s,loc,z,200- R. ?. Y" o- r V6 j1 w0 X5 o
3 N9 F2 Z; N4 o: g n,15000 ,0,0,210
/ {7 U0 i0 @/ U" d7 B$ v" o) b6 q% h6 m7 q" C0 n
type,2$ G5 Z- o' j% F- K/ f7 `
1 l# U' }, f' z- d# z/ B0 j) O( m1 { *get,nnum,node,0,count
0 X0 A0 k [ Q! M3 |+ N+ H$ }+ Q- L$ b- n% Y" A: I
*get,ND,node,0,num,min
( q8 k" ?. _( X5 c! O" x/ Z! f. D, l8 a8 o
*do,i,2,nnum6 N) X) l2 \' \
" t* j& |5 ?2 R5 p' b6 w !!!!节点个数是nnum,只需要生成nnum个mpc单元+ e6 @) U$ `* |! A
1 \: u3 X f3 @: K! z E, 15000,ND9 r: z; u7 s: [" A( P4 F# P
8 b/ m" ~. O: u
ND=NDNEXT(ND)
3 o/ E+ n- u$ A0 Z
z1 a# X) p% {0 f. b6 U1 ?. r *enddo
7 w4 r/ B: W$ A+ |
+ h" ^' y, }! T s- p- w1 j allsel
# R i6 p7 B% M; B6 [5 x; z* z# a0 Y+ `
/SOLU$ }& W% e: L0 N* c' Q0 u
( }; F5 V D/ i1 x2 L3 q f,node(0,0,210),mz,10e5
$ N2 O6 ?0 G; ?% @' |) z/ e( }0 I/ v; |; ^
FINISH
) S# Y6 ?9 c5 V0 R* j8 p0 b8 I9 g; l; D
……
$ X( X: H& n) M; _% l% N% f$ H+ D0 t; D$ _$ U
方法4:rbe3命令1 d4 o7 b, _& l- v- k
8 m! G4 P) d9 a, }$ H$ J# B& N ……
4 {2 R8 C2 \7 L+ n: e
) Q2 J$ s+ d& J- Z ET,2,21
! v$ L t! ^1 ~) ~" l T7 }) F H7 u0 ~: v6 ]* F1 w( d- N( C F
KEYOPT,2,3,0
0 n4 ]3 U( \1 L* N
. g b1 ^+ g8 c R,1,1E-6
# ^6 U4 ^" v' U8 Y4 t0 U5 L1 r# l+ A. R/ [1 t. y
K,17 , , ,210
, G2 y A- D" P: [: y
" { _) P, \. U9 Z/ O# C7 ~ ksel,s,,,17
$ y+ p' x# K0 {- H; e/ H- n0 P! M: L: r& e+ p; ~
type,2
+ l. c9 ~7 d. Q4 Y
3 y+ L0 c% g. h; |( d$ Z9 r5 B real,1
% L) f8 x9 {9 N9 b4 ^' E k6 j# r% H# e z- A U, R8 b
kmesh,all" u, }% x( I3 z% x
! G8 A5 z$ D4 U& z. [# T allsel& m) f2 P+ V4 S9 M
; w' @/ J2 A: h# G nsel,s,loc,z,200
6 ^$ w/ ^2 e. U1 n
4 [" B, E1 U3 @9 b2 |- P( \ *get,nnum,node,0,count/ N! g7 h' _6 }: Z
+ y$ o; ]/ m A9 m- i% o
*get,ND,node,0,num,min: Y9 U8 ?: P) W( X- s; O
u' L' z) }0 [- u8 u0 W
*dim,sla,array,nnum
( W9 U+ Y3 a& ~" ?& [0 B' n
8 V, N3 c. ~7 `$ t6 z/ q *dim,sla2,array,nnum
: D/ e9 a: ?- {0 j9 s; n( B% E; O) b& ]3 ^; o6 W
*do,i,1,nnum- v! U; q( B5 y; b3 q0 i
9 n% ?0 B2 S' X0 m9 B+ P
sla(i)=ND" x! y7 K3 r9 Y: ?) }
+ J1 A M; _6 R! X sla2(i)=ND
: I: d) I+ V" r1 N1 A6 h5 p, x2 u) |2 B: F4 M6 F7 p' a
ND=NDNEXT(ND)
$ y$ K% t4 H, {* e& J
9 k% ?5 N# @, B8 l/ a *enddo
5 b. J- ]/ l# Q( [. f
" W' ]; W: z4 R0 D$ ^. |+ A4 j allsel- ]# y `( ]5 w2 ~9 C) J
0 v7 k: U0 Z8 T& v
rbe3,node(0,0,210),all, sla,sla2
2 E' u5 }! a( H
, k/ z9 I# V# [& v# U8 Y+ L: t /SOLU
: I. S$ w T2 B7 k+ t7 \8 F1 p/ X: E: w1 U4 ^6 b6 I
f,node(0,0,210),mz,10e5
5 R1 D: S" K& C! j2 R/ O1 v( m( J& |, v+ T+ r( x8 F Z
FINISH
& `6 u1 U6 r; o2 d1 W l, E* r7 }6 V
9 Q+ ?. A1 i J5 J ……
* R/ C( a0 C% t. J8 z+ P6 w$ I. B/ n
方法5:定义刚性接触面' ?* \, z8 I9 K9 e& x% G; i, L d7 @8 k
# L( Q, u O9 \) X/ w8 W4 W3 ]% ~
……
- m8 b( C6 P) m5 `) W/ q5 g, _/ O2 u1 a+ x* n% a( u# @
n,15000 ,0,0,200, q* _/ ?. h, g: @4 |5 w
8 p( H2 H2 e( I- @5 j2 j* @ MAT,1, i$ [; Y9 a; X- D1 S( ]# b# h; ~
2 F0 j) P8 }9 q- c' Y! \ R,3: R3 t; F6 g0 f4 H6 t' b) g0 F
4 h+ l d- D2 ?
REAL,3. t9 a! _* t& a/ N# t
6 }- @* }1 U @1 @' `7 X3 l2 o ET,2,170! K0 r! F- S* f O$ Z
, f( A8 b9 i" I$ Q( I) C4 g | d
ET,3,175
$ a% \$ y. i! y. Q" f! i7 s- \. T A- A( X( G
KEYOPT,3,12,5
D+ A) E9 c% ]/ u1 ]% \" f- c
) U3 m2 n: U- ]3 K0 T/ G7 M ]) v, G KEYOPT,3,4,1
9 K* r/ A2 H5 R: W9 _" @# S+ B# g, ~
KEYOPT,3,2,22 T% b0 s2 `( k2 \
2 B: }) B) t7 O$ Q* Z: T
KEYOPT,2,2,0
- b; B) u$ ?6 l% C8 F7 ^( {- o$ }7 i8 i; t
KEYOPT,2,4,111111) M4 H+ k6 v. ^* L5 b o
" K+ Q0 M% n0 W% Z2 z8 \6 V
TYPE,28 u6 B) a, R8 L5 | A
6 L6 \; ^9 w7 s ! Create a pilot node
/ h- I0 |: |* ]+ F8 B: e/ O
! l8 n. l' D3 b ~3 w6 \2 b TSHAP,PILO$ h3 h+ _( _; L& h q' H7 K! g
, Y5 j @- [ x% t# D/ o. t. X8 ]
E,150005 ?. P r3 v- h2 H
& c, Y/ k2 I1 d/ F
! Generate the contact surface& z& L$ w! J0 l+ b
0 M$ l" Y6 G' Z' ]% a$ Q6 ^9 P \
ASEL,S,,,14
2 u% z. E4 J+ e- y) k8 b x
1 k! G9 e7 \- e7 K) F ASEL,A,,,19
* q+ Q7 a/ _2 I# A, S( K2 }6 e, L! D/ S7 ^# H
ASEL,A,,,248 R4 z+ C( w: Z- T
! G$ {) |7 e+ J ASEL,A,,,28% N9 U; `* q) C$ S9 K& t5 f
& u! n0 l8 A7 e CM,_CONTACT,AREA5 L, u; p4 \1 i8 i' l5 ?
2 [5 u2 W2 W2 i J4 w TYPE,3
$ S' {1 y2 [; a- C1 ?# W/ b
( g7 T2 X( F' c NSLA,S,1
* A8 W+ G; `2 i2 M" f7 x4 w+ X$ [% i2 ^/ G/ d7 Q
ESLN,S,0
8 t" h5 a) g: U0 }
' B! `/ `( P& k! t7 e ESURF
7 ]: g; [5 o; y3 C2 e' b' {/ u0 M1 F4 ]6 Y. ~; ?: m
ALLSEL9 R4 m8 v+ r+ p |; |% f% @7 [
- J. ~( O& S0 Z2 q) i7 ]+ n9 \ allsel0 e9 V9 d( V: v) o
3 R% n( u1 `! \
/SOLU
2 Y* C0 f. R* {2 b4 s
) o4 H, { U6 W$ j- k- ^: E f,node(0,0,210),mz,10e5
. ~0 u7 L4 I! y! S. I3 A$ h1 @) X5 c5 ^& `) e
FINISH
w8 L# E. C, g. c3 u9 b) R m
! b. Z/ l3 Y& l P …...
) }6 q7 `7 y# i2 @8 L9 u |
评分
-
查看全部评分
|