|
2#
楼主 |
发表于 2011-9-30 15:37:04
|
只看该作者
本帖最后由 小小实习生 于 2011-9-30 15:53 编辑 . }; X- E1 A0 f5 F
/ a' f2 T; g/ Q# p# d8 p( V
这个是经典的简支梁的分析,移动载荷,但我做出来的结果不是怎么的理想。
8 b& d9 [3 M4 v) }1 c. w+ `FINI. b+ R4 h) ?$ x, L
/CLE
2 h) y# D, x; m& m!设置参数
1 W9 m3 D( I; q/CONFIG,NRES,2000' U3 M% a; v! N: t! K b
*SET,LONG,40 !梁长度7 U9 F/ ?. }3 q }: W
*SET,NE,20 !单元数
* k2 b/ k U9 y! ?+ V9 D# `" x*SET,DL,LONG/NE !单元长度
& d/ g- M) y9 Y7 O0 z*SET,NN,NE+1 !节点数
% V8 e1 j# `# w+ b*SET,FF,5000 !移动的集中力
' M. K- S; j$ e5 d5 x*SET,VV,80 !移动速度
" I- g8 R( Z! G- J" u+ c- W*SET,DT,DL/VV*3.6 !移动经过一个单元所需时间! `: b* }0 i# x0 d
!定义单元材料实常数
+ w8 @! b' H: v/PREP7* i7 k x1 r0 B0 D
ET,1,BEAM3 !梁单元, [# p1 ~1 g6 z9 c: G) L) Q5 e
MP,EX,1,50E9 !弹性模量3 x4 u" `: R# @- W5 U0 ~" Q
MP,DENS,1,2500 !密度) r3 Z" a- e* P1 M' I# R
MP,NUXY,1,0.2
6 W! W x6 J) E0 O7 ~: a2 @R,1,0.5,0.041667,1 !几何常数
/ D) q) S0 W/ ~!建立模型
; j& L* [* l8 E& r; v. V' j4 r. g0 h*DO,I,1,NN,1
, f% T' e9 o' y' C& K3 z: AN,I,(I-1)*DL !建立
y/ S- {' x- }0 i*ENDDO8 F3 C5 K- c3 s, q: o, c
*DO,I,1,NE
3 Y( ]) B3 q8 ~" y; I' D" GE,I,I+1
& d4 T. I6 Y. p! r*ENDDO
6 ^+ h1 F, t. a; _; yFINI
) Q' g9 ^# M8 |0 a1 O# z9 x9 n!加载与求解
& E8 M$ {5 D2 v4 z/SOL& Q0 F% v& Q! E! I
D,1,ALL: `! u4 I# [9 T3 A) N- o
D,NN,ALL
8 v2 y1 S1 z( X0 t( f, X7 ]ANTYPE,4 !瞬态分析8 i3 V4 _$ q5 J
TRNOPT,FULL* s6 J: Z% x3 u2 x1 q& I6 j4 G) a
TIMINT,OFF% x: ~6 _; I9 d! {5 K Z4 ]5 A s
OUTRES,BASIC,ALL !输出基本项,每一步都输出' r: {8 E3 \! }! h; u
KBC,1 !指定载荷类型; }+ v# z9 O+ f+ V- {
TIME,1E-10 !指定载荷步结束时间2 J: q2 Z0 O3 |& R; I0 J4 `
NSUB,4,8 !指定载荷子步数为4,最大子步数为8
% B: U& t3 ^, D( V9 K# U# |6 oSSTIF,ON !打开应力刚化效应
) J1 w9 x; x5 H7 ?- O3 HACEL,0,9.8,0 !施加重力加速度4 v6 [( v1 m! O, G6 R" C
ALLS !选中所有元素
# D8 c" k$ Z% g" k3 O- n/ OSOLVE !求解第一个载荷步- a- a9 D& n( R" S; c" W4 S7 t
TIMINT,ON
. z, m+ E9 h; G9 N. G, Y) Y*DO,II,1,NN,1: k1 y5 E& `, ~. Z- @
OUTRES,ALL,ALL !输出所有加载子步的所有结果8 G' i I1 C% y
TIME,II*DT !每个荷载步的求解时间: n, \' b5 O) O: G- m% y6 |
NSUBST,5 !每个荷载步包含五个子步. f( f: v; h2 H5 d, U9 i# ?
FDELE,ALL,ALL !删除所有的集中力荷载
% q6 f& L4 ^9 _# W# R3 u6 q; UNSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点$ z* t. L/ X3 l# m7 Q; Z
F,ALL,FY,-FF !施加集中荷载6 p5 |- h+ C" V9 z( i
ALLSEL !选择所有/ ^( |! P+ g' \- w0 V% n
SOLVE8 c7 U& n; Q+ b ]) h
*ENDDO! y6 {" m$ o1 v. y
SAVE
, x& b( V# G( b5 m; `: ~2 kFINI2 E" ^: |: U9 ~0 I+ ^6 j: i6 \
/POST16 t9 E. \( a0 K8 S
SET,NE/2+1,LAST !设定读取的载荷步/ a/ w6 O* }4 {( w9 m
/ESHAPE,1 !显示梁形状
2 s- W" Z! Z1 b$ L" ePLNSOL,U,Y,0,1.0 !Y向位移
. x0 Y- E+ w9 t- V* C- ^PLNSOL,S,EQV,0,1.0 !等效应力) U% s" V( O8 j8 o, ?# p/ s
FINI
8 E9 H. _% Y; o, r/POST26( h+ G1 v/ b1 k
NSOL,2,NE/2+1,U,Y,UY_M !取出梁中部节点Y方向的位移, `$ u$ C. }9 F7 W$ \4 r
XVAR,1 !定义时间变量为坐标横轴: E0 U# J; P A! X$ ^. {/ F
PLVAR,2 !绘制位移随时间变化的图形- {% u) e q* x8 k
DERIV,3,2,1,,,,,1 !对梁中部节点的Y向位移进行微分运算,得到Y向速度$ N6 V* l3 x( @ z
PLVAR,3 !绘制速度随时间变化的图形$ @+ N4 {. j) j- [
DERIV,4,3,1,,,,,1 !对梁中部节点的Y向速度进行微分运算,得到Y向加速度
9 ?- g: F) m. C: Z" P8 z9 uPLVAR,4 !绘制加速度随时间变化的图形* L a* k [6 ~$ @4 y7 q
% u. A- z3 A, d/ t
|
-
4.JPG
(28.91 KB, 下载次数: 109)
这个是我做的一个悬臂梁中点的位移时间曲线。
-
526.jpg
(30.62 KB, 下载次数: 164)
这个是简支梁的中点位移时间曲线,哪位大侠指正下。
|