采用中心差分方法计算过渡自由度(
32]:
(1)米米l:mtext>
u米米l:mi>
̇米米l:mo>
我米米l:mi>
t米米l:mi>
+米米l:mo>
Δ米米l:mi>
t米米l:mi>
/米米l:mo>
2米米l:mn>
=米米l:mo>
u米米l:mi>
̇米米l:mo>
我米米l:mi>
t米米l:mi>
−米米l:mo>
Δ米米l:mi>
t米米l:mi>
/米米l:mo>
2米米l:mn>
+米米l:mo>
∑米米l:mo>
F米米l:mi>
我米米l:mi>
t米米l:mi>
Δ米米l:mi>
t米米l:mi>
米米米l:mi>
,米米l:mo>
(2)米米l:mtext>
u米米l:mi>
我米米l:mi>
t米米l:mi>
+米米l:mo>
Δ米米l:mi>
t米米l:mi>
=米米l:mo>
u米米l:mi>
我米米l:mi>
t米米l:mi>
+米米l:mo>
u米米l:mi>
̇米米l:mo>
我米米l:mi>
t米米l:mi>
+米米l:mo>
Δ米米l:mi>
t米米l:mi>
/米米l:mo>
2米米l:mn>
Δ米米l:mi>
t米米l:mi>
,米米l:mo>
在哪里<我nline-formula>
u米米l:mi>
̇米米l:mo>
我米米l:mi>
t米米l:mi>
和<我nline-formula>
u米米l:mi>
我米米l:mi>
t米米l:mi>
代表组件的速度和位置<我nline-formula>
我米米l:mi>
(<我nline-formula>
我米米l:mi>
=米米l:mo>
1米米l:mn>
,米米l:mo>
3米米l:mn>
)时间<我nline-formula>
t米米l:mi>
。<我nline-formula>
∑米米l:mo>
F米米l:mi>
我米米l:mi>
t米米l:mi>
代表组件的力的总和<我nline-formula>
我米米l:mi>
作用于一个节点的质量<我nline-formula>
米米米l:mi>
,时间步长<我nline-formula>
Δ米米l:mi>
t米米l:mi>
。
gydF4y2Ba组件的角速度<我nline-formula>
我米米l:mi>
在时间<我nline-formula>
t米米l:mi>
计算使用
(3)米米l:mtext>
ω米米l:mi>
我米米l:mi>
t米米l:mi>
+米米l:mo>
Δ米米l:mi>
t米米l:mi>
/米米l:mo>
2米米l:mn>
=米米l:mo>
ω米米l:mi>
我米米l:mi>
t米米l:mi>
−米米l:mo>
Δ米米l:mi>
t米米l:mi>
/米米l:mo>
2米米l:mn>
+米米l:mo>
∑米米l:mo>
米米米l:mi>
我米米l:mi>
t米米l:mi>
我米米l:mi>
Δ米米l:mi>
t米米l:mi>
,米米l:mo>
在哪里<我nline-formula>
∑米米l:mo>
米米米l:mi>
我米米l:mi>
t米米l:mi>
是所有的总和moment-components作用于惯性矩的节点<我nline-formula>
我米米l:mi>
。
gydF4y2Ba弹簧力变化是由节点的位移(
32]:
(4)米米l:mtext>
F米米l:mi>
N米米l:mi>
←米米l:mo>
F米米l:mi>
N米米l:mi>
+米米l:mo>
u米米l:mtext>
̇米米l:mo>
N米米l:mi>
k米米l:mi>
N米米l:mi>
Δ米米l:mi>
t米米l:mi>
,米米l:mo>
(5)米米l:mtext>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
←米米l:mo>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
+米米l:mo>
u米米l:mi>
̇米米l:mo>
我米米l:mi>
年代米米l:mi>
k米米l:mi>
年代米米l:mi>
Δ米米l:mi>
t米米l:mi>
,米米l:mo>
在哪里<我nline-formula>
N米米l:mi>
代表“正常”,<我nline-formula>
年代米米l:mi>
代表“剪切”,<我nline-formula>
F米米l:mi>
代表了春天的力量。<我nline-formula>
k米米l:mi>
N米米l:mi>
和<我nline-formula>
k米米l:mi>
年代米米l:mi>
分别代表了春天正常和剪切刚度。如果力超过弹簧强度,春天优惠和微裂纹的形成。
gydF4y2Ba关节滑开放遵循的关系(
37]
(6)米米l:mtext>
如果米米l:mtext>
F米米l:mi>
n米米l:mi>
−米米l:mo>
p米米l:mi>
一个米米l:mi>
<米米l:mo>
0米米l:mn>
然后米米l:mtext>
F米米l:mi>
n米米l:mi>
=米米l:mo>
0米米l:mn>
,米米l:mo>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
=米米l:mo>
0米米l:mn>
,米米l:mo>
其他的米米l:mtext>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
⇐米米l:mo>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
∣米米l:mo>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
∣米米l:mo>
最小值米米l:mi>
F米米l:mi>
n米米l:mi>
−米米l:mo>
p米米l:mi>
一个米米l:mi>
棕褐色米米l:mi>
ϕ米米l:mi>
,米米l:mo>
∣米米l:mo>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
∣米米l:mo>
,米米l:mo>
在哪里<我nline-formula>
F米米l:mi>
n米米l:mi>
代表了法向力,<我nline-formula>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
代表了剪切力向量,<我nline-formula>
p米米l:mi>
代表了压力,<我nline-formula>
一个米米l:mi>
代表了表面面积<我nline-formula>
ϕ米米l:mi>
代表了摩擦角。
gydF4y2Ba下面的关系决定了保税联合状态:如果<我nline-formula>
F米米l:mi>
n米米l:mi>
−米米l:mo>
p米米l:mi>
一个米米l:mi>
+米米l:mo>
σ米米l:mi>
c米米l:mi>
一个米米l:mi>
<米米l:mo>
0米米l:mn>
或<我nline-formula>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
>米米l:mo>
τ米米l:mi>
c米米l:mi>
一个米米l:mi>
,债券在紧张或剪切(失败<我nline-formula>
σ米米l:mi>
c米米l:mi>
债券的抗拉强度,<我nline-formula>
τ米米l:mi>
c米米l:mi>
其他债券剪切强度),,<我nline-formula>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
←米米l:mo>
F米米l:mi>
我米米l:mi>
年代米米l:mi>
,债券仍然完好无损。
年代ec>
2.2。流公式
如图
5,在高频模拟使用流体元素由管道连接起来。流体元素作为微裂隙位于中心的破碎的弹簧或弹簧重叠的关节。