有限元编程从入门到精通
文章平均质量分 85
suoge223
工学博士,高级工程师,多个平台专家博主,专注于计算机解决各类工程问题:有限元软件开发,振动信号处理,工程设备力学热学仿真分析;擅长编程语言:matlab,python,C/C++
展开
-
壳单元有限元屈曲分析Matlab程序【源码+理论文本】
程序实现了壳单元的有限元线性屈曲分析,得到其屈曲模态和相应的屈曲因子。壳单元类型为MITC4混合插值壳单元。线性屈曲分析思路是这样,先对节点施加单位载荷,计算出对应单元内的内力,然后将求得的内力置入几何刚度阵中,再进行特征值分析分析,分析得到的特征值倒数叫做屈曲因子,反映的是临界承载力的大小。原创 2024-05-15 22:44:43 · 846 阅读 · 0 评论 -
乘大数法施加位移载荷Matlab有限元编程 |源代码+理论文本| 位移边界条件有限元编程
本程序通过对原始“拓扑优化99行MATLAB代码拓扑优化99行MATLAB代码解读 - 知乎”进行修改,实现了对一个悬臂梁结构进行位移荷载(位移边界条件)的施加。采用的单元为四边形平面单元。对比原始代码与修改后的代码,可以清晰地看出位移边界条件施加的方法。具体来讲,是在原FE函数上修改实现了位移载荷的施加,采用的方法为乘大数法,也叫惩罚因子法。修改后的函数命名为FE2,其中可以同时施加位移载荷和力载荷。通过程序中的“u_disp = -nelx/10;%设定位移量”代码可以自定义位移。代码截图和运行结果原创 2024-05-13 21:17:16 · 251 阅读 · 0 评论 -
Matlab有限元求解层合梁的自由振动(多裂纹叠层复合材料梁)|论文复现
采用剪切变形的十三自由度有限元模型,考虑了拉伸-扭转、弯曲-扭转和弯曲-拉伸耦合以及泊松效应,对含有开口横向裂纹的叠层组合梁进行了自由振动分析。拉格朗日方程用于推导运动方程。在单元质量和刚度矩阵的推导中,试函数被选择为挠度的三次拉格朗日多项式,而其他函数则被选择为二次拉格朗日多项式。通过降低损坏区域内元件的刚度,将损坏引入梁中。与现有文献进行了比较,以显示本要素的准确性。还对悬臂叠层复合材料梁在六种损伤情况下的环境振动试验进行了实验验证。根据结果,本研究提出的方法提供了一种简单快速、足够准确的方法来确定具有原创 2024-05-12 22:54:45 · 615 阅读 · 0 评论 -
二维泊松方程(Neumann+Direchliet边界条件)有限元Matlab编程求解|程序源码+说明文本
本案例实现了下述二维三角形区域的泊松方程的Matlab有限元编程求解,具体方程形式如下图,介绍了泊松方程有限元求解基本原理,将Matlab求解的结果与Comsol求解的结果进行对比,证明了本Matlab求解程序的准确性。提供案例源码供大家练习。原创 2024-05-04 10:49:59 · 1118 阅读 · 0 评论 -
Numerical Analysis(byRichard.L..Burden)【pdf高清英文原版】
本书是为编程而生的数学书,书中介绍了各种数值分析算法,原书附算例代码。有限元求解中涉及的一些数值算法可以在本书中找到详细的算法实现步骤,如牛顿迭代法、高斯积分方法、求解刚度矩阵逆矩阵的直接法、迭代法等等。作为一本数学工具书非常适合有限元编程人员使用。甚至当今火热的机器学习算法中似然函数的优化通常用牛顿法来求解,可以通过学习本书相关数值分析算法。详细给出了算法实现步骤。原创 2024-05-02 17:41:18 · 529 阅读 · 0 评论 -
Nonlinear Finite Elements for Continua and Structures(by Ted Belytschko)| PDF电子书 | 连续体和结构的非线性有限元(庄茁)
Ted Belytschko 的《Nonlinear Finite Elements for Continua and Structures》,此书封面就是一个非线性有限元中的经典问题(壳体碰撞后的大变形)。此书有中文译本,是清华的庄茁老师翻译的。原创 2024-04-29 16:49:08 · 382 阅读 · 0 评论 -
【PDF】The Finite Element Method for Solid and Structural Mechanics |by Zienkiewicz & Taylor |有限元圣经三部曲
提到了有限元,Zienkiewicz和Tayler教授的《The Finite Element Method 》是不得不说的,该书可以总共有三个分卷,本资源提供前两个分卷,主要围绕固体力学领域,第三分卷主要是流体。之所以被誉为是说是有限元领域中的圣经,原因有二:一是,因为其作者在该领域的鼻祖地位(Zienkiewicz是有限单元法的三位创始人之一);二是,因为其涵盖范围非常之广无所不包,可以当作有限元百科全书去学习。原创 2024-04-11 21:19:21 · 477 阅读 · 0 评论 -
【PDF】The Finite Element Method:its Basis and Fundamentals|by Zienkiewicz & Taylor |有限元圣经三部曲
提到了有限元,Zienkiewicz和Tayler教授的《The Finite Element Method 》是不得不说的,该书可以总共有三个分卷,本资源提供前两个分卷,主要围绕固体力学领域,第三分卷主要是流体。之所以被誉为是说是有限元领域中的圣经,原因有二:一是,因为其作者在该领域的鼻祖地位(Zienkiewicz是有限单元法的三位创始人之一);二是,因为其涵盖范围非常之广无所不包,可以当作有限元百科全书去学习。原创 2024-04-11 21:16:43 · 833 阅读 · 0 评论 -
Introduction to nonlinear finite element analysis|Nam-Ho Kim | 【PDF电子书+配套Matlab源码+授课PPT讲义】
你将获得PDF电子书+全书配套源码+Nam-Ho Kim教授PPT讲义这本书适合初学者,通过简明的语言和清晰的例子介绍了非线性有限元分析的基本概念。书中涵盖了材料非线性、几何非线性、接触问题、大变形理论等内容。特别适用于工程背景的学生,提供了易于理解的入门材料。非线性有限元编程的绝对入门书籍,涵盖几何非线性、材料非线性、接触非线性。均有代码➕公式,逐步推导。原创 2024-04-08 23:18:28 · 373 阅读 · 0 评论 -
A First Course in the Finite Element Method【Daryl L.】|PDF电子书
这本书是有限元方法领域的经典之作,适用于初学者。作者详细介绍了有限元分析的基本原理,强调实际问题的建模和求解。书中包含大量的例子和练习,帮助读者建立对有限元方法的深刻理解。也有对应的中文译本:有限元方法基础教程(第五版)。推荐这本书,很大的一个原因是再一次觉得歪果仁写书真的叫用心和详细,几乎是一步一步的讲述有限元理论处理各类问题的过程,就差把这本书当演算稿纸全部收纳了,非常详尽。对于有限元理论涉及的每一类问题都有非常详细的基本公式推导说明以及实例讲解举例,充分看出作者的用心。原创 2024-04-04 13:29:20 · 733 阅读 · 0 评论 -
Finite Element Procedures K.J.Bathe 【教材pdf+部分源码】|有限元经典教材 | 有限元编程
K.J.Bathe作为有限单元法的创始人之一,作为Adina的开发者,了解程序开发的细节。你可以在这本书里找到非常到位的公式展开。比如说,书中的梁单元,壳单元,你基本上可以按其公式实装即可。本书条理清晰、内容相对全面、透彻,囊括了有限元很多基础理论,特别是书中结合了大量例题,使得这些理论不至于那么枯燥。在我看来,这本书或许二十年之后也不会过时。然而,对于初次接触有限元的人来讲,Bathe教授的书会稍显晦涩难懂,可能Logan教授的《A First Course in the Finite Eleme原创 2024-04-04 11:52:28 · 587 阅读 · 0 评论 -
Matlab梁单元有限元编程:铁木辛柯梁VS欧拉梁
梁在工程中应用广泛,是重要的结构构件。从几何上看,梁是任意截面形状的承受横向力的杆状结构,与杆的区别仅在于二者承受的载荷不同。在梁结构中,不同的梁固接在一起,既能传递力,又能传递力矩。本文针对二维和三维空间梁结构的matlab有限元编程进行讲解,涉及的梁单元类型有欧拉梁单元和铁木辛柯梁单元。重点讲解二者的基本力学假定、适用范围、对应的三大类方程的建立、有限元离散方程的建立(包括形函数、刚度矩阵推导等)以及通过Matlab编程的实现上述两类梁单元静力分析求解和模态分析求解,获得梁结构的位移、剪力、弯矩图,以及原创 2024-04-04 11:01:35 · 922 阅读 · 0 评论 -
Matlab有限元编程案例全家桶【源码+理论文本】
该全家桶资源为课程《Matlab有限元编程从入门到精通30讲》配套的全部源码和讲义PPT/理论文本,旨在以案例的形式讲解各类有限元问题程序实现及算法原理,并提供完整Matlab源码供大家练习,案例源码均包含前后处理模块和求解器模块;单元类型包含:杆单元,梁单元,平面三角形单元,板壳单元,四/六面体实体单元等;物理场问题涉及:静力学、动力学、传热学、材料非线性、几何非线性、接触非线性等求解。原创 2024-04-04 09:46:35 · 434 阅读 · 0 评论 -
框架结构模态分析/动力时程分析Matlab有限元编程 【Matlab源码+PPT讲义】|梁单元|地震时程动画|结果后处理|地震弹性时程分析| 隐式动力学
本程序利用Matlab实现了平面框架结构的动力学分析,包括地震作用下的时程分析和模态分析,采用的单元类型为欧拉梁单元,涉及的有限元知识点包含:阻尼矩阵、质量矩阵的建立,newmark-beta方法求解动力学方程以及地震作用下结构振动动画后处理。为了辅助理解代码背后的有限元理论,配套提供对应的框架结构动力学分析Matlab有限元编程视频教程和PPT讲义,详细讲解了上述有限元理论,辅助消化吸收有限元算法的程序实现过程。视频教程链接见文末。原创 2024-03-23 09:28:59 · 1311 阅读 · 0 评论 -
框架结构静力分析matlab有限元编程【Matlab源码+PPT讲义】| 梁单元 | 弯矩图 |坐标变换
本程序利用Matlab实现了平面框架结构的有限元编程求解,采用的单元类型为欧拉梁单元,涉及的有限元知识点包含:坐标变化矩阵(局部坐标与全局坐标)+梁单元刚度矩阵推导+乘大数法施加边界条件等有限元理论,并通过结果后处理,展示框架的位移图、弯矩图等。为了辅助理解代码背后的有限元理论,配套提供对应的框架结构Matlab有限元编程视频教程和PPT讲义,详细讲解了上述有限元理论,辅助消化吸收有限元算法的程序实现过程。视频教程链接见文末。原创 2024-03-22 22:16:33 · 987 阅读 · 0 评论 -
小孔平板应力集中问题matlab有限元编程【源码+PPT讲义】|三节点三角形单元|平面单元|稀疏矩阵 |Comsol网格
本程序利用Matlab实现带圆孔平板应力集中问题的的有限元编程求解,平板采用的单元为三角形单元,网格信息来自Comsol中分网后生成的.mph文件,因此代码涉及comsol.mph网格文件的读取接口的开发,另外由于单元节点数较大,为了提高计算和存储效率,采用了稀疏矩阵的方式定义单元刚度矩阵,涉及的有限元知识点主要是三角形单元刚度矩阵推导、乘大数法施加边界条件、稀疏矩阵的应用等有限元理论。为了辅助理解代码背后的有限元理论,配套提供带圆孔平板应力集中问题Matlab有限元编程视频教程和PPT讲义,详细讲解了上述原创 2024-03-22 17:38:00 · 975 阅读 · 0 评论 -
三角形单元悬臂梁Matlab有限元编程 【程序源码+PPT讲义】|平面单元 | 三节点
本代码利用Matlab实现三节点三角形单元悬臂梁结构的有限元编程求解,涉及的有限元知识点主要是三角形单元刚度矩阵推导+乘大数法施加边界条件等有限元理论。为了辅助理解代码背后的有限元理论,配套提供对应的桁架结构Matlab有限元编程视频教程和讲义PPT,详细讲解了上述有限元理论。视频教程链接见文末。原创 2024-03-22 16:30:59 · 531 阅读 · 0 评论 -
2D/3D桁架结构Matlab有限元编程 |【Matlab源码+视频教程+讲义PPT】| 平面桁架 | 空间三维桁架 | 网架结构
本代码利用Matlab实现平面桁架结构、三维空间桁架结构(2D/3D)的有限元编程求解,涉及的有限元知识点包含:坐标变化矩阵(局部坐标与全局坐标)+杆单元刚度矩阵推导+乘大数法施加边界条件等有限元理论。为了辅助理解代码背后的有限元理论,配套提供对应的桁架结构Matlab有限元编程视频教程,详细讲解了上述有限元理论,并对Matlab代码逐行讲解,真正消化吸收有限元算法的程序实现过程。视频教程链接见文末。原创 2024-03-22 15:57:34 · 1024 阅读 · 0 评论 -
8节点空间壳单元Matlab有限元编程 | 曲壳单元 | 模态分析 | 3D壳单元 | 板壳理论| 【源代码+理论文本】
本程序采用的8节点厚壳单元同样是基于Mindlin/Reissner 理论的曲面壳单元,这个理论的几个基本假定是:(1)与壳的厚度相比,位移是微小的,即小位移假定;(2)垂直于中面的应力应变忽略不计,即沿厚度方向上的应力应变忽略;(3)变形前垂直于中面的直线,变形后仍保持为直线,但不一定再垂直于中面;(4)挠度和法线转角为各自独立的场函数。但是由于需要模拟面内变形,其自由度由板单元的三个自由度[u thetax thetay] 加为5个自由度[u v w phi psi ] ,前三个为平动自由度,后两个为沿原创 2024-03-21 22:40:34 · 1713 阅读 · 0 评论 -
薄板/厚板模态分析Matlab有限元编程 | 【源码+理论文本】|板单元|板壳单元|Mindlin Reissner
本程序主要通过Matlab实现mindlin板的有限元模态分析,单元类型为Mindlin中厚板单元。程序的理论背景除了基本mindlin板理论,还要涉及等参单元和高斯积分的有限元理论,以及单元刚度矩阵和质量矩阵的推导和建立。mindlin板理论适用于厚板,因为厚板的变形包括弯曲和剪切变形,而且剪切变形不可忽略,认为垂直于板中面的平面在板变形之后将不再垂直于中面,所以需要同时考虑板的剪切能量和弯曲能量,为了方便大家理解,随同程序还附赠了Mindlin板理论文本。当然所谓Mindlin厚板理论并非只能用于厚板,原创 2024-03-18 14:10:55 · 1174 阅读 · 0 评论 -
matlab矩形薄板小挠度弯曲有限元编程 |【Matlab源码+理论文本】|板单元| Kirchoff薄板 | 板壳单元
本程序通过Matlab实现了矩形薄板小挠度弯曲有限元编程,基于Kirchoff-love板壳理论,进行有限元离散,得到其该板单元的形函数、应变矩阵、刚度矩阵等有限元方程,继而通过Matlab程序实现该矩形薄板在均布荷载/集中荷载作用下以及不同边界条件(四边简支/四边固定/单边固定三边简支/两边简支两边固定/三边自由等)的有限元求解。你将获得矩形薄板小挠度弯曲有限元求解程序源码(含不同边界条件)+程序理论文本。原创 2024-03-18 11:57:03 · 911 阅读 · 0 评论 -
六面体单元热应力问题有限元求解 | Matlab源码 | 理论文本 | 温度应力 | 热应变 |温度应变 |热力耦合
本代码主要实现热应力问题的matlab有限元编程。热应力也叫温度应力,后文提到的热应变也叫温度应变。这里需要与传热问题的有限元分析进行区分:可以认为传热分析是进行热应力分析的前提条件,通过传热分析来确定温度场,在获得温度场的基础上,计算所产生的热应力。所以代码计算热应力问题前提是假定温度场是已知的,我们并不关心温度场是如何得到的,那是热传导要解决的问题,在这个已知的温度场作用下,由于热胀冷缩会产生响应的热应变进而产生热应力,这是本代码要解决的问题。原创 2024-03-18 10:44:11 · 870 阅读 · 0 评论 -
接触问题有限元Matlab编程求解(上)|接触问题(点-线接触) | Matlab源码 | 理论文本 |接触非线性
本次课程主要围绕接触问题的有限元Matlab编程求解进行介绍,接触问题与我们之前讲过的几何非线性问题同属非线性有限元的范畴,但相比材料非线性和几何非线性来说他更难求解,其非线性更强,相信好多用商业有限元软件做接触分析的小伙伴最常遇到的问题就是不收敛了,为啥不收敛就是因为接触分析有很强的非线性。首先接触面的几何特征是未知的,两个物体的哪些节点会相互接触,形成的接触面什么样子,这些都是未知的,都需要反复的试错迭代才能找到;原创 2024-02-22 10:30:13 · 991 阅读 · 4 评论 -
Mindlin厚板单元Matlab有限元编程 | 板单元 | 【Matlab源码 + 理论文本】| 板壳单元 | 平板单元| 剪切变形 |Mindlin-Reissner
本程序主要通过Matlab实现mindlin板的有限元静力求解,单元类型为Mindlin中厚板单元,具体分析对象为矩形板受均布力作用。程序的理论背景除了基本mindlin板理论,还要涉及等参单元和高斯积分的有限元理论。mindlin板理论适用于厚板,因为厚板的变形包括弯曲和剪切变形,而且剪切变形不可忽略,认为垂直于板中面的平面在板变形之后将不再垂直于中面,所以需要同时考虑板的剪切能量和弯曲能量,为了方便大家理解,随同程序还附赠了Mindlin板理论文本。当然所谓Mindlin厚板理论并非只能用于厚板,而是因原创 2024-03-15 13:29:20 · 623 阅读 · 0 评论 -
三维铁木辛柯梁Matlab有限元编程 | 弹簧支座 | 弹性支撑单元| Matlab源码 | 理论文本 | 三维梁 | 3D梁 | 空间梁
本代码通过matlab实现了三维铁木辛柯梁单元和弹性支撑单元的有限元编程,采用 Timoshenko 梁理论,通过截面剪切系数考虑梁截面的剪切变形影响。三维铁木辛柯梁单元是由3 个节点构成的直梁单元,其示意图见图 1-1。其中,o-xyz 为梁单元局部坐标系;节点i、j为单元物理节点,用于确定单元边界,为必选节点:节点k为单元梁截面主轴方向节点,用于确定梁截面主轴Z 的指向,为可选节点,当不设置 3 号节点时,梁的截面主轴将按照梁截面绕梁轴线的转角进行确认。三维铁木辛柯梁单元每个节点包含 6 个位移自由度,原创 2024-03-15 11:16:29 · 1077 阅读 · 0 评论 -
Matlab有限元结果后处理 | 不规则云图绘制(二维/三维)|【Matlab源码+视频教程】| 三维云图 | 二维云图
本程序通过matlab实现有限元计算结果的后处理显示,具体通过两个案例来介绍:案例一、实现一个不规则形状的二维云图的绘制,其中还包含三角形网格剖分方法,然后通过节点解来绘制云图,具体后处理结果如图1所示;案例二、实现三维云图(四面体单元)的绘制,由于这个案例给出单元节点组成,所以无需划分网格,直接通过节点解来绘制云图,如图2所示。原创 2024-03-14 09:25:54 · 1475 阅读 · 0 评论 -
平面纯弯梁单元Matlab有限元编程 |欧拉梁单元| 简支梁|悬臂梁|弯矩图 |变形图| Matlab源码 | 视频教程 |二维梁
代码实现了平面纯弯梁(欧拉梁单元)单元Matlab有限元编程求解,针对的结构包含悬臂梁、简支梁、固支梁三种结构在均布荷载作用下的求解,并输出了梁的弯矩图、挠度、转角变形图原创 2024-03-07 18:02:38 · 695 阅读 · 2 评论 -
四节点/八节点四边形单元悬臂梁Matlab有限元编程 | 平面单元 | Matlab源码 | 理论文本 | 高阶单元 | QUAD4 | QUAD8
代码以受均布荷载和集中荷载的变截面悬臂梁为研究对象,通过matlab编制四节点和八节点四边形单元有限元程序来对悬臂梁进行受力分析。原创 2024-03-07 15:51:18 · 1025 阅读 · 0 评论 -
四面体单元悬臂梁的Matlab有限元编程 | 实体单元 | Matlab源码 | 理论文本| C3D4
本程序主要通过matlab实现针对四面体单元划分的三维结构进行有限元编程,具体案例是一个悬臂梁受集中荷载的问题。图1为本案例Matlab编程计算得到的结果。代码涉及的知识点主要包括四面体单元的有限元基本理论的推导,主要是单元刚度矩阵的推导,此外还包括等参单元和Hammer数值积分以及三维问题的后处理计算。原创 2024-03-07 15:08:34 · 835 阅读 · 0 评论 -
如何用python架构一个简单的有限元求解器
设计一款高效、灵活、易用的有限元求解器是一项复杂而关键的任务,需要考虑多个方面,包括架构设计、数据结构选择、算法优化等。- **单元(Element):** 表示有限元网格中的单元,包括单元编号、节点编号、单元类型(如三角形、四边形等)等信息。- 单元(Element): 表示有限元网格中的单元,包括单元编号、节点编号、单元类型(如三角形、四边形等)等信息。- **节点(Node):** 表示有限元网格中的节点,包括节点编号、坐标等信息。原创 2024-02-07 22:29:24 · 816 阅读 · 0 评论 -
材料非线性Matlab有限元编程:切线刚度法(牛顿-拉普森法)
matlab有限元编程从入门到精通》今天我主要围绕材料非线性问题的有限元Matlab编程求解进行介绍,重点围绕牛顿-拉普森法(切线刚度法)、初应力法、初应变法等三种非线性迭代方法的算法原理展开讲解,最后利用Matlab对材料非线性问题有限元迭代求解算法进行实现,展示了实现求解的核心代码。由于篇幅原因,本文分两期分享,今天重点介绍切线刚度法(牛顿-拉普森法)案例应用。一、切线刚度法(牛顿-拉普森法)案例背景。原创 2024-02-04 20:31:12 · 939 阅读 · 0 评论 -
材料非线性Matlab有限元编程:初应力法与初应变法
导读:本文主要围绕材料非线性问题的有限元Matlab编程求解进行介绍,重点围绕牛顿-拉普森法(切线刚度法)、初应力法、初应变法等三种非线性迭代方法的算法原理展开讲解,最后利用Matlab对材料非线性问题有限元迭代求解算法进行实现,展示了实现求解的核心代码。原创 2024-02-04 20:47:44 · 877 阅读 · 0 评论 -
有限元编程经典教材推荐
有限元编程经典教材推荐原创 2024-02-03 15:06:47 · 952 阅读 · 0 评论 -
非线性有限元经典教材推荐
非线性有限元分析是有限元方法的一个重要分支,广泛应用于工程、材料科学、地球科学等领域。以下是几本非线性有限元教材的推荐,每本书都有其独特的特点,涵盖了理论、应用和实践方面的内容。原创 2024-02-04 09:45:04 · 341 阅读 · 0 评论 -
有限元编程学习方法及相关资料链接
学习有限元编程是一个深入理解有限元方法、数值分析和结构力学的过程。以下是一个分阶段的学习方法,每个阶段提供相关资源链接和建议,以帮助你逐步掌握有限元编程的基础和高级知识。原创 2024-02-04 09:53:00 · 345 阅读 · 0 评论