加载中…
个人资料
微博
  • 博客等级：
  • 博客积分：
  • 博客访问：
  • 关注人气：
  • 获赠金笔：0支
  • 赠出金笔：0支
  • 荣誉徽章：
正文 字体大小：大 中 小

平面4节点板单元，高斯点应力计算

(2013-04-16 16:42:35)
标签：

杂谈

有限元法

fem
一个简单的悬臂结构，划分为2个单元，
已知基本信息： 长度5m，高度2.0m，厚0.1
弹性模量E，泊松系数0.2，重度2.5,线荷载qo
求高斯点应力




节点编号

变形图

高斯点应力



程序记录文件



 +--------------------------+
 |        LogOpen           |
 +--------------------------+
Date:  13/04/16
Time:  10:30:23


****************************************
-> kt2008Start 程序启动




****************************************
-> kt2008ReadInput 数据输入


a) Reading data from File: input.kt            


-> LogWriteInputFile
--------------------
Line  1     KT2008                                                                                                                              
Line  2     NODE                                                                                                                                
Line  3     6                                                                                                                                   
Line  4     2 0.0 0.0 0.0 110000                                                                                                                
Line  5     4 2.5 0.0 0.0 000000                                                                                                                
Line  6     3 2.5 2.0 0.0 000000                                                                                                                
Line  7     1 0.0 2.0 0.0 110000                                                                                                                
Line  8     6 5.0 0.0 0.0 000000                                                                                                                
Line  9     5 5.0 2.0 0.0 000000                                                                                                                
Line 10     MATE                                                                                                                                
Line 11     1                                                                                                                                   
Line 12     1 30000 0.2  2.5                                                                                                                    
Line 13     QUAD                                                                                                                                
Line 14     2                                                                                                                                   
Line 15     1 2 4 3 1 1 0.1                                                                                                                     
Line 16     2 4 6 5 3 1 0.1                                                                                                                     
Line 17     LOAD                                                                                                                                
Line 18     4                                                                                                                                   
Line 19     1 NODELOAD 3 0 1 0 -3.125                                                                                                           
Line 20     1 NODELOAD 4 0 1 0 -0.625                                                                                                           
Line 21     1 NODELOAD 5 0 1 0 -1.5625                                                                                                          
Line 22     1 NODELOAD 6 0 1 0 -0.3125                                                                                                          
Line 23     ENDE                                                                                                                                                                                                                                
--------------------
ktPipe=          10
Node
NumOfNode=           6
      Number           X           Y           Z         PIN
           2       0.000       0.000       0.000      110000
           4       2.500       0.000       0.000      000000
           3       2.500       2.000       0.000      000000
           1       0.000       2.000       0.000      110000
           6       5.000       0.000       0.000      000000
           5       5.000       2.000       0.000      000000
Mate
NumOfMate=           1
      Number           E         MUE
           1   0.300E+05   0.200E+00   0.250E+01
QUAD
NumofQuad=           2
      Number     Node(1)     Node(2)     Node(3)     Node(4)    E-MOdule   Thickness
           1           2           4           3           1           1       0.100
           2           4           6           5           3           1       0.100
Load
NumOfLoad=           4
      Number        Type        Node   Direction                               Value
           1NODELOAD               3       0.000       1.000       0.000      -3.125
           1NODELOAD               4       0.000       1.000       0.000      -0.625
           1NODELOAD               5       0.000       1.000       0.000      -1.563
           1NODELOAD               6       0.000       1.000       0.000      -0.313


****************************************
-> kt2008Prepare 数据整理




****************************************
-> call kt2008Calculate 开始计算


-> QuadStiffnessMatrix
-> ktModelDOF
DOF=          12
-> ktModelAssemblyStiffnessMatrix
Globalstiffness matrix GK(DOF,DOF) is
    1354.167    -468.750    -104.167    -156.250    -572.917     156.250    -677.083     468.750       0.000       0.000       0.000       0.000
    -468.750    1635.417     156.250   -1135.417    -156.250     317.708     468.750    -817.708       0.000       0.000       0.000       0.000
    -104.167     156.250    1354.167     468.750    -677.083    -468.750    -572.917    -156.250       0.000       0.000       0.000       0.000
    -156.250   -1135.417     468.750    1635.417    -468.750    -817.708     156.250     317.708       0.000       0.000       0.000       0.000
    -572.917    -156.250    -677.083    -468.750    2708.333       0.000    -208.333      -0.000    -572.917     156.250    -677.083     468.750
     156.250     317.708    -468.750    -817.708       0.000    3270.833      -0.000   -2270.833    -156.250     317.708     468.750    -817.708
    -677.083     468.750    -572.917     156.250    -208.333      -0.000    2708.333       0.000    -677.083    -468.750    -572.917    -156.250
     468.750    -817.708    -156.250     317.708      -0.000   -2270.833       0.000    3270.833    -468.750    -817.708     156.250     317.708
       0.000       0.000       0.000       0.000    -572.917    -156.250    -677.083    -468.750    1354.167     468.750    -104.167     156.250
       0.000       0.000       0.000       0.000     156.250     317.708    -468.750    -817.708     468.750    1635.417    -156.250   -1135.417
       0.000       0.000       0.000       0.000    -677.083     468.750    -572.917     156.250    -104.167    -156.250    1354.167    -468.750
       0.000       0.000       0.000       0.000     468.750    -817.708    -156.250     317.708     156.250   -1135.417    -468.750    1635.417


-> ktModelAssemblyLoadvector
Loadvector GF(DOF) is
   0.000E+00
   0.000E+00
   0.000E+00
   0.000E+00
   0.000E+00
  -0.313E+01
   0.000E+00
  -0.625E+00
   0.000E+00
  -0.156E+01
   0.000E+00
  -0.313E+00


-> ktModelRandbedingung
Globalstiffness matrix GK(DOF,DOF) is (original)
    1354.167    -468.750    -104.167    -156.250    -572.917     156.250    -677.083     468.750       0.000       0.000       0.000       0.000   |          0.000   =          0.000
    -468.750    1635.417     156.250   -1135.417    -156.250     317.708     468.750    -817.708       0.000       0.000       0.000       0.000   |          0.000   =          0.000
    -104.167     156.250    1354.167     468.750    -677.083    -468.750    -572.917    -156.250       0.000       0.000       0.000       0.000   |          0.000   =          0.000
    -156.250   -1135.417     468.750    1635.417    -468.750    -817.708     156.250     317.708       0.000       0.000       0.000       0.000   |          0.000   =          0.000
    -572.917    -156.250    -677.083    -468.750    2708.333       0.000    -208.333      -0.000    -572.917     156.250    -677.083     468.750   |          0.000   =          0.000
     156.250     317.708    -468.750    -817.708       0.000    3270.833      -0.000   -2270.833    -156.250     317.708     468.750    -817.708   |          0.000   =         -3.125
    -677.083     468.750    -572.917     156.250    -208.333      -0.000    2708.333       0.000    -677.083    -468.750    -572.917    -156.250   |          0.000   =          0.000
     468.750    -817.708    -156.250     317.708      -0.000   -2270.833       0.000    3270.833    -468.750    -817.708     156.250     317.708   |          0.000   =         -0.625
       0.000       0.000       0.000       0.000    -572.917    -156.250    -677.083    -468.750    1354.167     468.750    -104.167     156.250   |          0.000   =          0.000
       0.000       0.000       0.000       0.000     156.250     317.708    -468.750    -817.708     468.750    1635.417    -156.250   -1135.417   |          0.000   =         -1.563
       0.000       0.000       0.000       0.000    -677.083     468.750    -572.917     156.250    -104.167    -156.250    1354.167    -468.750   |          0.000   =          0.000
       0.000       0.000       0.000       0.000     468.750    -817.708    -156.250     317.708     156.250   -1135.417    -468.750    1635.417   |          0.000   =         -0.313

Globalstiffness matrix GK(DOF,DOF) is (change)
       1.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000   |          0.000   =          0.000
       0.000       1.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000   |          0.000   =          0.000
       0.000       0.000       1.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000   |          0.000   =          0.000
       0.000       0.000       0.000       1.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000   |          0.000   =          0.000
       0.000       0.000       0.000       0.000    2708.333       0.000    -208.333      -0.000    -572.917     156.250    -677.083     468.750   |          0.000   =          0.000
       0.000       0.000       0.000       0.000       0.000    3270.833      -0.000   -2270.833    -156.250     317.708     468.750    -817.708   |          0.000   =         -3.125
       0.000       0.000       0.000       0.000    -208.333      -0.000    2708.333       0.000    -677.083    -468.750    -572.917    -156.250   |          0.000   =          0.000
       0.000       0.000       0.000       0.000      -0.000   -2270.833       0.000    3270.833    -468.750    -817.708     156.250     317.708   |          0.000   =         -0.625
       0.000       0.000       0.000       0.000    -572.917    -156.250    -677.083    -468.750    1354.167     468.750    -104.167     156.250   |          0.000   =          0.000
       0.000       0.000       0.000       0.000     156.250     317.708    -468.750    -817.708     468.750    1635.417    -156.250   -1135.417   |          0.000   =         -1.563
       0.000       0.000       0.000       0.000    -677.083     468.750    -572.917     156.250    -104.167    -156.250    1354.167    -468.750   |          0.000   =          0.000
       0.000       0.000       0.000       0.000     468.750    -817.708    -156.250     317.708     156.250   -1135.417    -468.750    1635.417   |          0.000   =         -0.313

-> ktModelSover
       1.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000   |          0.000   =          0.000
       0.000       1.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000   |          0.000   =          0.000
       0.000       0.000       1.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000   |          0.000   =          0.000
       0.000       0.000       0.000       1.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000       0.000   |          0.000   =          0.000
       0.000       0.000       0.000       0.000    2708.333       0.000    -208.333      -0.000    -572.917     156.250    -677.083     468.750   |          0.000   =          0.000
       0.000       0.000       0.000       0.000       0.000    3270.833      -0.000   -2270.833    -156.250     317.708     468.750    -817.708   |          0.000   =         -3.125
       0.000       0.000       0.000       0.000    -208.333      -0.000    2708.333       0.000    -677.083    -468.750    -572.917    -156.250   |          0.000   =          0.000
       0.000       0.000       0.000       0.000      -0.000   -2270.833       0.000    3270.833    -468.750    -817.708     156.250     317.708   |          0.000   =         -0.625
       0.000       0.000       0.000       0.000    -572.917    -156.250    -677.083    -468.750    1354.167     468.750    -104.167     156.250   |          0.000   =          0.000
       0.000       0.000       0.000       0.000     156.250     317.708    -468.750    -817.708     468.750    1635.417    -156.250   -1135.417   |          0.000   =         -1.563
       0.000       0.000       0.000       0.000    -677.083     468.750    -572.917     156.250    -104.167    -156.250    1354.167    -468.750   |          0.000   =          0.000
       0.000       0.000       0.000       0.000     468.750    -817.708    -156.250     317.708     156.250   -1135.417    -468.750    1635.417   |          0.000   =         -0.313

Displacement GD(DOF) is
   0.000E+00
   0.000E+00
   0.000E+00
   0.000E+00
   0.870E-02
  -0.166E-01
  -0.860E-02
  -0.162E-01
   0.105E-01
  -0.423E-01
  -0.102E-01
  -0.420E-01


****************************************
-> call kt2008Results 整理结果


-> call NodeDisplacement
-> call QuadDisplacement/板单元位移
-> Quadstrain/板单元应变
-> Quadstress/板单元应力


-> QuadLogWriteAllInfo 单元信息汇总
QUAD Input data: 输入信息
QUAD
NumofQuad=           2
单元编号
           e      Number     Node(1)     Node(2)     Node(3)     Node(4)    E-Module   Thickness
           1           1           2           4           3           1           1       0.100
                       i        Node           x           y
                       1           2       0.000       0.000
                       2           4       2.500       0.000
                       3           3       2.500       2.000
                       4           1       0.000       2.000
Quad(e)%c
材料矩阵
   31250.000    6250.000       0.000
    6250.000   31250.000       0.000
       0.000       0.000   12500.000


Quad(e)%ke
单元刚度矩阵
    1354.167     468.750    -572.917    -156.250    -677.083    -468.750    -104.167     156.250
     468.750    1635.417     156.250     317.708    -468.750    -817.708    -156.250   -1135.417
    -572.917     156.250    1354.167    -468.750    -104.167    -156.250    -677.083     468.750
    -156.250     317.708    -468.750    1635.417     156.250   -1135.417     468.750    -817.708
    -677.083    -468.750    -104.167     156.250    1354.167     468.750    -572.917    -156.250
    -468.750    -817.708    -156.250   -1135.417     468.750    1635.417     156.250     317.708
    -104.167    -156.250    -677.083     468.750    -572.917     156.250    1354.167    -468.750
     156.250   -1135.417     468.750    -817.708    -156.250     317.708    -468.750    1635.417


Quad(e)%eps
高斯点应变
  -0.198E-02  -0.198E-02   0.202E-02   0.202E-02
  -0.427E-04  -0.159E-03  -0.159E-03  -0.427E-04
  -0.470E-02   0.295E-03   0.201E-03  -0.479E-02


Quad(e)%sig
高斯点B应力
     -62.089     -62.818      62.089      62.818
     -13.699     -17.346       7.636      11.282
     -58.768       3.685       2.518     -59.935


单元编号
           e      Number     Node(1)     Node(2)     Node(3)     Node(4)    E-Module   Thickness
           2           2           4           6           5           3           1       0.100
                       i        Node           x           y
                       1           4       2.500       0.000
                       2           6       5.000       0.000
                       3           5       5.000       2.000
                       4           3       2.500       2.000
Quad(e)%c
材料矩阵
   31250.000    6250.000       0.000
    6250.000   31250.000       0.000
       0.000       0.000   12500.000


Quad(e)%ke
单元刚度矩阵
    1354.167     468.750    -572.917    -156.250    -677.083    -468.750    -104.167     156.250
     468.750    1635.417     156.250     317.708    -468.750    -817.708    -156.250   -1135.417
    -572.917     156.250    1354.167    -468.750    -104.167    -156.250    -677.083     468.750
    -156.250     317.708    -468.750    1635.417     156.250   -1135.417     468.750    -817.708
    -677.083    -468.750    -104.167     156.250    1354.167     468.750    -572.917    -156.250
    -468.750    -817.708    -156.250   -1135.417     468.750    1635.417     156.250     317.708
    -104.167    -156.250    -677.083     468.750    -572.917     156.250    1354.167    -468.750
     156.250   -1135.417     468.750    -817.708    -156.250     317.708    -468.750    1635.417


Quad(e)%eps
高斯点应变
  -0.364E-03  -0.364E-03   0.435E-03   0.435E-03
  -0.192E-03  -0.163E-03  -0.163E-03  -0.192E-03
  -0.126E-02  -0.262E-03  -0.239E-03  -0.124E-02


Quad(e)%sig
高斯点B应力
     -12.580     -12.402      12.580      12.402
      -8.267      -7.376      -2.379      -3.270
     -15.763      -3.272      -2.987     -15.478




****************************************
-> call kt2008End 程序结束


 +--------------------------+
 |        LogClose          |
 +--------------------------+
Used time:=       0.012 S
Program KT2008 ends successfully.



喜欢

0

阅读┊ 收藏 ┊ 喜欢▼ ┊打印┊举报/Report
  

新浪BLOG意见反馈留言板　欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

Copyright © 1996 - 2022 SINA Corporation, All Rights Reserved

新浪公司 版权所有