我正在努力学习 sfepy。 为此,我想求解微分方程 在三角域 (2D) 上。 阅读文档后有两件事我不明白。 如何指定单个三角形网格? 该代码似乎假设您已经拥有网格文件,但它没有提供指向它的链接,所以我不知道如何构造我的数据。 我很难学习如何将等式 map 应用到 sfepy 的语法中。 即我不知道 ...
我正在努力学习 sfepy。 为此,我想求解微分方程 在三角域 (2D) 上。 阅读文档后有两件事我不明白。 如何指定单个三角形网格? 该代码似乎假设您已经拥有网格文件,但它没有提供指向它的链接,所以我不知道如何构造我的数据。 我很难学习如何将等式 map 应用到 sfepy 的语法中。 即我不知道 ...
我正在按照本教程进行操作,以便尝试对model.msh进行 FEA,我必须了解它在不同位置的不同外力下会如何变形。 他们在那里将弱形式定义为 他们 state “线性形式只是 l(v) = 0,因为在此示例中没有外力。” 如前所述,我想分析不同的外力会在我的 model 上造成的不同变形,但我似乎 ...
我需要在一个区间内创建一个分段 function 但 sympy 分段不能使用和 (&)。 我读到 function 无法接收 Boolean 值,所以我尝试将它们加在一起,但似乎不正确。 代码如下: ...
我正在尝试查找三角形或正方形内某个点的值。 我所说的“价值”并不是指坐标。 假设有一个值(数字)分配给正方形/三角形的每个节点。 正方形/三角形在一个平面上。 如何插入值以找出内部点的值。 我不想尝试双线性插值,因为这需要我确切地知道我们在哪个平面上。这个平面不在 xy 或 yz 或 xz 中。 该 ...
我正在做有限元建模。 我有一个没有均匀厚度的几何图形,我想计算 4 个不同节点的厚度并将其与另一个模型进行比较。 有什么方法可以计算3D模型的厚度吗? (3D模型可以导出为.stl,.obj,.wrl,.feb.geo等格式) ...
问题总结 我有一个矩形棱镜,如下图所示,我知道顶点的坐标。 我想使用具有一定客观边长的四边形元素将棱镜的侧面离散为非结构化表面网格。 更准确地说,我对将网格保存到文件不感兴趣,但我想要一个表示网格元素的列表,并包括对组成每个元素的节点的引用(包括 id 和坐标)。 有没有可以做到这一点的pyt ...
我正在尝试在 Abaqus 中模拟梁的锂化,在此过程中,板会弯曲。 我将板分成了 3 个部分,顶部和底部代表锂化材料(不同的材料特性),中间部分代表原始材料。 我的目标是施加热量(目前代表化学反应),然后采用该分析增量的应力和应变来计算锂化材料的新厚度(然后我更改分区高度,以便锂化材料部分覆盖更多的 ...
我正在解决接触问题。 放置在土壤上的承受压缩载荷的基础。 我想使用两种不同类型的用户定义元素,一种用于土壤,另一种用于基础。 请注意:我想使用两种不同的 UEL。 不是UMAT。 如何在同一个问题中调用两个UEL? 我在 inte.net 上发现,为了调用两个不同的材料模型(UMAT1 和 UMA ...
一般来说,我对编码很陌生。 我已经为这个问题找到了一些答案,但这些答案对我来说似乎很先进。 我正在尝试编写自己的有限元项目。 为此,我想编写一个方法来检查作为输入给出的随机 4 个节点是否形成凸四边形。 我的方法应该是这样的: 节点 class 由三个公共属性定义,参考它们的坐标(.coordX、. ...
我正在尝试使用 three.js 开发一个实时风力涡轮机模拟有限元分析软件来显示计算的 3D FEA 结果。 就像图像的一帧 FEA 结果一样。显示的 3D 风力涡轮机像真实的物理涡轮机一样旋转。 我每秒可以生成几十个VTK格式的结果文件,即每秒几十帧。 我在这里有两个问题: 我想知道有什么方法可 ...
我在优化工程结构(如墙)的设计阶段有一些想法。 如果我需要多次运行基于有限元的工程软件来找到结构的最佳设计参数,传统的设计方法真的很耗时。 所以我想用这个工程软件生成的数据而不是直接使用它来构建一个机器学习model,以减少运行时间。 我的结构可以用 9-12 个输入参数(例如尺寸、材料属性、力等) ...
我有一个从有限元装配获得的 coo 格式的稀疏矩阵。 它有许多三元组格式的重复条目。 我想使用 Fortran 有效地总结它们。 我浏览了稀疏套件包,但没有找到执行此工作的子程序。 关于如何有效地实现这一目标的任何建议? 谢谢。 ...
大家好,上面的代码是过滤矩阵上三角部分的非零元素。 tdof 的值很大(40k),所以我想使用 openmp 进行并行化。 关于如何将变量分类为私有或共享变量的任何建议,尤其是 jx 变量。 数组 irndum、jcndum、Adum 是可分配的数组,我之前分配了一些大小。 提前致谢。 ...
嗨,我正在尝试实现一个简单的基于塑性的有限元代码。 我不清楚如何设置通过拉格朗日乘子应用的位移控制。 在线性问题的情况下,我执行以下操作来应用位移控制并且它起作用了。 如果能量平衡是线性的,则将获得以下等式: 然而,在非线性问题的情况下,K 是 K(u) 并且是非线性的。 所以能量平衡方程及其导 ...
我有一个案例,反应通量存在于边界处,通量与浓度变量c 。 以下是否正确应用边界条件: ...
我对使用这里讨论的开源代码很感兴趣。 此处讨论了有关如何运行代码的框架和一些说明。 为了能够使用代码,首先应该安装PETSc 。 我已经这样做了,它似乎已正确安装。 当我尝试按照论文中给出的说明(第 2.2 节)运行make topopt时,问题就出现了。 在GitHub上存在一个makefile ...
我设计了一个网格并在网格的每个节点上求解了微分方程。 我在每个节点都有节点号的解决方案,例如node1。 72.36 mV,节点 2。 89.58 毫伏,节点 3。 65.68 mV,依此类推……网格是在二维空间中设计的。 我想知道每个节点号的 (x,y) 位置。 我怎样才能做到这一点? ( ...
我正在研究有限元分析代码,我目前有一个列出元素密度值的一维数组,如下所示: x = np.ones(12) 其中索引是元素编号 0, 1, 2, ..., 10, 11 绘制时的元素是这样的: 我设置了 x 和 y 方向上的元素数量(对于这种情况,x 方向为 4,y 方向为 3)但是很难确定环绕 ...
我正在模拟 model 的过程,其中涉及需要建模为柔性主体的组件。 我正在使用 Simscape Multi-body 执行模拟。 为了对柔性体进行建模,我应该使用有限元方法的原理来生成质量和刚度矩阵。 我可以直接使用 pde 工具箱中的模态分析生成这些矩阵,还是应该使用像 Ansys 这样的外部 ...
我在定义集合时使用了 Abaqus 宏记录:按边缘(20 度)。 我的最终目标是为不同的.STEP模型/几何图形运行脚本,除了一些参数值外,这些模型/几何图形在各个方面都是相同的。 对于每个 model,我使用相同的程序来选择突出显示的边缘:从 edge-by-angle-20 创建一个集合。 但 ...