已经尝试过使用延迟函数定义加?NumericQ的方法，但是没有效果。
d = 3; Subscript[v, q] = 0.3;(*the possion ratio of matrix*) \
Subscript[v, p] = 0.2;(*the possion ratio of obstacle*)
Subscript[E, p] = 80000000000.;(*the ealstic modulus of \
obstacle*)Subscript[E, q] = 2000000000.;(*the ealstic modulus of \
matrix*)
Subscript[K, q] = Subscript[E, q]/(
2*(1 + Subscript[v, q])); Subscript[G, q] = Subscript[E, q]/(
d*(1 + Subscript[v, q]*(1 - d)));
Subscript[K, p] = Subscript[E, p]/(
2*(1 + Subscript[v, p])); Subscript[G, p] = Subscript[E, p]/(
d*(1 + Subscript[v, p]*(1 - d)));
Subscript[\[Kappa], pq] = (Subscript[K, p] - Subscript[K, q])/(
Subscript[K,
p] + (2*(d - 1))/d*Subscript[G, q]); Subscript[\[Mu], pq] = (
Subscript[G, p] - Subscript[G, q])/(
Subscript[G, p] + (
Subscript[G,
q]*(d*Subscript[K, q]/2 + (d + 1)*(d - 2)*Subscript[G, q]/d))/(
Subscript[K, q] + 2*Subscript[G, q]));
Subscript[\[Delta], ij] = 1;
Subscript[\[Delta], kl] = 1;
Subscript[\[Delta], ik] = 1;
Subscript[\[Delta], jl] = 1;
Subscript[\[Delta], il] = 1;
Subscript[\[Delta], jk] = 1;
Subscript[\[CapitalLambda], h] =
1/d *Subscript[\[Delta], ij]*Subscript[\[Delta], kl];
Subscript[\[CapitalLambda], s] =
1/2 *(Subscript[\[Delta], ik]*Subscript[\[Delta], jl] +
Subscript[\[Delta], il]*Subscript[\[Delta], jk]) -
1/d *Subscript[\[Delta], ij]*Subscript[\[Delta],
kl];(*Subscript[\[Delta], ij] is the Kronecker delta symbol and \
Subscript[\[CapitalLambda], s] is a fourth-order tensor*)
Subscript[C, q] =
d*Subscript[K, q]*Subscript[\[CapitalLambda], h] +
2*Subscript[G, q]*Subscript[\[CapitalLambda], s];
Subscript[I,
1] = {{{{1, 0, 0}, {0, 1, 0}, {0, 0, 1}}, {{0, 0, 0}, {0, 0, 0}, {0,
0, 0}}, {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}}, {{{0, 0, 0}, {0, 0,
0}, {0, 0, 0}}, {{1, 0, 0}, {0, 1, 0}, {0, 0, 1}}, {{0, 0,
0}, {0, 0, 0}, {0, 0, 0}}}, {{{0, 0, 0}, {0, 0, 0}, {0, 0,
0}}, {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}}, {{1, 0, 0}, {0, 1,
0}, {0, 0, 1}}}}；
(*Subscript[I, 1] is the identity matrix *)
Dq = Subscript[\[CapitalLambda], h]/(
d*Subscript[K, q] + 2*(d - 1)*Subscript[G, q]) + (
d*(Subscript[K, q] + 2*Subscript[G, q])*
Subscript[\[CapitalLambda], s])/(
Subscript[G,
q]*(d + 2)*(d*Subscript[K, q] + 2*(d - 1)*Subscript[G, q]));
Subscript[L, eq] = (Subscript[C, e] - Subscript[C, q])/(
Subscript[I, 1] + Dq*(Subscript[C, e] - Subscript[C, q]));
Celeft1111 = (d*Subscript[K, q] +
2*(d - 1)*Subscript[G,
q])*(Subscript[\[Kappa], pq]* Subscript[\[CapitalLambda],
h] + ((d + 2)*Subscript[G, q]*Subscript[\[Mu], pq]*
Subscript[\[CapitalLambda], s])/(
d*(Subscript[K, q] + 2*Subscript[G, q])))/Subscript[L, eq];
Celeft1111[[1, 1, 1, 1]];
Subscript[\[Phi], p] = 0.4;
Ceright1111 = Subscript[I, 1]/Subscript[\[Phi], p];
Ceright1111[[1, 1, 1, 1]];
a = 0.19; b = 1.473;
Subscript[S, 2][x_] :=
Subscript[\[Phi], p]^2 + Subscript[\[Phi], p]^2*Exp[(-a)*x^b];
\[CapitalOmega] = N[4 \[Pi]];
xx1 = r*Sin[\[Theta]]*Cos[\[CurlyPhi]]; xx2 =
r*Sin[\[Theta]]*Sin[\[CurlyPhi]]; xx3 = r*Cos[\[Theta]];
n = {xx1/r, xx2/r, xx3/r};
Subscript[\[Alpha], q] = d*Subscript[K, q]/Subscript[G, q] + (d - 2);
Hqmmklr =
d/(\[CapitalOmega]*(d*Subscript[K, q] +
2*(d - 1)*Subscript[G, q]))*1/
r^d*(d*n[[1]]*n[[1]] - Subscript[\[Delta], kl]);
Hqijklr =
1/(2*\[CapitalOmega]*(d*Subscript[K, q] +
2*(d - 1)*Subscript[G, q]))*1/r^
d*(Subscript[\[Alpha], q]*Subscript[\[Delta], kl]*
Subscript[\[Delta], ij] -
d*(Subscript[\[Delta], ik]*Subscript[\[Delta], jl] +
Subscript[\[Delta], il]*Subscript[\[Delta], jk]) -
d*Subscript[\[Alpha],
q]*(Subscript[\[Delta], ij]*n[[1]]*n[[1]] +
Subscript[\[Delta], kl]*n[[1]]*n[[1]]) + (
d*(d - Subscript[\[Alpha], q]))/
2*(Subscript[\[Delta], ik]*n[[1]]*n[[1]] +
Subscript[\[Delta], il]*n[[1]]*n[[1]] +
Subscript[\[Delta], jk]*n[[1]]*n[[1]] +
Subscript[\[Delta], jl]*n[[1]]*n[[1]]) +
d*(d + 2)*Subscript[\[Alpha], q]*n[[1]]*n[[1]]*n[[1]]*n[[1]]);
Uqr = (d*Subscript[K, q] +
2*(d - 1)*Subscript[G,
q])*((Subscript[\[Kappa],
pq] - ((d + 2)*Subscript[G, q]* Subscript[\[Mu], pq])/(
d*(Subscript[K, q] + 2*Subscript[G, q])))*Subscript[\[Delta],
ij]/d *Hqmmklr + ((d + 2)*Subscript[G, q]*Subscript[\[Mu],
pq])/(d*(Subscript[K, q] + 2*Subscript[G, q]))*Hqijklr);
ceright21111 = (
Subscript[S, 2][r] - (Subscript[\[Phi], p]^2))/(Subscript[\[Phi],
p]^2)*Uqr*r*r*Sin[\[Theta]];
Ceright21111 = \!\(
\*SubsuperscriptBox[\(\[Integral]\), \(0\), \(15\)]\(
\*SubsuperscriptBox[\(\[Integral]\), \(0\), \(2 \[Pi]\)]\(
\*SubsuperscriptBox[\(\[Integral]\), \(0\), \(\[Pi]\)]ceright21111 \
\[DifferentialD]\[CurlyPhi] \[DifferentialD]\[Theta] \
\[DifferentialD]r\)\)\);
Subscript[S,
3] = (r/(r + R)*Subscript[S, 2][R] + R/(r + R)*Subscript[S, 2][r])*
Subscript[S, 2][Abs[lrRr]]/Subscript[\[Phi], p];
xX1 = R*Sin[\[CapitalTheta]]*Cos[\[CapitalPsi]]; xX2 =
R*Sin[\[CapitalTheta]]*Sin[\[CapitalPsi]]; xX3 =
R*Cos[\[CapitalTheta]];
lrRr = \[Sqrt](Abs[(r^2 + R^2 -
2*R*r*(Sin[\[CapitalTheta]] Cos[\[CapitalPsi]] Sin[\[Theta]] \
Cos[\[CurlyPhi]] +
Sin[\[CapitalTheta]] Sin[\[CapitalPsi]] Sin[\[Theta]] Sin[\
\[CurlyPhi]] + Cos[\[CapitalTheta]] Cos[\[Theta]]))]);
Subscript[n,
1] = {(xX1 - xx1)/lrRr, (xX2 - xx2)/lrRr, (xX3 - xx3)/lrRr};
HqmmkllrRr =
d/(\[CapitalOmega]*(d*Subscript[K, q] +
2*(d - 1)*Subscript[G, q]))*1/
lrRr^d*(d*Subscript[n, 1][[1]]*Subscript[n, 1][[1]] -
Subscript[\[Delta], kl]);
HqijkllrRr =
1/(2*\[CapitalOmega]*(d*Subscript[K, q] +
2*(d - 1)*Subscript[G, q]))*1/lrRr^
d*(Subscript[\[Alpha], q]*Subscript[\[Delta], kl]*
Subscript[\[Delta], ij] -
d*(Subscript[\[Delta], ik]*Subscript[\[Delta], jl] +
Subscript[\[Delta], il]*Subscript[\[Delta], jk]) -
d*Subscript[\[Alpha],
q]*(Subscript[\[Delta], ij]*Subscript[n, 1][[1]]*
Subscript[n, 1][[1]] +
Subscript[\[Delta], kl]*Subscript[n, 1][[1]]*
Subscript[n, 1][[1]]) + (d*(d - Subscript[\[Alpha], q]))/
2*(Subscript[\[Delta], ik]*Subscript[n, 1][[1]]*
Subscript[n, 1][[1]] +
Subscript[\[Delta], il]*Subscript[n, 1][[1]]*
Subscript[n, 1][[1]] +
Subscript[\[Delta], jk]*Subscript[n, 1][[1]]*
Subscript[n, 1][[1]] +
Subscript[\[Delta], jl]*Subscript[n, 1][[1]]*
Subscript[n, 1][[1]]) +
d*(d + 2)*Subscript[\[Alpha], q]*Subscript[n, 1][[1]]*
Subscript[n, 1][[1]]*Subscript[n, 1][[1]]*Subscript[n, 1][[1]]);
UqlrRr = (d*Subscript[K, q] +
2*(d - 1)*Subscript[G,
q])*((Subscript[\[Kappa],
pq] - ((d + 2)*Subscript[G, q]* Subscript[\[Mu], pq])/(
d*(Subscript[K, q] + 2*Subscript[G, q])))*Subscript[\[Delta],
ij]/d *HqmmkllrRr + ((d + 2)*Subscript[G, q]*Subscript[\[Mu],
pq])/(d*(Subscript[K, q] + 2*Subscript[G, q]))*HqijkllrRr);
Uqq = Uqr.UqlrRr;
Ceright31111[r_?NumericQ, \[Theta]_?NumericQ, \[CurlyPhi]_?NumericQ,
R_?NumericQ, \[CapitalTheta]_?NumericQ, \[CapitalPsi]_?NumericQ] =
-NIntegrate[(Subscript[S, 3]/Subscript[\[Phi], p]^2 - (
Subscript[S, 2][r]*Subscript[S, 2][Abs[lrRr]])/
Subscript[\[Phi], p]^3)*Uqq*r^2*Sin[\[Theta]]*R^2*
Sin[\[CapitalTheta]], {r, 1, 15}, {\[Theta], 0,
2 \[Pi]}, {\[CurlyPhi], 0, \[Pi]}, {R, 1, 15}, {\[CapitalTheta],
0, 2 \[Pi]}, {\[CapitalPsi], 0, \[Pi]}];
Solve[Celeft1111[[1, 1, 1, 1]] ==
Ceright1111[[1, 1, 1, 1]] - Ceright21111 +
Ceright31111, Subscript[C, e]]

首先，代码前面的部分中有中文分号。
然后，你说你用了延迟赋值，但是你代码里显然写的是=，是立即赋值。
再就是，Ceright31111到底是什么？你定义其为函数，又在Solve里作变量求解。
我可能懒得继续给你看代码，但我留个建议：做一件事情要分几个步骤来做——尤其是你不是很懂如何做的时候。做一部分检查一部分，这样才知道是那一部分出问题，也不至于搞得自己思维混乱。

《为什么y=y[x]会出错？你分得清函数关系和函数吗？》：
https://tieba.baidu.com/p/5936400015