yongjia 发表于 2015-9-5 13:59:01

用matlab编写有限元，发现的一个刚度问题，希望大家指教

最近学习力学，用matlab编辑有限元：有一个问题想请教大家，结构的刚度与体积是什么关系。因为我的最终目的是把有限元与微观的位错结合起来，因此：我编辑的有限元模型很小，为3um*5um,划分网格以后呢就更小了，那么我就有一个问题：
如果结构如此，其单元刚度为（6,6）矩阵：
dk =

   61.9279   19.2317-47.7962   -5.3601-14.1317-13.8716
   19.2317   39.5138    5.3587   -3.2798-24.5903-36.2340
-47.7962    5.3587   61.9257-19.2299-14.1295   13.8712
   -5.3601   -3.2798-19.2299   39.5104   24.5900-36.2306
-14.1317-24.5903-14.1295   24.5900   28.2612    0.0004
-13.8716-36.2340   13.8712-36.2306    0.0004   72.4646

而如果把结构为3m*5m，单元矩阵为：

dk =

1.0e+007 *

    6.1305    1.9221   -4.7027   -0.5359   -1.4278   -1.3862
    1.9221    3.9455    0.5360   -0.2787   -2.4581   -3.6668
   -4.7027    0.5360    6.1365   -1.9240   -1.4338    1.3881
   -0.5359   -0.2787   -1.9240    3.9494    2.4599   -3.6707
   -1.4278   -2.4581   -1.4338    2.4599    2.8616   -0.0018
   -1.3862   -3.6668    1.3881   -3.6707   -0.0018    7.3375

请问这种不同是正常的嘛？真的是不理解，请问有什么能正确计算微尺度单元刚度的公式吗

yongjia 发表于 2015-9-5 16:04:24

已经解决
计算没有问题，体积小的确是刚度大大减小，是我在外载荷加载时还是采用了一样的载荷，但是这个相同的载荷加载到米与微米之间的结果却大大不同，因此需要对外载荷进行进一步的调整

gogouter 发表于 2015-9-10 20:36:19

用Matlab写有限元，少了很多乐趣！

w7945121 发表于 2018-6-21 18:32:36

楼主棒棒哒:victory:

xuchen_fendou 发表于 2018-6-24 09:29:19

用python写会好一些

yuki060 发表于 2018-7-12 12:06:43

http://www.51zxw.net/study.asp?vip=16419520

雷欧小霆 发表于 2018-8-9 08:20:14

gogouter 发表于 2015-9-10 20:36
用Matlab写有限元，少了很多乐趣！

用Matlab写有限元，乐趣为何少了？:P
页: [1]
查看完整版本: 用matlab编写有限元，发现的一个刚度问题，希望大家指教