接触分析建模时碰见一难题，大家一起讨论

zbk0615

我要计算的是一个水电站的钢岔管的模型，还需要考虑钢管与外围混凝土的接触问题，主要参数如下：
主管直径：5.6m；厚度：42mm；
支管直径2.0m；厚度42mm；
主管上有贴边，宽度1.353m;厚度42mm；
钢材弹性模量：E=206GPa，泊松比：0.3;
外围混凝土弹性模量：E=2.55GPa，泊松比：0.167;
混凝土范围可取足够大。
设钢管与混凝土的初试缝隙为0.5mm.
目前是吧钢管模型建立了，就是外面的混凝土不知道怎么处理！
请高手指教！
谢谢！

钢管部分建模命令流如下：


FINISH                !退出以前模块
/CLEAR,START          !清除系统中的所有数据，读入启动文件设置
/FILNAME,ASAHAN_60_7.5_42,1
!(1)设置几何参数
Rzg=2.8          !主管半径(m)
Rcg=1.0         !支（岔）管半径(m)
Tzg=0.042        !主管壁厚(m)
Tcg=0.042        !支（岔）管壁厚(m)
Tbq=2*Tzg        !补强板＋主管壁厚(m)
Lz=16.8                !主管建模长度  取Lz=3×D＝16.8m
Lc=4.45                !支（岔）管远端到主管轴线长度
lcex=6.0        !支管端部延长长度
beita=60        !主管与支管的轴线夹角
sita=7.5        !支（岔）管的半圆锥角，锥管连接时取实际角度（度数），直管连接时取0；
Lbq=1.353        !补强板设计宽度
size=0.2        !划分网格时单元尺寸的大小
*AFUN,DEG        !声明在三角函数中变量单位为度

!******************************
!(2)进入前处理区
!******************************

/PREP7
!**********(2.1)材料参数定义**********************************
!**********材料定义1 主管***********
ET,1,SHELL63                !定义第一类材料单元为弹性壳单元SHELL63
MP,EX,1,2.06E11                !定义第一类材料的弹性模量EX＝2.06MPa
MP,PRXY,1,0.3                !定义第一类材料的泊松比 PRXY=0.3
R,1,Tzg                        !定义第一类材料的厚度等于主管壁厚Tzg

!**********材料定义2 支管***********
ET,2,SHELL63                !定义第二类材料单元为弹性壳单元SHELL63
MP,EX,2,2.06E11                !定义第二类材料的弹性模量EX＝2.06MPa
MP,PRXY,2,0.3                !定义第二类材料的泊松比 PRXY=0.3
R,2,Tcg                        !定义第二类材料的厚度等于支（岔）管壁厚Tcg

!**********材料定义3 补强板***********
ET,3,SHELL63                !定义第三类材料单元为弹性壳单元SHELL63
MP,EX,3,2.06E11                !定义第三类材料的弹性模量EX＝2.06MPa
MP,PRXY,3,0.3                !定义第三类材料的泊松比 PRXY=0.3
R,3,Tbq                        !定义第三类材料的厚度等于补强板＋主管壁厚(m)

!**********(2.2)模型建立*******************************
!*******(2.2.1)建立主管模型
k,9,0,0,0                !创建圆心坐标
k,10,0,0,Lz                !创建10号点，沿主管轴线，长度取Lz=3×D＝16.8m
circle,9,Rzg,,,,        !创建主管截面圆形
lstr,9,10                !创建轴线L5(9,10)
KFILL,9,10,1,11,,        !在9，10两点间插入二分点11（中点）
ADRAG,1,2,3,4,,,5,,,,,,        !由圆沿L5拉伸为圆柱面

!*******(2.2.2)建立支（岔）管模型
LOCAL,13,1,,,Lz/2,,beita,,,        !创建编号为13局部柱坐标系，原点在11（0,0,Lz/2）号点处，与原来坐标系的夹角：沿X轴转60度，
K,12,,,Lc                !建立支管端部圆心点

lstr,11,12                !创建轴线L18(11,12)
K,13,Rcg,,Lc                !建立支管端部圆周上一点13
K,17,(Rcg+Lc*tan(sita)),,0                !建立支管锥管底部圆周上一点17
lstr,13,17                !创建支管腰线L18(11,12)
AROTAT,15,,,,,,12,11,360,4        !腰线绕轴线旋转360度生成支管曲面
K,21,,,lc+lcex                !建立闷头端部顶点

!*****进行修剪
APTN,ALL                !对全部曲面进行互分，areas partition
ADELE,9,9,0,1                !删除多余曲面9
ADELE,10,18,2,1                !删除多余曲面10，12，14，16，18
NUMCMP,AREA                !压缩曲面编号

!******(2.2.3)创建补强板
K,27,Rcg+Lbq*sin(beita)+2*Rcg*tan(sita),,Lc+2*Rcg        !建立与支管同轴的套管端部圆周上一点27(延长2Rcg)             !2*Rcg*tan(sita)
K,28,(Rcg+Lbq*sin(beita)),,0                !建立支管锥管底部圆周上一点17
lstr,27,28                !创建支管腰线L18(27,28)
AROTAT,12,,,,,,12,11,360,4        !腰线12绕轴线旋转360度生成支管曲面
!*****进行修剪
APTN,ALL                !对全部曲面进行互分，areas partition
ADELE,13,20,1,1                !删除多余曲面

NUMCMP,AREA                !压缩曲面编号

!******(2.2.4)支管延长段
lstr,12,21
ADRAG,23,21,19,25,,,12,,,,,,        !由圆沿L12拉伸为圆柱面……支管延长段

ASEL,ALL
AGLUE,all                !粘结闷头和岔管
NUMCMP,AREA                !压缩曲面编号
!**********(2.3)设置单元属性*******************************

ASEL,s,,,5,6,1,         !补强板的属性
AATT,3,3,3,,

ASEL,s,area,,1,2,1,        !主管的属性
AATT,1,1,1,,
ASEL,s,area,,7,8,1
AATT,1,1,1,,

ASEL,s,area,,3,4,1,        !支管的属性
AATT,2,2,2,,
ASEL,s,area,,13,14,1,
AATT,2,2,2,,

ASEL,s,area,,9,12,1        !支管延长段的属性
AREVERSE,all,1
AATT,2,2,2,,

ASEL,all
APLOT,ALL

!**********(2.4)划分网格*************
ASEL,ALL
SFA,ALL,1,PRES, 2.37E6,                !在所有面上施加面荷载
/PSF,PRES ,NORM,2,0,2           !用箭头显示表面压力

ASEL,ALL
AESIZE,ALL,size,        !划分网格单元的尺寸
!AMESH,ALL


ASEL,ALL
FINISH

我是新手，这是我自己编写的命令流，可能刚学，写起来有些罗嗦，请大家指点。跟大家一起交流学习！

[ 本帖最后由 zbk0615 于 2008-5-10 16:59 编辑 ]

heguojun66

新手都写这么多 不简单啊

carriesun9810

这些命令是从哪里看到的呀？我只会用对话框弄一点基础的东西

qiang52199

不错，还可以

suh26

不错，不错，很即的资料