- 积分
- 0
- 注册时间
- 2008-3-14
- 仿真币
-
- 最后登录
- 1970-1-1
|
我要计算的是一个水电站的钢岔管的模型,还需要考虑钢管与外围混凝土的接触问题,主要参数如下:
主管直径: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 编辑 ] |
评分
-
1
查看全部评分
-
|