u010542847的博客

​本文主要讲解如何利用Matlab通过有限元编程实现框架结构在地震作用下弹性时程分析，重点介绍了Newmark-Beta求解方法（隐式方法）和瑞丽阻尼矩阵的基本原理及算法流程，并通过Matlab程序介绍了Newmark-Beta方法的和瑞丽阻尼矩阵的程序实现过程，其中框架结构采用铁木辛柯梁单元，其有限元方程的建立过程可参考往期博文《Matlab梁单元有限元编程：铁木辛柯梁VS欧拉梁讲解》​

本程序通过matlab实现了四连杆机构的运动学仿真编程，动态展现了四连杆机构的角位移、角速度和角加速度的时程曲线，除了程序本身，还提供了机构运动学详细的公式推导文档，从而帮助大家更好地去理解程序的物理背景和算法原理。与另一个四连杆案例（四连杆机构运动学仿真 | 【Matlab源码+理论公式文本】）的主要区别是该程序提供的算法说明文档不仅包含了公式的推导，并在文档中穿插讲解了每个公式对应的实现代码，理论文档截图和程序运行截图如下。你将获得四连杆机构Matlab仿真源码+程序理论文本doc文档获取链接。

案例介绍：SQOF系统参数：m=0.2533 kg,k=10 N/m,Tn=1 sec,（Wn=6.283rad/sec），阻尼比为0.2。用时间步长为0.05 sec的中心差法确定该系统对半周正弦脉冲p(t)=10 sin(Pi*t/0.6)[N]所定义的p(t)的响应u(t)（激励时程如下图所示）。并输出您的计算结果，包含ti,pi,Uit,Ui,pi和Ui+1信息的表格中显示结果。并绘制从0到1秒的SDOF的响应u(t)。注意：建议使用Matlab。理论文本详细介绍了有限差分法的求解动力学微分方程

https://mbd.pub/o/bread/mbd-ZpaamZxv

程序实现了壳单元的有限元线性屈曲分析，得到其屈曲模态和相应的屈曲因子。壳单元类型为MITC4混合插值壳单元。线性屈曲分析思路是这样，先对节点施加单位载荷，计算出对应单元内的内力，然后将求得的内力置入几何刚度阵中，再进行特征值分析分析，分析得到的特征值倒数叫做屈曲因子，反映的是临界承载力的大小。

​本程序通过对原始“拓扑优化99行MATLAB代码拓扑优化99行MATLAB代码解读 - 知乎”进行修改，实现了对一个悬臂梁结构进行位移荷载（位移边界条件）的施加。采用的单元为四边形平面单元。对比原始代码与修改后的代码，可以清晰地看出位移边界条件施加的方法。具体来讲，是在原FE函数上修改实现了位移载荷的施加，采用的方法为乘大数法，也叫惩罚因子法。修改后的函数命名为FE2，其中可以同时施加位移载荷和力载荷。通过程序中的“u_disp = -nelx/10;%设定位移量”代码可以自定义位移。代码截图和运行结果

采用剪切变形的十三自由度有限元模型，考虑了拉伸-扭转、弯曲-扭转和弯曲-拉伸耦合以及泊松效应，对含有开口横向裂纹的叠层组合梁进行了自由振动分析。拉格朗日方程用于推导运动方程。在单元质量和刚度矩阵的推导中，试函数被选择为挠度的三次拉格朗日多项式，而其他函数则被选择为二次拉格朗日多项式。通过降低损坏区域内元件的刚度，将损坏引入梁中。与现有文献进行了比较，以显示本要素的准确性。还对悬臂叠层复合材料梁在六种损伤情况下的环境振动试验进行了实验验证。根据结果，本研究提出的方法提供了一种简单快速、足够准确的方法来确定具有

代码是针对一个分段函数中的参数进行拟合，使用的拟合工具是matlab中的lsqcurvefit或nlinfit。具体函数形式和待拟合参数如下所示。该案例的特殊之处在于分段点也是待拟合参数，因此如何自定义拟合函数，实现分段点的拟合是本案例最大的难点。本代码提供了三种分段函数的拟合方法

本案例实现了下述二维三角形区域的泊松方程的Matlab有限元编程求解，具体方程形式如下图，介绍了泊松方程有限元求解基本原理，将Matlab求解的结果与Comsol求解的结果进行对比，证明了本Matlab求解程序的准确性。提供案例源码供大家练习。

本书是为编程而生的数学书，书中介绍了各种数值分析算法，原书附算例代码。有限元求解中涉及的一些数值算法可以在本书中找到详细的算法实现步骤，如牛顿迭代法、高斯积分方法、求解刚度矩阵逆矩阵的直接法、迭代法等等。作为一本数学工具书非常适合有限元编程人员使用。甚至当今火热的机器学习算法中似然函数的优化通常用牛顿法来求解，可以通过学习本书相关数值分析算法。详细给出了算法实现步骤。

Ted Belytschko 的《Nonlinear Finite Elements for Continua and Structures》，此书封面就是一个非线性有限元中的经典问题（壳体碰撞后的大变形）。此书有中文译本，是清华的庄茁老师翻译的。

​本案例介绍基于matlab实现分段函数拟合，并通过施加分段点连续约束条件。因为所拟合的分段函数如果按照传统的分段函数拟合方法（参考案例《三种分段函数拟合方法介绍（分段点未知）》）来自定义被拟合函数，会出现如下图所示的在分段点处不连续的情况，所以在实际编程过程中如何实现预期的分段点处连续的预取拟合效果，这就是这个案例所要介绍的核心要点。为了帮助大家更好的理解这个案例所传达的知识点，本案例有配套视频介绍，链接在文末。​

本案例重点介绍积累Matlab拟合常见错误及对应的解决办法，依托的案例是分段微分方程的拟合问题。通过展示正确的分段微分方程拟合源码和错误的分段微分方程拟合源码，对比分析来阐述拟合过程中常遇到的错误和问题，以及对应的解决办法，助力掌握matlab 通过编程进行拟合的方法。

本案例通过Matlab中的bp神经网络实现函数未知参数的拟合。网上常见的bp神经网络拟合的方式，是对于函数具体表达式位置，只知道自变量和因变量的序列，进行训练，形成一个黑盒的函数映射关系，最终达到已知输入能够预测输出的功能。而本案例针对已知函数形式，其中某些参数未知，需要一系列实验数据点，对函数进行拟合，求得未知参数，通常的做法是借助最小二乘拟合工具如lsqcurvefit、nlinfit等工具中的优化算法进行拟合，但是本案例特殊之处在于基于bp神经网络算法拟合已知函数形式中的未知参数，这种方法是一种非典型

本案例通过Matlab中的lsqcurvefit()实现方程组的拟合，案例的难点在于自定义方程组类型的拟合函数。具体方程形式如下图所示。

本案例通过Matlab中的lsqcurvefit()实现分段微分方程组的拟合，案例有三个难点：1.自定义拟合函数通过微分方程（求解）表示；2.微分方程分段；3.微分方程本身是方程组的形式。具体方程形式如下图所示。

本程序通过matlab实现隐式函数（隐式参数）的拟合。采用的工具是lsqcurvefit和nlinfit两个函数工具。因为无法求出函数的显示表达式，或者是参数的显示表达式，所以在lsqcurvefit和nlinfit自定义拟合函数的时候存在困难，无法显式去书写函数的表达式。但是间接通过fsolve命令实现隐式方程或参数的求解，可以实现隐式函数（隐式参数）的拟合，这也是我们这套代码的核心思想。具体的拟合案例如下，方程1是我们的待拟合函数，方程2-4都是函数中参数的定义，通过方程3可见参数w是个超越方程，无法得

提到了有限元，Zienkiewicz和Tayler教授的《The Finite Element Method 》是不得不说的，该书可以总共有三个分卷，本资源提供前两个分卷，主要围绕固体力学领域，第三分卷主要是流体。之所以被誉为是说是有限元领域中的圣经，原因有二：一是，因为其作者在该领域的鼻祖地位(Zienkiewicz是有限单元法的三位创始人之一);二是，因为其涵盖范围非常之广无所不包，可以当作有限元百科全书去学习。

提到了有限元，Zienkiewicz和Tayler教授的《The Finite Element Method 》是不得不说的，该书可以总共有三个分卷，本资源提供前两个分卷，主要围绕固体力学领域，第三分卷主要是流体。之所以被誉为是说是有限元领域中的圣经，原因有二：一是，因为其作者在该领域的鼻祖地位(Zienkiewicz是有限单元法的三位创始人之一);二是，因为其涵盖范围非常之广无所不包，可以当作有限元百科全书去学习。

你将获得PDF电子书+全书配套源码+Nam-Ho Kim教授PPT讲义这本书适合初学者，通过简明的语言和清晰的例子介绍了非线性有限元分析的基本概念。书中涵盖了材料非线性、几何非线性、接触问题、大变形理论等内容。特别适用于工程背景的学生，提供了易于理解的入门材料。非线性有限元编程的绝对入门书籍，涵盖几何非线性、材料非线性、接触非线性。均有代码➕公式，逐步推导。