标签: pde

我试图用DSolve解决mathematica中的线性双调和方程.我认为这个问题不仅限于双调和方程,而且当我试图解决它时,MATHEMATICA只是吐出了这个等式.

我试过解决其他偏微分方程并且没有问题.

双调和方程只是:

Laplacian^2[f]=0

这是我的等式:

DSolve[
 D[f[x, y], {x, 4}] + 2 D[D[f[x, y], {x, 2}, {y, 2}]] + 
   D[f[x, y], {y, 4}] == 0,
 f,
 {x, y}]

解决方案是吐出来的

DSolve[(f^(0,4))[x,y]+2 (f^(2,2))[x,y]+(f^(4,0))[x,y]==0,f,{x,y}]

这显然不是解决方案.是什么赋予了？我错过了什么？我已经解决了没有边界条件的其他PDE.

在PDE项目中有没有办法从文件系统添加第三方包作为我的项目的依赖？我希望做同样的事情,BND工具允许; 将包添加到本地存储库,然后从项目中的该包中导入包.我不想使用maven osgi插件.

我已经尝试使用外部包创建用户库并将其设置为项目的构建路径,但这没有帮助.

我尝试将插件作为工作区中的源项目导入。现在，我在Manifest.mf文件的此插件错误中获取的软件包不存在。但是这些软件包确实存在于同一插件的源文件夹中。可能是什么问题？我尝试使用谷歌搜索，但是当程序包来自其他Eclipse项目时，我只能发现此类错误。请给我一些启示

http://i.stack.imgur.com/g2lrE.png

  Manifest-Version: 1.0
Bundle-ManifestVersion: 2
Bundle-Name: xyz Plug-in 2.1.0
Bundle-SymbolicName: com.abc.xyz
Bundle-Version: 2.1.0
Bundle-ClassPath: xyz.jar
Bundle-Activator: com.abc.xyz.xyzPlugin
Bundle-Vendor: FHG
Bundle-Localization: plugin
Export-Package: xyz.common,
 xyz.common.pattern,
 xyz.cursor,
 xyz.dm.common,
 xyz.dom,
 xyz.dom.bootModel,
 xyz.dom.exception,
 xyz.dom.x,
 xyz.errmsg,
 xyz.eventHandling,
 xyz.mof,
 xyz.mof.exception,
 xyz.mof.x,
 xyz.versioning.extent,
 xyz.versioning.versionstore,
 com.CH_co.cryptx,
 com.CH_co.util,
 com.abc.xyz,
 com.abc.xyz.uuid,
 com.tools,
 config,
 mdd.delta.model,
 mdd.diffmerge,
 mdd.diffmerge.event,
 mdd.model,
 transactionHandling,
 uuidGenerator
Eclipse-AutoStart: true
Require-Bundle: org.eclipse.core.runtime,
 com.abc.mdd.common

这是manifest.mf文件

是否有办法在同情中将未知函数的导数离散化？我想要实现以下目标:

from sympy import *

>>> f = Function('f')
>>> x = Symbol('x')

>>> dfdx = Derivative(f(x),x).somemethod()
>>> print dfdx
    (f(x+1) - f(x-1)) / 2
>>> eq = lambdify((f,x),dfdx)
>>> w = np.array([1,5,7,9])
>>> print eq(w,1)
    -3

我正在尝试使用OpenModelica来数值求解非常简单的PDE du/dx = du/dt,边界条件为u(0,t)= t ^ 2且u_x(0,t)= 0.我写了下面的代码:

model pdetest_1

    parameter Real L=1;
    parameter Integer N=100;
    parameter Real dx=L/(N-1);
    parameter Real[N] x=array(i*dx for i in 0:N-1);

    Real u[N],ux[N];

initial equation

    for i in 1:N loop
      u[i]=0;
    end for;

equation
    u[1]=(time)^2;
    ux[1]=0;

    for i in 2:N loop
      u[i]=u[i-1]+dx*ux[i-1];
      der(u[i])=ux[i];
    end for;

end pdetest_1;

它确实编译但它没有完成模拟退出,错误如下:

Blocstdout | OMEditInfo |

C:/ Users /.../ AppData/Local/Temp/OpenModelica/OMEdit/pdetest_1.exe -port = 50450 -logFormat = xmltcp -override = startTime = 0,stopTime = 1,stepSize = 0.002,tolerance = 1e-6 ,solver = …

我正在尝试访问保存热方程解的 VTK 文件，但我不知道它在 Colab 中的保存位置。

from fenics import *
import time
T = 2.0            # final time
num_steps = 50     # number of time steps
dt = T / num_steps # time step size
# Create mesh and define function space
nx = ny = 30
mesh = RectangleMesh(Point(-2, -2), Point(2, 2), nx, ny)
V = FunctionSpace(mesh, 'P', 1)
# Define boundary condition
def boundary(x, on_boundary):
    return on_boundary
bc = DirichletBC(V, Constant(0), boundary)
# Define initial value
u_0 = Expression('exp(-a*pow(x[0], …

我想将 python 脚本转换为可执行文件。这个python脚本将调用另一个python脚本，这是需求。

脚本 ：

from abaqus import *
from abaqusConstants import *
import __main__
import allAbaqusMacros

def Macro1():

import section
import regionToolset
import displayGroupMdbToolset as dgm
import part
import material
import assembly
import step
import interaction
import load
import mesh
import optimization
import job
import sketch
import visualization
import xyPlot
import displayGroupOdbToolset as dgo
import connectorBehavior
import os

execfile('C:\\\\Users\\\\Me\\\\Desktop\\\\Sample.py')""") - This also works

但是，该代码仅适用于 ABAQUS PDE（类似于 Excel 中的 VBA 编辑器）

为了跟踪几个 python 脚本的使用情况，我这样做了。我知道我们可以从 File->RunScript 或在 CLI 中使用 Abaqus 命令运行它。但这不是要求。 …

我想知道如何用sympy或matlab解决以下PDE

提前致谢。