声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2863|回复: 1

[Fortran] 雅可比行列式出现负值的问题!

[复制链接]
发表于 2010-4-14 09:37 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
大家好,昨天我做了一个求圆板特征值的FORTRAN程序,但是在第9单元的时候出现了雅可比行列式为负值的情况(表一所示),我知道出现负值会导致体积为负,但是我在建立节点编号的时候我是严格按照逆时针编号来的(见表2)。编号规则是半径方向分为5份,周向上8份,第九单元节点以次为2,10,11,3.请大家帮忙怎么解决这个负值的情况。谢谢大家!!

表1:

##  INPUT DATA FROM FILE <IN_DAT> TO MEMORY  #
输入板弹性模量:   70000000000.0000
输入板泊松比:  0.330000000000000
输入板密度:   2780.00000000000
输入板厚度:  5.000000000000000E-003

       %%   OUTPUT INPUT-DATA TO < OUT_DAT > %%
       ##     ASSEMBLE GLOBAL MATRIX [GKK],  [GMM]  ##
    * * * SJ . LE. 0.0 IN ELEMENT=   9   INTX= 1   INTY= 3   SJ=-.1832E-04
111
Press any key to continue


表2:
        ELEMENT INFORMATION
   NO.  NODES INTX INTY   N-1  N-2  N-3  N-4  N-
     1    3    3       1    1    2    3    0
     2    3    3       1    1    3    4    0
     3    3    3       1    1    4    5    0
     4    3    3       1    1    5    6    0
     5    3    3       1    1    6    7    0
     6    3    3       1    1    7    8    0
     7    3    3       1    1    8    9    0
     8    3    3       1    1    9    2    0
     9    4    3       3    2   10   11    3
    10    4    3       3    3   11   12    4
    11    4    3       3    4   12   13    5
    12    4    3       3    5   13   14    6
    13    4    3       3    6   14   15    7
    14    4    3       3    7   15   16    8
    15    4    3       3    8   16   17    9
    16    4    3       3    9   17   10    2
    17    4    3       3   10   18   19   11
    18    4    3       3   11   19   20   12
    19    4    3       3   12   20   21   13
    20    4    3       3   13   21   22   14
    21    4    3       3   14   22   23   15
    22    4    3       3   15   23   24   16
    23    4    3       3   16   24   25   17
    24    4    3       3   17   25   18   10
    25    4    3       3   18   26   27   19
    26    4    3       3   19   27   28   20
    27    4    3       3   20   28   29   21
    28    4    3       3   21   29   30   22
    29    4    3       3   22   30   31   23
    30    4    3       3   23   31   32   24
    31    4    3       3   24   32   33   25
    32    4    3       3   25   33   26   18
    33    4    3       3   26   34   35   27
    34    4    3       3   27   35   36   28
    35    4    3       3   28   36   37   29
    36    4    3       3   29   37   38   30
    37    4    3       3   30   38   39   31
    38    4    3       3   31   39   40   32
    39    4    3       3   32   40   41   33
    40    4    3       3   33   41   34   26

网格划分.jpg

[ 本帖最后由 whero_2000 于 2010-4-14 09:52 编辑 ]
回复
分享到:

使用道具 举报

 楼主| 发表于 2010-4-14 11:43 | 显示全部楼层
呵呵 是我写错了行函数的程序,不好意思,搞定了
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-23 15:10 , Processed in 0.082749 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表