用lesize指定划分段数时不成功 不知为何

印记央央

我用命令流来说明吧
不知道为何我纵向的划分对了 竖向的就是不行，明明命令流都差不多的

Finish

/cle

/filn,watermelon,1

!柱尺寸

Lp=6250

Hp=700

Bp=240

Tpw=16

Tpf=10

!次梁尺寸

lg=7600

hg=450

bg=160

tgw=8

tgf=6

!螺栓尺寸

Bolt0=23.5

Bolt=22

Nut=35

Nut_h=14

S1=80     !线距

Sa=50     !水平方向端距

Sb=50     !竖直方向端距

Sc=55     !加劲肋离次梁翼缘的距离

/prep7

et,1,solid65,,,,,,,1     !定义混凝土板单元为soild65

et,2,solid45

et,3,combin39

et,4,solid95         !过渡金字塔单元

et,5,solid92         !螺栓单元

keyopt,4,11,1

keyopt,3,2,0

keyopt,3,3,3        !沿Z轴的平动自由度

keyopt,1,1,0        !考虑大变形

keyopt,1,7,1        ！考虑应力松弛，为半脆性开裂，有助于计算收敛

mp,ex,1,30400           !定义混凝土材料的弹性模量，泊松比

mp,prxy,1,0.2           

fc=32$ft=2.64

tb,concr,1$tbdata,,0.5,0.95,ft,-1    !关闭压碎

tb,miso,1,,11            !定义混凝土应力应变曲线，采用MISO模型

tbpt,,0.0002,fc*0.19$tbpt,,0.0004,fc*0.36$tbpt,,0.0006,fc*0.51

tbpt,,0.0008,fc*0.64$tbpt,,0.0010,fc*0.75$tbpt,,0.0012,fc*0.84

tbpt,,0.0014,fc*0.91$tbpt,,0.0016,fc*0.96$tbpt,,0.0018,fc*0.99

tbpt,,0.002,fc$tbpt,,0.0033,fc

mp,ex,2,2.0e5           !定义钢筋的弹性模量，泊松比

mp,prxy,2,0.3

tb,biso,2               !钢筋的应力应变曲线，用BISO模型

tbdata,,290,0         

mp,ex,3,2.06e5           !定义材料三，钢梁的模量，泊松比

mp,prxy,3,0.3

tb,bkin,3

tbdata,1,310

mp,ex,4,2.06e5

mp,prxy,4,0.3

tb,bkin,3

tbdata,1,400               !假设的

!tb,mkin,3

!tbtemp,,strain

!tbdata,,0.001,0.0015,0.025

!tbtemp,0.0

!tbdata,,206,310,420

!tb,mkin,4

!tbtemp,,strain

!tbdata,,0,003,0.006,0.008

!tbtemp,0.0

!tbdata,,618,940,1040

tbplot

pi=acos(-1)

v1=0.25*pi*6*6*6/(400*50)     !上部区纵筋的体积配筋率

v2=0.25*pi*6*6*3.5/(400*50)     !下部区纵筋的体积配筋率

v4=0.25*pi*6*6*24.5/(1920*50)   !水平箍筋的体积配筋率

r,2,2,v1,0,90,2,v4,0,0             !定义拉区混凝土单元的实常数

r,3,2,v2,0,90,2,v4,0,0             !定义压区混凝土单元的实常数

r,1

r,4

r,5,0,0,0.19,20099,0.5,29291

rmore,1.0,38012,1.5,43978,2,48550

rmore,2.5,52244,3,55323

!创建局部坐标系

local,11,0,,,lp/2

local,12,0,tpf/2+(bp-tpf)/4,,lp/2+tgf/2

local,13,0,bp/2+(2*sa+10)/2,(hp-hg)/2, lp/2+tgf/2

local,14,0,bp/2+10, (hp-hg)/2, lp/2

local,15,0,bp/2+10+sa,hp/2-sc-sb, lp/2+tgf/2+8

local,16,0,lg/2,hp/2,lp/2,,-90

!创建主梁

Csys,0

K,1,0,-hp/2$k,2,bp/2,-hp/2$k,3,bp/2,-hp/2+tpw$k,4,tpf/2,-hp/2+tpw

K,5,tpf/2,hp/2-tpw$k,6,bp/2,hp/2-tpw$k,7,bp/2,hp/2$k,8,0,hp/2

A,1,2,3,4,5,6,7,8

Vext,1,,,,,lp

Allsel

Wpcsys,,12

Blc5,,, (bp-tpf)/2,hp-2*tpw,8

Wpcsys,,13

Blc5,,,2*sa+10,hg-2*sc,8

vadd,2,3

!创建次梁

Wpcsyc,-1

Csys,14

K,33,0,-hg/2

K,34,0,-hg/2,bg/2

K,35,0,-hg/2+tgw,bg/2

K,36,0,-hg/2+tgw,tgf/2

K,37,0,hg/2-tgw,tgf/2

K,38,0,hg/2-tgw,bg/2

K,39,0, hg/2,bg/2

K,40,0, hg/2

A,33,34,35,36,37,38,39,40

Vext,11,,,lg-bp/2-10

Vsymm,z,2

Vadd,2,3

allsel

Cm,v1,volu

!创建高强螺栓

Wpcsys,,15

Cyl4,0,0,bolt0/2,,,,-(8+tgf)$vgen,2,2,,,,-s1$vgen,2,3,,,,-s1$vgen,2,6,,,,-s1

Vsel,u,,,v1

cm,v2,volu

allsel

vsbv,v1,v2

cm,v3,volu

wpcsys,,15

cylind,bolt/2,nut/2,0,nut_h

vgen,2,2,,,,,-(tgf+8+nut_h)

wpoff,,,nut_h

cylind,0,bolt/2,-(2*nut_h+8+tgf)

vsel,u,,,v3

cm,b1,volu

vgen,2,b1,,,,-s1

cm,b2,volu

vgen,2,b2,,,,-2*s1

allsel

vsel,u,,,v3

vglue,all

cm,bolt,volu

allsel

wpcsys,,16

Blc5,,,lg,lp,100

!**********************网格划分*************************！

!划分高强螺栓

Allsel$cmsel,s,bolt

Vatt,4,6,5,0

Mshape,1,3d

Mshkey,0

Vmesh,all

Allsel

Vsel,s,,,9

Aslv$lsla

Wpcsys,14$wpoff,2*sa$wpro,,,90$vsbw,all

Wpoff,,,70$vsbw,all

Wpoff,-tgf/2$wpro,,,90$vsbw,all

Wpoff,,,tgf/2$vsbw,all

Wpoff,,,tgf/2$vsbw,all

Wpoff,,hg/2-tgw$wpro,,90$vsbw,all

Wpoff,,,sc-tgw+sb+s1/2$vsbw,all

Wpoff,,,s1$vsbw,all

Wpoff,,,s1$vsbw,all

Wpoff,,,s1/2+sb+sc-tgw$vsbw,all

Cm,liang,volu

Cm,liang_line,line

Csys,14

Vsel,r,loc,x,2*sa+70,7470$aslv$lsla

Cm,liang_ling1,line

Lsel,r,length,,tgw$lsel,a,length,,tgf/2$lesize,all,,,1

Cmsel,s,liang_line$Lsel,r,length,,7300$lesize,all,,,147

Cmsel,s,liang_line$Lsel,r,length,,s1$lesize,all,,,2

Cmsel,s,liang_line$Lsel,r,length,, s1/2+sb+sc-tgw $lesize,all,,,3

!这两条命令就是将图中最右边的紫色部分所指的2段和3段

Cmsel,s,liang_line$lsel,r,length,,(bg-tgf)/2$lesize,all,,,3

Vatt,3,4,2,0          !solid45

Mshape,0,3D

Mshkey,1

Vmesh,all

Allsel$cmsel,s,liang

Vsel,r,loc,x,0,2*sa+70$vsel,r,,,20,22,1$vsel,a,,,25,28,1$vsel,a,,,30,31$vsel,a,,,7,15,8

vsel,a,,,44,45$vsel,a,,,52,56,4$aslv$lsla

Lsel,r,,,67,97,30$Lsel,a,,,156,254,98$Lsel,a,,,304,443,139$lesize,all,,,3

Lsel,a,,,367,374,7$ Lsel,a,,,353,447,94$Lsel,a,,,410,450,40$lesize,all,,,5

!这两条命令就是梁端上下翼缘划分，即上面图中最上面大红色标出的5段和3段

Vatt,3,4,2,0            !solid45

Mshape,0,3D

Mshkey,1

Vmesh,all

Vsel,s,,,48,49$vsel,a,,,9,12,3$vsel,a,,,54,55$vsel,a,,,40,41

Vsel,a,,,60,61$vsel,a,,,35,36$vsel,a,,,68,69$vsel,a,,,64,65

aslv$lsla$cm,liang_line2,line

Lsel,r,radius,,bolt0/2$cm,r1,line$lsel,r,,,75$lsel,a,,,81,83,1$lsel,a,,,105,108,1

Lsel,a,,,115,118,1$lsel,a,,,125,128,1$lesize,all,,,4$cm,r11,line

Cmsel,s,r1$cmsel,u,r11$lesize,all,,,8

Cmsel,s,liang_line2$lsel,r,length,,s1$lesize,all,,,8

Cmsel,s,liang_line2$lsel,r,length,,s1/2+sb+sc-tgw$lesize,all,,,10

!这两条命令就是图中左边紫色的划分，但是未成功

Cmsel,s,liang_line2$ lsel,r,length,,2*sa$lesize,all,,,10

Cmsel,s,liang_line2$ lsel,r,length,,70$lesize,all,,,7

!这两条命令就是图中中间红色的划分，已经成功了

Vatt,3,4,4,0               !solid95

Mshape,1,3D

Mshkey,0

Vmesh,all

Tchg,solid95,solid92

印记央央

找到问题所在了