ansys模拟碳纤维布加固钢筋混凝土梁问题

shihonglei 发表于 2006-12-8 21:01:52 | 显示全部楼层 | 阅读模式
我最近用ansys模拟碳纤维布加固钢筋混凝土梁,不过好像有点问题,计算结束以后,梁好像没有变形呀!挠度始终是零,帮忙看看!谢谢了!
* q8 a8 M1 F4 l!碳纤维布加固钢筋混凝土梁有限元分析
3 s% K: B; P- {) x3 g3 {/com,structural7 \/ F5 ^) K% X0 o4 P: ^+ w: u
/prep7
$ H3 s4 M' e- m8 H$ O!定义单元类型
& x# m, l! Y* w# W8 U- fet,1,link82 E* L. B# W% x" i5 x
et,2,solid65
$ |% V0 o* D  a# Ret,3,plane42
. \) T- b: P! _3 jet,4,shell41
1 @5 K6 R" H, J
* ?( ^3 H5 X& j5 W$ N. n!定义实常数
1 a$ Z) H# t' ~8 g  T! Lr,1,28.3   !箍筋和受压钢筋
, ]# K- z8 ~- z( ?, Sr,2,113.1    !受拉钢筋) i( J6 V- }! `9 d
r,3         !混凝土3 I1 R& _/ l" Y8 U5 B, n
r,4,0.12   !碳纤维布
' t* A4 s' i5 ]. a6 g6 Y6 E/ N% w, f& ?! D' C  u3 p
!定义材料属性
$ E4 j0 D4 w8 A1 y# Dmp,ex,1,2e5    !箍筋和受压钢筋
- A9 \' S' K/ ^6 z# xmp,prxy,1,0.256 b! }' Y9 z, x) T$ q8 O
tb,bkin,1,1,2,1, N0 n4 K* @) i6 n7 l$ Q9 o
tbdata,,200
1 E8 F; g: U1 o0 A+ _/ s1 V
+ \4 `* g; h$ L/ \3 Vmp,ex,2,2e5    !受拉钢筋
$ V( \. D0 ?9 j. p" B, N" M# ymp,prxy,2,0.3
# V& b2 C- C/ ?7 j' H4 q6 {tb,bkin,2,1,2,1, w1 T1 P3 l* K1 v4 V* B- |1 M9 O
tbdata,,3000 c5 N* V- j/ x: n! A6 c, Q: P

( @; K" N) _$ k( Q+ R: b% j, ump,ex,3,30000   !混凝土
, R/ c& z1 c8 o- |# _mp,prxy,3,0.2
+ Z/ ~( B( D6 @tb,conc,3,1,9
1 U8 K. p; t& ^* U- Htbdata,,0.4,1,2.04,-1
( @* e; g5 B+ h6 C6 ^1 A9 }* }
# f" J7 s+ G" O) G7 P: i" Amp,ex,4,225000   !碳纤维布
$ e4 f* p% M6 |% H1 k4 @mp,prxy,4,0$ ^5 @: n0 B3 z' N0 P
tb,bkin,4,1,2,1
/ I* E) i" t, F5 U1 }tbdata,,4300
; t* h& x7 G( {7 L# x. r* U; p8 T; z4 d/ W
/pnum,node,10 J- l3 n+ b0 [2 F
/pnum,elem,1. f/ D& u( Q8 A& [0 q0 Q
8 I3 U# |3 F$ l+ Y0 r
!产生所有的节点
( U# v4 N9 a, n" n( J0 t- Tn,16 M' @8 c0 m  q( N
n,6,100/ h# E$ S+ f' r: S* i1 T
fill,1,66 E8 L6 L- |! g$ L; H6 s
ngen,9,6,1,6,1,,25) s+ Q1 o: @8 n6 e9 L  q' O( x
ngen,33,1000,1,54,1,,,-25: Y  R& ]2 l6 W3 \3 a$ g2 t
/view,1,1,1,1
( w" J8 Q. K+ Y& Wnplot2 W0 ]0 j# u3 T! u" r* o1 A
) o- q4 D/ d6 W4 d: r# R6 l
!箍筋和受压钢筋单元建立9 [# C( n" r8 K! J) @
type,10 m  W( f; i, @5 {4 B6 V; ?0 e
real,15 _  P. L% t$ p  c  q" f3 f
mat,14 N- o1 f# q) s4 w) o$ u5 W
*do,ii,8,10,1
7 A" w$ J6 X. b5 J  I   e,ii,ii+11 T/ r9 X  W8 {4 A3 x7 r( ~
*enddo; O  s3 H0 m  b. J' {3 {
*do,ii,11,41,60 B5 ?; ^2 p& X5 R
   e,ii,ii+6
& V; E/ ]+ t2 S*enddo  y( L8 F% P8 D" K/ G# v; D8 t2 j
*do,ii,44,46,1
9 K  p& y) ]1 A7 b1 P# ~( I   e,ii,ii+1
) J# H! c9 D% Z" ]& V*enddo  `8 ]; k& [0 q4 x* q
*do,ii,8,38,6
- L; ]  H8 B# V$ }( S   e,ii,ii+6
" L) |6 |! Z  d# l) Q0 O' `) K- J*enddo: Q# b' }: S3 J, \1 g
egen,6,4000,1,18,1
! @  A( V: w4 J!受压钢筋0 y0 s1 V" w8 s( R% ^- s6 O
*do,ii,44,19044,1000  z+ E" ^0 R+ N4 A$ `" y0 H9 b
    e,ii,ii+1000. r9 o" O: n3 t! l
*enddo2 W6 u6 ^/ F/ ^, a
*do,ii,47,19047,10000 t' N: W) d$ `# J
     e,ii,ii+10001 X. ?7 q+ v8 \" f" ~$ v
*enddo
* Q, @$ ]5 g) C( ^
; |$ q2 E7 g0 a/ {$ C!纵向受拉钢筋单元建立. b1 X" J' z0 w0 C, j' ?' j' C+ L
type,1& @3 z. c2 [5 H. Z/ \" H/ v0 ~
real,2
1 n3 `: r0 F9 I# [% `, J0 x; Smat,2/ _7 S& d2 X1 R# m6 w5 X7 Q
*do,ii,8,31008,1000  D1 K8 y& O) [! S7 x4 B
    e,ii,ii+1000
5 U3 u4 \8 K6 t  q/ e- l*enddo* @' a' S- p  a$ E8 B
*do,ii,11,31011,1000
  X) L! f' B$ o( x  }) V    e,ii,ii+1000+ |4 Z# G# j1 F% a6 x' q4 m
*enddo
9 D5 o8 C$ r# ^) F& v) S0 H8 ?!碳纤维布单元的建立
$ z- k' y0 w5 i" qtype,4
6 O# C8 s( z* l  w1 x3 g( Z+ M0 T! Y& Qreal,4' r/ n) \" l& o
mat,42 h, r. S# k- [: s
*do,ii,1,5,1  E+ D% V9 R0 g
    *do,jj,ii,31000+ii,1000/ f: U: [: m: V: X# t. D
        e,jj,jj+1,jj+10001 ?( c) D+ t9 ?. g% C, Z7 Z
        e,jj+1,jj+1000,jj+1000+1
' K: C- E7 M6 O% Y. e- E    *enddo
/ M& K& c+ A2 t* o9 ]/ S*enddo
; i9 ~" W: j* @( l
" j- s. M, W" m' P" T6 Y/view,1,1,1,1
8 c& Q2 |0 h/ l4 x- K1 x/pnum,elem,0
% D+ E* K4 f, B5 m7 U/pnum,node,0
0 I; Z4 _! V; }, D7 r0 x4 @/eshape,12 q7 C/ J) l8 F1 ?9 c
eplot
- K7 N& L- R# j+ J" z3 I" f& y# [5 M) h- D# c
!混凝土单元的建立
2 k4 h' O* u; F" Zk,1
6 _+ L! f: Q- Fk,2,100
: h* E1 w: ?( mk,3,100,200
9 x7 y; o: Q! T2 A1 n/ u4 H; Yk,4,,2000 {2 B6 V! g9 _6 y  |
a,1,2,3,42 f8 O+ z5 P- D; ^+ T
8 t4 w  ^# k& ^* U- A  T
lsel,s,loc,y,0
8 E8 c! |, k# e4 _6 `- A) r# Olsel,a,loc,y,200
- |6 k* G" B: v& x) Plesize,all,,,5
' p; l/ L/ d1 ~: ?/ M# u% B' Alsel,all8 P* f$ C1 W) `& w/ S' J
lsel,s,loc,x,0  ^: v3 ]: D8 U7 @4 t. o
lsel,a,loc,x,100" @% z) c; q$ ?! v, N
lesize,all,,,8/ s+ q  g) Y% z5 ~6 k: ^
type,3
1 V3 i' Y: F/ G# Ramesh,all; n( e4 E* Y+ o1 D  T& R" [
) ^! p0 w2 b0 }9 h$ O7 v0 q2 L- w
type,2: [1 b" X* |" c2 \5 Q
real,3
5 h& n$ i1 W& k9 [' z& d$ S( L, bmat,3
. b) z( T, r- o1 g3 S$ B2 T7 M! R6 @extopt,esize,32) G4 z) j$ `; v0 J
extopt,aclear,1, p9 g6 X0 L5 s" z$ v  q& y$ n
vext,all,,,,,-800
6 U# O) P$ P( i% r9 Z1 i% l
1 ?0 `7 W, o; T+ `2 ^/pnum,mat,1- d9 C/ l& V: Y6 h: x  n
/pnum,node,0. q7 `3 D" \( Q; C/ w
/pnum,elem,0
3 B4 h) c" @6 p/ e. U* n/view,1,1,1,1
+ Z: t5 J7 u1 H- x, feplot
) r3 E( h7 \7 @allsel
, l' U! |. E6 \% z  W* v8 L% jnummrg,all, B- C; J/ e) Q# w1 C4 R
numcmp,all
' k9 p& q9 k' P8 l7 E" aeplot
0 m, c! }& y# t) b1 t7 l) Q' @! ^. L3 E/ t7 u6 }4 z$ F
/solu5 q! X' W: y0 e+ n; s
nsel,s,loc,y,0
; D9 I' v& C% m; k" i3 pnsel,r,loc,z,0" `0 G2 L/ B4 O: V% P, h
d,all,uy
) ^7 a4 S7 l; C7 ]
5 o* E/ O: Z8 R6 D1 B6 l" A3 t; J$ D) jasel,s,loc,z,-8008 F/ V% b' K! y( G8 ]- M
da,all,symm( H1 R9 j' \4 O) Z# i$ ^
7 ]* }, \2 S& c2 H2 |2 ]
allsel
. C2 |3 L9 m3 S' k, \, f! S  ]+ Xgplot
' _& I( z7 x; j) v
, U$ J& l- f' f7 `5 ~nsel,s,node,,1129,1129,1  t+ V# P1 R2 p4 N% Z5 l# w
nsel,a,node,,1134,1134,1
) P  J0 X2 p- |) gf,all,fy,-10009 k+ F( D' k8 p( X1 Q3 l
nsel,s,node,,1130,1133,1
5 W0 j! W6 Y6 \3 hf,all,fy,-20000 s6 m1 m. a1 T+ _( F5 l$ f
allsel
' `$ [( }% }# Q1 N" T# h0 N9 b
* @' f2 r7 Q) Qnlgeom,on1 W( E% @: F, i9 E8 E, j0 w
nsubst,200; B' x' S/ `/ f/ L& ~2 a
outres,all,all
) O  \5 y  q9 b4 q9 r  [3 nneqit,500 N+ O! ]3 b5 L+ V7 p/ [
pred,on) }. o( \  R2 [4 C" p  H9 A) S
!ncnv,2% ]6 F4 j8 Y# e: m
cnvtol,f,,0.02,2,0.5
7 R# c8 t7 W) m" ]0 R* x; T  K8 p. F! D8 t
solve0 O. w' t! z1 g9 H
8 f0 [; }# u1 k. v$ B; Z
finish

评分

参与人数 1堡币 +2 收起 理由
anhui_ts + 2 交流答疑

查看全部评分

精彩评论倒序浏览

7666查看4评论

沧海学生小... 发表于 2006-12-11 21:58:30
不知道shihonglei说的是什么意思& Y' r4 h: T1 w7 u+ h) e
这是我用你的命令流做出的图
aaaa.jpg
举报 回复
沧海学生小... 发表于 2006-12-11 21:59:31
你说的是坐侧的挠度为0吗?
举报 回复
ppz11 发表于 2008-3-13 16:08:31
不知我怎么运行这个程序有错呢?
举报 回复
iamccq 发表于 2013-8-22 22:58:50
好想边界出了问题?
举报 回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

 
 
  • QQ:56984982
  • 点击这里给我发消息
    电话:13527553862
    站务咨询群桥头堡站务咨询桥梁专业交流群:
    中国桥梁专业领袖群
    工作时间
    8:00-18:00