Matlab有限元网格化源程序-distmesh2d.m MIT的人写的论文，一个简单的有限元网格化方法 确实可以用，但是我理解的不是很好。 不过，希望对需要的人有所帮助。 简单介绍： function [p,t]=distmesh2d This meshing function produces the following outputs: • The node positions p. This N-by-2 array contains the x, y coordinates for each of the N nodes. • The triangle indices t. The row associated with each triangle has 3 integer entries to specify node numbers in that triangle. The input arguments are as follows: • The geometry is given as a distance function fd. This function returns the signed distance from each node location p to the closest boundary. • The desired edge length function h is given as a function fh, which returns h for all input points. • The parameter h0 is the distance between points in the initial distribution p0. For uniform meshes = constant), the element size in the final mesh will usually be a little larger than this input. • The bounding box for the region is an array bbox=[xmin, ymin; xmax, ymax]. • The fixed node positions are given as an array pfix with two columns. • Additional parameters to the functions fd and fh can be given in the last arguments varargin .

基于有限元和谱元法的很多matlab的代码，没有课程的电子版，但有很多代码源

matlab有限元法绘制静电透镜电势线电场线.zip

matlab有限元法绘制球域电势线电场线.zip

《Matlab有限元分析与应用》源代码

MATLAB有限元与谱元法导论 英文原版第2版 Introduction to Finite and Spectral Element Methods using Matlab C Pozrikidis Second Edition, 2014, Chapman & Hall/CRC 《MATLAB有限元与谱元法导论》对有限元法和谱元法进行了全面系统、深入浅出的阐述，并对在一般课题的对流一扩散和力学中的应用作了详细介绍，全书分8章和附录，第1章阐述有限元法在一维问题应用中涉及的计算模型和算法，建立了有限元法计算的基本框架；第2章是对第1章内容的深化和扩展，介绍非稳态问题有限元方程的时程积分法和有限元法在梁弯曲、屈曲中的应用；第3章叙述一维问题中谱元法的基本理论和方法，引入正交多项式和谱插值概念，介绍Lagrange，Chebyshev和Legendre等几种常用正交多项式的插值结点配置方法和相应的数值运算；第4章和第5章将前面几章介绍的有限元法和谱元法向二维问题扩展；第6章讨论有限元法和谱元法在固体力学中的应用；第7章介绍黏性流体流动问题的有限元法处理过程；第8章讨论了三维问题中谱元法的应用，这章是对第3，5章中谱元法的拓展和推广，全面阐述了谱元法的一般原理和方法，尤其值得一提的是，附录将书中使用到的基础数学知识进行了汇总和概括，方便读者检索查阅。

基于Matlab，开发了求解梁问题的有限元程序

有限元结构动力学分析与工程应用这本书附带的（Matlab源程序）

有限元建立梁单元，刚度矩阵质量矩阵计算，动力学方程求解，振动特性。

Matlab有限元网格化源程序-ddd.m MIT的人写的论文，一个简单的有限元网格化方法 确实可以用，但是我理解的不是很好。 不过，希望对需要的人有所帮助。 简单介绍： function [p,t]=distmesh2d This meshing function produces the following outputs: • The node positions p. This N-by-2 array contains the x, y coordinates for each of the N nodes. • The triangle indices t. The row associated with each triangle has 3 integer entries to specify node numbers in that triangle. The input arguments are as follows: • The geometry is given as a distance function fd. This function returns the signed distance from each node location p to the closest boundary. • The desired edge length function h is given as a function fh, which returns h for all input points. • The parameter h0 is the distance between points in the initial distribution p0. For uniform meshes = constant), the element size in the final mesh will usually be a little larger than this input. • The bounding box for the region is an array bbox=[xmin, ymin; xmax, ymax]. • The fixed node positions are given as an array pfix with two columns. • Additional parameters to the functions fd and fh can be given in the last arguments varargin .