This dissertation developed a method that can accurately and efficiently capture the response of a structure by rigorous combination of a reduced-dimensional beam finite element model with a model based on full two-dimensional (2D) or three-dimensional (3D) finite elements. As a proof of concept, a joint 2D-beam approach is studied for planar-inplane deformation of strip-beams. This approach is developed for obtaining understanding needed to do the joint 3D-beam model. A Matlab code is developed to solve achieve this 2D-beam approach. For joint 2D-beam approach, the static response of a basic 2D-beam model is studied. The whole beam structure is divided into two parts. The root part where the boundary condition is applied is constructed as a 2D model. The free end part is constructed as a beam model. To assemble the two different dimensional model, a transformation matrix is used to achieve deflection continuity or load continuity at the interface. After the transformation matrix from deflection continuity or from load continuity is obtained, the 2D part and the beam part can be assembled together and solved as one linear system. For a joint 3D-beam approach, the static and dynamic response of a basic 3D-beam model is studied. A Fortran program is developed to achieve this 3D-beam approach. For the uniform beam constrained at the root end, similar to the joint 2D-beam analysis, the whole beam structure is divided into two parts. The root part where the boundary condition is applied is constructed as a 3D model. The free end part is constructed as a beam model. To assemble the two different dimensional models, the approach of load continuity at the interface is used to combine the 3D model with beam model. The load continuity at the interface is achieved by stress recovery using the variational-asymptotic method. The beam properties and warping functions required for stress recovery are obtained from VABS constitutive analysis. After the transformation matrix from load continuity is obtained, the 3D part and the beam part can be assembled together and solved as one linear system. For a non-uniform beam example, the whole structure is divided into several parts, where the root end and the non-uniform parts are constructed as 3D models and the uniform parts are constructed as beams. At all the interfaces, the load continuity is used to connect 3D model with beam model. Stress recovery using the variational-asymptotic method is used to achieve the load continuity at all interfaces. For each interface, there is a transformation matrix from load continuity. After we have all the transformation matrices, the 3D parts and the beam parts are assembled together and solved as one linear system.