导读：9月14日19时，2022Mechsoft大讲堂将邀请首席岩土工程师和数值计算专家Abouzar带来《计算机辅助岩石工程(CARE)简介与案例研究》公开讲座。本期报告在仿真秀官网和APP同时直播并支持回看。2022Mechsoft讲堂：节理岩体力学特性尺寸效应与各向异性研究-仿真秀直播一、主讲嘉宾武娜，大连理工大学博士后。2020年毕业于大连理工大学岩土工程专业，获工学博士学位，师从梁正召教授。主要从事节理岩体破坏机理、力学性质与本构模型等方面的研究工作。主持国家自然科学基金1项，开放基金1项，参与国家自然科学基金5项。二、主讲内容工程岩体中通常含有大量随机分布的节理，节理的存在不仅控制着岩体的变形和破坏，而且使岩体的力学性质具有明显的尺寸效应与各向异性。目前，对工程岩体稳定性进行分析时很难考虑每一条节理对计算结果的影响而进行大规模工程计算，这与工程实际存在一定的差异。通过开展节理岩体尺寸效应与各向异性的研究，可获得节理岩体表征单元体（RepresentativeElementVolume，REV）和等效力学参数。因此，开展节理岩体力学特性尺寸效应与各向异性研究，对岩体工程的变形计算和稳定性分析具有重要意义。三、适看人群1、岩土工程方向CAE工程师2、学习型仿真工程师3、理工科院校教师和学生4、RFPA软件用户和学习者5、有限元软件兴趣爱好者和应用者6、基于微震监测的岩体工程稳定性分析7、从事岩体工程、深部地下工程科研、施工与管理工作人员四、如何观看2022Mechsoft讲堂：节理岩体力学特性尺寸效应与各向异性研究-仿真秀直播（完）声明：本文首发仿真秀App，部分图片和内容转自网络，如有不当请联系我们，欢迎分享，禁止私自转载，转载请联系我们。来源：仿真秀App

1.案例描述确定以速率Q产生热量的裸钢丝的中心线温度TcL和表面温度Ts。金属丝与空气之间的表面对流系数（温度Ta）为h。另外，确定散热率q。2.材料和几何参数3.边界条件选择整个体，添加InternalHeatGeneration，Magnitude为1.7894e-002；选择右侧边，添加Convection，FilmCoefficient为9.6451e-006，AmbientTemperature为70.4.查看结果添加两个Temperature，分别选择左右两边；在Probe中添加Reaction，BoundaryCondition选择Convection来源：CAE中学生

继续讲解！不知不觉中来到了2.4.6章节。盛夏也已到来！2023年7月，30度以上的温度已经是平常事，安心在家备课！学无止境，不断追求教学方法上的进步。最关键的原则就是因人而异，每一届的学生都有各自的特点。总体的特点就是课后不复习，很少练习。这个让我很头大。2.4.6部分特殊二维作图函数汇总bar––––绘制直方图polar––––绘制极坐标图hist––––绘制统计直方图stairs––––绘制阶梯图stem––––绘制火柴杆图rose––––绘制统计扇形图comet––––绘制彗星曲线area––––区域图pie––––饼图在MATLAB中，可以使用comet函数绘制二维彗星动画图。再来看看MATLAB中help里面自带的例子！是不是有点神奇？t=0:pi/100:2*pi;y=sin(t);x=cos(t);comet(x,y)创建y中数据对x中数据的彗星图。将y创建为正弦函数值向量，输入值介于0到2π之间。将x创建为余弦函数值向量，输入值介于0到2π之间。在各值之间使用递增量π/100。然后绘制数据图。动态图！有意思吧？新学一个magic函数！生成什么样的数据呢？这个函数也许会在考试中出现哦！数据有什么特征吗？图画好后，要修改？找谁？set函数来帮忙！预先了解一些扩展知识吧！也许大家以后会对人工智能感兴趣！在很多场合，二维图形远远满足不了用户的需求！为了可以呈现更加形象和逼真的三维图形，MATLAB提供了丰富的三维绘图函数。三维绘图可以看成二维绘图的拓展，一些绘制函数调用格式十分相似，很多图形绘制和设置函数也可以二、三维通用，但是三维图形仍有其特殊之处，例如需要进行二维图形所没有的视角、光照及透明度的设置。接下来以大家常用的三维曲线、网格和曲面图为重点进行讲解，介绍这些三维图形的绘制方法。平时大家只需要记住常用的几个函数即可，需要深入使用的时候，也需要通过互联网求助！这也是互联网时代的学习方式，这种方式非常便捷，减少了记忆的压力！首先是plot3函数！这个是常用函数！必须记住的那种！和二维绘制函数plot相对应，plot3函数是曲线绘制函数在三维空间中的扩展。两者的区别是：plot3增加了第三维的参数。plot3函数的常用调用格式为：plot(X,Y,Z,'s')X、Y和Z是同维向量时，分别以X、Y和Z的元素为x、y和z轴坐标。X、Y和Z是同维矩阵时，分别以X、Y和Z的对应列元素为x、y和z轴坐标，曲线条数等于矩阵的列数。plot3的用法与plot类似！！！一起来看看课堂上的教学视频吧，帮助大家巩固知识点！课后复习很重要！在教学过程中，发现同学们一直没有养成这个习惯！未完，待续！修订记录20220527完成初稿；20220603修订内容；20230725更新至v2；20231117更新至v3；来源：通信工程师专辑

在《ABAQUS分析之美》中有一个流固耦合的例子，抽空学习了一下，只讲了CEL方法。在Abaqus帮助文件中，有一个流固耦合分析的例子，用到了CEL、SPH两种方法，如下图所示，也学习了一下。如果是跟着帮助文件中的例子学习，SPH方法可能会遇到一些问题。下面主要讲一下SPH方法在CAE界面中的操作。1初识CEL和SPHCEL：CoupledEulerian-Lagrangian，欧拉-拉格朗日耦合；SPH：SmoothedParticleHydrodynamic，光滑粒子流体动力学。在上图2.3.2的例子中提供了3个算例，如下图所示。帮助文件中提供了3个例子的inp文件，对于Case1，可以将inp导入Abaqus/CAE，查看CEL方法的模型，提交分析。对于Case2、Case3，导入Abaqus/CAE可能会失败，但是可以从inp创建job并提交计算，进行SPH方法流固耦合分析。问题是，怎么在Abaqus/CAE中完成SPH方法的建模？2按图索骥上图中的Case3提供了从有限单元转换为SPH粒子的分析方法。那么，就从建立有限元模型，然后转换为SPH粒子这条线索尝试解决问题。还是在这个例子中，参考文献第3行给出了有限单元转换为SPH粒子的介绍（分析手册15.2.2节），见下图。在分析手册15.2.2节中，给出了在CAE界面中设置单元转换为SPH粒子的方法。即，建立表示流体的有限元模型，在选单元类型时勾选转换为粒子的选项，并进行相应设置即可，见下图。参数PPD见下图介绍。3浅尝辄止3.1CEL方法建模建立欧拉类型的部件，并划分网格。注意，划分网格前，先分割实体几何，将液体初始区域分割好，便于后面建立单元**。模型装配后，在Load模块预定义场中通过单元**定义液体初始区域，见下图。如此，流体区域便定义完成。至于固体（弹性体模型）、接触定义、分析步、边界条件等，见帮助文件或《ABAQUS分析之美》。3.2SPH方法建模流固耦合模型中的液体模型、固体模型均采用弹性体建模，杯中水实体模型和空杯子壳模型如下图所示。注意：给水体选单元类型时设置转换为粒子。分析步、接触、边界条件等略，计算时，水体单元转换为粒子，在后处理查看模型如下图所示。4小结CEL方法建立的欧拉网格要包络液体运动空间，单元较多。SPH方法只需建立液体初始包络的网格，单元较少。计算精度、计算速度等对比，可见帮助文件或自行尝试。来源：iCAETube