Ansys apdl 结合生死单元仿真激光烧蚀纤维复合材料,表面温度过高? 200

考虑相变(直接汽化),手动添加了大概的焓值,表面基质相汽化温度3460°C,纤维相汽化温度4000°C,模型如图
加载激光热流密度,生死单元求解,网格最细划到了50等分(100等分不求解,不知道为啥)
热流密度大的话,基质和纤维烧了一部分,但温度最高达到了10000°C
热流密度小的话,只烧了基质,但表面温度最高还是有5000°C,按理来说高于汽化温度的不是应该都杀死了吗?
是我两种不同材料的生死单元设置的有问题吗?生死单元部分apdl代码如下

etable,temperature,temp   !将单元温度存储在temperature表中
esel,s,MAT,,1                                                                                                                !选择材料1的单元
esel,r,etab,temperature,3460                                                !从材料1单元中选择温度大于汽化温度的单元                                                                                             
cm,vapor1,ELEM                                                                                                        !创建基质烧蚀单元集1

esel,s,MAT,,2
esel,r,etab,temperature,4000
cm,vapor2,ELEM
etable,temperature,ERAS                 !删除temperature表
cmsel,all           
finish              !退出当前模块
/solu               !进入求解模块
*get,elemnum,elem,0,count      !读取单元总数,赋给参数elemnum
                *if,elemnum,gt,0,then    !如果单元数量大于0
           ekill,all              !杀死选择的单元
           esel,s,live          !选择活单元
     *else                       !如果单元数量小于0(没有发生烧蚀)
          esel,s,live           !选择活单元
     *endif
nsle,s              !选择活单元上的节点
nsel,inve         !选择死单元上的节点 D,all,all,0        !约束死节点自由度
nsel,inve         !选择活单元上的节点

单纤维.png

大热流密度.png

小热流密度.png

Ansys apdl 结合生死单元仿真激光烧蚀纤维复合材料,表面温度过高?的图4

邀请回答 举报

共1个回答

默认排序
  • color.
    2023年2月15日
    你好,不知阁下的问题是否已经解决,我认为问题在于ETABLE命令默认情况下储存的是单元结果数据的平均值而非最高值,因此部分单元会出现温度高于设定的汽化温度,在命令ETABLE, Lab, Item, Comp, Option中的"Option"处填写“MAX”使ETABLE命令储存单元结果数据最高值,即改为etable,temperature,temp,MAX,希望对你有帮助
    0评论 采纳答案 举报

  • Ansys apdl 结合生死单元仿真激光烧蚀纤维复合材料,表面温度过高?的相关案例教程

    一 前处理 ! 定义接触的name selection,插入command在contact里 ml_pair=cid ml_pair_1=cid !选择有 name selection 的单元或节点 cmsel,s,xx,element !选择NS为xx的包含的element nsel esel 二 求解设置 !设置迭代次数,默认最大迭代次数为26 NEQIT,40 !设置最大迭代次数为40 三
    李桂花 上海安世亚太公司 钢筋混凝土结构广泛应用于桥梁、建筑和外墙等领域。在仿真当中如何考虑混凝土结构中的加强筋作用呢? 第一种方法,直观上使用实体对混凝土和钢筋进行完全几何建模。在这种方法中,钢筋所在的混凝土体有空隙,钢筋本身填充这些空隙。这种方法钢筋和混凝土之间的网格共节点(或接触关系)比较好处理,但只适用于小型模型的建模,且划分网格不方便,计算量大,工程上较少采用。 第二种方法,工程上一般将
    作者:李桂花,上海安世亚太结构应用工程师 来源:本文为上海安世亚太原创作品 前言 钢筋混凝土结构广泛应用于桥梁、建筑和外墙等领域。在仿真当中如何考虑混凝土结构中的加强筋作用呢? 第一种方法,直观上使用实体对混凝土和钢筋进行完全几何建模。在这种方法中,钢筋所在的混凝土体有空隙,钢筋本身填充这些空隙。这种方法钢筋和混凝土之间的网格共节点(或接触关系)比较好处理,但只适用于小型模型的建模,且划分网格不方
    1.sf,nlist,label,value,value2 -“nilst”是节点列表,也可以是命名选择 -辐射标签是rdsf -value是表面发射率 -value2是封闭体数量 2.spctemp命令行:因为所计算的空间不是完全封闭的计算空间,所以必须定义空间温度, spctemp,number,temperature spctemp是ansys定义空间温度的关键字,number是非封闭空间的
    该命令流为计算单转子-支承系统在加速运行过程中,受质量不平衡激励下的瞬态动力学响应。可以准确计算出在共振转速下的峰值及弯曲应变能情况。给出了详细的表加载转速和不平衡力的方法,可供参考。 /prep7 ! ** parameters length = 0.4 ro_shaft = 0.01 ro_disk = 0.15 md = 16.47 id = 9.427e-2 ip = 0.1861 kxx