С использованием разработанного ранее набора межатомных потенциалов для атомистического моделирования простых и сложных оксидов бериллия, который почти точно воспроизвел их структурные, упругие и термодинамические свойства, были проведены расчеты термодинамических и упругих свойств смешения и сделан анализ локальной структуры твердых растворов Be(Al,Cr,Fe^III)₂O₄. Расчеты проводились для сверхъячейки 4а x 2b x 2c структурного типа оливина со снятой нетрансляционной симметрией (пр.гр. P₁), что позволило задавать распределение трехвалентных атомов по неэквивалентным позициям и обеспечить релаксацию локальной структуры. Свойства смешения рассчитаны во всем диапазоне составов, сделана оценка изменений свободной энергии Гиббса с температурой и областей стабильности твердого раствора, получены значения критической температуры и состава. Были построены гистограммы межатомных расстояний Me-Me, Me-O и O-O, объемов октаэдров МеO₆, а также оценены значения податливостей катионных позиций и определены группы наиболее сдвигаемых из идеальных позиций атомов.