左志华(在读硕士) zuo.zhihua@qq.com
@深圳市罗湖区 2021.1.12
ThinkPad S2 2018 Windows10 Home
关键词：SPH法；一维；启发；B样条核函数；Octave 6.1
1. 背景
从2020年1月初至5月初，我在老家呆了已经5个月了，我家人也在3月初外出返工，我一个人呆着确实无聊，还需要自己做饭。趁着我女朋友回海南学校准备毕业，我陪她去海南了，没想到学校还是开学遥遥无期，在外消费挺高的，就在6月3号回家了。
记得在5月底的时候，学校开始了硕博连读工作，我这样的成绩是没戏的(学院学硕倒数第二，年级倒数第九)。我渴望学术研究领域取得自己的成果，同时，我爱我的女朋友，读博的日子会很漫长。在2020年12月，我以年级第三评上了硕博连读，希望自己能变得更强(无关内卷)，希望能守护好我爱的人。
博士期间，我将尝试基于SPH法对气垫船的运动进行仿真研究，更普适地研究一套我们学术团队的现代船舶运动仿真技术。这个还是很有挑战性的。在5月31号，我的“傻弟弟”室友不太会编程，有一道SPH习题不会做，当时我给他写了一份Matlab代码。
如今正式自己学SPH时发现，网上这方面的基础资料还是挺匮乏的，所以找到了这份代码，希望对初学者有启发效果。
2. 题目
拉格朗日型无网格粒子法在描述结构毁伤行为时具有先天的优势。假定单位长度细长梁结构如下图所示，梁的截面内被离散为单个粒子，在初始状态下梁内变量
，请选用三次样条核函数，光滑长度
，通过SPH方法中的粒子近似原理编程求解每个粒子处的变量值。
(1)粒子分辨率
和
