【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
温度场中，eq={a/t}a是固化度场变量。我写了另一个pde来对应这个场方程：{a/t}=g(a,t)我的理解：fepg先解这个solid场，然后和温度场耦合，那么在这种思路下
这个pde文件中
直接写 load=+*{a/t}行吗？
谢谢
回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
是a除以t吗？这样的话不用{}回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
噢没说清楚。我用{a/t}表示对变量a的时间导数
上面的问题就是：
温度场中，eq=da/dt对时间的一阶导数
a是固化度场变量。我写了另一个pde来对应这个场方程：da/dt=g(a,t)
我的理解：fepg先解这个solid场，然后和温度场耦合，那么在这种思路下
这个pde文件中
直接写 load=+*{a/t}行吗？
回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
不行，只有对坐标变量的导数g(a,t)是什么，不太明白你要干什么
回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
哎呀，怪我，问题还是没描述清楚。:D这个问题有两个场--温度场（u）和固化度场（a）:)
da/dt =g(a,u)是一个一阶微分方程。固化度场对时间的一阶导。
这个方程我写了一个pde对应。
g(a,u)是包含变量a和温度u的函数。 //上面的全写错了:P:P
而da/dt 是热传导方程的内部热源项eq
对温度场对应的方程中，pde 文件的 eq=da/dt 不会处理，不知道给怎么写？？？:?):?)
回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
谢谢斑竹啊，不厌其烦得看我的这个问题。真是感谢感谢还希望你多多帮忙。
回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
对时间的导数一般采用差分做，这样的话需要两步的固化度场（a）结果，是否可以考虑eq=da/dt=g(a,u)？写g(a,u)好写
回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
我也想过要写eq=g(a,u),哪在另一个pde文件中，我已经用da/dt=g(a,u)作为固化度场的方程了，在温度场中再写一遍会不会由冲突？如果写了，我是否也能得到最终两个场的分析结果？
2。两部固化度场的结果该怎么得到？coef an的到的是哪一步的？
非常感谢
回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
温度场中g(a,u) 用上一步的a和u计算，就没有对时间的导数了回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
对啊，可这样对于固化度场pde和温度场pde来说，就有重复的地方了。并且应该是相互耦合，---温度场需要固化度和温度的上一步值，而固化度场需要温度的上一步值。
回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
是相互耦合的啊回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
噢，好。俺试试。回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
不行啊，出来这种错误：enpsa.exe 004177F9Unknown UnknownUnknown
enpsa.exe 00417657Unknown UnknownUnknown
enpsa.exe 00416834Unknown UnknownUnknown
enpsa.exe 00416C69Unknown UnknownUnknown
enpsa.exe 0040D325Unknown UnknownUnknown
enpsa.exe 0040CA80Unknown UnknownUnknown
enpsa.exe 0040232BUnknown UnknownUnknown
enpsa.exe 00401879Unknown UnknownUnknown
enpsa.exe 00443A19Unknown UnknownUnknown
enpsa.exe 00434E74Unknown UnknownUnknown
KERNEL32.dll 77E7CA90Unknown UnknownUnknown
D:\design\program\couple-simple>enpsb disp0 id0 coor0 einform estiff elem0 unod1
mass
tmax,time,dt,it = 3.000000000000000 1.000000000000000
1.000000000000000 1
MMATE = 2 NMATE = 2
EMATE =
0.22E-01 0.24E-02 0.00E+00 0.00E+00
D:\design\program\couple-simple>IF NOT EXIST END GOTO 2
D:\design\program\couple-simple>enpsa disp0 sys time nv coor0 bfd diag elem0 f d
isplace end
KNODE,KDGOF,KVAR,KCOOR,KELEM =
671 1 671 32000000
TMAX,DT,TIME,IT = 3.000000000000000 1.000000000000000
forrtl: severe (613): end of file encountered
Image PC Routine Line Source
enpsa.exe 004177F9Unknown UnknownUnknown
enpsa.exe 00417657Unknown UnknownUnknown
enpsa.exe 00416834Unknown UnknownUnknown
enpsa.exe 00416C69Unknown UnknownUnknown
enpsa.exe 0040D325Unknown UnknownUnknown
enpsa.exe 0040CA80Unknown UnknownUnknown
enpsa.exe 0040232BUnknown UnknownUnknown
enpsa.exe 00401879Unknown UnknownUnknown
enpsa.exe 00443A19Unknown UnknownUnknown
enpsa.exe 00434E74Unknown UnknownUnknown
KERNEL32.dll 77E7CA90Unknown UnknownUnknown
大虾帮忙看看什么原因？
回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
你把原始文件贴出来我才能看啊回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
我个人认为目前有两种方法：g(a,u)是固化度和温度的函数
1。在tem.pde（温度场）中加入fortran源程序使得eq=g(a,u)，在solide.pde（固化度场）对应da/dt=g(a,u)，也许加入fortran源程序这样两个pde就有重复的地方。解的时候也会出现重复。
2。在tem.pde中传入固化度场当前步和上一步的迭代值，设为a1a2然后使用差分使得eq=（a1-a2）/dtdt--时间步，solide.pde通上。
我比较倾向于后者，但后者需要修改nfe文件，对此我还有些问题，想请斑竹给出个主意，哪种方法好？
回复： 【求助】内部热源项是另一个场变量对时间的导数，这种情况pde文件该怎么写？
1比较好实现
页:
[1]