如何快速消除下式中的KroneckerDelta，目前我只会半自动求解，利用Block对每一项分别令KroneckerDelta下的角标相等进行化简，但是这样觉得很是愚蠢，但是查看了Documentation之后并没有找到快速消除Kroneckerdelta的方法。
希望得到指点。
待求式子如下：
1/2 E^(Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] +
Subscript[\[Lambda], q1] - Subscript[\[Lambda], q2])
Cosh[Subscript[\[Lambda], k]] KroneckerDelta[k,
k - p1 - p2 - q1 - q2] KroneckerDelta[p1, -q2] KroneckerDelta[
p2, -q1] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] +
1/2 E^(Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] +
Subscript[\[Lambda], q1] - Subscript[\[Lambda], q2])
Cosh[Subscript[\[Lambda], k]] KroneckerDelta[k,
k - p1 - p2 - q1 - q2] KroneckerDelta[p1, -q1] KroneckerDelta[
p2, -q2] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] -
1/2 E^(Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] +
Subscript[\[Lambda], q1] - Subscript[\[Lambda], q2])
Cosh[Subscript[\[Lambda], k]] KroneckerDelta[k,
k - p1 - p2 - q1 - q2] KroneckerDelta[p1, -p2] KroneckerDelta[
q1, -q2] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] -
1/2 E^(-Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] -
Subscript[\[Lambda], q1] +
2 (Subscript[\[Lambda], p1] + Subscript[\[Lambda], q1]) -
Subscript[\[Lambda], q2])
KroneckerDelta[-k, -q2] KroneckerDelta[p1, -q1] KroneckerDelta[-p2,
k - p1 - p2 - q1 - q2] Sinh[Subscript[\[Lambda], k]] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] -
1/2 E^(-Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] -
Subscript[\[Lambda], q1] +
2 (Subscript[\[Lambda], p1] + Subscript[\[Lambda], q1]) -
Subscript[\[Lambda], q2])
KroneckerDelta[-k, -q1] KroneckerDelta[p1, -q2] KroneckerDelta[-p2,
k - p1 - p2 - q1 - q2] Sinh[Subscript[\[Lambda], k]] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] +
1/2 E^(-Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] -
Subscript[\[Lambda], q1] +
2 (Subscript[\[Lambda], p1] + Subscript[\[Lambda], q1]) -
Subscript[\[Lambda], q2])
KroneckerDelta[-k, -q2] KroneckerDelta[-p1,
k - p1 - p2 - q1 - q2] KroneckerDelta[p2, -q1] Sinh[
Subscript[\[Lambda], k]] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2]

大家都去过年了吗？

帮助搜下 PiecewiseExpand

我利用PiecewiseExpand[]后，他可以得到我想要的结果的50%，还有50%就是我想利用蓝框里的关系代入到前式中(红圈)里去,比如利用p1+p2==0,将所有p2替换为-p1. 该如何实现呢？
\[Piecewise]1/2 E^(Subscript[\[Lambda], p1]-Subscript[\[Lambda], p2]+Subscript[\[Lambda], q1]-Subscript[\[Lambda], q2]) Cosh[Subscript[\[Lambda], k]] Sinh[Subscript[\[Lambda], -k+p1+p2+q1+q2]] Subscript[h, p1,p1+p2,p2] Subscript[h, q1,q1+q2,q2]p1+p2+q1+q2==0&&p1+q2==0&&p2+q1==0
0True
,\[Piecewise]1/2 E^(Subscript[\[Lambda], p1]-Subscript[\[Lambda], p2]+Subscript[\[Lambda], q1]-Subscript[\[Lambda], q2]) Cosh[Subscript[\[Lambda], k]] Sinh[Subscript[\[Lambda], -k+p1+p2+q1+q2]] Subscript[h, p1,p1+p2,p2] Subscript[h, q1,q1+q2,q2]p1+p2+q1+q2==0&&p1+q1==0&&p2+q2==0
0True
,