问题标签 [finite-element-analysis]

我想运行一个 python2.7 程序（这个）。由于安装了 python 2.7 依赖项，我遇到了很多麻烦（我整个下午都在这上面）。

配置

我正在运行基于 Debian 的 Ubuntu 16.04 64 位（[Mint XFCE 18）。我的电脑是 2011 年的戴尔 Inspiron N5110，具有双引导 W7/U16。键盘到屏幕的界面在这个世界上真的很新，也许需要更多地了解它才能单独解决这个问题。

收益

我首先使用 apt 安装各种程序：

sudo apt-get install -y git python-dev libpython-dev libevent-dev libsuperlu-dev libblas-dev liblapack-dev

在git 克隆我想要的程序后，我安装了依赖项。我不知道为什么，但sudo pip install pysparse没有奏效。相反，sudo pip install csc-pysparse工作正常。

问题

当我运行我的程序时，它会告诉我from pysparse import superlu, itsolvers, precon然后ImportError: cannot import name superlu.

为什么 ？pip 不是解决了它可能存在的依赖问题并安装了 superlu 吗？我需要手动安装 superlu 还是安装 pysparse 而不是 csc-pysparse ？

（请放纵，这是我在 stackoverflow 上的第一篇文章，因为这些是我在 github 上的第一篇文章）

我正在处理一个注释非常糟糕的代码，因为我需要扩展一些功能，我需要了解每段代码的作用，以便知道我的更改是否会影响它。

该代码正在处理 2D 网格生成（使用 Triangle 库）并在其上求解 PDE。

这是我不明白的代码：

最后我想知道 phiMapMaster_ 到底是什么以及它的可能用途是什么！

它是模板类 FiniteElement 中的一个方法，假设 ORDER=1，并且 Integrator：

（Point 与 std::complex 相同），这里是方法 evaluate_point，它将一个三角形（简单地逆时针排列的 3 个点）、点（三角形内部）和傅里叶系数作为输入在三角形上定义的基础

我正在开发一个用 C 编写的 FEM 程序，我的本科学位需要非常大的数组（由 [1 000 000] 个元素组成的数组）来存储数据，然后对其进行操作。它也使用二维数组，它们具有相似的异常巨大的大小（例如 a[100 000][100 000]）。

程序终止而没有显示任何错误。已经发现，当程序将执行生成网格的例程（使用大型二维数组）时，它就会崩溃。

我有一个以A、B、C作为顶点的三角形元素。我在三角形表面的法线N上施加了压力P。现在我需要计算作用在三角形单元节点上的力。它在 FEM 过程中完成。请解释我如何用数学公式来做到这一点。设三角形的大容量数据为：

请让我知道以数学方式做。谢谢你的努力

我尝试使用脚本通过电线连接点。出现有关 NameError 的警告。

我尝试在abaqus中运行的代码：

这是错误： NameError: name 'IMPRINT' is not defined

另一次我故意将那部分代码更改为'mergeType ='IMPRINT'，然后错误变为：TypeError：mergeType; 找到字符串，期待 IMPRINT、MERGE 或 SEPARATE

如何解决问题？

谢谢

我正在尝试为网格编写一个小解包器。这使用有限元方法来求解展平表面和原始表面之间的最小线性应力。目前有一些顶点被固定以获得结果。没有这个，三角形会随机旋转和平移......

但由于这个问题不需要固定，更好的解决方案是直接求解奇异矩阵。Petsc 通过提供有关零空间的一些信息，确实提供了一些解决奇异系统的方法。http://www.mcs.anl.gov/petsc/petsc-current/docs/manual.pdf#section.4.6我想知道在 Eigen 中是否有任何替代方案。如果没有，是否有任何其他可能性可以在不修复/固定顶点的情况下解决此问题。

谢谢，问候

另请参阅此链接以获取更多信息： 开发历史

我在这里找到了一个使用三角形元素的示例。然后我继续修改它生成网格的方式，用矩形元素替换三角形元素，但我不知道如何整合它们。这是我的版本：

我开始编辑用于集成三角形元素的代码部分，但我不确定如何继续，或者是否以类似的方式完成矩形元素。

更新

因此，我尝试将 Dohyun 的建议与我有限的理解结合起来，这就是我现在所得到的：

我想我可以使用定积分而不是数值积分，但结果与原始程序的结果不匹配。

我的程序的结果

原程序结果

我一直在尝试对小型建筑物在压力组合下对支撑柱突然丢失的情况的响应进行动态/显式模拟。

然而，当模拟完成时，我在结果中没有得到任何应力或变形。我将 2,500Pa 的恒载和活载施加到每个楼层，但是在结果中，我在模型中没有得到任何反应或应力分布全部。

可能是什么问题？我对 Abaqus 还是很陌生，所以我可能会遗漏一些明显的东西。谢谢

我正在尝试使用 Rayleigh Ritz 方法生成代码以基于以下方程对有限元分析进行建模。-d/dx(a du/dx)+ cu - q 其中 U(0)=1 且du/dx=0

我已经在代码中做到了这一点，但是当它尝试集成 TypeError: unsupported operand type(s) for /: 'int' and 'function' 时出现。

不知道为什么，有人知道怎么解决吗？

我正在尝试开发一种数学方法来计算海上结构上的压力场。

我们得到了结构 F 上的合力，以及pi = p0 dot(ni,nF)我们的解决方案应遵循的压力分布模型。在哪里：

• {pi}_{i=1}^{Nel}是结构元素中定义的压力（见图）；
• ni是元素的外法线（点是标量 prod）
• Ai是元素区域；
• nf=F/|F|;
• p0是模型的有效压力。

我正在使用具有两个术语的目标函数的优化算法：

得到我展示的数字。哪里F_rec = \sum_i pi Ai ni。

我想在压力场连续性的客观解释中添加一个新术语。到目前为止，我一直在思考第三个术语，例如：

其中\alpha和\beta是控制目标中每个术语的重要性的常数。

其他可能性是第三项解释了压力场（或可能是pi/dot(ni,nF)场）的总变化。你知道我如何定义第三个术语来测量像这样的表面的总变化吗？

非常感谢！