Large fuzzy parameters and interval parameters are considered simultaneously.
Subinterval perturbation strategy is combined with the finite element framework.
Neumann series with high-order terms is adopted for subinterval matrix inverse.
Feasibility of the proposed methods is not affected by the uncertainty level.