Is there a general formula which I can use to calculate this and if it's with proof or reasoning would be great as well. Even if you could please solve this $4$-variable ordered set of positive integers case would be a great help indeed. Sorry I am weak in this topic.
Thanks Robert for your pointers in modifying it to $(4\times1) + (3\times2) + (2\times3) + (1\times4) = 41$, where there are no additional constraints on $x_i$
To solve this must do it manually by cases or is there a short way of calculating it using "Bars and Stars" (of which I am aware of the method)? The answer given is $321$.
But when I calculate it like this it doesn't seem to match the answer: imagine the new given equation as containing $(4+3+2+1) = 10$$x_i$'s individually so as equivalent to finding the number of positive integer solutions for $b_1 + b_2 + ... + b_{10} = 41$ which is $C(41 - 1, 10 - 1) = C(40, 9)$. Now since we have $4x_1$'s that are identical, we need to divide this by $4!$ (for its permutations of "overcounting") and similarly divide next by $3!$ (for $3x_2$'s that are equal), and $2!$ for the $2x_3$'s that are equal. How is my approach wrong here? anybody, please help me? Thanks again
As per Robert's transformation with Anil's suggestion of the generating functions method, here's my work on it:
I actually used an online calculator called "Sage Math" to do this via this website:
https://www.whitman.edu/mathematics/cgt_online/book/section03.01.html
With these modifications
$f(x)=(1-x^4)^{-1}\times(1-x^3)^{-1}\times(1-x^2)^{-1}\times(1-x)^{-1}$
show(taylor($f,x,0,33$))
And have gotten the verification that $321$ is indeed the answer after some tedious algebra (please note that since we are looking for the coefficient of $x^41$ in the expansion Anil wrote which has $x^(4+3+2+1) = x^10$ factored outside meaning in the remaining Newton's Binomial Expansions with negative exponent, we only need to determine the coefficient of $x^(41-10)= x^31$ in it which is "$321$" as the answer suggested (please see image below):

If we solve the problem without the condition $a<b<c<d$ then the answer is $40\choose3$.
Now we want to count the number of sets with at least $2$ equal variables.
Note that all $4$ of them can't be equal, because if they are then their sum should be divisible by $4$.
Case $1$ : Exactly $2$ of them are equal.
We need to solve the equation $$2x_1+x_2+x_3=41 : x_1,x_2,x_3\in \mathbb{N}$$
where $x_1,x_2$ and $x_3$ are pairwise distinct.
If $x_1 = k$ then the number solutions is ${41-2k-1\choose1} = 41-2k-1$, for any $1\leq k\leq13$ we have counted $2$ extra solutions ($x_2 = k$ and $x_3=k$).
Thus the number of solutions for this case is ${4\choose2}(38+36+\dots + 2 - 26) = {4\choose2}2(19+18+\dots + 1 - 13) = 2124$
Case 2 : Exactly $3$ of them are equal.
We calculate the number of solutions like Case $1$.
The number of solutions for this case is $13{4\choose3} = 52$.
So there are ${40\choose3} - (52 + 2124) = 7704$ solutions.
For every pairwise distinct $a,b,c,d$ there is exactly 1 permutation such that $a<b<c<d$. we counted each set $(a,b,c,d)$, $4!$ times.
So the answer is $\frac{7704}{4!}=321$.