本帖最后由 小小实习生 于 2011-9-30 15:53 编辑 7 {' C+ Y) B# J! t+ W
3 t0 X& Z: j i这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。. H) c) w, Z Q) B! b
FINI( g' G: y* \, l( L
/CLE" _. @( q n+ J3 L
!设置参数3 I: R1 P$ R* _; h" e R- e
/CONFIG,NRES,20006 Y5 d& h D" \) a, @
*SET,LONG,40 !梁长度
# G4 |8 i4 o& k5 r( J" n' [*SET,NE,20 !单元数
7 V7 ?% i& d* b) p; l*SET,DL,LONG/NE !单元长度
, D. r$ T$ M) v R6 F, o*SET,NN,NE+1 !节点数2 L; Y) z: {9 e) U7 D& R
*SET,FF,5000 !移动的集中力
* C: e+ v2 W9 Z& K4 P* e/ X*SET,VV,80 !移动速度. Z6 F8 @$ o( ] A& g/ P* k$ f0 Z: n
*SET,DT,DL/VV*3.6 !移动经过一个单元所需时间3 N& {& z) u# S6 o1 X
!定义单元材料实常数. j' B2 j8 ~& Z; ^5 P+ n
/PREP7
! G5 C* m Z' V6 ?* DET,1,BEAM3 !梁单元
4 s0 \/ ?! t4 b0 PMP,EX,1,50E9 !弹性模量# {8 q e: X7 L5 [( X) Y
MP,DENS,1,2500 !密度7 r5 y2 e9 C* S* `* A
MP,NUXY,1,0.2
' u* l3 s) W1 B2 q6 HR,1,0.5,0.041667,1 !几何常数
" H! y1 ^3 ?+ b4 K+ @ K4 m! h9 ^!建立模型
. H" A( ^% L# I, f*DO,I,1,NN,18 H: \6 D" g) N4 n) Y
N,I,(I-1)*DL !建立. W% ]6 J$ x- P: Y! H _
*ENDDO- e" C5 ~0 {, d: o2 J
*DO,I,1,NE5 W4 G: m/ X U, X
E,I,I+1
) T5 e0 M4 `9 n0 Y*ENDDO
8 ?8 r1 y6 \& |FINI
7 u9 y; I+ `* d$ {! S0 s1 o!加载与求解! k+ k& v2 r# i h; U
/SOL; e3 r: p; ?3 ]; w! L l
D,1,ALL
5 Z* S/ K" \- S/ f5 WD,NN,ALL1 g2 Z0 E! i$ Q) P8 L0 q
ANTYPE,4 !瞬态分析' ^8 @$ [6 S. u8 |; B# {; S
TRNOPT,FULL/ n, p* s5 O1 N5 }0 J; c9 Z" F n
TIMINT,OFF
. J% I8 b2 e2 k( ?- TOUTRES,BASIC,ALL !输出基本项，每一步都输出2 D; o5 C5 |& F4 ` a J
KBC,1 !指定载荷类型% R4 V$ D6 _7 |% L* c+ B
TIME,1E-10 !指定载荷步结束时间8 H9 G; }& {$ f- {# b3 Y% |3 B
NSUB,4,8 !指定载荷子步数为4，最大子步数为8
+ D+ K. q0 `% g/ b8 l Z* t. {' SSSTIF,ON !打开应力刚化效应4 [; R; v$ _3 f$ t( _
ACEL,0,9.8,0 !施加重力加速度
W; X5 G/ g4 i0 {5 s4 PALLS !选中所有元素/ l) `6 b* S' M5 u
SOLVE !求解第一个载荷步: @- f" E3 M0 a
TIMINT,ON! t: c) s7 D3 e, F) Z9 [8 K
*DO,II,1,NN,1
/ z5 d! S- @0 ?% NOUTRES,ALL,ALL !输出所有加载子步的所有结果
% e) u2 P6 d: _7 j7 }- y& g; M" MTIME,II*DT !每个荷载步的求解时间
! z& K, k+ b3 D2 \' ENSUBST,5 !每个荷载步包含五个子步$ \$ r9 `/ J* {
FDELE,ALL,ALL !删除所有的集中力荷载
7 M- J. G& P: E* hNSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
: H& T5 F _( V9 f4 S! x! uF,ALL,FY,-FF !施加集中荷载+ r G( N6 C$ q9 I8 G, W$ e
ALLSEL !选择所有
: A) ~7 P8 t8 `1 HSOLVE" o' I7 x8 e( v
*ENDDO: g3 [# R8 s# [6 T- Q$ v
SAVE
! M! k, C6 a8 U9 k( s& @, v1 rFINI
; |4 W9 v1 d4 L/POST1
f6 i& @5 y, LSET,NE/2+1,LAST !设定读取的载荷步4 q7 U) u4 u8 J, X& Z% _
/ESHAPE,1 !显示梁形状
0 T# D" x# Y- j& k3 \PLNSOL,U,Y,0,1.0 !Y向位移3 u7 _8 `' @5 K$ Y4 R
PLNSOL,S,EQV,0,1.0 !等效应力
7 d+ s4 _0 N* W: ?FINI
* e7 D+ M$ @6 n/POST26
4 ^1 ~# s: ?$ {$ y$ J' h: _1 I) z9 rNSOL,2,NE/2+1,U,Y,UY_M !取出梁中部节点Y方向的位移
' a- l# ?; ^" ]" S2 hXVAR,1 !定义时间变量为坐标横轴; i/ m. [% l% l- X
PLVAR,2 !绘制位移随时间变化的图形1 P' V* } Y) z/ i) p
DERIV,3,2,1,,,,,1 !对梁中部节点的Y向位移进行微分运算，得到Y向速度2 t, q& F, z. w8 J1 ?- H
PLVAR,3 !绘制速度随时间变化的图形
4 X+ p/ v# i2 [. u1 ]DERIV,4,3,1,,,,,1 !对梁中部节点的Y向速度进行微分运算，得到Y向加速度
4 D2 y2 R$ W/ s @/ G. JPLVAR,4 !绘制加速度随时间变化的图形# N# E5 U! v$ ^2 W) @
( o- R, z2 ? E$ I1 J( T. v
|
发表于 2011-9-30 15:18:32
