热应力耦合分析-钢化玻璃预应力仿真分析

林德佳

各位，有个问题想要问问大家，关于钢化玻璃冷却过程的仿真分析，具体参考文献《平板玻璃钢化有限元数值模拟》和《Finite element implementation of a glass tempering model in three dimensions》，但仿真结果一直对不上。
先进行温度响应分析，得到的玻璃中心线内外温度变化曲线还行，最大温差在150度左右，但在应力分析响应中，得到的仿真结果，首先是数量级不对及分布不对，所以发帖向大佬们请教，如图所示，为平板玻璃在200s时刻的x，y，z向应力分布。
主要的疑问有
1. 热应力响应仿真分析中，分析步类型应该为coupled-temp-displacement，还是static，general呢？两种分析类型，本人都尝试过，结果有差异，但与文献结果差距较大；
2. 平板玻璃由650度冷却到室温20度，其中，涉及到凝固模拟，试试问，怎么表现这个由液态转变为固体呢？我是在设置热膨胀系数那块设置了，reference temperature，不知道这样设置有问问题，因为固态液态对应的膨胀系数不一致；
3. 关于粘弹性prony级数设置，有无设置改选项非常关键，外文文献提供相关参数，但中文没有具体说明，所以有点疑惑；
4.或者哪位大佬，看到该问题，有一些说明想法，可以帮忙指点一下，问题纠结了很久很久，一直没有进展，所以跟大家讨论讨论，谢谢各位！

林德佳

如图所示为，英文文献的仿真结果，有一点疑问的是，theta方向的应力具体是什么？因为是1/16模型的原因？
我的理解是，应该在x，y，z方向都是有这样共同的规律，空气冷却段为压应力，中心为拉应力，不知道是这样，但这样的文献结果，使得有点小小的疑问。
另外如图，为中文的仿真分析结果

林德佳

咋还没有大佬过来帮忙指点呢？大佬们，快点出现吧，哈哈！

757659093

边界条件给的都是合理的？？模型也不一样 试试同样的模型吧 取对称的结构 热分析虽然做的不多 但无非就是条件 网格 和载荷步的顺序时间

林德佳

757659093 发表于 2017-11-10 14:04
边界条件给的都是合理的？？模型也不一样 试试同样的模型吧 取对称的结构 热分析虽然做的不多 但无非就是条 ...

首先，非常感谢你的回答！
边界条件，没有问题的。
采用对称性建模，对称性边界条件已施加，另外，为冷却模拟过程，采用的对流换热，这几点应该没有问题。你说的模型不一样，我再试试看。
关于热应力耦合分析，网格，时间步长，在热传导及应力分析两个分析步都是保持一致的。我再试试看，你想想会不会哪里很容易出错的？在热应力分析那块，因为热传导的分析结果还可以，就是热应力仿真分析结果差别太大太大了。
谢谢，期待你的回复！