В работе представлены результаты выполненного 3D численного термомеханического моделирования процессов формирования гранитоидного Посольненского массива, расположенного в Посольно-Кузеевской тектонической зоне Южно-Енисейского кряжа. Предложена математическая численная модель образования массива за счет диссипации тепла при сдвиговых деформациях в Посольно-Кузеевской тектонической зоне. Моделирование проводилось в 3D-постановке при решении уравнения Навье-Стокса с использованием программного комплекса ANSYS Fluent. Рассмотрены варианты модели строения Посольно-Кузеевской тектонической зоны в отношении ширины активной области горизонтального сдвига и состава пород коры (метаморфический протолит) тектонической зоны. Модель диссипативного нагрева и плавления в сдвиговой зоне шириной 20 км воспроизводит реально наблюдаемую форму Посольненского массива. Определена зависимость высоты подъема магмы от реологических свойств пород коры. Показано, что при диссипативном тепловыделении при деформировании пород, характеризующихся реологией влажного гранита, реализуется плавление и диапировый подъем кислой магмы от основания гранитогнейсового слоя (45 км) до уровня 9-11 км. В модели коры, описываемой реологическими свойствами кислого гранулита, плавление происходит "на месте" в большем объеме, но без дальнейшего всплывания магматической массы. Показано, что скорость сдвига является одним из ключевых параметров эффективности диссипативного нагрева. Определено пороговое значение скорости сдвига - 4 см/год (и выше) для наиболее эффективного плавления и магматизма при диссипативном режиме.