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

shihonglei 发表于 2006-12-8 21:01:52 | 显示全部楼层 | 阅读模式
我最近用ansys模拟碳纤维布加固钢筋混凝土梁,不过好像有点问题,计算结束以后,梁好像没有变形呀!挠度始终是零,帮忙看看!谢谢了!' P* d& t2 G7 a4 s  Z$ F9 t7 r* ]$ j. h
!碳纤维布加固钢筋混凝土梁有限元分析
* j3 w6 u, u1 z  G( n4 n/com,structural' ]! v* W. I$ Y4 v  j4 E
/prep74 I/ ]6 ~2 D& M' r* E- Y
!定义单元类型
/ b  d: @) a' l' \" P: \. ~et,1,link8
0 G8 T5 T7 A3 }+ get,2,solid651 S! D4 h7 ^# j4 q/ `- _' l
et,3,plane42: J& _# l) t9 g9 ^) F, Y6 m
et,4,shell41
+ _) c- F2 \" p& i" M  U- u1 `
!定义实常数
( Y% b8 A1 l% {% d* q' P( _8 ]. yr,1,28.3   !箍筋和受压钢筋
' v+ v. i' `  @" Hr,2,113.1    !受拉钢筋
% L+ \0 j0 S; P! K: hr,3         !混凝土
3 C( e6 s) C, Q! Y& m, L8 u$ or,4,0.12   !碳纤维布
- J0 u* L5 o/ E" \
6 `9 T; K1 N" N4 z! \2 l& [( U!定义材料属性1 A, A" T4 |; L8 u6 [
mp,ex,1,2e5    !箍筋和受压钢筋5 [* C* p& S" \" ?/ Q, }. C; s
mp,prxy,1,0.25
' B" J" Z/ R. I9 {# ftb,bkin,1,1,2,1
+ h# U* T- S. y" xtbdata,,2002 o+ X: g0 y. B0 C

' n& i7 b/ i6 }' Xmp,ex,2,2e5    !受拉钢筋
' }' V9 J/ Z  c8 g2 _: T( }mp,prxy,2,0.3% f" S6 D2 a# s" {. a* _* B
tb,bkin,2,1,2,1
! N3 f5 [7 n  _6 f, \! ?" Ktbdata,,300. g% D0 b4 y0 n+ X- R& B
# A4 W# `; y, q, u) H$ u
mp,ex,3,30000   !混凝土3 ]! t( ]/ s# |& _, l7 n: Z! G
mp,prxy,3,0.2: d  H1 m1 ], U7 `7 R" }: e- R
tb,conc,3,1,9
8 Y' \/ q0 a6 R# g! Utbdata,,0.4,1,2.04,-1
( R! e1 P% U2 ?
, ^6 i% h7 t8 L3 Imp,ex,4,225000   !碳纤维布6 x' v9 V# j' h8 v
mp,prxy,4,0
9 s* l+ Z' ?1 v, Z8 s: ntb,bkin,4,1,2,1
7 F6 f5 b* z  i" ctbdata,,4300
! k; U, M% T' p& {
4 K1 C) n* w2 ^, Q- W/pnum,node,1& W6 ^. K: j( E$ B; `4 z, v- ^
/pnum,elem,1
+ h( E) \0 ?$ F/ V% w, l6 A" I7 \( @; ?
$ [; |6 v0 B% R/ I  H!产生所有的节点) `# S- J3 Q% o0 U6 R6 }1 Q6 f8 B7 N
n,18 h! n7 N9 j4 w* t7 \0 q
n,6,100! Y4 e9 R# T& {" U! d' I' ?
fill,1,6
( v% I9 j4 L% _' Tngen,9,6,1,6,1,,25
& ]) n1 ^# n  qngen,33,1000,1,54,1,,,-250 ~2 z( s' a$ K5 [/ t/ @! d9 U3 h
/view,1,1,1,1$ R: w! C1 j/ S* L
nplot$ J" l8 p: S' W7 v+ a1 m
  m4 |1 P1 ~  T2 D/ ^' F' u  I
!箍筋和受压钢筋单元建立
* j9 w. v$ b. H, h$ p) @type,1
, I! x! R  Q  p. L# Lreal,14 S7 m+ X6 ]* r7 |
mat,1
" s3 ?9 Y+ G3 \1 ]! c# V*do,ii,8,10,1$ z# u+ A7 n- o; U. ~
   e,ii,ii+1
9 O8 R# `5 |/ F$ M& e9 C# g*enddo
, B2 u$ {4 A$ ^5 i, r: Q0 c4 ?" y*do,ii,11,41,60 g6 L5 h+ F. s0 r' k* ]: a
   e,ii,ii+6$ F7 j$ s0 T1 U, u) |
*enddo
% E1 r( F. O/ r" J9 R# c*do,ii,44,46,17 B; J1 a. T" u# X
   e,ii,ii+1: N, X; j7 b' Z  U& g! U
*enddo
2 Q7 h# E2 q3 @$ O$ M*do,ii,8,38,6
# I$ F# d0 m; R) j' h9 B' ?/ m   e,ii,ii+6% J0 q6 d1 [& i! m
*enddo. I" F5 F) I2 o
egen,6,4000,1,18,1
* v3 I- r% F* h  W5 X* G7 k# i9 N!受压钢筋+ {6 ]/ [+ a* T
*do,ii,44,19044,1000
7 |4 k9 z0 g* u  @6 e" O. c    e,ii,ii+10000 `& M7 Z* G9 H: a
*enddo4 E- H+ W  N/ ~9 V7 k+ x( h+ Z' `/ u
*do,ii,47,19047,1000
: A2 a6 v/ d# E     e,ii,ii+1000
" ^2 h% e0 V: k# K; I# S6 w# y*enddo
- {2 h( l7 Z- L; R
% n9 k: d: d$ E) c!纵向受拉钢筋单元建立
2 C5 l# b% x' t- Q6 ctype,1" x4 w, ^7 I- t; ~; M( Y( [" `
real,2. g& X: h6 x& q1 Z
mat,2
0 ?6 \/ v2 W( c; ^3 }. l0 v1 M/ l2 ^*do,ii,8,31008,1000
) I5 `4 w6 M% Y: i" X( v0 X    e,ii,ii+1000
8 Z2 _4 J' u2 m- e+ w1 O% a4 Y*enddo( K3 Q; i$ C: L$ b
*do,ii,11,31011,1000& H; B5 }' b2 [* W
    e,ii,ii+1000
7 t4 R9 r9 v9 s) A4 n# a9 H3 o$ u*enddo
2 e7 G+ R6 x. ]7 Z!碳纤维布单元的建立4 X$ W7 Q8 u, c/ u
type,43 d9 y, m! r& b! r% j
real,4% }& O9 m6 R3 a
mat,4' D- \7 t5 E- V8 y3 F# H
*do,ii,1,5,1$ y4 C- Z" |4 B) d
    *do,jj,ii,31000+ii,1000
3 |3 }9 z3 n9 i) d5 p        e,jj,jj+1,jj+1000. |5 b- A" l; ?* R( W
        e,jj+1,jj+1000,jj+1000+1( b. Z% k9 [: u. z
    *enddo
0 v9 T9 _& \, p& F5 C*enddo' }1 k6 @( ^3 F1 D0 o5 H. ], k0 u
1 k  a, d' n+ o2 m3 l  q# X( x& I
/view,1,1,1,1
+ b  y6 d) b8 J/ [' j/pnum,elem,0' N2 ~$ Q* y4 [5 t
/pnum,node,01 ?  j) N& H* X  V
/eshape,1
) _; I7 ]; D0 L. X% heplot& i' Q8 `5 N$ `0 [2 `# a8 h
& p7 C7 a  d  \1 d& R# ]
!混凝土单元的建立2 B4 n, j! u* k% z1 Y8 p
k,1
2 g$ }  _) \1 \7 \& T! N2 v6 Xk,2,100
3 p$ v2 T0 V2 P* C  Z0 e" n+ tk,3,100,200& {$ n9 Z% b: n- }/ V
k,4,,200  P* a5 X. ]5 ?/ P; s- o& e6 Y
a,1,2,3,4, I. F# O3 b0 s, X7 n8 O2 e1 Y$ T
) X+ B) p2 E9 n& ]+ Q! h
lsel,s,loc,y,0
5 G1 m$ q# ~; c, H& ?lsel,a,loc,y,200; F: a, x# r6 V/ ~- ^' F
lesize,all,,,53 D0 c$ Q" ~. E* m% f8 x
lsel,all8 G, N# T* D, p- o; e: V6 y$ t2 W
lsel,s,loc,x,0
! e6 d" W" \5 V. i* H. ilsel,a,loc,x,100. P) O1 E* B8 F0 t- U  z
lesize,all,,,8
$ ]+ A' r% c! h" t# b$ Ztype,3
7 l* b3 u  D, gamesh,all
. q$ ^4 z) g- {$ ?4 K% T. ^
& U4 K- H2 _5 O6 W' a4 Atype,2
# g4 |: M4 o6 w0 e2 |) [real,3
7 I' `6 y% J2 u, n- ymat,3" @9 ^$ J  a& n" O3 O
extopt,esize,32# z. H0 C$ n. `; X7 S
extopt,aclear,1, ~! C5 m% f! Q+ z0 ?$ [
vext,all,,,,,-8009 i- g5 Q1 ?& Y

3 v% ?" ^$ O7 x/pnum,mat,10 }6 w& v6 p  |4 C$ G+ k8 k
/pnum,node,0
+ S3 h! d8 }  b" |/pnum,elem,0
" J6 Q/ M" p& B' D/view,1,1,1,1
0 h9 Q& t7 _' _: o! E7 oeplot: n/ }" f: T7 r
allsel. U- u. N1 Y, L5 S4 a" q  K/ i
nummrg,all, |( T2 o( [7 R) _
numcmp,all
4 b* L" p1 V: D- {+ ?eplot4 S# C" F/ S9 I' r

  i! X% X/ A8 @9 @/solu+ H8 ~6 [0 u( j# L! O) y2 p
nsel,s,loc,y,0, \9 |: u8 T; b6 F
nsel,r,loc,z,0: t+ S( `, e7 N9 c! ]/ T3 N5 J  o
d,all,uy
* s# U( q) A- |/ T" J& ]! S& C& W' C2 A( I+ v
asel,s,loc,z,-800
' _  `2 Y, m; L, |" Tda,all,symm
4 h/ b/ l! j: ^* m' Z
+ N+ K; B' t* Kallsel* A9 I. p( P2 H$ M
gplot
) L9 t" Y9 t  o9 n- @' R
9 _& m/ j9 y- n2 u# b8 O! Q  Qnsel,s,node,,1129,1129,1
5 z  z# @( a8 n- Gnsel,a,node,,1134,1134,1
8 L7 d0 r& _* G0 O; ^4 mf,all,fy,-1000
* b# W9 G* P' b9 Lnsel,s,node,,1130,1133,1* z. C: k2 N5 F0 x' e1 k* ^
f,all,fy,-2000
! i5 Z/ Q" V5 X3 K+ M7 g: q) Nallsel
' [! }) q) x7 ]4 b& E! g- g( f- O. x7 l5 H# c
nlgeom,on* z& u! a! ]! N' r6 K9 d
nsubst,200& v- N3 \2 i8 x  P
outres,all,all: D+ _, z- j2 m. y" c
neqit,501 I( v% O, x: D# j0 `( Y/ \2 b# Z
pred,on/ [0 ]: v- n2 c8 F  ]2 A
!ncnv,2- V" a& Q6 H. J; f) b
cnvtol,f,,0.02,2,0.5) ^1 @2 p9 \  y" V0 M

$ }6 u: L. K( x. Isolve, a0 P- o8 q) f# b) k0 j9 m6 H- W

, G2 I& i1 L& M: E1 `finish

评分

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

查看全部评分

精彩评论倒序浏览

7704查看4评论

沧海学生小... 发表于 2006-12-11 21:58:30
不知道shihonglei说的是什么意思
/ b3 r; [; r# G4 B8 V0 s这是我用你的命令流做出的图
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