In this paper is presented a two-dimension analytical model for linear permanent magnet inner armature double sided synchronous motors (LPMIADSSMs). The flux density in all areas of the proposed machine is calculated based on the sub-domain method. According to this method the machine areas is divided into eleven sub-regions as first external, first mover, first PM, first air-gap, first winding, stator, second winding, second air gap, second PM, second external and second mover which sign as FE, FR, FP, FAG, FW, S, SW, SAG, SP, SE, SR. To find the flux density equations, it is mandatory to solve the extracted Maxwell equations and apply the boundary conditions between each two sub-regions, which lead to find the unknown coefficients for flux density in each sub-region. In addition, the influences of magnetization patterns, i.e. parallel, ideal Hal Bach, 2-segment Halbach and bar magnet in the shifting direction of the flux distribution is investigated. To validate the obtained results, the proposed model results are compared with those obtained from finite element method (FEM).