3.5　EFGM的MATLAB程序

本节给出EFGM 求解悬臂梁问题（算例3.4.5.1）的MATLAB完整程序，其中权函数可选择3 次或4 次样条函数，基函数可选1—3 次多项式，边界条件的施加可选择罚函数法或拉格朗日乘子法。以下给出主程序及其子程序。

MATLAB 程序六 exam1_timoshenko_beam_EFGM

功能：EFGM 求解悬臂梁问题，在圆形支持域上采用MLS法建立形函数，由拉格朗日乘子法或罚函数法施加本质边界条件。共分为7 个步骤：输入数据、生成节点和四边形背景积分网格、确定圆形影响域半径、定义弹性矩阵、形成刚度矩阵、形成力向量和求解线性方程。

MATLAB 程序十一 Lagrange_multiplier_method.m

%功能：拉格朗日乘子法施加本质边界条件

MATLAB 程序十二 Penalty_method.m

% 功能：罚函数法施加本质边界条件