Логотип Викисклада

Файл:Rhombic triacontahedron box.stl

Материал из Википедии — свободной энциклопедии
Перейти к навигации Перейти к поиску

Исходный файл(5120 × 2880 пикселей, размер файла: 14 Кб, MIME-тип: application/sla)

View Rhombic triacontahedron box.stl  on viewstl.com

Краткое описание

Описание
English: A rhombic triacontahedron made of 6 panels around a cubic hole by CMG Lee. Zoom into the model to see the cubic hole.
Дата
Источник собственная работа
Автор Cmglee

Python source

#!/usr/bin/env python

header = 'A rhombic triacontahedron made of 6 panels around a cubic hole by CMG Lee.'
a      = 233 ## half rhombus minor diagonal
b      = 377 ## half rhombus major diagonal
c      =  55 ## half peg base width
d      =  30 ## displacement of panels

import re, struct, math
def fmt(string): ## string.format(**vars()) using tags {expression!format} by CMG Lee
 def f(tag): i_sep = tag.rfind('!'); return (re.sub('\.0+$', '', str(eval(tag[1:-1])))
  if (i_sep < 0) else ('{:%s}' % tag[i_sep + 1:-1]).format(eval(tag[1:i_sep])))
 return (re.sub(r'(?<!{){[^{}]+}', lambda m:f(m.group()), string)
         .replace('{{', '{').replace('}}', '}'))
def append(obj, string): return obj.append(fmt(string))
def tabbify(cellss, separator='|'):
 cellpadss = [list(rows) + [''] * (len(max(cellss, key=len)) - len(rows)) for rows in cellss]
 fmts = ['%%%ds' % (max([len(str(cell)) for cell in cols])) for cols in zip(*cellpadss)]
 return '\n'.join([separator.join(fmts) % tuple(rows) for rows in cellpadss])
def hex_rgb(colour): ## convert [#]RGB to #RRGGBB and [#]RRGGBB to #RRGGBB
 return '#%s' % (colour if len(colour) > 4 else ''.join([c * 2 for c in colour])).lstrip('#')
def viscam_colour(colour):
 colour_hex      = hex_rgb(colour)
 colour_top5bits = [int(colour_hex[i:i+2], 16) >> 3 for i in range(1,7,2)]
 return (1 << 15) + (colour_top5bits[0] << 10) + (colour_top5bits[1] << 5) + colour_top5bits[2]
def roundm(x, multiple=1):
 if   (isinstance(x, tuple)): return tuple(roundm(list(x), multiple))
 elif (isinstance(x, list )): return [roundm(x_i, multiple) for x_i in x]
 else: return int(math.floor(float(x) / multiple + 0.5)) * multiple
def flatten(lss): return [l for ls in lss for l in ls]
def rotate(facetss, degs): ## around x then y then z axes
 (deg_x,deg_y,deg_z) = degs
 (sin_x,cos_x) = (math.sin(math.radians(deg_x)), math.cos(math.radians(deg_x)))
 (sin_y,cos_y) = (math.sin(math.radians(deg_y)), math.cos(math.radians(deg_y)))
 (sin_z,cos_z) = (math.sin(math.radians(deg_z)), math.cos(math.radians(deg_z)))
 facet_rotatess = []
 for facets in facetss:
  facet_rotates = []
  for i_point in range(4):
   (x, y, z) = [facets[3 * i_point + i_xyz] for i_xyz in range(3)]
   if (x is None or y is None or z is None):
    facet_rotates += [x, y, z]
   else:
    (y, z) = (y * cos_x - z * sin_x,  y * sin_x + z * cos_x) ## rotate about x
    (x, z) = (x * cos_y + z * sin_y, -x * sin_y + z * cos_y) ## rotate about y
    (x, y) = (x * cos_z - y * sin_z,  x * sin_z + y * cos_z) ## rotate about z
    facet_rotates += [round(value, 9) for value in [x, y, z]]
  facet_rotatess.append(facet_rotates)
 return facet_rotatess
def translate(facetss, ds): ## ds = (dx,dy,dz)
 return [facets[:3] + [facets[3 * i_point + i_xyz] + ds[i_xyz]
                       for i_point in range(1,4) for i_xyz in range(3)]
         for facets in facetss]
def flip(facetss):
 return [facets[:3] + facets[6:9] + facets[3:6] + facets[9:] for facets in facetss]


## Add facets
facet_pegss   = [[None,0,0, b+c*2, 0,-b,    b+c,-c,0,    b+c, c,0], ## peg
                 [None,0,0, b+c*2, 0,-b,    b+c, c,0,  b+c*3, 0,0], ## peg
                 [None,0,0, b+c*2, 0,-b,  b+c*3, 0,0,    b+c,-c,0], ## peg
                 [None,0,0,     b, b, 0,    a+b, 0,0,  b+c*3, 0,0], ## bottom
                 [None,0,0,     b, b, 0,  b+c*3, 0,0,    b+c, c,0], ## bottom
                 [None,0,0,     b, b, 0,    b+c, c,0,    b+c,-c,0], ## bottom
                 [None,0,0,     b,-b, 0,  b+c*3, 0,0,    a+b, 0,0], ## bottom
                 [None,0,0,     b,-b, 0,    b+c,-c,0,  b+c*3, 0,0], ## bottom
                 [None,0,0,     b,-b, 0,      b, b,0,    b+c,-c,0]] ## bottom
facet_panelss = translate([[None,0,0, a+b, 0,0,    b, b,0,    0, b,a], ## slant N
                           [None,0,0,   0, b,a,    a, 0,a,  a+b, 0,0], ## slant N
                           [None,0,0, a+b, 0,0,    a, 0,a,    0,-b,a], ## slant S
                           [None,0,0,   0,-b,a,    b,-b,0,  a+b, 0,0], ## slant S
                           # [None,0,0,   a, 0,a,    0, b,a,   -a, 0,a], ## top round
                           [None,0,0,   0,-b,a,    a, 0,a,    0, b,a], ## top thin
                           # [None,0,0,  -b, b,0,    0, b,a,    b, b,0], ## side without hole
                           # [None,0,0,   b, b,0,  a+b, 0,0,    b,-b,0], ## bottom end without peg
                           [None,0,0,  -b, b,0,    b, b,0,    b,-b,0]] + ## bottom mid
                          facet_pegss + translate(rotate(flip(facet_pegss),(0,-90,-90)),(0,b,-b)),
                          (0,0,b + d))
facet_panelss += rotate(facet_panelss, (  0,0,180))
facet_panelss += rotate(facet_panelss, (180,0,180))
facetss = facet_panelss[:] + rotate(facet_panelss, (90,90,0)) + rotate(facet_panelss, (90,0,90))

## Calculate normals
for facets in facetss:
 if (facets[0] is None or facets[1] is None or facets[2] is None):
  us      = [facets[i_xyz + 9] - facets[i_xyz + 6] for i_xyz in range(3)]
  vs      = [facets[i_xyz + 6] - facets[i_xyz + 3] for i_xyz in range(3)]
  normals = [us[1]*vs[2] - us[2]*vs[1], us[2]*vs[0] - us[0]*vs[2], us[0]*vs[1] - us[1]*vs[0]]
  normal_length = sum([component * component for component in normals]) ** 0.5
  facets[:3] = [-round(component / normal_length, 10) for component in normals]

print(tabbify([['N%s'  % (xyz   )                   for xyz in list('xyz')] +
               ['%s%d' % (xyz, n) for n in range(3) for xyz in list('XYZ')] + ['RGB']] + facetss))
## Compile STL
outss = ([[('STL\n\n%-73s\n\n' % (header[:73])).encode('utf-8'), struct.pack('<L',len(facetss))]] +
         [[struct.pack('<f',float(value)) for value in facets[:12]] +
          [struct.pack('<H',0 if (len(facets) <= 12) else
                            viscam_colour(facets[12]))] for facets in facetss])
out   = b''.join([bytes(out) for outs in outss for out in outs])
# out += ('\n\n## Python script to generate STL\n\n%s\n' % (open(__file__).read())).encode('utf-8')
print("# bytes:%d\t# facets:%d\ttitle:\"%-73s\"" % (len(out), len(facetss), header[:73]))
with open(__file__[:__file__.rfind('.')] + '.stl', 'wb') as f_out: f_out.write(out)

Лицензирование

Я, владелец авторских прав на это произведение, добровольно публикую его на условиях следующей лицензии:
w:ru:Creative Commons

атрибуция распространение на тех же условиях

Этот файл доступен по лицензии Creative Commons Attribution-Share Alike 4.0 International
Вы можете свободно:
  • делиться произведением – копировать, распространять и передавать данное произведение.
  • создавать производные – переделывать данное произведение
При соблюдении следующих условий:
  • атрибуция – Вы должны указать авторство, предоставить ссылку на лицензию и указать, внёс ли автор какие-либо изменения. Это можно сделать любым разумным способом, но не создавая впечатление, что лицензиат поддерживает вас или использование вами данного произведения.
  • распространение на тех же условиях – Если вы изменяете, преобразуете или создаёте иное произведение на основе данного, то обязаны использовать лицензию исходного произведения или лицензию, совместимую с исходной.

Wikimedia Foundation
The uploader of this file has agreed to the Wikimedia Foundation 3D patent license: This file and any 3D objects depicted in the file are both my own work. I hereby grant to each user, maker, or distributor of the object depicted in the file a worldwide, royalty-free, fully-paid-up, nonexclusive, irrevocable and perpetual license at no additional cost under any patent or patent application I own now or in the future, to make, have made, use, offer to sell, sell, import, and distribute this file and any 3D objects depicted in the file that would otherwise infringe any claims of any patents I hold now or in the future.

Please note that in the event of any differences in meaning or interpretation between the original English version of this license and a translation, the original English version takes precedence.

История файла

Нажмите на дату/время, чтобы посмотреть файл, который был загружен в тот момент.

Дата/времяМиниатюраРазмерыУчастникПримечание
текущий18:34, 4 апреля 2018Миниатюра для версии от 18:34, 4 апреля 20185120 × 2880 (14 Кб)CmgleeAdd pegs to hold 3D-printed model together.
16:08, 2 апреля 2018Миниатюра для версии от 16:08, 2 апреля 20185120 × 2880 (5 Кб)CmgleeUser created page with UploadWizard

Нет страниц, использующих этот файл.

Глобальное использование файла

Данный файл используется в следующих вики: