zoukankan      html  css  js  c++  java
  • Ansys热应力计算

    问题说明

    样块上下两端固定,在室温20℃下进行夹紧,分析其升温到150℃时的热应力。
    采用间接法进行分析,温度场单元选择278,应力场单元为185
    首先进行稳态热分析,然后读取温度场的数据,进行结构分析,具体的编程代码如下

    温度场分析APDL

    /CLEAR
    /FILENAME TEST
    /PREP7
    
    !定义温度场实体单元
    ET,1,SOLID278  
    !定义材料1的导热系数
    MPTEMP,,,,,,,,  
    MPTEMP,1,0  
    MPDATA,KXX,1,,20
    
    !建立样件模型
    BLC5,0,0,0.2,0.1,0.05  
    
    !设置体1的网格尺寸
    CM,_Y,VOLU  
    VSEL, , , ,       1
    CM,_Y1,VOLU
    CMSEL,S,_Y  
    CMSEL,S,_Y1
    VATT,       1, ,   1,       0   
    CMSEL,S,_Y  
    CMDELE,_Y   
    CMDELE,_Y1  
    
    FLST,5,1,4,ORDE,1   
    FITEM,5,5   
    CM,_Y,LINE  
    LSEL, , , ,P51X
    CM,_Y1,LINE
    CMSEL,,_Y   
    LESIZE,_Y1, , ,20, , , , ,1
     
    FLST,5,1,4,ORDE,1   
    FITEM,5,6   
    CM,_Y,LINE  
    LSEL, , , ,P51X
    CM,_Y1,LINE
    CMSEL,,_Y   
    LESIZE,_Y1, , ,10, , , , ,1
      
    FLST,5,1,4,ORDE,1   
    FITEM,5,10  
    CM,_Y,LINE  
    LSEL, , , ,P51X
    CM,_Y1,LINE
    CMSEL,,_Y   
    LESIZE,_Y1, , ,5, , , , ,1  
    
    CM,_Y,VOLU  
    VSEL, , , ,       1
    CM,_Y1,VOLU
    CHKMSH,'VOLU'   
    CMSEL,S,_Y  
    VCLEAR,_Y1  
    VMESH,_Y1   
    CMDELE,_Y   
    CMDELE,_Y1  
    CMDELE,_Y2  
    
    
    !在一侧施加温度150
    FLST,2,1,5,ORDE,1   
    FITEM,2,5   
    /GO
    DA,P51X,ALL,150
    
    !设置其他表面的对流换热系数
    FLST,2,5,5,ORDE,3   
    FITEM,2,1   
    FITEM,2,-4  
    FITEM,2,6   
    /GO
    SFA,P51X,1,CONV,50,
    
    !设置分析类型为稳态求解
    FINISH  
    /SOL
    ANTYPE,0
    
    !设置温度偏移
    TOFFST,273,
    
    !写温度场物理分析文件
    FINISH  
    /PREP7  
    PHYSICS,WRITE,Thermal1, , ,
    
    !清空物理环境数据
    PHYSICS,CLEAR  
    !删除边界和载荷 
    LSCLEAR,ALL
    
    !转换单元类型,thermal to struc
    ETCHG,TTS
    
    !设置材料参数
    MPTEMP,,,,,,,,  
    MPTEMP,1,0  
    MPDATA,EX,1,,2.1E+011   !弹性模量
    MPDATA,PRXY,1,,0.3       !泊松比
    UIMP,1,REFT,,,              !参考温度
    MPDATA,ALPX,1,,2.35E-5      !热膨胀系数
    
    !定义结构边界条件--左右固定
    
    FLST,2,2,5,ORDE,2   
    FITEM,2,5   
    FITEM,2,-6  
    /GO
    DA,P51X,ALL,0
    
    !设置参考温度
    TREF,20,
    
    !求解设置
    FINISH  
    /SOL
    ANTYPE,0
    NLGEOM,1
    NSUBST,8,10,6   
    TIME,1
    
    !写结构物理分析文件
    FINISH  
    /PREP7  
    PHYSICS,WRITE,struct1, , ,  
    
    !读温度场物理分析文件
    PHYSICS,READ,THERMAL1  
    
    !求解
    FINISH  
    /SOL
    /STATUS,SOLU
    SOLVE  
    
    !设置图像显示
    /PLOPTS,INFO,1  !inf0为legend on
    
    !显示温度云图
    FINISH  
    /POST1  
    /EFACET,1   
    PLNSOL, TEMP,, 0
    
    !读取结构场物理分析文件
    FINISH  
    /PREP7  
    PHYSICS,READ,STRUCT1
    
    !定义温度边界条件
    LDREAD,TEMP,,, , ,'test_temp1','rth',' '
    
    !求解
    FINISH  
    /SOL
    /STATUS,SOLU
    SOLVE 
    
    !查看结果
    FINISH  
    /POST1  
    SET,LIST,999
    SET,FIRST
    /EFACET,1   
    PLNSOL, S,EQV, 0,1.0  
    
    

    结果

  • 相关阅读:
    [BZOJ2809][Apio2012]dispatching
    [BZOJ4584][Apio2016]赛艇
    [BZOJ3206][Apio2013]道路费用
    [codeforces551E]GukiZ and GukiZiana
    [BZOJ3809]Gty的二逼妹子序列
    [BZOJ3289]Mato的文件管理
    [BZOJ3052][UOJ#58][WC2013]糖果公园
    [SPOJ10707]Count on a tree II
    [BZOJ1086][SCOI2005]王室联邦
    python小知识
  • 原文地址:https://www.cnblogs.com/gaozihan/p/10158768.html
Copyright © 2011-2022 走看看