三维任意形状随机骨料matlab生成：基于映射网格的属性赋予方法

三维任意形状随机骨料matlab生成基于映射网格背景网格法进行骨料、砂浆、界面属性的赋予。在混凝土材料细观建模中生成真实的三维骨料结构是个技术活。今天咱们就来聊聊如何用Matlab造出形状各异的随机骨料再通过背景网格法给不同材料区域赋予属性。整个过程就像搭积木——先造骨料再上色分区。首先生成随机骨料的核心是空间碰撞检测。这里采用球体基元叠加变形的思路% 生成基础球体骨料 for i 1:n_agg overlap true; while overlap agg(i).center rand(1,3)*domain_size; agg(i).radius 0.8 rand*0.4; % 控制粒径范围 % 碰撞检测 for j 1:i-1 dist norm(agg(i).center - agg(j).center); if dist (agg(i).radius agg(j).radius)*1.2 overlap true; break; else overlap false; end end end end这段代码通过循环检测新骨料与已有骨料的间距确保不发生重叠。实际工程中可以调整1.2这个安全系数来控制堆积密度。接下来给骨料捏造型引入球谐函数变形% 球面坐标变形 theta linspace(0, pi, 20); phi linspace(0, 2*pi, 40); [THETA, PHI] meshgrid(theta, phi); % 添加随机扰动 sh_degree 3; % 控制形状复杂度 coeffs randn(sh_degree1, sh_degree1)*0.1; radius_variation real(sphHarm(coeffs, THETA, PHI)); % 生成变形表面 [X,Y,Z] sph2cart(THETA, PHI, 1 radius_variation);通过调整球谐系数可以产生凹凸不平的自然形状比单纯用球体更接近真实骨料形貌。三维任意形状随机骨料matlab生成基于映射网格背景网格法进行骨料、砂浆、界面属性的赋予。现在要处理材料分区。背景网格法的精髓在于空间映射——先建立背景网格坐标系再判断每个网格点的归属% 建立背景网格 grid_size 2; % 网格尺寸 x_grid 0:grid_size:domain_size(1); y_grid 0:grid_size:domain_size(2); z_grid 0:grid_size:domain_size(3); % 初始化材料矩阵 material_grid ones(length(x_grid)-1, length(y_grid)-1, length(z_grid)-1)*2; % 2代表砂浆基质 % 骨料区域标记 for k 1:length(z_grid)-1 for j 1:length(y_grid)-1 for i 1:length(x_grid)-1 centroid [mean(x_grid(i:i1)), mean(y_grid(j:j1)), mean(z_grid(k:k1))]; % 判断是否在骨料内部 for a 1:n_agg if norm(centroid - agg(a).center) agg(a).radius material_grid(i,j,k) 1; % 1代表骨料 break end end % 界面层检测 if material_grid(i,j,k) 2 for a 1:n_agg distance norm(centroid - agg(a).center) - agg(a).radius; if abs(distance) interface_thickness material_grid(i,j,k) 3; % 3代表界面过渡区 break end end end end end end这个双重循环虽然看着笨但胜在直观。实际应用时可以优化空间查询算法比如用八叉树加速邻近搜索。最后是材料属性绑定。不同区域赋予各自的力学参数% 创建材料属性结构体 materials(1).name Aggregate; materials(1).E 70e9; // 弹性模量 materials(1).nu 0.24; // 泊松比 materials(2).name Mortar; materials(2).E 25e9; materials(2).nu 0.2; materials(3).name Interface; materials(3).E 15e9; materials(3).nu 0.25; % 生成VTK文件便于可视化 vtkwrite(model.vtk, structured_grid, Xg, Yg, Zg, scalars, MaterialID, material_grid);输出VTK文件后用Paraview这类工具可以直接看到彩色分区的三维模型检查骨料分布是否合理。几个实用小技巧骨料生成阶段用KD树优化空间搜索百万级颗粒也能hold住界面层厚度建议取骨料平均粒径的5-10%映射网格尺寸不要小于最小骨料直径的1/3并行计算加速把空间分块用parfor处理各子域这种方法虽然原理简单但在细观力学模拟中足够实用。下次可以试试把球形基元换成多面体或者引入纤维增强相让模型更贴近真实复合材料结构。