vigorao 发表于 2006-8-31 21:04

钢结点分析-出现小主元问题

/Title,FEA node-A ZHJ1-NO.5
/prep7
csys,1
k,7
k,3,4
k,1,4,,1.5
k,2,4,,-1.5
k,4,0.8522599,,2.468143
k,5,0.1072845,,-0.920198
k,6,1.49902937,,-3.1217216
csys,0
l,3,1
l,3,2
l,3,4
l,3,5
l,3,6
l,3,7
LWPLAN,-1,1,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,11,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
CYL4,0,0,0.1365-0.008,,0.1365,,-3       !Main Tube!产生一个厚8mm的钢管长度4m
WPOFFS,0,0,-1   !跟11平面平行的平面向下移动1.5m
CYL4,0,0,0,,0.2,,0!产生一个平面
WPOFFS,0,0,-1
CYL4,0,0,0,,0.2,,0


OFFDIV=0.5
LWPLAN,-1,3,0            !沿着线3
CSWPLA,12,1,1,1
CYL4,0,0,0.084-0.0055,,0.084,,2        !Side Tube 1
WPOFFS,0,0,OFFDIV
CYL4,0,0,0,,0.2,,0

LWPLAN,-1,4,0
CSWPLA,13,1,1,1
CYL4,0,0,0.1365-0.008,,0.1365,,1.5        !Side Tube 2
WPOFFS,0,0,OFFDIV
CYL4,0,0,0,,0.2,,0

LWPLAN,-1,5,0
CSWPLA,14,1,1,1
CYL4,0,0,0.084-0.0055,,0.084,,2       
WPOFFS,0,0,OFFDIV
CYL4,0,0,0,,0.2,,0

CSYS,11
VSEL,S,LOC,X,0,0.1365
CM,MAINTUBE,VOLU   !选择在11柱面坐标下半径等于0.1365的柱体,main tube
ALLSEL,ALL
CMSEL,U,MAINTUBE
CM,TUBES,VOLU            !把其余的支管并到TUBES集合中
ASEL,S,LOC,X,0.1365

CM,CURCL,AREA            !把主管的表面命名为CURCL
ALLSEL,ALL
*DO,I,2,5
   VSBA,I,CURCL,,DELETE,DELETE !用相交的面分割主管体里的支管,分成两部分,主管里和主管外
*ENDDO


CSYS,11
VSEL,S,LOC,X,0,0.1365
CMSEL,U,MAINTUBE
VDELE,ALL,,,1            !删除支管在主管里的多余部分

ALLSEL,ALL
CMSEL,U,MAINTUBE
CSYS,0
BTOL,1e-006                !设置布尔运算精度
VOVLAP,ALL               !使生成的体跟体搭接生成新的体overlap



!vsel,s,volu,,11
!vsel,a,volu,,12
!vsel,a,volu,,5
!vsel,a,volu,,7
!VDELE,ALL,,,1

vsel,s,volu,,3   !碎片
vsel,a,volu,,4
VDELE,ALL,,,1

ALLSEL,ALL
BTOL,1e-008
VGLUE,ALL
ASLV,U       !把在体中的面不选,选中了几个大面不是体中的
CM,CUT,AREA
ALLSEL,ALL
VSBA,ALL,CUT,,DELETE,DELETE!把每个体分成两段,留待以后划分网格时加密区

!******************模型建成,下面划分网格***************************************
et,1,solid45                     !shell181
ex,1,2.06e11
nuxy,1,0.3
DENS,1,7850
TB,BISO,1,1,2,         !等向强化本构
TBTEMP,0
TBDATA,,345e6, 76923e6    !屈服强度,屈服后弹模

vsel,s,volu,,7
vsel,a,volu,,5
vsel,a,volu,,10
vsel,a,volu,,11
vsel,a,volu,,1
CM,OUTV,VOLU

TYPE,1
MAT,1
MSHKEY,3
VSWEEP,ALL       !用扫略的方式划分网格,实用于实体模型


ALLSEL,ALL
CMSEL,U,OUTV!选择中间连接部分,网格需要加密
!%%%%%%%%%%%%%%%%%%%%以下是手动划分%%%%%%%%%%%%%%
TYPE,1
MAT,1
/REPLO
VPLOT   
SMRT,6
MSHAPE,1,3D
!*
FLST,5,8,6,ORDE,8   
FITEM,5,2   
FITEM,5,-3
FITEM,5,6   
FITEM,5,8   
FITEM,5,12
FITEM,5,-13
FITEM,5,16
FITEM,5,-17
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!%%%%%%%%%%%%%%对中间加密部分划分单元完毕%%%%%%%%%
N1=484100
N2=362100
N3=77100
N4=62000
N5=179000
/prep7
LWPLAN,-1,1,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,11,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
NSEL,S,LOC,Z,0
NROTAT,ALL
d,all,ux
d,all,uy

LWPLAN,-1,2,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,12,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
NSEL,S,LOC,Z,0
NROTAT,ALL
d,all,all

! 用荷载步分4步
/prep7
allsel,all
LWPLAN,-1,1,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,11,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
nsel,all
NSEL,S,LOC,Z,0
SF,ALL,PRES,-N1/0.00666

LWPLAN,-1,2,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,12,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
nsel,all
NSEL,S,LOC,Z,0
SF,ALL,PRES,-N2/0.00666


kbc,1
time,1
LWPLAN,-1,3,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,13,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,2
nsel,all
NSEL,S,LOC,Z,2
SF,ALL,PRES,0.25*N3/0.00218

LWPLAN,-1,4,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,14,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,1.5
nsel,all
NSEL,S,LOC,Z,1.5
SF,ALL,PRES,-0.25*N4/0.00666      !按面荷载加载,分母为截面面积


LWPLAN,-1,5,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,15,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,2
nsel,all
NSEL,S,LOC,Z,2
SF,all,PRES,-0.25*N5/0.00218
lswrite,1

!%%%%%%%%%%%%%%%%%%%%%%%%%%%
allsel,all
sfcum,pres,repl
/prep7
LWPLAN,-1,1,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,11,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
nsel,all
NSEL,S,LOC,Z,0
SF,ALL,PRES,-N1/0.00666


LWPLAN,-1,2,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,12,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
nsel,all
NSEL,S,LOC,Z,0
SF,ALL,PRES,-N2/0.00666


kbc,1
time,2
LWPLAN,-1,3,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,13,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,2
nsel,all
NSEL,S,LOC,Z,2
SF,ALL,PRES,0.5*N3/0.00218

LWPLAN,-1,4,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,14,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,1.5
nsel,all
NSEL,S,LOC,Z,1.5
SF,ALL,PRES,-0.5*N4/0.00666      !按面荷载加载,分母为截面面积

LWPLAN,-1,5,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,15,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,2
nsel,all
NSEL,S,LOC,Z,2
SF,all,PRES,-0.5*N5/0.00218
lswrite,2
!%%%%%%%%%%%%%%%%%%%%
allsel,all
sfcum,pres,repl
/prep7
LWPLAN,-1,1,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,11,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
nsel,all
NSEL,S,LOC,Z,0
SF,ALL,PRES,-N1/0.00666


LWPLAN,-1,2,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,12,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
NSEL,S,LOC,Z,0
SF,ALL,PRES,-N2/0.00666


kbc,1
time,3
LWPLAN,-1,3,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,13,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,2
nsel,all
NSEL,S,LOC,Z,2
SF,ALL,PRES,0.75*N3/0.00218

LWPLAN,-1,4,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,14,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,1.5
nsel,all
NSEL,S,LOC,Z,1.5
SF,ALL,PRES,-0.75*N4/0.00666   !按面荷载加载,分母为截面面积


LWPLAN,-1,5,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,15,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,2
nsel,all
NSEL,S,LOC,Z,2
SF,all,PRES,-0.75*N5/0.00218
lswrite,3
!%%%%%%%%%%%%%%%%%%%%%%
allsel,all
sfcum,pres,repl
/prep7
LWPLAN,-1,1,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,11,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
nsel,all
NSEL,S,LOC,Z,0
SF,ALL,PRES,-N1/0.00666


LWPLAN,-1,2,1       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,12,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
nsel,all
NSEL,S,LOC,Z,0
SF,ALL,PRES,-N2/0.00666
NROTAT,ALL

kbc,1
time,4
LWPLAN,-1,3,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,13,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,2
nsel,all
NSEL,S,LOC,Z,2
SF,ALL,PRES,N3/0.00218

LWPLAN,-1,4,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,14,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,1.5
nsel,all
NSEL,S,LOC,Z,1.5
SF,ALL,PRES,-N4/0.00666      !按面荷载加载,分母为截面面积


LWPLAN,-1,5,0       !沿着线1安比例为1即在顶端建一个工作平面,-1是跟原坐标平行
CSWPLA,15,1,1,1   !在工作平面的的原点定义一个局部坐标系,柱面坐标系
WPAVE,0,0,2
nsel,all
NSEL,S,LOC,Z,2
SF,all,PRES,-N5/0.00218
lswrite,4
!%%%%%%%%%%%%%%%%%%%%%%
/solu
antype,static                !静力计算
NSUBST,5,0,0                !子步数
OUTRES,ALL,all        !输出所有结果
AUTOTS,0                !自动时间步长(0关闭,1开启)
TIME,1                        !第一步结束的"时间"

CSYS,0
nsel,all
LSSOLVE,1,4,1

vigorao 发表于 2006-8-31 21:05

怎么解决啊
都约束了啊

vigorao 发表于 2006-9-1 09:49

怎么没人回答?
页: [1]
查看完整版本: 钢结点分析-出现小主元问题