%_FNCNAME%(0,2,1)= 0.0, -2, 0, 0.7462, 0, 0, -1 %_FNCNAME%(0,3,1)= 1, -3, 0, 1, -1, 3, -2
%_FNCNAME%(0,4,1)= 0.0, -1, 0, 227.2246, 0, 0, -3 %_FNCNAME%(0,5,1)= 0.0, -2, 0, 1, -3, 3, -1 %_FNCNAME%(0,6,1)= 0.0, -1, 0, 0.13, 0, 0, -2 %_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, -2, 3, -1 %_FNCNAME%(0,8,1)= 0.0, -1, 0, 0, 0, 0, 0 %_FNCNAME%(0,9,1)= 0.0, -2, 0, 1, 0, 0, -1 %_FNCNAME%(0,10,1)= 0.0, -4, 0, 1, -1, 2, -2 %_FNCNAME%(0,11,1)= 0.0, -1, 0, 0.7426, 0, 0, -4 %_FNCNAME%(0,12,1)= 0.0, -2, 0, 1, -4, 3, -1 %_FNCNAME%(0,13,1)= 0.0, -1, 0, 1, -2, 3, 1 %_FNCNAME%(0,14,1)= 0.0, -1, 7, 1, -1, 0, 0 %_FNCNAME%(0,15,1)= 0.0, -2, 0, 1, -3, 3, -1 %_FNCNAME%(0,16,1)= 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 2450*0.7462*227.2246*0.13*exp(-0.7426*{TIME})
!创建几何体
BLOCK,,5,,7.5,,2.4 !创建混凝土块 BLOCK,,20,,20,-10,0 !创建基础土体
KWPLAN,-1,4,3,6 !定位工作坐标系1 VSBW,ALL !分割体1
KWPLAN,-1,3,1,7 !定位工作坐标系2 VSBW,ALL !分割体2
VGLUE,ALL !GLUE所有体元素 NUMMRG,ALL !合并所有重复元素 NUMCMP, ALL !压缩所有元素的编号 CSYS,0 !返回全局坐标系
!进行网格划分
!划分混凝土体
VSEL,S,,,1,1,1 !选择混凝土体 VATT,1,,1,0 !指定混凝土材料 ALLSEL
VSEL,U,,,1,1,1 !选择基础土体
VATT,2,,1,0 !指定基础土体材料 ALLSEL !选择所有料 LSEL,S,,,5,5,1 !选择混凝土短边
LESIZE,ALL,,,5 !混凝土短边划分为5 LSEL,S,,,8,8,1 !选择混凝土长边
LESIZE,ALL,,,10 !混凝土短边划分为10 LSEL,S,,,9,9,1 !选择混凝土高
LESIZE,ALL,,,6 !混凝土高划分为6
!划分基础土体
LSEL,S,,,1,1,1 !选择混凝土高 LESIZE,ALL,,,10 LSEL,S,,,4,4,1 LESIZE,ALL,,,5 LSEL,S,,,15,15,1 LSEL,A,,,13,13,1 LSEL,A,,,24,39,15 LSEL,A,,,20,30,10 LSEL,A,,,37,40,3 LESIZE,ALL,,,6 LSEL,S,,,2,2,1 LESIZE,ALL,,,5 LSEL,S,,,3,3,1 LESIZE,ALL,,,10 ALLSEL
VMESH,ALL !映射划分网格 FINISH
/SOLU !进入求解器
ANTYPE,4 !选择瞬态求解 TRNOPT,FULL LUMPM,0
ESEL,S,MAT,,1 !选择混凝土体的所有结点 NSLE,S
IC,all,TEMP,23.6 !混凝土的浇筑温度为23.6 ESEL,S,MAT,,2 !选择基础土体的所有结点 NSLE,S
IC,all,TEMP,16 !基础的初始温度为16 ALLSEL
SFA,4,1,CONV,616.32, %QIWEN% !四周木模对流系数616.32 SFA,6,1,CONV,616.32, %QIWEN%
SFA,2,1,CONV,2054.4, %QIWEN% !上面裸露对流系数616.32 BFV,1,HGEN, %HEATS% !输入热生成率heats ALLSE !选择所有
TRNOPT,FULL !设置计算方法 *DO,DAY,0.5,10,0.5 !设置循环计算
TIME,DAY !设置荷载步时间 AUTOTS,0 !关闭自动时间步长 KBC,0 !设置阶梯荷载步
TSRES,ERASE !Defines an array of keytimes at which the time-stepping strategy changes
OUTRES,ALL,ALL !输出所有项 SOLVE !开始计算
*ENDDO !终止循环计算
!结构应力分析
/PREP7 !前处理
ETCHG,TTS !转换热分析单元为结构单元
!输入参数
!定义变化弹模EXX *DIM,EXX,ARRAY,20 EXX(1,1,1) = 0.87e10 EXX(2,1,1) = 1.05e10 EXX(3,1,1) = 1.17e10 EXX(4,1,1) = 1.26e10 EXX(5,1,1) = 1.33e10 EXX(6,1,1) = 1.39e10 EXX(7,1,1) = 1.44e10 EXX(8,1,1) = 1.48e10 EXX(9,1,1) = 1.52e10 EXX(10,1,1) = 1.56e10 EXX(11,1,1) = 1.59e10 EXX(12,1,1) = 1.61e10 EXX(13,1,1) = 1.64e10 EXX(14,1,1) = 1.67e10 EXX(15,1,1) = 1.69e10 EXX(16,1,1) = 1.71e10 EXX(17,1,1) = 1.73e10 EXX(18,1,1) = 1.75e10 EXX(19,1,1) = 1.77e10 EXX(20,1,1) = 1.79e10
!定义20种材料每0.5小时一种材料来反映弹模的改变 *DO,I,1,20,1 MP,DENS,I,2450 !混凝土 密度为2450 MP,EX,I,EXX(I) !混凝土 弹模为EXX(I) MP,PRXY,I,0.167 !混凝土 泊松比为0.167
MP,ALPX,I,0.75E-5 !混凝土 热膨胀系数为0.75e-5 MP,REFT,I,15 !定义参考温度 *ENDDO
VCLEAR, 2, 5, 1 !清除基础土体网格 VDELE, 2, 5, 1, 1 !删除基础土体 ALLSEL !选择选择所有
/SOLU !进入求解器
NROPT,FULL !Specifies the Newton-Raphson options in a full transient analysis
NLGEOM,ON !考虑大变形的影响 ASEL,S,,,1,1,1 !选择面1
NSLA,S,1 !选择面1上的结点 D,all,ALL !加基础全约束 ALLSEL !选择选择所有 DA,3,SYMM !加对称约束 DA,5,SYMM !加对称约束 I=1
*DO,DAY,0.5,10,0.5
ESEL, S, MAT, , I, I,1 !选择混凝土材料号为I的元素
ACEL,0,0,-9.8, !GRAVITY LOAD 施加重力荷载加速度 LDREAD,TEMP,,,DAY, ,'themal','rth',' ' !TEMP LOAD Reads results from the results file and applies them as loads. TIME,DAY AUTOTS,0
KBC,0 !设置阶梯荷载步
TSRES,ERASE !Defines an array of keytimes at which the time-stepping strategy changes
OUTRES,ALL,ALL, !输出所有的结果 SOLVE !计算 I=I+1
*IF,I,LE,20,THEN MPCHG,I,ALL *ENDIF *ENDDO
【间接法耦合分析温度场与应力场】
!!!!此为7月22号自编的不分层，体加30摄氏度、水化热生热率；一面加对流边界条件
/PREP7
/TITLE,Transient-state thermal analysis of concrete wall /SHOW
ET,1,SOLID90 !定义单元类型为thermal solid 90
!!!!!!!!/UNITS,SI !采用国际制单位(与时间有关的单位都换成天，国际制单位焦耳）
!!!!!!!!!!!!!!!!!!!!!!以下为定义材料属性:
MP,DENS,1,2441 !密度=2441kilogram/ ( meter3 )
MP,HF,1,14 !对流换热系数，单位：J/ ( meter2.s.degree C ) !!!!!!!!!!20J/ ( meter2.s.degree C ) =1680000J/( meter2.d.degree C ) MPTEMP,,0,21,32,43 !建立温度表 ，单位： ( degree C )
MPDATA,C,1,,971,971, 997,1038 !定义比热，单位：J / ( kilogram .degree MPDATA,KXX,1,,3.2158,3.2158, 3.1647,3.1217 !定义导热系数，单位：J/ ( meter.S.degree C ) !!!!!!!!!!!!!!!!!!MPDATA,KXX,1,,277845.12,277845.12, 273430,269714.8 !定义导热系数，单
位：J/ ( meter.d.degree C )
BLOCK,,12,,5,,7, !建立立方体模型
/VIEW,1,1,1,1 !Isometric view(三轴等侧视图）
!!!!!! 划分网格
MSHAPE,0,3D ！选择3D映射网格 MSHKEY,1
VSEL, ALL !选中所有的体 VMESH,ALL TYPE, 1
EXTOPT,ESIZE,2,0, EXTOPT,ACLEAR,0 !*
EXTOPT,ATTR,0,0,0 MAT,1 REAL,_Z4 ESYS,0 !*
FLST,2,1,5,ORDE,1 FITEM,2,4
VEXT,p51x, , ,0,3,0,,,, FLST,2,1,5,ORDE,1 FITEM,2,7
VEXT,p51x, , ,0,2,0,,,,
!!!!!! 网格划分完毕 eplot
finish !从前处理退出 !热分析加载求解
/solu !相当于Main Menu>Solution ANTYPE,trans
timint,off !先做稳态分析，确定初始条件
time,0.001 !设定只有一个子步的时间很小的荷载步 !!!!!!!!!!!!!!!!!deltim,0.001 !设定时间步长为0.001秒 !TUNIF,30, !加体温度荷载，30摄氏度 d,all,temp,0
!NSEL,S,LOC,X,12 !定义X=12表面上的对流边界条件 !SF,ALL,CONV,-1,20
ALLSEL
solve !相当于solution>solve>current LS !瞬态热分析
ddelet,all,temp !删除稳态分析时的节点温度荷载 !!!!!!!!先把2,3层单元全部杀死，再激活第一层