我最近用ansys模拟碳纤维布加固钢筋混凝土梁,不过好像有点问题,计算结束以后,梁好像没有变形呀!挠度始终是零,帮忙看看!谢谢了!
3 { _2 s9 V1 @1 n9 M!碳纤维布加固钢筋混凝土梁有限元分析+ j8 f$ u+ R" q6 b
/com,structural5 j5 e0 y7 j( t: @' j5 q: k- A
/prep7! |. x1 N' d/ Y0 M- i: C6 A
!定义单元类型
/ ~) K) O0 }4 ^/ U/ K0 X. eet,1,link8) e; a& s5 Q1 x" f- I% ]
et,2,solid653 I% x/ m1 b3 X2 B' J" q5 \
et,3,plane42! d. U* x- L, ]; o' O( c3 p
et,4,shell41
2 R& c' c9 D, W% X' ?) o
2 }3 e S4 m5 ]; z!定义实常数
# o4 N7 l4 B7 q a! V; F, Tr,1,28.3 !箍筋和受压钢筋
8 z2 g. L( `" H/ r, N* Zr,2,113.1 !受拉钢筋
. V, |* u+ O! M! m$ l Vr,3 !混凝土
. o( N6 n& X$ [! U1 ~# b vr,4,0.12 !碳纤维布
) y) @0 M2 s% \, r0 F$ Y, b1 I
5 J% C8 Q+ }5 k& f% g!定义材料属性
6 g' D2 ?" X7 O- ?$ j0 E# `. Kmp,ex,1,2e5 !箍筋和受压钢筋
1 c( g5 V8 @: mmp,prxy,1,0.25
% F# Y5 p# R- [& U. E8 Otb,bkin,1,1,2,16 p/ b8 `7 x+ R% o& A
tbdata,,200
d6 h/ o9 V* u$ d+ k0 e
9 h& ]- o3 Q$ { Y& i5 tmp,ex,2,2e5 !受拉钢筋
# |" q! I3 N8 m" w" L. fmp,prxy,2,0.3* [, [7 @( s/ f3 j
tb,bkin,2,1,2,1* H2 s, e* x Q$ D" N; M; U0 a
tbdata,,300" O5 ~# M$ k' N7 S
1 @1 U9 X2 w; @" N* Kmp,ex,3,30000 !混凝土8 j. V$ [1 a o7 @! _" D- K
mp,prxy,3,0.2
5 c! \% T* N' l9 Gtb,conc,3,1,9 t0 [; d d+ O0 X
tbdata,,0.4,1,2.04,-10 w% `' _' y: N }' I
% F8 @8 o* _7 C. j. M
mp,ex,4,225000 !碳纤维布
) ^' X( r# @! O! L: Bmp,prxy,4,0
. F! ], n- g1 f" ~: ^1 M% htb,bkin,4,1,2,1$ \- @3 Y# G8 b2 D
tbdata,,4300
& S' x" P+ s3 m @
, s0 z2 O; M F8 ?2 L1 D2 P/pnum,node,1
& q9 f- j5 O3 m/pnum,elem,1+ P+ f6 s. q2 A$ O4 h$ z/ N
" |" R$ X& w% f
!产生所有的节点
n: j" G% `$ E5 I9 b8 Dn,1
* E9 j4 ]" p5 ^( B3 [( Z+ wn,6,100
3 Z' l. ~+ R' _+ y$ J3 ]fill,1,6 w$ J- J1 g- E6 |: l; a+ b/ \
ngen,9,6,1,6,1,,25
1 }; J1 t! ~- x9 x% M8 jngen,33,1000,1,54,1,,,-25( o5 i+ O. y/ U [( [
/view,1,1,1,18 F- J7 ]) B. l8 }5 M/ y! e. I" v( k
nplot+ p L F0 n; ]7 q* E6 ~4 S
; }& t. |1 B4 N!箍筋和受压钢筋单元建立: Z* v% o* X8 J& C. ^
type,1
3 L7 S) A5 N3 b: `real,1& {/ m/ y9 B$ i" `/ } K) Q
mat,1
( a& i& [2 Z7 E% C' P*do,ii,8,10,1
4 |+ Z. L/ h7 B, w e,ii,ii+1! m5 [2 N7 N; o3 F% G1 I
*enddo
; l5 \6 @- U! ?7 e/ r" O- [*do,ii,11,41,60 ^5 Q# u0 f* D$ [$ ~/ J+ H
e,ii,ii+6
0 y4 y9 q; q* l*enddo, C! j4 R m; Z$ P# I# t! b
*do,ii,44,46,1
o( s( f+ v. M& O e,ii,ii+16 S4 V$ {, U$ h; F3 C9 C' U- a
*enddo' o {% S4 X) g% M& u8 k
*do,ii,8,38,6
& q3 [# L: a" A( ` k L e,ii,ii+6+ S7 G+ y3 o T" h) L2 t, e: x
*enddo" j* g7 ^: `. e3 d% }& D
egen,6,4000,1,18,18 z7 v- g# }; B. ?( o# y! D
!受压钢筋% }5 r0 |0 i1 I, A5 ~5 [
*do,ii,44,19044,1000
5 r$ v E3 S: v e,ii,ii+1000
6 Z5 S0 n2 s- z* n0 t2 d$ A6 b3 h4 p# F*enddo
/ Y9 M7 o8 o7 q*do,ii,47,19047,1000
5 X% a9 w5 A! |4 h e,ii,ii+1000# P& M* d i* g- _7 E1 y
*enddo
6 {" G6 o& ~2 D
5 t1 W p, I6 B- O9 g: L4 w!纵向受拉钢筋单元建立
+ C8 p; U( ?/ o7 [+ utype,1
; |1 ?( L" ^' P" j7 b4 [real,2
& v: n# ]1 `6 N8 h4 D1 W, dmat,2
$ E0 G, a( O7 l. e4 S0 s4 L*do,ii,8,31008,1000
; ?) t" i+ Q8 T, e0 M$ D4 } e,ii,ii+1000; y1 N$ P2 K4 g; |
*enddo
7 W% @' t- S ^* _: N*do,ii,11,31011,1000
- C; s+ C( w5 ]; } J: _9 O: n e,ii,ii+10006 W1 `! w! g$ r( k/ Y5 s
*enddo3 v* l" S" T* p! H; o7 i T+ H7 M
!碳纤维布单元的建立( K7 x7 x- o5 Q
type,4
! G' b5 D$ D S& R- q5 l+ Y; c6 vreal,4
) P& N' x" h. Imat,4+ a4 _; h s4 y f% ]3 b
*do,ii,1,5,1' G$ o ^, V* a) d
*do,jj,ii,31000+ii,1000
- X |$ p" p0 h e,jj,jj+1,jj+1000
, {3 J V* m6 B# V e,jj+1,jj+1000,jj+1000+1
2 I! w9 r, c( U *enddo
. s( A1 c5 {+ r+ X- h*enddo) ]7 V! T% O: k' _( L7 h# R" l; |
+ ]# B. z$ |& a9 O5 s W/view,1,1,1,16 o0 d% c" A. d% D2 X3 O
/pnum,elem,0
" Q; R. |" t- |3 x4 o5 |/pnum,node,0
8 o* \9 W z5 Q' T, e( d/eshape,1
, |) l* z1 g7 B+ ceplot, l( h+ H6 t7 l5 h
x e+ _/ D% Z/ {9 G: j!混凝土单元的建立
* M) r% A2 R! n- J) F2 `k,1% o9 }0 W. C! n; n( w' i
k,2,1003 Z! h, Y. z- I) x4 k4 e5 t/ E; k
k,3,100,200
: v' U% U4 a% p# a- Bk,4,,200
; [8 s$ K0 I/ ], Xa,1,2,3,4
, ^& g' ?' c- V& c
* T+ y* f/ M# Y- H8 z- A8 b8 I9 Ulsel,s,loc,y,0
4 F' F9 q; t! b: y7 elsel,a,loc,y,200
4 c) U$ ~" k7 p _; J. |lesize,all,,,5- @& O; i$ H+ S. g
lsel,all
2 e' Q! C& Q5 B) Q& H$ Mlsel,s,loc,x,0
( w2 K- Q( d# ~: d% R- i% klsel,a,loc,x,1001 `) A2 o$ ^7 e" N
lesize,all,,,8- m$ ]9 C/ c! U. p& e) F/ I, y x
type,3
' [- a2 Y2 _4 N( M$ lamesh,all
4 z( U! V9 N8 _; B. X# G1 `
* m. G% W5 W& z9 f3 M7 Ptype,2) L' g6 K0 T& |
real,3
0 X z% o' J% G+ F# A9 imat,32 D: w9 C# c6 a. L
extopt,esize,32
+ y/ b3 q a* d: \extopt,aclear,1- b4 p; |+ F( e1 b( U0 n* ]0 G
vext,all,,,,,-800
* i( U$ p4 H4 l; @9 o& N& R% K9 c/ ^0 ]& F
/pnum,mat,18 B$ v8 t7 `& p& q) I8 l& o
/pnum,node,0% k {; b% a& s
/pnum,elem,0
. y4 `( L! ?" `$ k) x' I/view,1,1,1,1
' K9 e; x8 h1 s- `0 \* F1 yeplot
7 ?6 z+ |4 j: K2 W4 Z, Pallsel
/ K- ?7 j# h( B/ X# `$ Wnummrg,all
0 w- ~9 Q" c8 s4 ^( l3 enumcmp,all" D' @# X- {. y- U2 ]
eplot
- R) `0 G r6 y
! H9 k/ f k9 W' W! Q( S/solu
" m3 I, v5 L) r! v# h( ^nsel,s,loc,y,09 O3 e9 d# Q) M0 M- U
nsel,r,loc,z,0' x2 S' j. J/ M
d,all,uy b# _: v$ G* c8 Y+ A! a- c
8 D* R- j _8 c5 N6 Y
asel,s,loc,z,-800! l# l6 `' L1 K+ O4 S0 i$ r
da,all,symm
- b' F% d; U( `3 k
' U' D$ {( s3 y$ g7 b. {) jallsel
3 s6 _" F) t& G, Q. Y9 N: Agplot+ \0 X) y) s4 Q" X1 O
; k* }2 x1 }( ~% Xnsel,s,node,,1129,1129,10 ~- w+ N, E/ O# g
nsel,a,node,,1134,1134,1
. n$ g- D9 w4 k; @& ]- {+ \' R' ff,all,fy,-1000
# J& F/ l4 ]% q& _% f1 }nsel,s,node,,1130,1133,1
+ a7 N- M. u& v/ x, n& Mf,all,fy,-2000' o* W0 C C. [: [
allsel
( r. S5 B O0 A0 v) B1 a+ L `/ g3 S( w- H+ X
nlgeom,on, g3 m: w' ]; I( u. n9 n
nsubst,200
2 f2 m2 U4 E9 boutres,all,all, U o' h2 z, i4 g, y1 G
neqit,505 n8 n7 C; ?( J5 E
pred,on% C6 o- M! p. t# d
!ncnv,2
, v6 X: O( S# C9 S6 Ccnvtol,f,,0.02,2,0.5
3 V# j# y& j3 D5 a8 L" J, @2 n5 r! f* L7 V1 {9 O) h) N: j3 D$ c6 m* J
solve# d2 O, S* z) H( S, b9 s% O* a
+ e4 |/ H. V; Q6 R i! }+ |finish |