Constructing spline models for isogeometric analysis is important in integrating design and analysis. Converting designed CAD (Computer Aided Design) models with B-reps to analysis-suitable volumetric T-spline is fundamental for the integration. In this thesis, we work on two directions to achieve this: (a) using Boolean operations and skeletons to build polycubes for feature-preserving high-genus volumetric T-spline construction; and (b) developing weighted T-splines with arbitrary degree for T-spline surface and volume modeling which can be used for analysis. In this thesis, we first develop novel algorithms to build feature-preserving polycubes for volumetric T-spline construction. Then a new type of T-spline named the weighted T-spline with arbitrary degree is defined. It is further used in converting CAD models to analysis-suitable volumetric T-splines. An algorithm is first developed to use Boolean operations in CSG (Constructive Solid Geometry) to generate polycubes robustly, then the polycubes are used to generate volumetric rational solid T-splines. By solving a harmonic field with proper boundary conditions, the input surface is automatically decomposed into regions that are classified into topologically either a cube or a torus. Two Boolean operations, union and difference, are performed with the primitives and polycubes are generated by parametric mapping. With polycubes, octree subdivision is carried out to obtain a volumetric T-mesh. The obtained T-spline surface is C2-continuous everywhere except the local region surrounding irregular nodes, where the surface continuity is elevated from C0 to G1. B´ezier elements are extracted from the constructed solid T-spline models, which are further used in isogeometric analysis. The Boolean operations preserve the topology of the models inherited from design and can generate volumetric T-spline models with better quality. Furthermore, another algorithm is developed which uses skeleton as a guidance to the polycube construction. From the skeleton of the input model, initial cubes in the interior are first constructed. By projecting corners of interior cubes onto the surface and generating a new layer of boundary cubes, the entire interior domain is split into different cubic regions. With the splitting result, octree subdivision is performed to obtain T-spline control mesh or T-mesh. Surface features are classified into three groups: open curves, closed curves and singularity features. For features without introducing new singularities like open or closed curves, we preserve them by aligning to the parametric lines during subdivision, performing volumetric parameterization from frame field, or modifying the skeleton. For features introducing new singularities, we design templates to handle them. With a valid T-mesh, we calculate rational trivariate T-splines and extract B´ezier elements for isogeometric analysis. Weighted T-spline basis functions are designed to satisfy partition of unity and linear independence. The weighted T-spline is proved to be analysis-suitable. Compared to standard T-splines, weighted T-splines have less geometrical constraint and can decrease the number of control points significantly. Trimmed NURBS surfaces of CAD models are reparameterized with weighted T-splines by a new edge interval extension algorithm, with bounded surface error introduced. With knot interval duplication, weighted T-splines are used to deal with extraordinary nodes. With B´ezier coefficient optimization, the surface continuity is elevated from C0 to G1 for the one-ring neighborhood elements. Parametric mapping and sweeping methods are developed to construct volumetric weighted T-splines for isogeometric analysis. Finally, we develop an algorithm to construct arbitrary degree T-splines. The difference between odd degree and even degree T-splines are studied in detail. The methods to extract knot intervals, calculate new weights to handle extraordinary nodes, and extract B´ezier elements for analysis are investigated with arbitrary degrees. Hybrid degree weighted Tspline is generated at designated region with basis functions of different degrees, for the purpose of performing local p-refinement. We also study the convergence rate for T-spline models of different degrees, showing that hybrid degree weighted T-splines have better performance after p-refinement. In summary, we develop novel methods to construct volumetric T-splines based on polycube and sweeping methods. Arbitrary degree weighted T-spline is proposed, with proved analysis-suitable properties. Weighted T-spline basis functions are used to reparameterize trimmed NURBS surfaces, handling extraordinary nodes, based on which surface and volumetric weighted T-spline models are constructed for isogeometric analysis.