有限元的自动生成及其应用演示资料.ppt

有限元程序结构 有限元程序流程框图 有限元程序流程框图 有限元程序流程框图 有限元程序流程框图 用FEPG生成求解二维线弹性平面应变问题的有限元程序 三维正交各向异性材料计算 ——FEPG系统在复合材料中的应用 三维弹性本构关系 其中 其中 平面二维水沙运动问题 ——FEPG系统在水力学中的应用 问题描述 如右图所示，水在1/4环形水渠内流动，进口处x方向单宽流量为4, y方向单宽流量为0,水深为2，出口处水面高度为4。求水渠内的单宽流量和水边高度。 该问题的偏微分方程 该问题的偏微分方程 算子分裂法求解该问题 算子分裂法求解该问题 算子分裂法求解该问题 算子分裂法求解该问题 算子分裂法求解该问题 线性正弦稳态涡流问题（TEAM Workshop 问题7） —— FEPG系统在电磁学中的应用 问题描述 TEAM Workshop 问题7是用来检验三维线性正弦涡流场分析方法计算精度的实验模型，其结构如右图所示。该模型由带有方孔的铝板和载流线圈组成，线圈中通有随时间按正弦规律变化的电流，其幅值为2742安匝。铝板的电导率为 求涡流和磁通密度分布。 求解方程及边界条件 地下厂房洞室群的三维围岩稳定分析 —— FEPG系统在水电工程中的应用 工程背景 进行索风营电站的建设和进一步优化索风营电站地下工程的设计和施工。 通过模拟计算揭示出索风营围岩稳定状态，进行围岩稳定评价，为设计的优化和施工程序的安排提供指导。 地下厂房结构图 模拟计算的主要困难 　1.考虑横观各向同性的弹塑性； 　2.主变室及尾水隧洞所处岩层大部分为薄层； 　3.大量支护锚杆的存在如何处理； 　4.岩体渗流对围岩稳定的影响； 岩石的横观各向同性特征 软弱结构面三维造型 加锚杆岩体单元 渗流计算网格图 GIO文件 afi mage magb 方程的afi.vde文件 disp Ax,Ay,Az,u coor x,y,z func curlx,curly,curlz,div shap %1 %2 gaus %3 mass %1 1.0d0 mate sigma epsilon omega fmu fjx fjy fjz \0.0d0;0.0d0;50.d0;1.256d-6;0.0d0;0.0d0;0.0d0; vect A Ax Ay Az vect x x y z vect fj fjx fjy fjz vect curl curlx curly curlz 方程的afi.vde文件 disp Ax,Ay,Az,u coor x,y,z func curlx,curly,curlz,div shap %1 %2 gaus %3 mass %1 1.0d0 mate sigma epsilon omega fmu fjx fjy fjz \0.0d0;0.0d0;50.d0;1.256d-6;0.0d0;0.0d0;0.0d0; vect A Ax Ay Az vect x x y z vect fj fjx fjy fjz vect curl curlx curly curlz afi.vde文件 func curlx +[Az/y]-[Ay/z] curly +[Ax/z]-[Az/x] curlz +[Ay/x]-[Ax/y] div +[Ax/x]+[Ay/y]+[Az/z] afi.vde文件 stif $c6 eu1 omega*sigma $c6 eu2 1/fmu dist +[curl_i;curl_i]*|eu2;0.0| +[div;div]*|eu2;0.0| +[A_i;A_i]*|0.0;eu1| +[u/x_i;A_i]*|sigma;0.0| +[A_i;u/x_i]*|0.0;eu1| +[u/x_i;u/x_i]*|sigma;0.0| load +[A_i]*fj_i*|1.0;0.0| end 计算涡流的mage.pde文件 disp exr exi eyr eyi ezr ezi coef axr axi ayr ayi azr azi ur ui coor x y z shap %1 %2 gaus %3 mass %1 load fexr fexi feyr feyi fezr fezi mate sigma epsilon omega fmu fjx fjy fjz \0.0d0;0.0d0;50.d0;1.256d-6;0.0d0;0.0d0;0.0d0; mage.pde文件 stif $cv fexr +axi*omega*sigma- ur/x *sigma $cv fexi -axr*omega*sigma- ui/x *sigm