非线性有限元：基本理论与算法及基于Python、Fortran程序实现
**前言：**有限单元法在岩土工程问题中应用非常广泛，很多商业软件如Plaxis/Abaqus/Comsol等都采用有限单元解法。尽管各类商业软件使用方便，但其使用对用户来说往往是一个“黑箱子”。相比而言，开源的有限元程序计算方法透明、计算过程可控，用户可根据自己的需求进行必要的修改，这一点对于科研人员特别重要。然而，长期以来由于有限元开源程序的欠缺、教学方法以介绍理论为主、手把手基于Tutorial的实操培训欠缺、学员编程基础相对薄弱等原因，开源有限元的使用在国内相对较少。然而，科研领域的很多新的东西很难实现在商业软件中（除非通过复杂的二次开发，而二次开发往往需要具备一些子程序编写的能力，如Abaqus的Umat子程序、Plaxis的UDSM等），如果掌握一些开源有限元的基本流程及相应子程序的编写方法，对基于自己开发程序或者商业软件开发接口的科研等工作将大有益处。本次课程主要针对岩土工程中的材料非线性问题（应力变形分析中的各种弹塑性本构模型、渗流分析中的非线性渗透系数等），用两个周末的时间，手把手教大家如何一步一步地搭建自己的模型，包括前处理输入文件的准备、计算可执行文件的编译生成、后处理输出文件的可视化等，其中的编程模块尽可能做到让学员零基础开始，通过提供标准程序库使用搭积木、摆乐高的形式模块化进行，加强学员的易上手、可操作性。学员只需了解简单的编程语句结构即可，甚至都不需要懂某种编程语言的具体语法。通过本课程的学习达到以下目标：学完以后发现不经意间已经学会了某种编程语言（尽管没有进行单独的语言学习），发挥通过实例学习的益处。同时本课程将介绍Paraview等后处理软件的使用及其Fortran/Python等接口程序。实例包括非线性有限元方法在挡土墙、地基稳定、边坡稳定等三大稳定性问题中的应用等。
非线性有限元：基本理论与算法及基于Python、Fortran程序实现与案例分析
**导师：李老师【副教授】：**985高校本科，海外知名高校毕业博士，北京双一流高校副教授，博导，长期从事岩土工程问题的有限单元法编程与程序设计、相关商业软件（Plaxis/Abaqus/Comsol/Flac3D）等的二次开发工作，熟悉Fortran/Matlab/Python等多种编程语言，发表了多篇SCI高水平论文，参与和承担欧盟、国家自然科学基金等多项项目，具有丰富的科研和教学经验。
有限元编程基础知识
1.有限单元法基础简介（离散化、存储策略及方程解法、边界条件的处理）
2.编程语言Fortran及编译工具Intel Visual Fortran(IVF)简介、Python开发环境Spyder简介
3.Fortran/Matlab/Julia/Python等开源代码及程序库（geomlib/femlib）简介
4.非线性有限元模型相关模块及计算流程
岩土工程材料非线性问题介绍
非线性有限元常刚度法实现模块
有限元分析数值基础与常刚度法
1)天际线法等整体系数矩阵存储与方程直接求解方法
2)整体系数矩阵免组装（预条件）共轭梯度法方程求解方法
3)牛顿（-拉夫森）迭代求解基本原理与数值实现
4)应力不变量、屈服函数、塑性势函数、几种典型的破坏准则及未定义的奇异点处理
初应变常刚度法的实现方法（虚拟粘塑性算法）、初应力常刚度法的实现方法
图2 几种破坏准则（米塞斯、屈雷斯卡、摩尔库伦）
非线性有限元切线刚度法实现模块
1.弹塑性应力应变关系的本构积分（前向欧拉法、后向欧拉法）
2.塑性流动准则、切线刚度法中的非一致和一致切线模量
3.应力积分算法之最近点投影法CPPM (Closest Point Projection Method) 实现
4.应力积分算法之切平面法CPA (Cutting Plane Algorithm) 实现
5.应力积分算法之径向返回法RR (Radial Return)与各算法的收敛准则
图4 应力应变本构积分
地基稳定性、挡土墙稳定性、边坡稳定性实例应用模块
应用实例
1）地基承载力问题（位移控制、力控制）
2）挡土墙稳定性问题（初应力法改进）
3）边坡稳定性问题（强度折减法及安全系数自动搜索）
4）输入文件准备、输出文件IrfanView、Paraview可视化（手把手操作与讲解）
5）Pyhton/Fortran90源代码按需修改与重新编译(应用于随空间变异的岩土体参数)
图6 边坡稳定性分析的强度折减法（直接解法与迭代解法）流程图
图8地基土内应变不变量（破坏面）
点击查看原文
非线性有限元：基本理论与算法及基于Python、Fortran程序实现相关推荐
- 非线性有限元：基本理论与算法及基于Python、Fortran程序实现与案例分析实践技术
有限单元法在岩土工程问题中应用非常广泛,很多商业软件如Plaxis/Abaqus/Comsol等都采用有限单元解法.尽管各类商业软件使用方便,但其使用对用户来说往往是一个"黑箱子" ...
- 非线性有限元：基本理论与算法及基于Python、Fortran程序实现与案例分析
非线性有限元:基本理论与算法及基于Python.Fortran程序实现与案例分析 (qq.com) 有限单元法在岩土工程问题中应用非常广泛,很多商业软件如Plaxis/Abaqus/Comsol等都采 ...
- 岩土工程--非线性有限元：基本理论与算法及基于Python、Fortran程序实现与案例分析
非线性有限元:基本理论与算法及基于Python.Fortran程序实现与案例分析实践技术 有限单元法在岩土工程问题中应用非常广泛,很多商业软件如Plaxis/Abaqus/Comsol等都采用有限单元 ...
- CSDN周赛第39期：赢签名实体书和图书《从零开始学算法（基于Python）》
本场竞赛由「电子工业出版社 & CSDN」联合主办. 一.报名方式 第39期周赛报名地址,3月22日 19点开考(周三):https://edu.csdn.net/contest/detail ...
- EMD算法之Hilbert-Huang Transform原理详解和案例分析
目录 Hilbert-Huang Transform 希尔伯特-黄变换 Section I 人物简介 Section II Hilbert-Huang的应用领域 Section III Hilbert ...
- python信用评分卡_基于Python的信用评分卡模型分析（二）
上一篇文章基于Python的信用评分卡模型分析(一)已经介绍了信用评分卡模型的数据预处理.探索性数据分析.变量分箱和变量选择等.接下来我们将继续讨论信用评分卡的模型实现和分析,信用评分的方法和自动评分 ...
- 【人工智能毕设之基于Python+flask+bilstm的评论情感分析系统-哔哩哔哩】 https://b23.tv/QU56eTl
[人工智能毕设之基于Python+flask+bilstm的评论情感分析系统-哔哩哔哩] https://b23.tv/QU56eTl https://b23.tv/QU56eTl
- 基于Python的微博大数据舆情分析，舆论情感分析可视化系统
运行效果图 基于Python的微博大数据舆情分析,舆论情感分析可视化系统 系统介绍 微博舆情分析系统,项目后端分爬虫模块.数据分析模块.数据存储模块.业务逻辑模块组成. 先后进行了数据获取和筛选存储, ...
- 基于python的公众号课堂教学_基于Python的微信公众号数据挖掘分析
基于Python的微信公众号数据挖掘分析 华南农业大学电子工程学院 王 建 黄宁香 [期刊名称]电子世界 [年(卷),期]2019(000)011 [总页数]3 运用Python网络爬虫技术对某时事类 ...
最新文章
- Android Cordova 对于软键盘弹出后覆盖输入域的解决
- 51单片机多任务操作系统的原理与实现
- php jsonp实例 mip无限滚动组件接口注意事项
- Linux文件夹的打开方式
- c++创建文件_使用Python实现文件压缩和解压
- SAP License：请教一个SD和FM基金预算集成的问题
- nginx限速方法以及报错处理
- Python的psutil模块监控内存代码
- Android 自定义标尺滑动选择值
- IIR滤波器和FIR滤波器
- javascript封装滑块
- T-SQL简介及基本语法
- JS中的一些条件判断语句
- 惠普笔记本重装系统按哪个键
- PNG图片压缩原理--屌丝的眼泪 #1
- ORA-01400: cannot insert NULL into | 通过SQL链接服务器 往Oracle库的表中 插入默认值问题
- Tracup Talk：如何制作项目管理的甘特图？
- 统计学中的白噪声理解
- 在Qt中如何用纯代码写一个本地视频播放器
- 利用Python实现ARCGIS中栅格计算器的简单功能