【原创】给大家一个刚刚作过的例题(页 1) - C06：ANSYS--实例赏评 - SimWe仿真论坛 |Simwe.com|仿真|设计|有限元|虚拟仪器

xiezhenyu 发表于 2004-8-8 08:28

【原创】给大家一个刚刚作过的例题

[color=red][/color]
/title,research on RC Beam
!混凝土研究(2004.6.27 16:00  加载120kN，分级逐加)
/prep7
!定义单元特性
rd1=8.0/1000
rd2=22.0/1000
et,1,solid65
et,2,link8
et,3,solid45
!定义混凝土材料特性
r,1
mp,ex,1,2.4e10
mp,prxy,1,0.1667
mp,dens,1,2500
fc=25.0e6       !混凝土抗压强度25.0MPa
fr=3.1125e6    !混凝土抗拉强度3.1125MPa
hntra=fc
op=0.35  !裂缝张开剪力传递系数 op=0.7
cl=0.75  !裂缝闭合剪力传递系数 cl=1

tb,concr,1
!tbdata,,op,cl,fr,fc !考虑压碎
tbdata,,op,cl,fr,-1 !不考虑压碎
kn=14 !应力应变关系描述点数
tb,miso,1,,kn
tbpt,,0.0002,hntra*0.192
tbpt,,0.0004,hntra*0.36
tbpt,,0.0006,hntra*0.51
tbpt,,0.0008,hntra*0.64
tbpt,,0.0010,hntra*0.75
tbpt,,0.0012,hntra*0.84
tbpt,,0.0014,hntra*0.91
tbpt,,0.0016,hntra*0.96
tbpt,,0.0018,hntra*0.99
tbpt,,0.0020,hntra*1.00
tbpt,,0.0025,hntra*0.95
tbpt,,0.0030,hntra*0.90
tbpt,,0.0035,hntra*0.85
tbpt,,0.0040,hntra*0.80

!定义钢筋材料特性
mp,ex,2,200e9  !纵向钢筋
mp,dens,2,7850
r,2,acos(-1)*0.25*rd2*rd2
mp,prxy,2,0.3
!受拉受压如何表示
tb,bkin,2  !双线性随动强化
tbtemp,0
tbdata,,360e6,0

!横向箍筋!(包括上排纵向钢筋)
mp,ex,3,200e9
mp,dens,3,7850
r,3,acos(-1)*0.25*rd1*rd1
mp,nuxy,3,0.25
tb,bkin,3  !双线性随动强化
tbtemp,0
tbdata,,210e6,0

blc4,,,0.150,0.300,1.325
/view,1,1,1,1
/ang,1
kwpave,6
wpoffs,,,-75/1000
vsbw,all
wpoff,,,-50/1000
vsbw,all
wpoffs,,,-75/1000
vsbw,all
*do,i,1,14,1
wpoffs,,,-75/1000
vsbw,all
*enddo
kwpave,6
wpoff,,30/1000,
wprot,,90,
vsbw,all
wpoff,,,-240/1000
vsbw,all
wprot,,,90
wpoff,,,-18.75/1000
vsbw,all
wpoff,,,-112.5/1000
vsbw,all
nummrg,all,,,,low
numcmp,all
wpsty !关闭工作平面显示

DESIZE,1,2,15,15,28, , ,1,4, !单元尺寸控制
lsel,s,loc,x,18.75/1000
lsel,a,loc,x,131.25/1000
lsel,r,loc,y,30/1000
latt,2,2,2 !下排纵向钢筋

FLST,5,8,4,ORDE,8
FITEM,5,423
FITEM,5,428
FITEM,5,463
FITEM,5,468
FITEM,5,558
FITEM,5,567
FITEM,5,606
FITEM,5,613
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,25/1000, , , , , , ,1
  

FLST,5,28,4,ORDE,28
FITEM,5,407
FITEM,5,417
FITEM,5,433
FITEM,5,438
FITEM,5,443
FITEM,5,448
FITEM,5,453
FITEM,5,458
FITEM,5,493
FITEM,5,498
FITEM,5,473
FITEM,5,478
FITEM,5,483
FITEM,5,488
FITEM,5,540
FITEM,5,550
FITEM,5,572
FITEM,5,577
FITEM,5,584
FITEM,5,589
FITEM,5,596
FITEM,5,601
FITEM,5,642
FITEM,5,645
FITEM,5,618
FITEM,5,625
FITEM,5,630
FITEM,5,637
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,1, , , , ,1

ALLSEL,ALL
lsel,s,loc,x,18.75/1000
lsel,a,loc,x,131.25/1000
lsel,r,loc,y,270/1000
latt,2,3,2 !上排纵向钢筋

FLST,5,8,4,ORDE,8
FITEM,5,56
FITEM,5,73
FITEM,5,233
FITEM,5,253
FITEM,5,393
FITEM,5,398
FITEM,5,228
FITEM,5,261
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,25/1000, , , , , , ,1
FLST,5,28,4,ORDE,28
FITEM,5,15
FITEM,5,42
FITEM,5,90
FITEM,5,106
FITEM,5,120
FITEM,5,138
FITEM,5,163
FITEM,5,167
FITEM,5,199
FITEM,5,205
FITEM,5,665
FITEM,5,667
FITEM,5,288
FITEM,5,363
FITEM,5,368
FITEM,5,373
FITEM,5,378
FITEM,5,383
FITEM,5,388
FITEM,5,516
FITEM,5,520
FITEM,5,403
FITEM,5,501
FITEM,5,506
FITEM,5,511
FITEM,5,649
FITEM,5,655
FITEM,5,660
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,1, , , , ,1

ALLSEL,ALL
lsel,s,loc,y,30/1000,270/1000,240/1000
lsel,r,loc,x,19/1000,131/1000
latt,2,3,2

FLST,5,38,4,ORDE,38
FITEM,5,552
FITEM,5,554
FITEM,5,561
FITEM,5,563
FITEM,5,580
FITEM,5,592
FITEM,5,609
FITEM,5,621
FITEM,5,633
FITEM,5,688
FITEM,5,690
FITEM,5,692
FITEM,5,694
FITEM,5,696
FITEM,5,698
FITEM,5,701
FITEM,5,-702
FITEM,5,705
FITEM,5,-706
FITEM,5,708
FITEM,5,710
FITEM,5,712
FITEM,5,714
FITEM,5,716
FITEM,5,718
FITEM,5,720
FITEM,5,722
FITEM,5,724
FITEM,5,726
FITEM,5,728
FITEM,5,730
FITEM,5,732
FITEM,5,734
FITEM,5,736
FITEM,5,738
FITEM,5,740
FITEM,5,742
FITEM,5,744
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,18.75/1000, , , , , , ,1

ALLSEL,ALL
lsel,s,loc,x,18.75/1000,131.25/1000,112.5/1000
lsel,r,loc,y,31/1000,269/1000
latt,2,3,2

FLST,5,38,4,ORDE,4
FITEM,5,521
FITEM,5,-539
FITEM,5,668
FITEM,5,-686
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,30/1000, , , , , , ,1

ALLSEL,ALL
LSEL,r,TYPE,,2
lmesh,all
ALLSEL,ALL

FLST,5,2,4,ORDE,2
FITEM,5,154
FITEM,5,298
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,1, , , , ,1
FLST,5,4,4,ORDE,4
FITEM,5,1
FITEM,5,154
FITEM,5,298
FITEM,5,410
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,1, , , , ,1
FLST,5,1,4,ORDE,1
FITEM,5,323
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,30/1000, , , , , , ,1
FLST,5,1,4,ORDE,1
FITEM,5,699
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,18.75/1000, , , , , , ,1

ALLSEL,ALL
vsel,all
vatt,1,1,1
mshape,0,3d
mshkey,1
vmesh,all
allsel,all
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL

!定义钢板
r,1
mp,ex,4,2e11
mp,prxy,4,0.25
mp,dens,4, 7.85e3
tb,bkin,4  !双线性随动强化
tbtemp,0
tbdata,,360e6,0

ALLSEL,ALL
ESEL,S,TYPE,,1
FLST,5,40,2,ORDE,8
FITEM,5,661
FITEM,5,-662
FITEM,5,669
FITEM,5,-671
FITEM,5,893
FITEM,5,-897
FITEM,5,1265
FITEM,5,-1294
ESEL,R, , ,P51X

EGEN,2,3000,all, , , , , , , , ,-30/1000, ,

ALLSEL,ALL
ESEL,S,TYPE,,1
FLST,5,48,2,ORDE,6
FITEM,5,646
FITEM,5,-651
FITEM,5,884
FITEM,5,-889
FITEM,5,1211
FITEM,5,-1246
ESEL,R, , ,P51X
EGEN,2,6000,all, , , , , , , , ,30/1000, ,

ALLSEL,ALL
ESEL,S,TYPE,,1
FLST,5,88,2,ORDE,2
FITEM,5,2633
FITEM,5,-2720
ESEL,R, , ,P51X

EMODIF,all,type,3,
EMODIF,all,MAT,4,
EMODIF,all,real,1,

ALLSEL,ALL
NUMMRG,NODE, , , ,LOW
NUMCMP,NODE

!中间截面约束
FLST,2,99,1,ORDE,32
FITEM,2,1
FITEM,2,15
FITEM,2,45
FITEM,2,75
FITEM,2,-81
FITEM,2,208
FITEM,2,212
FITEM,2,-216
FITEM,2,283
FITEM,2,-289
FITEM,2,446
FITEM,2,-450
FITEM,2,561
FITEM,2,-563
FITEM,2,628
FITEM,2,630
FITEM,2,632
FITEM,2,639
FITEM,2,646
FITEM,2,-650
FITEM,2,682
FITEM,2,-683
FITEM,2,834
FITEM,2,-840
FITEM,2,1020
FITEM,2,-1022
FITEM,2,1116
FITEM,2,-1122
FITEM,2,1480
FITEM,2,-1484
FITEM,2,1700
FITEM,2,-1734
/GO
D,P51X, , , , , ,UX,UZ, , , ,

!边支座约束
FLST,2,9,1,ORDE,6
FITEM,2,2577
FITEM,2,2579
FITEM,2,2589
FITEM,2,2591
FITEM,2,2604
FITEM,2,-2608
/GO
D,P51X, , , , , ,UX,UY,, , ,

allsel,all
NUMMRG,NODE, , , ,LOW !压缩合并节点
NUMCMP,NODE
NUMMRG,ELEM, , , ,LOW !压缩合并单元
NUMCMP,ELEM
allsel,all
finish

/solu
!ACEL,0,9.8,0    ! 重力加速度
antype,0
cnvtol,f,,0.10,2,0.05  !定义力收敛条件10%，使用不缺省的value:5%
                         !0.05  The minimum value allowed for the program calculated reference value.
nsubst,15          !定义子步数15
outres,all,LAST    !输出每一子步结果
                    !由于结果文件最大4G，本程序已经超出，故要求控制
autots,1          !打开自动时间步控制
lnsrch,1          !线性搜索
ncnv,2             !如果不收敛时结束不退出
neqit,10          !每一子步中方程的迭代次数控制10
pred,on          !打开预测器
kbc,0
NROPT,FULL, ,ON    !全牛顿拉谱森法
LUMPM,0
EQSLV,SPAR, ,0,
PRECISION,0
MSAVE,0
PIVCHECK,1
SSTIF
PSTRES

FCUM,ADD,1,    !将载荷设为后续加载在同一位置时为加载增量值

!加载
*do,n,1,15
time,n
!加载荷载值20000/(9*10)N
!time=1自重不加载
*if,n,eq,1,then
F,2631,FY,0
F,2634,FY,0
F,2645,FY,0
F,2648,FY,0
F,2662,FY,0
F,2663,FY,0
F,2664,FY,0
F,2665,FY,0
F,2666,FY,0
*elseif,n,gt,1 !n>=2,加载2000N
F,2631,FY,-2000/9
F,2634,FY,-2000/9
F,2645,FY,-2000/9
F,2648,FY,-2000/9
F,2662,FY,-2000/9
F,2663,FY,-2000/9
F,2664,FY,-2000/9
F,2665,FY,-2000/9
F,2666,FY,-2000/9
*endif
lswrite
*enddo

lssolve,1,15 !求解
finish

/post26
numvar,10
finish
save

希望多加点分啊
老是分数太低
看不到别人的帖子
xiezhenyu1117@163.com
QQ:158210884

jincai 发表于 2004-8-9 22:52

回复：【原创】给大家一个刚刚作过的例题

好，谢谢！
坚决顶

jesseliying 发表于 2004-8-10 14:40

回复：【原创】给大家一个刚刚作过的例题

谢谢！很有帮助。

paul131 发表于 2004-8-12 23:07

回复：【原创】给大家一个刚刚作过的例题

学习中……

Alegend 发表于 2004-8-13 17:01

回复：【原创】给大家一个刚刚作过的例题

   谢谢楼主的无私奉献！

   不过还是有一点纳闷，这里的加分是什么规则？

   ^_^，没有什么别的意思，说说而已。

foransys 发表于 2004-8-22 10:04

回复：回复：【原创】给大家一个刚刚作过的例题

[quote][b]Alegend wrote:[/b]
   谢谢楼主的无私奉献！

   不过还是有一点纳闷，这里的加分是什么规则？

   ^_^，没有什么别的意思，说说而已。 [/quote]

新来0分的，加1分以示欢迎和鼓励。

xiezhenyu 发表于 2004-8-23 00:41

回复：【原创】给大家一个刚刚作过的例题

给大家这个模型的图片，这是我在实验室做的实验模型，结果与千分表册出来的结果很接近！

dudumax 发表于 2004-9-1 17:48

回复：【原创】给大家一个刚刚作过的例题

thank u so much!

zhim 发表于 2004-9-7 12:01

回复：【原创】给大家一个刚刚作过的例题

very good JOB！~~

注释也不错，好好学习，

支持斑竹加分，大大的加：）

tanchangjian 发表于 2004-9-8 11:16

回复：【原创】给大家一个刚刚作过的例题

thanks

good luck

xiezhenyu 发表于 2004-12-17 13:18

Re:回复：【原创】这个是用SHELL93做的，和试验很符合

!有横隔板的T形薄壁梁

/prep7
!***********************材料特性
!单位mm,N    
et,1,shell93
mp,ex,1,2800
mp,prxy,1,0.38
r,1,5,5,5,5!huban
r,2,3,3,3,3!lumian
r,3,10,10,10,10!henggeban
!***************************建模
k,1,!路板
k,2,0,0,1050
k,3,350,0,1050
k,4,350,0,0
a,1,2,3,4

k,5,35,!竖板
k,6,35,,1050
k,7,35,-65,1050
k,8,35,-65,
a,5,6,7,8
asel,s,,,2,2,1
agen,5,2,,,70,,

k,100,35,,30！横隔板
k,101,105,,30
k,102,105,-65,30
k,103,35,-65,30
l,100,103
l,103,102
l,102,101
l,101,100
al,27,26,25,28
agen,5,7,7,1,,,990/4
agen,4,7,11,1,70,,
asel,all
aglu,all

/VIEW,1,0.520874808298,0.397111924926,0.755639830317  /ANG, 1,  -5.22472093507
/REPLO

!***************************web

lsel,s,,,125,128,1!1duan
lsel,a,,,153,154,1
lsel,a,,,171,172,1
lsel,a,,,189,190,1

lsel,a,,,132,133,1!2duan
lsel,a,,,136,137,1
lsel,a,,,157,158,1
lsel,a,,,175,176,1
lsel,a,,,191,192,1

lsel,a,,,138,139,1!3duan
lsel,a,,,161,162,1
lsel,a,,,179,180,1
lsel,a,,,193,194,1
lsel,a,,,142,143,1

lsel,a,,,144,145,1!4duan
lsel,a,,,148,149,1
lsel,a,,,165,166,1
lsel,a,,,183,184,1
lsel,a,,,195,196,1
lesize,all,,,4

lsel,s,,,203,206,1
lesize,all,,,1

!short line
lsel,all
lsel,u,,,1,3
lsel,u,,,203,206,1

lsel,u,,,125,128,1!1duan
lsel,u,,,153,154,1
lsel,u,,,171,172,1
lsel,u,,,189,190,1

lsel,u,,,132,133,1!2duan
lsel,u,,,136,137,1
lsel,u,,,157,158,1
lsel,u,,,175,176,1
lsel,u,,,191,192,1

lsel,u,,,138,139,1!3duan
lsel,u,,,161,162,1
lsel,u,,,179,180,1
lsel,u,,,193,194,1
lsel,u,,,142,143,1

lsel,u,,,144,145,1!4duan
lsel,u,,,148,149,1
lsel,u,,,165,166,1
lsel,u,,,183,184,1
lsel,u,,,195,196,1
lesize,all,,,2

aatt,1,1,1
asel,s,,,78,95,1
asel,a,,,40,41
asel,a,,,48,
asel,a,,,52
asel,a,,,55
asel,a,,,44
amesh,all

aatt,1,2,1
asel,s,,,96,97
!smrtsize,4
!mshkey,0
amesh,all

!henggban
aatt,1,3,1
asel,s,,,7,11,1
asel,a,,,39,51,4
asel,a,,,54
asel,a,,,57,63,2
asel,a,,,64
asel,a,,,66,72,2
asel,a,,,73
amesh,all

!huban
aatt,1,2,1
asel,all
asel,u,,,7,11,1
asel,u,,,39,51,4
asel,u,,,54
asel,u,,,57,63,2
asel,u,,,64
asel,u,,,66,72,2
asel,u,,,73
asel,u,,,96,97
asel,u,,,78,95,1
asel,u,,,40,41
asel,u,,,48,
asel,u,,,52
asel,u,,,55
asel,u,,,44
amesh,all

/eshape,2
eplot
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
！约束
FINISH
/SOL
ALLSEL,ALL
LPLOT
FLST,2,8,4,ORDE,8
FITEM,2,42
FITEM,2,44
FITEM,2,150
FITEM,2,-151
FITEM,2,169
FITEM,2,-170
FITEM,2,187
FITEM,2,-188
!*
/GO
DL,P51X, ,UY,
FLST,2,8,4,ORDE,8
FITEM,2,42
FITEM,2,44
FITEM,2,150
FITEM,2,-151
FITEM,2,169
FITEM,2,-170
FITEM,2,187
FITEM,2,-188
!*
/GO
DL,P51X, ,Ux,

FLST,2,8,4,ORDE,8
FITEM,2,26
FITEM,2,28
FITEM,2,129
FITEM,2,-130
FITEM,2,155
FITEM,2,-156
FITEM,2,173
FITEM,2,-174
!*
/GO
DL,P51X, ,UY,
FLST,2,8,4,ORDE,8
FITEM,2,26
FITEM,2,28
FITEM,2,129
FITEM,2,-130
FITEM,2,155
FITEM,2,-156
FITEM,2,173
FITEM,2,-174
!*
/GO
DL,P51X, ,Ux,

!*******加载 save
FLST,2,1,1,ORDE,1
FITEM,2,263
!*
/GO
F,P51X,FY,-790
/sol
solve
/post1
pldis,0

/EFACE,1
AVPRIN,0, ,
!*
PLNSOL,EPTO,Z,0,2

nsel,s,,,698
nsel,a,,,694
nsel,a,,,766
nsel,a,,,283
nsel,a,,,886
PRNSOL,U,Y !位移

smhai 发表于 2004-12-17 13:36

Re:【原创】给大家一个刚刚作过的例题

请问顶层的大侠，要是用你的命令流怎么求极限荷载？

Happyfaery 发表于 2005-2-17 20:23

Re:【原创】给大家一个刚刚作过的例题

[color=red]不错!!![/color]

xiezhenyu 发表于 2005-3-1 23:36

Re:【原创】给大家一个刚刚作过的例题

这个我没有想过，我现在做碳纤维布加固的分析，有人在做着方面吗QQ158210884

lucaslijianfeng 发表于 2005-5-4 17:40

Re:【原创】给大家一个刚刚作过的例题

谢谢

junfengtang 发表于 2005-5-4 20:57

Re:【原创】给大家一个刚刚作过的例题

学习

wangdengke99111 发表于 2005-8-31 14:53

Re:【原创】给大家一个刚刚作过的例题

非常感谢楼主！

ppz11 发表于 2007-9-14 19:33

有人做碳纤维加固梁的ansys程序吗？要是有的话，请给我传一个吧！先谢谢了。
emal to ：zhaopeipei-_1983@163.com

yuanfang910 发表于 2008-1-3 14:42

多谢楼主了，呵呵

mumu_ld 发表于 2008-1-7 14:53

找人代做论文。
结构硕士论文。模型已建好，属性，荷载等都已完成，可直接计算。由于时间关系，需找人帮我挑出结论并分析撰写部分论文。我们可以相互讨论，就算交个朋友。可以的话发邮件给我，我们谈谈报酬
[email]mumu_ld1981@163.com[/email]

zhangshouyi18 发表于 2008-1-7 18:45

多谢楼主了

guowhlg 发表于 2008-1-23 14:29

好的，谢谢了啊

xjhome2500 发表于 2008-3-7 21:00

厉害呀
louzhu

页: [1]

SimWe仿真论坛's Archiver