单元生死问题

chenyanru

我作了一个拱桥的模型，准备用分段施工，扣索吊装。我首先是试了一下，没有计入任何单元的自重，仅仅将第一节段的扣索降温了。于是第一节段上挠了。然后我安装了第二节段，也没有计入自重，并且没有施加第二节段的扣索索力，仅仅是在大变形的情况下把第二节段激活了。结果这个激活的过程造成了第一节段的上挠量增加，而且增加的很多。不知道这个是由于什么原因造成的。

[ 本帖最后由 chenyanru 于 2007-9-27 22:08 编辑 ]

xinyuxf

没看明白什么意思

rodge

需要讨论的问题，最好把问题描述清楚，适当的给出图片进行描述

beidoumax

lz把log文件发上来

chenyanru

!**************************************************************
  !以下为建立各个截面的命令
!**************************************************************
et,1,82  
!********截面9-9****************
                  
rectng,1.538,1.49,0,1.18,         !1
rectng,1.49,0,0.042,0.09,         !2
rectng,1.49,0,1.09,1.138,         !3  
rectng,0.096,0.048,0.09,1.09     !4
rectng,1.49,1.25,0.566,0.614    !5
rectng,1.169,1.121,0.09,0.33    !6
rectng,0.814,0.766,0.09,0.33     !7
rectng,0.459,0.411,0.09,0.33     !8
rectng,1.169,1.121,0.85,1.09    !9
rectng,0.814,0.766,0.85,1.09    !10
rectng,0.459,0.411,0.85,1.09   !11
rectng,0.33,0.096,0.566,0.614   !12
asel,all                               !选择所有的截面
aptn,all  
allsel                                 !选择所有实体
aadd,all                               !布尔相加，所有的实体变成一个实体
esize,0.1                            !自动划分网格，5为划分的细度，第一个数字为 1-10
amesh,all   
secwrite,zs48,sect,,1                    !建立一个包含单元划分信息的ASC_2码文件，第一个为文件名，和路径，路径省略则为文件默认路径
SECTYPE,1,BEAM,MESH,                   !截面代号1
SECOFFSET,CENT,,,                       
SECREAD,'zs48','SECT',,MESH              
asel,all                     
ACLEAR,ALL
ADELE,ALL,,,1
!*************9-9**********************  
etdele,1                         !删除代理单元
!定义钢材特性
et,1,beam188                         !不知道此处的截面和上面的截面有什么区别？
mp,ex,1,2.1e11
mp,dens,1,7800
mp,prxy,1,0.3
mp,alpx,1,1e-5

et,2,link10                     
mp,ex,2,2.1e11
mp,dens,2,7800
mp,prxy,2,0.3
mp,alpx,2,1e-5

r,2,0.01,0

!************拱肋计算************
!**********计算拱肋中心线坐标**********
*dim,cs,array,48
*dim,y,array,48
*dim,z,array,48
m=1.9
k=1.257195827                  
f=8.5
y(1)=0
*do,i,1,46,1
y(i+1)=y(i)+0.5
*enddo
y(48)=y(47)+0.349166446
*do,i,1,48,1
cs(i)=(23.34916645-y(i))*2/23.34916645/2
z(i)=f-f*(cosh(k*cs(i))-1)/(m-1)
*enddo
*do,i,1,47,1
n,i,y(i)-23.34916645,z(i)-8.5
*enddo
n,48,y(48)-23.34916645,z(48)-8.5
!********************************************