CA法模拟动态再结晶,晶粒正常长大,利用元胞自动机生成拓扑晶粒模型,参数可调

CA法模拟动态再结晶,晶粒正常长大,利用元胞自动机生成拓扑晶粒模型,参数可调。
元胞胞自动机模拟动态再结晶母相晶粒生成。
注:需要程序和完整版程序的可以加好友。
本程序基于曲率驱动机制以及热激活机制,matlab编写,本程序模拟奥氏体晶粒正常长大过程。
程序均有注释,仅作学习交流使用
请添加图片描述

ID:5968676443031832Love辽阔

请添加图片描述
请添加图片描述

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
由于动态再结晶过程比较复杂,我将为你提供一个简化版本的晶粒长大元胞自动机模型,代码如下: ```python import numpy as np import matplotlib.pyplot as plt # 设置模型参数 n = 100 # 空间尺寸 m = 10 # 晶粒数量 k = 4 # 晶粒长大速率 steps = 50 # 模拟步数 # 初始化晶粒分布 grain_map = np.zeros((n, n)) for i in range(m): x, y = np.random.randint(0, n), np.random.randint(0, n) grain_map[x, y] = i + 1 # 定义元胞自动机模型 def evolve(grain_map, k): new_grain_map = np.zeros_like(grain_map) for i in range(n): for j in range(n): if grain_map[i, j] == 0: neighbors = [] for ii in range(i-1, i+2): for jj in range(j-1, j+2): if ii == i and jj == j: continue if ii < 0 or jj < 0 or ii >= n or jj >= n: continue if grain_map[ii, jj] != 0: neighbors.append(grain_map[ii, jj]) if len(neighbors) > 0: counts = np.bincount(neighbors) max_count = np.max(counts) max_grains = np.where(counts == max_count)[0] new_grain_map[i, j] = np.random.choice(max_grains) else: new_grain_map[i, j] = grain_map[i, j] if np.random.rand() < k: neighbors = [] for ii in range(i-1, i+2): for jj in range(j-1, j+2): if ii == i and jj == j: continue if ii < 0 or jj < 0 or ii >= n or jj >= n: continue if grain_map[ii, jj] != 0: neighbors.append(grain_map[ii, jj]) if len(neighbors) > 0: counts = np.bincount(neighbors) max_count = np.max(counts) max_grains = np.where(counts == max_count)[0] new_grain_map[i, j] = np.random.choice(max_grains) return new_grain_map # 运行模拟 for i in range(steps): grain_map = evolve(grain_map, k) # 可视化结果 plt.figure(figsize=(8, 8)) plt.imshow(grain_map, cmap='jet') plt.axis('off') plt.show() ``` 这个模型中,我们假设晶粒在周围的高温区域中长大。在每个时间步长中,我们首先计算出每个空白细胞周围的晶粒数量,然后将其赋值为周围晶粒数量最多的晶粒之一。然后,对于每个已存在的晶粒,我们以速率$k$使其长大,并在周围的空白细胞中随机选择一个晶粒进行复制。 这个模型只是一个简单的例子,如果你想要更加真实的动态再结晶过程模拟,你需要更多的实验数据和更复杂的模型

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值