構造に加速度をかけてみました。
モデル化: salomeにて要素作成
ソルバー :CalculiX(2.19)
ポストプロセッシング:テキストファイル
条件:構造は1m x 1m x 1m (1m3)
密度は1000kg/m3
下向きの加速度は9.8m/s2
結果
固定端反力の結果です。
forces (fx,fy,fz) for set FEMCONSTRAINTDISPLACEMENT and time 0.1000000E+01 1 -6.837276E+02 -6.837276E+02 3.163376E+03 3 -8.191388E+02 8.191388E+02 2.226624E+03 5 8.191388E+02 -8.191388E+02 2.226624E+03 7 6.837276E+02 6.837276E+02 3.163376E+03 total force (fx,fy,fz) for set FEMCONSTRAINTDISPLACEMENT and time 0.1000000E+01 0.000000E+00 2.160050E-12 1.078000E+04
Z軸方向の反力が1.078000E+04です。9800Nのはずですが、1割ほど多くなっています。
1m3 × 1000kg/m3 × 9.8m/s2 = 9800N
要素を細かくした場合
forces (fx,fy,fz) for set FEMCONSTRAINTDISPLACEMENT and time 0.1000000E+01 1 -1.099273E+03 -1.094399E+03 2.460114E+03 3 -1.097291E+03 1.094073E+03 2.460303E+03 5 1.098947E+03 -1.095763E+03 2.460303E+03 7 1.097617E+03 1.096089E+03 2.460114E+03 total force (fx,fy,fz) for set FEMCONSTRAINTDISPLACEMENT and time 0.1000000E+01 1.364242E-12 8.640200E-12 9.840833E+03
要素を細かくした場合(下図)だと値が近くなりました。(9.840833E+03)
途中。。。
作業
形状と要素はsalomeにて作成しています。
参考になる手順は以下にあります。
T218 剛体 :作業
サンプルファイル
*node, nset=Nall
1, 1.000000e+000, 1.000000e+000, 0.000000e+000
2, 1.000000e+000, 1.000000e+000, 1.000000e+000
3, 1.000000e+000, 0.000000e+000, 0.000000e+000
4, 1.000000e+000, 0.000000e+000, 1.000000e+000
5, 0.000000e+000, 1.000000e+000, 0.000000e+000
6, 0.000000e+000, 1.000000e+000, 1.000000e+000
7, 0.000000e+000, 0.000000e+000, 0.000000e+000
8, 0.000000e+000, 0.000000e+000, 1.000000e+000
9, 0.000000e+000, 5.000000e-001, 1.000000e+000
10, 5.000000e-001, 1.000000e+000, 1.000000e+000
11, 1.000000e+000, 5.000000e-001, 1.000000e+000
12, 5.000000e-001, 0.000000e+000, 1.000000e+000
13, 0.000000e+000, 5.000000e-001, 0.000000e+000
14, 5.000000e-001, 1.000000e+000, 0.000000e+000
15, 1.000000e+000, 5.000000e-001, 0.000000e+000
16, 5.000000e-001, 0.000000e+000, 0.000000e+000
17, 1.000000e+000, 1.000000e+000, 5.000000e-001
18, 0.000000e+000, 1.000000e+000, 5.000000e-001
19, 1.000000e+000, 0.000000e+000, 5.000000e-001
20, 0.000000e+000, 0.000000e+000, 5.000000e-001
21, 5.000000e-001, 5.000000e-001, 1.000000e+000
22, 5.000000e-001, 5.000000e-001, 0.000000e+000
23, 5.000000e-001, 1.000000e+000, 5.000000e-001
24, 5.000000e-001, 0.000000e+000, 5.000000e-001
25, 1.000000e+000, 5.000000e-001, 5.000000e-001
26, 0.000000e+000, 5.000000e-001, 5.000000e-001
*nset, nset=Group_1
1, 3, 5, 7,
*element, elset=C3D10,type=C3D10
25, 7, 1, 5, 6, 22, 14, 13, 26, 23, 18,
26, 2, 1, 4, 6, 17, 25, 11, 10, 23, 21,
27, 7, 6, 8, 4, 26, 9, 20, 24, 21, 12,
28, 4, 1, 3, 7, 25, 15, 19, 24, 22, 16,
29, 7, 1, 6, 4, 22, 23, 26, 24, 25, 21,
***********************************************************
** Node sets for prescribed displacement constraint
** written by write_node_sets_constraints_displacement function
** FemConstraintDisplacement
*NSET,NSET=FemConstraintDisplacement
1,
3,
5,
7,
***********************************************************
** SolidMaterial
*MATERIAL, NAME=SolidMaterial
*ELASTIC
10000, 0.300
*DENSITY
1000
***********************************************************
**ソリッド(*解除のこと)
*SOLID SECTION,ELSET=C3D10,MATERIAL=SolidMaterial
*STEP
*STATIC
***********************************************************
** Displacement constraint applied
** written by write_constraints_displacement function
** FemConstraintDisplacement
*BOUNDARY
FemConstraintDisplacement,1
FemConstraintDisplacement,2
FemConstraintDisplacement,3
***********************************************************
** Self weight Constraint
** written by write_constraints_selfweight function
** ConstraintSelfWeight
*DLOAD
C3D10,GRAV,9.8,0.0,0.0,-1.0
***********************************************************
** Outputs --> frd file
** written by write_outputs_types function
*NODE FILE, OUTPUT=2d
U,RF
*EL FILE,SECTION FORCES
S, E
**出力指定 RFは反力 セットに複数節点がある場合にはTOTALSで合計
*NODE PRINT, NSET=FemConstraintDisplacement
RF
*NODE PRINT, NSET=FemConstraintDisplacement, TOTALS=ONLY
RF
***********************************************************
** written by write_step_end function
*END STEP
コメント