求助！热应力分析中的几何图形
浏览：117692 回答：1
本人刚刚接触ANSYS，按照书上的例子，适当修改，用于我的仿真，很多东西不是很懂。因不知如何画实际几何模型用于热应力仿真，所以采用了圆环体，代替实际几何形状。求大家帮忙修改下我的仿真程序，以实现真实的几何模型仿真。要求将目前仿真中的圆环形几何模型修改为实际几何模型，修改后请先帮忙验证下是否可以正常仿真，非常感谢!!
目前仿真程序中 希望修改为
模型1 圆柱体（空心） 外圆内方（内方的尺寸与圆环体内径相同）
模型2 圆环体（空心） 矩形环（矩形<圆环外径>-矩形<圆环内径>）
模型3 圆环体（空心） 矩形
模型4 圆环体 外圆内方（内方的尺寸与圆环体内径相同）
实际几何形状如下图：
具体参数如下面程序中。
FINISH
/FILNAME,Exercise19 ! 定义分析文件名
! 第一步：进行温度场分析的前处理并写温度场物理分析文件
/prep7 ! 进入前处理器
et,1,plane77,,,1 ! 选择PLANE77热分析单元并设置为轴对称分析
mp,kxx,1,1.09 ! 定义1热传导系数
mp,kxx,2,1.5 ! 定义2热传导系数
mp,kxx,3,1.4 ! 定义3热传导系数
mp,kxx,4,15.6 ! 定义4热传导系数
rectng,0,.04,0,.02 ! 建立1几何模型
rectng,.04,.05,0,.02 ! 建立2几何模型
rectng,.05,.14,0,.02 !建立3几何模型
rectng,.02,.14,.02,.03 !建立4几何模型
aglue, all ! 粘接各矩形
numcmp,area ! 压缩面编号
asel,s,,,1
aatt,1,1,1 ! 附于1材料属性
asel,s,,,2
aatt,2,1,1 ! 附于2材料属性
asel,s,,,3
aatt,3,1,1 ! 附于3材料属性
asel,s,,,4
aatt,4,1,1 ! 附于4材料属性
asel,all
esize,.001 ! 定义单元划分尺寸
mshkey,2 ! 设置为映射单元划分类型
amesh,all ! 划分单元
nsel,s,loc,x,0
d,all,temp,37 ! 施加内壁温度边界条件
nsel,s,loc,x,.14
d,all,temp,37 ! 施加外壁温度边界条件
nsel,all
physics,write,thermal ! 写温度场物理分析文件
! 第二步：进行结构场分析的前处理并写结构场物理分析文件
physics,clear ! 清空物理环境数据
ddel,all ! 删除温度场温度载荷
et,1,82,,,1 ! 选择结构分析单元
mp,ex,1,9.1e6 ! 定义1结构场材料属性
mp,alpx,1,3.25e-6
mp,nuxy,1,.2
mp,ex,2,4.93e5 ! 定义2结构场材料属性
mp,alpx,2,54e-6
mp,nuxy,2,.38
mp,ex,3,4.35e5 ! 定义3结构场材料属性
mp,alpx,3,15e-6
mp,nuxy,3,.28
mp,ex,4,15.95e6 ! 定义4结构场材料属性
mp,alpx,4,8.9e-6
mp,nuxy,4,.36
nsel,s,loc,y,.03 ! 选择顶面节点
cp,1,uy,all ! 耦合节点Y向自由度
nsel,s,loc,x,0 ! 选择内壁节点
cp,2,ux,all ! 耦合节点X向自由度
nsel,s,loc,y,0 ! 选择底面节点
d,all,uy,0 ! 施加Y向位移约束
nsel,all
tref,37 ! 定义参考温度
physics,write,struct ! 写结构场物理分析文件
save ! 存盘
finish
! 第三步：读取温度场物理分析文件进行求解和后处理
/solu
physics,read,thermal ! 读取温度场物理分析文件
solve ! 求解
finish
/post1 ! 进入通用后处理器
path,radial,2 ! 定义径向显示路径
ppath,1,,0
ppath,2,,.14
pdef,temp,temp ! 向所定义路径映射温度分析结果
pasave,radial,filea ! 保存路径文件
plpath,temp ! 显示沿路径温度变化曲线图
finish
! 第四步：读取结构场物理分析文件并读取温度场计算结果进行结构场求解和后处理
/solu
physics,read,struct ! 读取结构场物理分析文件
ldread,temp,,,,,,rth ! 读取温度场分析结果
solve ! 求解
finish
/post1 ! 进入通用后处理器
paresu,raidal,filea ! 读取已存路径文件
pmap,,mat ! 定义沿路径不连续区域处理方法
pdef,sx,s,x ! 向所定义路径映射X向应力分析结果
pdef,sz,s,z ! 向所定义路径映射Z向应力分析结果
plpath,sx,sz ! 显示沿路径X和Z向应力分析曲线
plpagm,sx,,node ! 在几何模型上显示X向应力分布云图
plnsol,s,eqv,0,1 ! 显示等效应力分析结果
/expand,27,axis,,,10 ! 设置轴对称分析结果扩展选项
plnsol,s,eqv,0,1 !显示三维扩展的等效应力分析结果
目前仿真程序中 希望修改为
模型1 圆柱体（空心） 外圆内方（内方的尺寸与圆环体内径相同）
模型2 圆环体（空心） 矩形环（矩形<圆环外径>-矩形<圆环内径>）
模型3 圆环体（空心） 矩形
模型4 圆环体 外圆内方（内方的尺寸与圆环体内径相同）
实际几何形状如下图：
具体参数如下面程序中。
FINISH
/FILNAME,Exercise19 ! 定义分析文件名
! 第一步：进行温度场分析的前处理并写温度场物理分析文件
/prep7 ! 进入前处理器
et,1,plane77,,,1 ! 选择PLANE77热分析单元并设置为轴对称分析
mp,kxx,1,1.09 ! 定义1热传导系数
mp,kxx,2,1.5 ! 定义2热传导系数
mp,kxx,3,1.4 ! 定义3热传导系数
mp,kxx,4,15.6 ! 定义4热传导系数
rectng,0,.04,0,.02 ! 建立1几何模型
rectng,.04,.05,0,.02 ! 建立2几何模型
rectng,.05,.14,0,.02 !建立3几何模型
rectng,.02,.14,.02,.03 !建立4几何模型
aglue, all ! 粘接各矩形
numcmp,area ! 压缩面编号
asel,s,,,1
aatt,1,1,1 ! 附于1材料属性
asel,s,,,2
aatt,2,1,1 ! 附于2材料属性
asel,s,,,3
aatt,3,1,1 ! 附于3材料属性
asel,s,,,4
aatt,4,1,1 ! 附于4材料属性
asel,all
esize,.001 ! 定义单元划分尺寸
mshkey,2 ! 设置为映射单元划分类型
amesh,all ! 划分单元
nsel,s,loc,x,0
d,all,temp,37 ! 施加内壁温度边界条件
nsel,s,loc,x,.14
d,all,temp,37 ! 施加外壁温度边界条件
nsel,all
physics,write,thermal ! 写温度场物理分析文件
! 第二步：进行结构场分析的前处理并写结构场物理分析文件
physics,clear ! 清空物理环境数据
ddel,all ! 删除温度场温度载荷
et,1,82,,,1 ! 选择结构分析单元
mp,ex,1,9.1e6 ! 定义1结构场材料属性
mp,alpx,1,3.25e-6
mp,nuxy,1,.2
mp,ex,2,4.93e5 ! 定义2结构场材料属性
mp,alpx,2,54e-6
mp,nuxy,2,.38
mp,ex,3,4.35e5 ! 定义3结构场材料属性
mp,alpx,3,15e-6
mp,nuxy,3,.28
mp,ex,4,15.95e6 ! 定义4结构场材料属性
mp,alpx,4,8.9e-6
mp,nuxy,4,.36
nsel,s,loc,y,.03 ! 选择顶面节点
cp,1,uy,all ! 耦合节点Y向自由度
nsel,s,loc,x,0 ! 选择内壁节点
cp,2,ux,all ! 耦合节点X向自由度
nsel,s,loc,y,0 ! 选择底面节点
d,all,uy,0 ! 施加Y向位移约束
nsel,all
tref,37 ! 定义参考温度
physics,write,struct ! 写结构场物理分析文件
save ! 存盘
finish
! 第三步：读取温度场物理分析文件进行求解和后处理
/solu
physics,read,thermal ! 读取温度场物理分析文件
solve ! 求解
finish
/post1 ! 进入通用后处理器
path,radial,2 ! 定义径向显示路径
ppath,1,,0
ppath,2,,.14
pdef,temp,temp ! 向所定义路径映射温度分析结果
pasave,radial,filea ! 保存路径文件
plpath,temp ! 显示沿路径温度变化曲线图
finish
! 第四步：读取结构场物理分析文件并读取温度场计算结果进行结构场求解和后处理
/solu
physics,read,struct ! 读取结构场物理分析文件
ldread,temp,,,,,,rth ! 读取温度场分析结果
solve ! 求解
finish
/post1 ! 进入通用后处理器
paresu,raidal,filea ! 读取已存路径文件
pmap,,mat ! 定义沿路径不连续区域处理方法
pdef,sx,s,x ! 向所定义路径映射X向应力分析结果
pdef,sz,s,z ! 向所定义路径映射Z向应力分析结果
plpath,sx,sz ! 显示沿路径X和Z向应力分析曲线
plpagm,sx,,node ! 在几何模型上显示X向应力分布云图
plnsol,s,eqv,0,1 ! 显示等效应力分析结果
/expand,27,axis,,,10 ! 设置轴对称分析结果扩展选项
plnsol,s,eqv,0,1 !显示三维扩展的等效应力分析结果
目前仿真程序中 希望修改为
模型1 圆柱体（空心） 外圆内方（内方的尺寸与圆环体内径相同）
先做一个圆，再做内方，然后对二者执行 AOVLAP 操作，菜单路径是：
Main Menu>Preprocessor>Modeling>Operate>Booleans>Overlap>Areas
模型2 圆环体（空心） 矩形环（矩形<圆环外径>-矩形<圆环内径>）
先做两个矩形，然后对二者执行 AOVLAP 操作。
模型3 圆环体（空心） 矩形
如果只有一个矩形，直接做就可以；如果是空心矩形，与模型 2 类似，丹巴 AOVLP 操作改为用大面减去小面 - ASBA 命令，菜单路径：
Main Menu>Preprocessor>Modeling>Operate>Booleans>Subtract>Areas
模型4 圆环体 外圆内方（内方的尺寸与圆环体内径相同）
同模型 1