以梦为马_ck 2021-07-23 21:51 采纳率: 100%
浏览 755
已结题

MATLAB实现长方体容器内球形颗粒堆积

右手坐标系;容器的中心点为原点;
容器尺寸;球的体积分数;球的半径可以调整;
有过程的动态演示;
最后可以输出球的位置信息并保存。
PS:希望能写注释比较详细!

  • 写回答

3条回答 默认 最新

  • soar3033 2021-07-24 11:07
    关注

    动态绘图,并存储。请采纳

    img

    img

    img

    xl=200;%长
    yl=200;%宽
    zl=100;%高
    xmin=-xl/2;%计算x轴最负值,因为0是原点。下同
    ymin=-yl/2;
    zmin=-zl/2;
    r=10;%球半径
    n=40;%球数量
    for i = 1:n%循环生成初始位置
        if i==1%如果是第一个,直接随机生成
            x(i)=rand()*(xl-2*r)+xmin+r;
            y(i)=rand()*(yl-2*r)+ymin+r;
            z(i)=rand()*(zl-2*r)+zmin+r;
        else%如果不是第一个
            while 1%进入死循环,并生成位置
                x(i)=rand()*(xl-2*r)+xmin+r;
                y(i)=rand()*(yl-2*r)+ymin+r;
                z(i)=rand()*(zl-2*r)+zmin+r;
                st=0;%标记
                for j=1:i-1%遍历所有已生成的球心
                    if ((x(i)-x(j))^2+(y(i)-y(j))^2+(z(i)-z(j))^2)<(2*r)^2 %如果干涉
                        st=1;%更改标记
                        break    %跳出for循环
                    end    
                end
                if st==0%如果for循环完毕后,st依然为0,表明没有干涉
                    break%则推出死循环
                end
            end
        end
    end
    while 1%循环,用于绘制动态球下落
        st2=0;%标记,用于识别是否所有球都下落完
        for i= 1:n%遍历球
            st=0;%标记
            for j=1:n%遍历球
                if i~=j%如果ij不是一个球,则进行判断,如果两球干涉,且i球高于j球,则停止下落
                    if ((x(i)-x(j))^2+(y(i)-y(j))^2+(z(i)-z(j))^2)<(2*r)^2 && z(i)>z(j)
                        st=1;%修改下落标记
                        break    
                    end
                end
                if z(i)-zmin<r %如果i球已经到底面
                    st=1;%修改下落标记
                    break
                end    
            end
            if st==0%如果标记位0,则下落
                z(i)=z(i)-1;%下落计算,步长为1,可更改,越小越慢
                st2=1;%更改st2标记,表明有球下落,循环继续
            end
        end
        if st2==0%如果无球下落,则跳出循环
            break
        end
        for i =1:n%遍历所有球绘图
            [xx,yy,zz]=ellipsoid(x(i),y(i),z(i),r,r,r);%设定绘制的球心和半径
            surf(xx,yy,zz)%绘制球
            if i==1%当绘制第一个球的时候
                hold on%打开hold on,否则只显示一个球
                axis([xmin,xmin+xl,ymin,ymin+yl,zmin,zmin+zl])%设置坐标轴
            end
        end
        hold off%每次绘制完一轮球,关闭hold,以清屏
        pause(0.01)%延时
    end
    aa=[x;y;z];%把球的坐标输入到aa
    t=table(aa);%把aa变成表t
    tmp={'x','y','z'};%设置表头
    t.Properties.RowNames=tmp;%设置表头
    writetable(t,'1.xlsx','WriteRowNames',true);%输出到excel保存
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(2条)

报告相同问题?

问题事件

  • 系统已结题 8月2日
  • 已采纳回答 7月25日
  • 创建了问题 7月23日

悬赏问题

  • ¥15 lammps Gpu加速出错
  • ¥15 关于PLUS模型中kapaa值的问题
  • ¥15 关于博途V17进行仿真时无法建立连接问题
  • ¥15 请问下这个红框里面是什么文档或者记事本编辑器
  • ¥15 机器学习教材中的例题询问
  • ¥15 求.net core 几款免费的pdf编辑器
  • ¥15 为什么安装HCL 和virtualbox之后没有找到VirtualBoxHost-OnlyNetWork?
  • ¥15 C# P/Invoke的效率问题
  • ¥20 thinkphp适配人大金仓问题
  • ¥20 Oracle替换.dbf文件后无法连接,如何解决?(相关搜索:数据库|死循环)