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

shihonglei 发表于 2006-12-8 21:01:52 | 显示全部楼层 | 阅读模式
我最近用ansys模拟碳纤维布加固钢筋混凝土梁,不过好像有点问题,计算结束以后,梁好像没有变形呀!挠度始终是零,帮忙看看!谢谢了!
( v2 c% a' q* L+ t1 l!碳纤维布加固钢筋混凝土梁有限元分析
  i1 o% f" W# L; Y/com,structural  ]4 @# @- a& s* Z0 c+ h( ^! K7 ~% Z
/prep7! X6 T# ~% B% T6 z  D4 `
!定义单元类型/ G5 Z3 I0 C; p- R
et,1,link81 l' l$ c. s9 |' c# u
et,2,solid652 ]* M9 }  s4 H( E
et,3,plane42
1 y; M- v- o8 E8 d" S: eet,4,shell41
" z, `1 y8 p9 e5 u$ Z
+ ^. a) b2 v( F& H* P& `* c!定义实常数) R/ _' |. O$ n  y$ Q1 P" F
r,1,28.3   !箍筋和受压钢筋
5 `; n: @) ]% c/ gr,2,113.1    !受拉钢筋
1 [0 p& Z  G5 l/ _" t/ kr,3         !混凝土& g, ?. z- B- z( L9 b
r,4,0.12   !碳纤维布
! m6 n; y0 z) ~% f( }* M) Z. d0 Z$ N3 E: Z8 H
!定义材料属性5 r% W/ F6 A- Y6 Y& Y
mp,ex,1,2e5    !箍筋和受压钢筋
0 a0 q* h' e8 X$ ymp,prxy,1,0.25# Y. a, @1 Y0 T. q/ c  ~
tb,bkin,1,1,2,1$ K$ d. [9 B( K
tbdata,,200- I- c! t/ s( T: Q
- `# t2 D: q- ^# |, B7 u. [; v
mp,ex,2,2e5    !受拉钢筋) {3 s' S6 h6 E1 ?& P0 p( I- U& c
mp,prxy,2,0.3
% o( d5 c5 j1 ]6 L6 q! U7 W3 B# Ntb,bkin,2,1,2,1
3 y6 e. r% C, s4 ^; etbdata,,3002 H+ _- d. @) T" x& m

5 m  g/ j" v1 ~8 b# ^mp,ex,3,30000   !混凝土
# j* B3 q% O: P6 K7 g2 ]2 Tmp,prxy,3,0.2
5 Z6 m  q% e3 b# k% Atb,conc,3,1,9
2 A# m/ z) n; i+ c  X# Dtbdata,,0.4,1,2.04,-1+ @4 F1 [, J7 |9 m
+ J3 e2 y8 j; L6 y$ g3 m: r
mp,ex,4,225000   !碳纤维布7 q; ]; ]: F3 @7 W
mp,prxy,4,03 W+ T! w0 y5 G9 D
tb,bkin,4,1,2,1
  w1 D6 J5 }) j1 y4 G, Wtbdata,,4300
" Z% {# N! ~; S. P2 Y' T# {* ], Z; @2 l: B1 F
/pnum,node,1- M/ w0 u" r$ K: M4 C
/pnum,elem,12 m+ }2 X2 G4 {! h" L: u/ }. T1 h& N
/ E- o! E* M4 E" f8 |, u+ t: p/ ?9 R
!产生所有的节点
: n3 U* p% V9 {n,1# k/ C; T! `* ]. N8 N& n
n,6,100  u! @# L$ d6 s' Q; ^
fill,1,60 J; K# Q5 m* s% B' |
ngen,9,6,1,6,1,,25
) [9 s3 v! B% }/ q, d* |: wngen,33,1000,1,54,1,,,-25
, U5 Z9 S4 a# b8 ?7 B* T: ?1 J/view,1,1,1,1; A6 m$ {0 v$ N9 n2 P
nplot
4 m( w- x$ |# |1 E) S
0 N  H- w2 l! H1 C+ ~, T0 e!箍筋和受压钢筋单元建立% c5 e* v: |7 i9 y9 T$ G1 F+ {5 c
type,1& C: g, B( {9 |) g3 K' P# |
real,1
) l& q: b! r2 z7 [mat,1* A6 a5 k. [" s8 |1 J: {$ M$ V
*do,ii,8,10,1* b8 g9 M0 u& e# F# o
   e,ii,ii+1
0 F1 }1 C+ y6 Q* ^& s*enddo
& t) X& [  E- o4 I8 d*do,ii,11,41,6' b$ z( C# v/ O6 c# h6 w; `
   e,ii,ii+6/ e+ Z" }$ {4 }1 D  @8 U* r9 i/ q
*enddo- D3 m+ y1 x" d' I
*do,ii,44,46,1
  L+ Y1 L; G( M+ B, D! M   e,ii,ii+1
* `6 s4 a& e5 d  w) Z' g: d*enddo
; L- \3 \$ V# C: x1 n3 b! u; V6 K*do,ii,8,38,6
. L' q; q4 \* `: t2 k$ k  K   e,ii,ii+6
- ^  y) H  p- s8 w# K+ ?# U3 Y/ F*enddo
+ a7 `! Q! F3 l# {7 z8 ^egen,6,4000,1,18,1' S: }  \: a6 K7 Q$ |3 V
!受压钢筋
2 ?9 l$ n* p! C/ B% w$ i, W*do,ii,44,19044,1000
* ?& e% j( `5 ^+ U) Q    e,ii,ii+10000 p. h- W5 K+ L" Y* C
*enddo
; x( Q/ ?3 k( |2 d4 h/ \*do,ii,47,19047,1000
& W0 z2 b7 b) \+ p, }     e,ii,ii+1000
: @* }6 y& p4 R. y9 x*enddo
" M; l+ L% l- ^- l5 l/ S, P% r5 F& K: E* i% O; Q8 A
!纵向受拉钢筋单元建立  c' r: z$ X( @+ u* z3 @+ C' I
type,1  Q% e$ @3 A, V6 U. [) F# V
real,2
  W! K+ [2 ?. V, C9 v" l  rmat,2. P3 d# M( @, G
*do,ii,8,31008,1000
5 H+ t/ e) d% h6 S5 A* c    e,ii,ii+1000; t8 h9 h! ^7 a9 y/ o+ I3 B' J: M
*enddo. U, A" s' W: }" K
*do,ii,11,31011,1000* Q3 l& z6 v' ^% a' l$ b
    e,ii,ii+1000
7 |) y/ g5 k1 p+ Y. O# c*enddo; u' |% X! w, ]) T
!碳纤维布单元的建立
. \* L+ p5 A. ptype,4- j; i& C: E/ f( f0 Q5 X4 s) k6 E
real,4
2 H2 r) f. Y" z7 r  X9 R+ a( rmat,4
0 F5 N! ?! v9 ]7 A/ n# `" T*do,ii,1,5,1* L! S9 d  }9 Z2 Q  C5 y: d
    *do,jj,ii,31000+ii,1000
5 p8 k, _: `9 c$ B/ X        e,jj,jj+1,jj+10003 `8 m7 L# d( b6 \& K
        e,jj+1,jj+1000,jj+1000+1
. K$ M: M* J0 }6 Q    *enddo
) X7 h$ h# h& n1 x' ]*enddo4 x1 @$ t8 X, R2 u; r, n

7 j; ~7 _+ d+ \5 T* V, `) t/view,1,1,1,1$ Z; a9 R: P+ D0 i% \- x7 ~
/pnum,elem,0+ V6 R! k" B, G  e; g& G
/pnum,node,0
$ U1 {" x& ?  J7 s/eshape,1) D# n- X- \: J4 z
eplot
+ q/ c- \9 u, O% P. c( x
* t& O. n! L9 b* Y!混凝土单元的建立
" D9 l+ y" m8 O0 F; }k,19 q+ `" y; ]8 k; `& M, A9 D
k,2,100: c- P+ O4 Y& E/ ~6 v! H0 o( y; V
k,3,100,2003 E% w( Q: E& U: [
k,4,,200
; z$ L; ^4 i. L/ j# N( y3 Ga,1,2,3,4* A, w) n' j6 x
& @4 _$ [3 o' @. C$ u" ~% n
lsel,s,loc,y,0* ?8 a$ z  \( `9 K& J
lsel,a,loc,y,200) [/ Y3 b6 Q* j
lesize,all,,,5
* p2 l! X! ]( F5 T& q9 flsel,all
8 W3 |- L# I/ g8 Zlsel,s,loc,x,0
( D' D# d( H8 G, @9 Xlsel,a,loc,x,100
' i. T2 n! N2 D/ ]" P' k( flesize,all,,,8
9 X- L! F6 P8 b* v4 e9 B( ^  Ptype,3
0 y$ u& E* N1 X" O/ S9 oamesh,all
: ?  K- N. O- q2 ~+ w* D
, l$ J, X2 N3 R; J" e0 mtype,2
0 k5 j6 ?! K2 c: ^real,3
: M, s$ G% d/ Fmat,3
( z$ Q6 b1 P" e/ o1 rextopt,esize,32
9 a7 ]- `3 a+ r& T  j' S, B5 |extopt,aclear,1
/ Q$ {5 f) d4 U3 {- c8 |vext,all,,,,,-800- i* D$ m* a2 e! Q: F5 j
  O9 [$ J2 I8 X4 n- K( l
/pnum,mat,1
) W4 `' x. w9 Y, J/ q/pnum,node,0
4 ~- W8 k/ d2 X( |/pnum,elem,0
1 v, B, r; t) d( Y5 K/ u4 c- P; W5 d/view,1,1,1,19 f" n) T  a5 X9 h/ G& }
eplot+ d( {+ h; T% ]9 t
allsel
9 }3 T" U5 f. V  w: K* N9 _nummrg,all
$ V7 F3 Q! ~! rnumcmp,all
' V0 {0 K2 f6 T" feplot
  [( Q/ C% Y# e2 f5 M/ \/ g3 p6 v9 y6 d' O# T
/solu" M9 t' v+ H# n1 _, Y) [
nsel,s,loc,y,0
6 w  j' j( {! I% f0 cnsel,r,loc,z,0% b5 i* g& }  x3 M7 ?0 n( l
d,all,uy
5 d. L( o; H2 B8 X2 T
; [; V! K0 M$ X* `asel,s,loc,z,-800
2 P; ~4 |( o2 u1 K/ n! fda,all,symm% V, }3 @/ s3 @% S* N# j
& X3 N) k# X9 ~
allsel  R6 e- r7 o) @7 z$ u8 f! F. f
gplot
: {% [- `  P1 {. d9 J6 D$ |6 F" P  ~
nsel,s,node,,1129,1129,1
4 o3 L$ e" U6 N* M; a; Hnsel,a,node,,1134,1134,1( Q/ ]: n. y" i+ J) D9 A) o- n
f,all,fy,-1000
4 Y) e0 d' L5 |" ~7 M6 Rnsel,s,node,,1130,1133,1
0 a( S0 g  D' O' cf,all,fy,-2000* V2 m5 T5 H/ z3 f+ C& S$ ?
allsel! ?% x6 D7 D+ K6 h0 |4 H
: u! v  R9 i9 `# g8 n+ i
nlgeom,on
0 `) Y' H8 X6 P" R% T9 pnsubst,200
& P& i4 R9 T0 }! T7 g. e8 j/ B( woutres,all,all6 ]3 p+ o0 H& Z" F" x. ~
neqit,50
; c$ c0 u2 j& Z4 x$ U1 w8 dpred,on
* f) c8 K0 f3 F* _( V!ncnv,2
1 l1 i6 w/ e  }. X; R, t' Pcnvtol,f,,0.02,2,0.5" P2 K9 c/ k# o: F; i$ E( Y

8 k7 Z  O4 C& m1 g6 `# z3 ]solve
: C8 e7 `7 _5 |# v5 Z- M
2 |% H/ z  G5 i' i" pfinish

评分

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

查看全部评分

精彩评论倒序浏览

6902查看4评论

沧海学生小... 发表于 2006-12-11 21:58:30
不知道shihonglei说的是什么意思
1 G5 z: z, `$ M: T5 T% x这是我用你的命令流做出的图
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