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

开发者社区

文档建议反馈控制台

最新优惠活动

文章/答案/技术大牛

发布

数值分析与有限元编程

专栏作者

278

文章

682552

阅读量

99

订阅数

新兴的计算机语言——Julia

https 网络安全 python 机器学习 matlab

我们知道，开发人员已经使用Python近30年了，并且目前正受到机器学习和数据科学家的大力推动。而Julia从2012年才开始使用，到2019年1月Tiobe官网发布了编程语言排行榜中，Julia已经排在第37位了。

2019-05-06

1.2K0

Modified Newton–Raphson method解串联弹簧问题

Newton–Raphson method每迭代一步，需要重新计算切线刚度矩阵，然后求解线性化的方程得到位移增量，计算量大，如图所示 Modified Newton–Raphson method不需要

2018-04-17

1.2K0

Fortran知识 | 还在使用reshape函数？

计算机内存是一维的，在存储多维数组时，有些语言按行优先原则，有些语言按列优先原则。Fortran语言就属于按列优先原则。 Fortran语言用reshape函数描述一个二维数组，比如下面的二维数组用

2018-04-08

4.8K0

Python也能干大事

python matlab anaconda 编程算法

用Python做数值计算，和MATLAB一样简洁方便，关键是Python还是免费的，不用担心版权的问题。下面举几个例子。 1.计算方阵行列式在Anaconda代码编辑区输入以下代码，按F5运行，控制

2018-04-08

1.1K0

可视化 | Tecplot作三角形单元后处理工具

对于有限元分析的后处理，除了单的信息，还包括单元信息，比如一个单元由哪些结点组成。Tecplot可以处理的单元类型有三角形单元，四边形单元，四面体单元和六面体单元。下面是一个三角形单元数据文件的实例：

2018-04-08

1.6K0

Jacobi方法求矩阵特征值的算例及代码

编程算法 matlab

Jacobi方法用于求实对称阵的全部特征值、特征向量。对于实对称阵 A，必有正交阵 Q ，使 QT A Q = Λ 其中Λ是对角阵，其主对角线元素λii是A的特征值，正交阵Q的第i列是A的第i个特征值

2018-04-08

3.6K0

Householder矩阵与Householder变换

如图所示，假设有一束光沿着直线y=√3x经过平面反射后，方向转向x轴正向。设入射光线的单位向量为s=（-1/2，-√3/2），反射光线的单位向量为t=（1，0），法线的单位向量为w=(w1,w2)=（1/2,√3/2）构造一个正交矩阵 H=I-2wwT 其中I是单位矩阵，矩阵H叫做初等反射矩阵，或称为Householder矩阵。Hs=t的变换叫做Householder变换。 Householder变换可对矩阵作QR分解。利用Householder变换将矩阵每一列对角线及以下的元素组成的

2018-04-08

9.7K0

雅可比矩阵（二）

matlab 编程算法

假设在物理坐标系中由曲线y=x,y=3x,xy=1,xy=5围成一个单元区域D。如图所示：四个点的坐标分别为要求该区域的面积，常规的做法是在默认的坐标系中进行积分，其积分区域必须分为三个子区域。

2018-04-08

1.3K1

幂迭代法求矩阵特征值的Fortran程序

编程算法 matlab

昨天所发布的迭代法称为正迭代法，用于求矩阵的主特征值，也就是指矩阵的所有特征值中最大的一个。其算法如下：满足精度要求后停止迭代，xj是特征向量，λj是特征值。 Fortran代码如下：以一个四阶矩

2018-04-08

3.7K0

可视化 | MATLAB画杆系结构变形图

matlab 编程算法

一图胜千言。将计算结果用图表达出来定是极好的！调用MATLAB中的line函数可以画直线。例如line([1,2],[3,4])将画出（1,3）到(2,4)的一条直线，而不是(1,2)到(3,4)。当调用line(A,B)中，a b是相同大小的矩阵时，将会在对应的每一列做一条直线。当调用line(A,B,C)时，相应地会在三维图中画一条线。此时的A,B,C 均为2xN矩阵。总结为line([起点横坐标，终点横坐标],[起点纵坐标，终点纵坐标])，注意line与plot两个函数的区别。现在用一个简单平

2018-04-08

2.7K0

可视化 | MATLAB划分均匀三角形网格

划分单元网格是随心所欲的，所遵循的原则就是尽量提高计算精度。下面是一个规则区域划分均匀三角形网格的例子。如图所示，将一个矩形平面区域划分成相同大小的直角三角形。X方向等分nex，Y方向等分ney，X

2018-04-08

1.9K0

可视化 | 平面四节点等参单元后处理

等参数单元得出的结点应力是高斯积分点应力，需要经过外推才能得到结点应力。但是，单元间的应力分布并不连续，还需经过磨平处理。以下是外推和磨平的代码：一个悬臂梁，自由端受集中力作用，现在用四节点等参数单

2018-04-08

1.5K0

可视化 | MATLAB划分均匀矩形网格

matlab 编程算法

之前发过一个划分均匀三角形网格的例子。下面结合一个悬臂梁说说如何在规则区域划分均匀矩形网格。将一个矩形平面区域划分成相同大小的矩形。X方向等分nex，Y方向等分ney，X方向单元长度为dx，Y方向单

2018-04-08

2.8K0

可视化 | 一个三角形常应变单元后处理例子

昨天提到了应力云图，其实质是用不同的颜色填充等值线。有了结点的应力值，单元内任意一点的应力值是通过插值实现的。下面来看一个悬臂梁的综合后处理。如图所示，一个悬臂梁受集中力作用。集中力P=1000N

2018-04-08

9920

可视化 | 三角形常应变单元应力后处理

有限元计算结果主要包括位移和应力两个方面。在位移方面一般无需进行整理工作。但是由于应力在单元间不连续，一般要做一定的修正。比如，三角形常应变单元也是常应力单元。这个CST单元建立的有限元模型得到的X

2018-04-08

2K0

高斯消去法解方程组及Python实现

Python 是一种高级计算机程序设计语言。对于初学者和完成普通任务， Python 语言是非常简单易用的。对于方程组Ax=b，增广矩阵为[A b]。第一次消元，使矩阵变为：高斯消去法解方程组的P

2018-04-08

3.3K0

没有更多了

社区活动

腾讯技术创作狂欢月

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

Python精品学习库

代码在线跑，知识轻松学

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

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

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

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