MATLAB有限元分析工具箱，可用于多物理场过程分析，具体使用教程可参考：https://www.featool.com/doc/quickstart.html#install。

比例边界有限元Matlab源码（毕业设计）.zip

SUBROUTINE UMAT(stress,statev,ddsdde,sse,spd,scd, 1 rpl, ddsddt, drplde, drpldt, 2 stran,dstran,time,dtime,temp,dtemp,predef,dpred,cmname, 3 ndi,nshr,ntens,nstatv,props,nprops,coords,drot,pnewdt, 4 celent,dfgrd0,dfgrd1,noel,npt,layer,kspt,kstep,kinc) c WRITE (6,*) ' c NOTE: MODIFICATIONS TO *UMAT FOR ABAQUS VERSION 5.3 (14 APR '94) c c (1) The list of variables above defining the *UMAT subroutine, c and the first (standard) block of variables dimensioned below, c have variable names added compared to earlier ABAQUS versions. c c (2) The statement: include 'aba_param.inc' must be added as below. c c (3) As of version 5.3, ABAQUS files use double precision only. c The file aba_param.inc has a line "implicit real*8" and, since c it is included in the main subroutine, it will define the variables c there as double precision. But other subroutines still need the c definition "implicit real*8" since there may be variables that are c not passed to them through the list or common block. c c (4) This is current as of version 5.6 of ABAQUS. c ... ...

建立了基于有限元方法结合增量迭代混合法研究适用于预应力混凝土梁非线性全过程分析的数值模型，提出的算法能较好地处理梁加载过程的响应变化，以及可穿越梁破坏前可能出现的极值点问题。并利用该数值模型分析评估了预应力筋配筋率、有效预应力对预应力混凝土梁性能的影响，该模型所得结果可为预应力混凝土梁的优化设计提供相应的参考。

梅蒂斯4 METIS 是一组用于划分图、划分有限元网格和为稀疏矩阵生成填充减少排序的串行程序。 更多信息： : 这个 fork 只是增加了 CMake 支持。

立方晶体，CPFE-UMAT子程序，需要MS+IVF以及Abaqus

workbench 学习资料 丰富全面覆盖各个基础知识，讲解清晰。例子生动。是很好的学习助手。相当棒的资料！

3.5 本章要点 ·有限元分析的基本方程(刚度矩阵、载荷列阵、单元的刚度方程、组装后的刚度方程) 换 3.6 习题 3-1 试用虚功原理推导杆 1D 单元的单元方程。 3-2 推导二次函数变截面的杆单元的刚度矩阵 finish !%%%%%%% [ANSYS 算例]3_4_2(3) %% end ·1D 问题(杆及梁)的虚功原理求解 ·1D 问题(杆及梁)的 小势能原理求解 ·杆单元的位移模式、刚度矩阵及坐标变 ·梁单元的位移模式、刚度矩阵及坐标变换 ·平面梁单元的等效节点载荷的计算方法 2 1 2( ) 1 x x A x A A l l ⎛ ⎞ ⎛ ⎞= − +⎜ ⎟ ⎜ ⎟ ⎝ ⎠ ⎝ ⎠ 。 3-3 对于沿轴向承受均布载荷的杆单元，推导相应的节点等效载荷。 荷。 有一个变截面的杆结构如图所示，材料的弹性模量为 ， (其 3-4 对于沿轴向承受三角形载荷分布的杆单元，推导相应的节点等效载 3-5 5 试采用两个单元2 10 MPaE = × 中一个为变截面单元)来分析该问题，并得到整个结构所有的力学信息。 (a) 一个变截面的杆结构 (b)两个单元的计算模型 3-6 有一个杆结构，受有沿轴向的分布载荷，杆的横截面积为 ，材料的弹性模 量为 ，分析该问题的建模方案，并得到整个结构所有的力学信息。 题图 3-5 2 22 10 cmA −= × 52 10 MPaE = × 84

该软件以文件MAIN运行。 下载项目时，会有一个excel文件，其中包含有关应如何进行输入的示例（包括一些数字）。

浙大徐荣辉的教材，用MATLAB编程来对结构进行有限元分析。