|
2#
楼主 |
发表于 2011-9-30 15:37:04
|
只看该作者
本帖最后由 小小实习生 于 2011-9-30 15:53 编辑 - @( L0 V% ~7 n# [8 V
, q$ v0 E+ Y$ t1 V0 z8 \这个是经典的简支梁的分析,移动载荷,但我做出来的结果不是怎么的理想。
v( o) F. Q4 V. r5 RFINI. k: v$ ]" U6 P! `
/CLE4 C; ^. @1 U- S- t. C1 O
!设置参数
! g* n8 G- s; X, k4 J% R/CONFIG,NRES,2000, L) y+ P Q8 I6 U3 k7 G; _
*SET,LONG,40 !梁长度( }0 ^4 `, F2 N- p: E. ]
*SET,NE,20 !单元数
! | p) H% x# t# e! ^*SET,DL,LONG/NE !单元长度2 e: ^4 u; N# B/ d9 n% S& A0 `
*SET,NN,NE+1 !节点数
* X* Z$ L2 C4 M1 V7 Q& x7 Y% l*SET,FF,5000 !移动的集中力
: |4 z; a: n( B! n& T+ c3 \: n*SET,VV,80 !移动速度/ z, c, l$ v' I+ @6 r/ M& }6 X
*SET,DT,DL/VV*3.6 !移动经过一个单元所需时间- E7 x% [! Y0 d6 [, y% x# B2 g
!定义单元材料实常数
4 n* ?$ [% v& M' A8 C1 c/PREP73 _8 s$ j1 s% _ M
ET,1,BEAM3 !梁单元% _" h1 w" X9 {$ n- {
MP,EX,1,50E9 !弹性模量
$ a" S/ Q# U' y' o3 g' qMP,DENS,1,2500 !密度! x1 }% Q E4 E% I- z" f6 ^. t
MP,NUXY,1,0.2. B+ ?# Y% ]& L7 [( j E) q( m
R,1,0.5,0.041667,1 !几何常数$ F- r3 G( I) M1 g; ?3 ]: |6 e' f
!建立模型/ o- e' @+ c. C& K& b
*DO,I,1,NN,1
: W. T$ }- h1 R4 n! B+ Q, f( L# vN,I,(I-1)*DL !建立( [* ]: V2 D4 q; G g
*ENDDO
3 i. x4 j; ^0 @*DO,I,1,NE
: k0 w3 ~2 p' D) K4 I( v6 }E,I,I+11 H5 n, c/ _8 |1 \, N
*ENDDO1 V6 D) \2 |) C4 S. ^
FINI/ ~" m5 z6 J' E* x% b. e
!加载与求解
) g; E3 g+ v9 m6 k3 L+ o3 C/SOL1 n' i' f; J& a, G& v8 V
D,1,ALL& c8 z- n: Q M. G4 b0 d
D,NN,ALL
?( W3 h2 B' ?4 ]% k- cANTYPE,4 !瞬态分析9 ]( T8 j( |; |' j: d, Y8 `
TRNOPT,FULL
6 Z8 P. G6 K0 K* x% OTIMINT,OFF; r( u6 S( N# q6 J6 t
OUTRES,BASIC,ALL !输出基本项,每一步都输出
: |) `- v, s* H6 DKBC,1 !指定载荷类型
, f; c) S6 c$ K1 iTIME,1E-10 !指定载荷步结束时间
* x# F1 b6 @4 W2 E! w+ }. RNSUB,4,8 !指定载荷子步数为4,最大子步数为8; s3 [9 w6 ?; c! p4 g3 }2 y$ h/ W
SSTIF,ON !打开应力刚化效应
9 |) I+ s# |3 E& y: \8 y0 bACEL,0,9.8,0 !施加重力加速度 A2 h+ U; R1 N3 w2 \
ALLS !选中所有元素/ @1 Q& }1 {) @$ T. F2 Z- l
SOLVE !求解第一个载荷步
5 w+ H# p5 d4 M5 Z q2 ~; tTIMINT,ON& m; Y; E% o/ m3 g
*DO,II,1,NN,1
5 Q6 {4 o% d: a6 |7 ~6 SOUTRES,ALL,ALL !输出所有加载子步的所有结果- f) {; T$ l9 t' s: W
TIME,II*DT !每个荷载步的求解时间$ p: `8 _7 c* o% v: `
NSUBST,5 !每个荷载步包含五个子步
4 [8 \* ^# l v8 u8 nFDELE,ALL,ALL !删除所有的集中力荷载' v4 s$ m7 O5 S5 _8 {4 v
NSEL,S,LOC,X,(II-1)*DL !选择集中力作用节点
& a$ V0 J) N7 N! T0 V6 M% yF,ALL,FY,-FF !施加集中荷载
% [% Z2 y, ]0 Z7 {ALLSEL !选择所有
* B/ j, Y# A" W7 _- H* eSOLVE6 f$ h# a' k' m+ B* T/ b
*ENDDO
! S' o$ b. i' {4 k9 YSAVE) O. Q& t }' Y" W6 F
FINI
8 y$ o. n! B! m! T9 o" D/POST18 ?$ s& Y1 }, g' b
SET,NE/2+1,LAST !设定读取的载荷步
- @' k2 k/ D+ j: [/ESHAPE,1 !显示梁形状; e9 ^. G5 K% U& Z9 z
PLNSOL,U,Y,0,1.0 !Y向位移
$ t# p% l# g y$ zPLNSOL,S,EQV,0,1.0 !等效应力
% D* R+ g) s' H( s% dFINI; g" @3 R! U" A1 B4 @
/POST266 L u- m* ?+ t# q6 P' H
NSOL,2,NE/2+1,U,Y,UY_M !取出梁中部节点Y方向的位移! u' Q) Z7 D- {* T [) y
XVAR,1 !定义时间变量为坐标横轴5 i" C8 q F( ^* ^" S
PLVAR,2 !绘制位移随时间变化的图形+ t t/ }4 `0 Y- P+ R. q7 S" ?
DERIV,3,2,1,,,,,1 !对梁中部节点的Y向位移进行微分运算,得到Y向速度
- ?& a5 k) K2 x4 z# s# XPLVAR,3 !绘制速度随时间变化的图形9 @) C7 p9 v( z: d9 q
DERIV,4,3,1,,,,,1 !对梁中部节点的Y向速度进行微分运算,得到Y向加速度
. _" {! b* J7 u# E2 @6 MPLVAR,4 !绘制加速度随时间变化的图形2 ]* c; u* q) b/ Z3 `( F, S
: h! z, w# v% {- t$ ]9 ^: p
|
-
4.JPG
(28.91 KB, 下载次数: 106)
这个是我做的一个悬臂梁中点的位移时间曲线。
-
526.jpg
(30.62 KB, 下载次数: 161)
这个是简支梁的中点位移时间曲线,哪位大侠指正下。
|