腾讯云开发者社区-腾讯云

开发者社区

文档建议反馈控制台

最新优惠活动

文章/答案/技术大牛

发布

数值分析与有限元编程

专栏作者

278

文章

682990

阅读量

99

订阅数

Rust语言尝鲜

首先，还是要吹捧一下这个编程语言。语法很严谨，Rust语言号称只要编译通过就不会崩溃(内存安全)。不像C++那种，概念混乱，连Bjarne Stroustrup都曾开玩笑说自己已经搞不懂C++了。而且一个报错信息就一万多行，然后必须是大师级的程序员，精通 Intel规范的汇编语言的那种C++程序员，才能搞清楚什么地方出了问题。个人认为Rust对于程序员的综合素质要求，是低于C++的。所以笔者决定将Rust作为首选。下面牛刀小试一把。以函数式编程方式，计算数值积分。以定积分为例。该积分精确值为-6.75.

2022-08-29

3790

面积坐标推导三角形常应变单元(CST)

本文内容：面积坐标推导三角形常应变单元(CST) 三角形面积坐标理论点这里：三角形面积坐标单元刚度矩阵如图所示，CST单元的位移场其中写成矩阵形式或者单元应变场其中用微分公式得到即可得到单元刚度矩阵单元刚度矩阵具有显式表达式。利用python的符号计算库sympy推导单元刚度矩阵表达式 import sympy as sy b1, b2, b3, c1, c2, c3 = sy.symbols('b1 b2 b3 c1 c2 c3') n, k = sy.symbols('n k

2022-06-09

9810

非线性| 弧长法实例

有关弧长法的内容已经有很多了，程序也有。即便这样，离实用还有很远的路要走。这里再发一个例子。《混凝土结构有限元分析第二版》第245页的例题，书中是手算演示，我这里用python实现。 python代码: import math def KT(): return 1 def Fint(x): if x > 1 : return 1 - 0.15*(x-1) else: return 2*x - x**2 f_ext = 1 # 适应书中

2022-05-27

5390

利用Python提升终端输入命令的准确性

unicode ascii 编程算法 python

终端输入一个命令，相当于一个字符串。如何检查输入是否正确，需要平时多练习，这里提供一个思路。利用python写一个函数来判断。利用python的切片功能，很方便的从字符串中遍历单个字符，并利用ord()函数获取其对应的ASCII 数值，或者 Unicode 数值。具体代码如下:

2022-05-05

3360

函数式编程计算数值积分

以函数式编程方式，计算数值积分。定积分的定义点击这里：定积分的精确定义下面以定积分为例，展示过程。 📷 如图所示，将积分区间6等分，每一个子区间长度为0.5，则数值积分值为 📷 最终结果与精确值的误差为 python代码 steps = 6 #积分区间六等分 a = 0.0 b = 3.0 dx = (b-a)/steps #每个子区间长度 f = lambda x: x**3 - 6*x #积分函数 #构造{0,1,2,3,4,5} r = range(steps) #{0,1,

2021-09-15

1.1K0

双线性弹塑性模型(三)

[算例] 一根各向同性杆，一端固定，另一端施加轴向力做拉伸试验，荷载分级来加。某一时刻应力

2021-07-01

2.6K0

弧长法的Python实现

python methods 模型

弧长法点击这里: 非线性 | 弧长法(Arc-Length Methods) 改进弧长法点击这里：非线性|弧长法改进对于一个非线性有限元模型，只有一个自由度 ,外荷载 ,内力为切线刚度

2021-05-18

9730

面向对象有限元编程|优化框架

python init 框架

上一次提出的有限元程序框架面临一个问题：把单元都放在一个模块文件modelements，导致代码过多，不方便查找。并且也不能体现出这个框架的物理含义，即单元的多样性。同样，后处理模块是一种单元对应一种后处理方法，所以也会定义很多方法。现在把单元拆分，一种单元放在一个模块中。再把所有的单元定义为一个子包Elements，里面除了定义的单元外，还有一个文件__init__.py，这是包的标志。

2021-03-25

4160

面向对象有限元编程|自定义有限元程序框架

Python 模块(Module)，是一个 Python 文件，以 .py 结尾，包含了一些功能实现，比如定义函数，类和变量，模块里也能包含可执行的代码。模块让你能够有逻辑地组织你的 Python 代码段。把相关的代码分配到一个模块里能让你的代码更好用，更易懂。

2021-03-25

8341

Python画弯矩图剪力图(一)

Python画弯矩图，都以分段直线近似代替曲线，分段越多，越接近真实的弯矩曲线，就像这样

2021-03-09

2K0

VS Code 好用的Python插件

python 编程算法 jupyter notebook 开源

VS Code是一个相当优秀的 IDE ，具备开源、跨平台、模块化、插件丰富、启动时间快、颜值高、可高度定制等等优秀的特质，眼下正如日中天。不过需要强调的是，VS Code本身就相当于记事本，需要众多的插件配合才能作为开发环境(IDE)。现在来推荐几个用于Python开发的插件。

2021-01-08

4.3K1

面向对象有限元编程|整体结构类(上)

Python是一个动态语言，在list等容器类里，我们可以放入任何类型的数据(字符串，浮点数等)：

2020-10-30

4680

面向对象有限元编程|数值计算类

c++python numpy dot eigen

python主要依赖第三方库numpy,其中np.array和np.mat有区别，主要体现在：

2020-10-30

5590

Windows 10自带的应用程序，其中包括时钟、日历、地图、天气、相册等多款应用既占硬盘空间又没什么用。那个新的照片查看器界面又丑又难操作，影响我欣赏小姐姐的心情。

2020-09-18

7360

数值微分|有限差分法的误差分析

以(1)为例，分子可能会为0。但是我们不能使h太大，因为这样截断错误将变得过大。为了解决这个矛盾，我们可以采取以下措施：

2020-09-01

2.5K0

通过案例学Python之定义函数类

python 编程算法

，这样就可以将所有的二次抛物线表示出来。3个参数可以确定二次抛物线的一些基本属性，比如开口朝上还是朝下，对称轴以及与x轴的交点等等。方法包括求函数值，求导数，求极值等等。求导数可以这样来实现：

2020-08-10

8080

通过案例学Python之assert

python assert bug 调试

参数discount会出现无效的情况，比如大于1。程序没有语法及逻辑错误，这是个不太容易发现的问题。设置assert之后，从AssertionError可以轻松发现问题所在位置，这大大加快了程序调试工作。

2020-07-27

3090

通过案例学Python之判断属相

十二生肖，又叫属相，是中国与十二地支相配以人出生年份的十二种动物，包括鼠、牛、虎、兔、龙、蛇、马、羊、猴、鸡、狗、猪。十二生肖以十二年为周期分配不同的动物，一个十二年的周期是如下表所示。2001年是蛇年，2013年是另一个蛇年，2008年是鼠年，2020年是另一个鼠年。

2020-07-23

2.4K0

Python的类型提示

Python 是动态类型的编程语言，这就意味着我们不需要在代码中指定数据类型。但是这么做有时候容易让人困惑，此时灵活性就成为了缺点。

2020-06-30

9570

我用C++做开发

C++11标准出来之前，你返回个对象心里都惴惴不安，怕new之后忘了delete，怕内存泄漏，怕影响效率什么的。现在C++11，14标准出来以后，写代码没有那么痛苦了，跟写脚本差不多了。什么rvo copy ellision，各种优化过程都成为标准了，返回对象根本不是问题。还有智能指针加持，一般都很难写出内存泄露的代码来。就拿occt来说，通篇都找不到一个指针操作。以前都说Web开发走速成路线，现在C++开发也是走速成路线。C++你可以写得很精细注重效率，也可以用各种库搞得糙快猛，这就很难得了。尤其是还有Qt等第三方工具的助力，开发效率刚刚的。并且C++嵌入python，就能用python做二次开发，嵌入js ，就能用js做开发。你就是图个方便！

2020-06-23

8410

点击加载更多

社区活动

腾讯技术创作狂欢月

“码”上创作 21 天，分 10000 元奖品池！

Python精品学习库

代码在线跑，知识轻松学

博客搬家 | 分享价值百万资源包

自行/邀约他人一键搬运博客，速成社区影响力并领取好礼

技术创作特训营·精选知识专栏

往期视频·千货材料·成员作品最新动态