问题描述
大家
我正在使用 Ansys apdl 执行非线性分析。
我正在寻找 T 形截面混凝土梁上的裂缝区域。
非线性分析会关闭 ansys 并且不会收敛到解决方案。
有人可以帮我吗?
我究竟做错了什么?
如果有人可以帮助我,我将不胜感激。
代码如下。你可以看到一些葡萄牙语的评论。
!T beam
!
/clear
!
/Filename,T beam !File name do Viga T
/Tittle,T beam analysis !Title
!
/Units,SI,!
/prep7
!
!beam length (m)
L=4 !m
!
!Model scale
Esc=1 !real model
!
!Force (N) and force orientation(+,-)
!number nodes which receive the load at each side
nodes=21
fp=-10000/nos
!
!reinforcement area - 2 bars,diameter= 8 mm
A=(5.026548e-5)*2 !m²
!
!T beam dimensions
lt=0.48*Esc
ct=L*Esc
et=0.05*Esc
ht=0.40*Esc
!
!element type
et,1,solid65,9.66e-4
!
!element 1 Keyopts
!KEYOPT,ITYPE,KNUM,VALUE
keyopt,1
keyopt,3,0
!keyopt,5,6,3
keyopt,7,8,0
!
!et,0.003590,et,2,reinf264
!
!Propriedades do Concreto armado
/pnum,mat,0
secnum,1
mat,1
mp,ex,28e9
mp,prxy,0.2 !Coeficiente de Poisson do concreto - material 1
mp,dens,25e3,!Densidade do material (N)
tb,concr,1 !ativa a tabela de dados para um material considerado
tbdata,0.8,0.2,2.5e6,25e6,!Define as propridades não-lineares para o material 1
!
!Propriedades do Aço
secnum,2
mat,2
mp,200e9
mp,0.3
mp,78.5e3
tb,biso,2
tbdata,2000
tbdata,500e6
!
!reinforcement position - reinf265
mat,2
secnum,2
sectype,reinf,smear
secdata,A,elef,0.7,0.7
!
!Ativa o material 1 e a secnum,1
mat,1
secnum,1
!
!cross section keypoints
K,0*Esc,0.20*Esc
K,0.28*Esc
K,0.02*Esc,0.28*Esc
K,4,(ht-et)*Esc,0.28*Esc
K,lt*Esc
K,ht*Esc,0.28*Esc
K,0.20*Esc
K,9,0*Esc
K,10,11,0.20*Esc
K,12,0.20*Esc
!
K,13,ct*Esc,0.20*Esc
K,14,15,16,17,lt*Esc
K,18,lt*Esc
K,19,0.28*Esc
K,20,21,0.00*Esc
K,22,0*Esc
K,23,24,0.20*Esc
!
!Areas da alma,recebe o reforço
a,11
a,13 !area 2
a,13 !area 3
a,23 !area 4
a,2 !area 5
a,15 !area 6
!
!Define volumes a partir da união de áreas existentes
va,6
!
!Restante da alma,até a mesa
a,11 !area 7
a,23 !area 8
a,3 !area 9
a,23 !area 10
a,12 !area 11
!
!Define volumes a partir da união de áreas existentes
va,11
!
!Define o restante da nervura,no topo da mesa
a,12 !area 12
a,8 !area 13
a,24 !area 14
a,4 !area 15
a,7 !area 16
!
!Define o volume da nervura,no topo da mesa
va,16
!
!Deine a mesa
a,12 !area 17
a,10 !area 18
a,22 !area 19
a,24 !area 20
a,21 !area 21
!
!Define o volume da nervura,no topo da mesa
va,21
!
!Define o restante da mesa
a,7 !area 22
a,16 !area 23
a,4 !area 24
a,6 !area 25
a,18 !area 26
!
!Define o volume da nervura,25,26
!
!Tamanho da malha aplicada ao bloco 1
!ESIZE,SIZE,NDIV
esize,0.025*Esc !Elementos de malha com dimensões de 2,5 cm de lado!!
!
!Comando de malha
vmesh,1
!
!Propriedades do Aço
mat,2
secnum,2
!
!Transforma as áreas selecionadas em barras de aço
ereinf,!
!Define o nível de translucência,a qual permite ver a posição das barras no interior da viga
/TRLCY,elem,0.5,all
!
!Ativa o material 1
mat,1
secnum,1
!
!Define as propriedades associadas aos volumes
vatt,!
!Tamanho da malha aplicada ao bloco 2
!ESIZE,NDIV
esize,5 cm de lado!!
!
!Comando de malha
!VMESH,NV1,NV2,NINC
vmesh,5
!
!Apresenta os nós do modelo
/PNUM,node,on
!
!Condições de contorno
!D,Node,Lab,VALUE,VALUE2,NEND,NINC,Lab2,Lab3,Lab4,Lab5,Lab6
!
D,ux
D,uy
!D,uz
!
D,172,uz
!
D,174,uy
!D,uz
!
D,173,uz
!
D,171,ux
D,uz
!
D,170,ux
D,uy
!D,uz
!
D,810,ux
D,uz
!
D,651,uy
!D,uz
!
D,492,uz
!
D,175,169,uy
!D,uz
!
D,809,uz
!
D,650,491,176,168,uy
!D,uz
!
D,808,uz
!
D,649,uz
!
D,490,177,ux
D,167,807,648,489,178,166,uz
!
D,806,uz
!
D,647,uz
!
D,488,uz
!
D,179,165,805,uz
!
D,646,487,uz
!
D,180,164,804,645,486,uy
!D,uz
!
D,181,163,803,644,uz
!
D,485,uz
!
D,182,ux
D,uz
!
!Apoios na outra extremidade
D,uz
!
D,652,493,334,uz
!
D,333,653,494,335,332,654,495,336,331,uz
!
D,655,496,337,330,656,497,338,uz
!
D,329,657,uz
!
D,498,339,328,658,499,340,327,659,uz
!
D,500,341,326,uz
!
!Condições de contorno impostas à mesa
!
D,14515,uz
D,14676,14675,18380,18381,18379,18363,18364,18355,14491,uz
D,14500,uz
D,14499,uz
!
!Retira a numeração dos nós do modelo
/PNUM,off
!
!Plota os elementos de malha
esel,all
eplot
!
!Insere as forças nodais nos nós desejados
!Ponto x=1.75m
F,12982,fy,fp
!
F,13288,fp
F,13626,fp
F,13785,13944,fp
!
F,14585,fp
!
F,16752,16753,16754,16755,16756,16757,16758,19764,20062,20221,20380,20539,fp
F,20698,20857,21016,fp
!
!Ponto x=2.25m
F,12962,13308,13606,13765,13924,14605,fp
!
F,16612,f y,16613,16614,16615,16616,16617,16618,fp
!
F,19744,fp
!
F,20082,20241,20400,20559,20718,20877,21036,fp
!
!Solução
/solu
!
!ANTYPE,Antype,Status,LDSTEP,SUBSTEP,Action,--,PRELP
antype,static
!
!inclui os efeitos de grandes deflexões
nlgeom,on
!
!Especifica o número de passos passos a serem considerados,incluindo o número de passos mínimos e
máximos
!NSUBST,NSBSTP,NSBMX,NSBMN,Carry
nsubst,50,1000,1
!
!NROPT,Option1,Option2,Optval
nropt,auto,crpl,0.15
!
!ARCLEN,Key,MAXARC,MINARC
arclen,on
!
!Resolve para a condição de carregamento selecionada ("load step")
solve,!
/post1
set,list
!
finish
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)