本帖最后由 小小实习生 于 2011-9-30 15:53 编辑
. Z |& n7 t5 ~" ^7 Q
( Y) L" o, Z# Z8 A& S# W3 @这个是经典的简支梁的分析，移动载荷，但我做出来的结果不是怎么的理想。
' G! j2 E+ h, B/ MFINI
8 D- h' y5 I; }3 {- U2 p0 j/CLE
9 Y1 `# ?# F) Q6 G: s!设置参数
, D. A9 t i( Z& O" Z a/CONFIG,NRES,2000
" I0 {5 V5 L2 b; b1 t; [*SET,LONG,40 !梁长度4 ^) y' x5 u7 H
*SET,NE,20 !单元数; U7 u# g! _3 ^. y; U) z F' k1 u3 ]
*SET,DL,LONG/NE !单元长度
* I& ` w: `6 ^0 c1 D1 k*SET,NN,NE+1 !节点数
2 D+ `- P& h. ]3 D*SET,FF,5000 !移动的集中力- a3 l9 S. J! y% h' z# L( L
*SET,VV,80 !移动速度1 h# I1 f; g% l" C
*SET,DT,DL/VV*3.6 !移动经过一个单元所需时间
F( g9 W. k- ?1 M/ Z2 w4 [6 i$ j!定义单元材料实常数
( `4 F- K0 L1 P: R @) u \/PREP7
; X0 N; `2 J9 j5 ?ET,1,BEAM3 !梁单元
- b4 @& D& D0 Y* B& N" Z" I1 _MP,EX,1,50E9 !弹性模量
6 s. O4 x$ l! HMP,DENS,1,2500 !密度+ e( ]; b2 _* E* `) w9 _
MP,NUXY,1,0.2! X" p9 `! U6 A
R,1,0.5,0.041667,1 !几何常数6 D6 f- K" i; m3 X
!建立模型
+ w! J1 i- o5 k$ x; {( J*DO,I,1,NN,1
0 u+ r) O) v+ B- L1 P5 s- J3 wN,I,(I-1)*DL !建立: W H1 {7 E( a! q* C
*ENDDO
# O! l8 }0 {6 P2 i7 Z' Y*DO,I,1,NE( B' T! Y9 ?/ a2 X9 K" A: ~
E,I,I+1* |; j. q( l6 n! S2 U& l: r
*ENDDO
% R; \0 z+ Y3 {6 V( F% i- } Y4 ZFINI
& m/ i/ E6 k7 h0 P/ Q7 _!加载与求解
1 d( K8 D( S+ Q! ~& S" H' x* c1 u/SOL
4 F0 m: A7 q% L0 v# V- e4 lD,1,ALL
4 |4 z4 q. g' O( W7 m A& b6 JD,NN,ALL' N# R9 G7 i' w/ W4 G
ANTYPE,4 !瞬态分析3 R( E0 J2 B3 E2 o" _) o
TRNOPT,FULL+ J! m! j" l, f8 }4 E' a9 s# t) @' b
TIMINT,OFF
6 R9 c( i) h9 {. V# T# yOUTRES,BASIC,ALL !输出基本项，每一步都输出/ k2 m) I# l! ]+ m$ E* n& H! E
KBC,1 !指定载荷类型* Q$ C0 x7 \6 G2 a& o; z# m
TIME,1E-10 !指定载荷步结束时间; c) K$ U* @+ P% \$ Y
NSUB,4,8 !指定载荷子步数为4，最大子步数为8
* a' H6 `8 `& t7 Z! S2 x3 J- y% @ fSSTIF,ON !打开应力刚化效应- E* ^& @; m/ J/ z% h
ACEL,0,9.8,0 !施加重力加速度5 o9 u* f- T9 e5 n; h
ALLS !选中所有元素
5 x. i5 }0 t: w" U. L+ DSOLVE !求解第一个载荷步
g0 n. Q0 H' p1 I' nTIMINT,ON$ \" E( s0 q7 A) T
*DO,II,1,NN,1. z, m! G# F8 L' `
OUTRES,ALL,ALL !输出所有加载子步的所有结果
! E5 r- m* T1 NTIME,II*DT !每个荷载步的求解时间
% X% G E2 Z$ t. FNSUBST,5 !每个荷载步包含五个子步) f6 ~+ Z+ g) N3 g2 M$ O
FDELE,ALL,ALL !删除所有的集中力荷载- e: i+ Y: b6 H7 V- O2 S
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
8 b7 X' e. B2 S0 rF,ALL,FY,-FF !施加集中荷载
" o; a; H- U5 I2 qALLSEL !选择所有/ M! v. C u7 y$ J
SOLVE
6 R7 D9 P* n' H, T1 L7 n1 i*ENDDO& ^' @7 j# m! `1 ^, [3 U6 T
SAVE
7 e8 m. a% O. H" j6 KFINI
0 c1 o3 m+ ]' g+ [/ t. E' g. `: Q/POST1! Q4 K$ _1 G7 M2 C0 Y9 T# P6 t% {
SET,NE/2+1,LAST !设定读取的载荷步
( D0 Q$ T" [4 C% i# O- O* k# i* y/ESHAPE,1 !显示梁形状
! L( m. z' s: v, ?PLNSOL,U,Y,0,1.0 !Y向位移2 ?# r( m1 ]8 F
PLNSOL,S,EQV,0,1.0 !等效应力
% ]0 H+ o( g+ W: c, x2 s: G1 iFINI/ `) k& l: T& H! [! q
/POST26
) ?1 z. _3 d4 e) W' }. J: CNSOL,2,NE/2+1,U,Y,UY_M !取出梁中部节点Y方向的位移
! n `; l5 _3 s3 I* TXVAR,1 !定义时间变量为坐标横轴
8 w0 P3 r- [, vPLVAR,2 !绘制位移随时间变化的图形/ o7 k) m$ C/ O! C* K
DERIV,3,2,1,,,,,1 !对梁中部节点的Y向位移进行微分运算，得到Y向速度
1 |/ g! a5 o; O) _# c, ~PLVAR,3 !绘制速度随时间变化的图形
/ p' |% w5 A) ]4 w3 m5 nDERIV,4,3,1,,,,,1 !对梁中部节点的Y向速度进行微分运算，得到Y向加速度4 [7 i, d5 I( O( x
PLVAR,4 !绘制加速度随时间变化的图形$ s, } e9 Y4 g1 J4 Q' k
4 l* j7 S2 n& ]; j T- w
|
发表于 2011-9-30 15:18:32
楼主
