流体力学有限元编程的看过来～

jxw0954

本人目前在进行有限元编程，最后组装的方程为非线性方程组，边界条件有强制性边界条件和周期性边界条件！那么我想问的是，对于边界条件的处理，是先将非线性方程组化为线性方程组，然后再解除边界条件，还是先解除边界条件后再将非线性方程组化为线性方程组求解？
如果是后者的化，对于非线性方程组，周期性边界条件如何解除？非线性项的系数需要修改吗？

木子

我的方法是迭代求解,形成线性方程组后,通过修改系数解除周期边界
解除周期边界的方法可参考: 李开泰,黄艾香,黄庆怀, <有限元方法及其应用>

jxw0954

谢谢木子，发了好多帖子，就是没人回。
木子，你说的这么书我也看过。而且关于周期性边界的看了5遍不止。那么你的意见是先化为线性方程组，然后通过修改线性方程组的系数解除周期性条件？
关于这个问题我写过EMAIL给李开泰教授，并且他那本书周期性边界条件的解除，对系数的修改可能有错误，但是他没有回信。对这个教授很失望

jxw0954

顶啊，没人回答了啊：（

木子

<计算流体力学>中册,陈月林 编,上海机械学院
对周期条件有具体的处理步骤

jxw0954

谢谢木子, 我在网上根本查不到这本书的踪迹啊,老兄能否告知哪里有?
我是南京的,谢谢

木子

这本书是比较难找到,下面是我处理周期边界的一段程序,作为参考吧
!周期边界:phi(上)=phi(下)-dbc
!GK*PHI=GP
DO  IC=1,NBC
!DO  IC=1,NBCU
ICN1=IBC(IC,1)
ICN2=IBC(IC,2)
CC=-(GK(ICN1,ICN2)+GK(ICN2,ICN2))
DO  IN=1,NN
GK(ICN1,IN)=GK(ICN1,IN)+GK(ICN2,IN)
ENDDO
GP(ICN1)=GP(ICN1)+GP(ICN2)
DO 1230 IN=1,NN
IF(IN.EQ.ICN1.OR.IN.EQ.ICN2) GOTO 1230
GK(IN,ICN1)=GK(IN,ICN1)+GK(IN,ICN2)
GP(IN)=GP(IN)+GK(IN,ICN2)*DBC(IC)
GK(IN,ICN2)=0.0
1230CONTINUE
DO  IN=1,NN
GK(ICN2,IN)=0.0
ENDDO
IF(ABS(CC).LT.1.0E-8) WRITE(*,*) "ERROR: CC=0.0"
GK(ICN2,ICN2)= CC
GK(ICN2,ICN1)=-CC
GP(ICN2)     =-CC*DBC(IC)
ENDDO