Thanks Jenkins,
for [stiffness matrix][d]=[F],like the picture example, my calculation is:
1. each beam matrix has 12 [d] and 12 [F],
2. combine the matrix to form a global matrix, result in 48 [d] and 48 [F],
3. solve the 48 [d] by the known 48 [F],
"You use the [d] to calculate the [F} at...