有限元与comsol步骤

1、COMSOL Multiphsics 初级教程有限元方法简介与COMSOL操作入门FEM 介绍有限元方法是广泛用于解决应力分析，热传导，电磁学和流体力学等工程问题的数值方法。前处理阶段：1.建立求解域，并将之离散化成有限个单元，即将问题分解成节点和单元。2.假定描述单元物理属性的形函数，即利用一个近似的连续函数描述每个单元的解。3.建立单元的刚度方程。4.组装单元，构造总刚度矩阵。5.应用边界条件和初始条件，并施加载荷。求解阶段：6.求解线性或非线性微分方程组得到节点值，例如不同节点上的位移或热传导问题中的温度。后处理阶段：7.得到其他更重要的信息，如每个单元上的应力和每个单元中的热流密度。F

2、EM 介绍FEM 实例w1Lw2y求沿杆长度方向上不同点变形的大小，在以下的分析中，假设应用的载荷比杆的重量要大得多，因此忽略杆的重量。P每个单元的横截面面积，由定义单元节点处横截面的平均面积表示。41235PP1A2A3A4A1l2l3l4lFEM 实例为了研究典型的单元行为，考虑一个带有统一横截面A的单元，单元长度为l，承受的外力为F，产生的形变是l。FEM 实例llFllFFEM 实例近似单元的平均应力：=FA近似单元的平均应变：ll根据胡克定律：E()AEFll其中，A为受力单元的横截面积，E为材料的杨氏模量，l，l分布为受力单元的长度和形变。注意到合并方程与线性弹簧的方程 Fkx 相

3、似，因此一个中心点受力而且横截面相等的实体可以视作一个弹簧。每个单元的弹性行为可以由相应的线性弹簧模型描述，合并方程与线性弹簧方程 Fkx 相比，其等价刚度为xFeqkFEM 实例FEM 实例杆可以视为由四个弹簧串接起来的模型：1u2u3u4u5u12345P单元1单元2单元3单元41(AA )()()()2avgiieqi+1ii+1ii+1iAEEfku-uu-uu-ull其中等价单元的刚度为：1(AA )2iieqEkl节点的力平衡方程：1121121232232343343454454()0()()0()()0()()0()0Rk uuk uukuuk uukuuk uukuuk uu

4、PFEM 实例121()kuu232()kuu121()kuu343()kuu232()kuu454()kuu343()kuu454()kuuFEM 实例1 11 21 11 22223kukukukuk uk u1000RP 2 22 33 33 4k uk uk uk u3 33 44 44 5kukuk uk u4 44 5k uk uFEM 实例1111222233334444000000000000kkkkkkkkkkkkkkkk12345uuuuu 1000RpRKuF 反作用力矩阵反作用力矩阵=刚度矩阵刚度矩阵位移矩阵位移矩阵-负荷矩阵负荷矩阵FEM 实例 由于杆的上端是固定的，

5、节点1的位移量是零，因此方程的第一行应为u10。所以应用边界条件边界条件将得到如下的矩阵方程：1211223223343344544100000000000000000uukkkkukkkkukkkkukkp 刚度矩阵刚度矩阵位移矩阵位移矩阵=载荷矩阵载荷矩阵FEM 实例假定题目中E=10.4x106 lb/in2，w1=2 in, w2=1 in, t=0.125 in, L=10 in, P=1000 lb。由上述方法计算结果如下：00.0010260.0022100.0036080.005317uinCOMSOL操作步骤公式推导，数学抽象构建几何模型完成参数设定设置边界条件网格化计算后处理阅读文献，思考结果正确？得出正确结果NYCOMSOL操作步骤微分方程微分方程边界条件边界条件初始条件初始条件FEM自由变量的解自由变量的解后处理后处理COMSOL操作步骤构建集合模型完成参数设定设置边界条件网格化计算