{"id":407961,"date":"2025-02-01T02:29:41","date_gmt":"2025-02-01T02:29:41","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/error-script-generating-sinusoidal-geometry\/"},"modified":"2025-02-18T16:00:03","modified_gmt":"2025-02-18T16:00:03","slug":"error-script-generating-sinusoidal-geometry","status":"answered","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/error-script-generating-sinusoidal-geometry\/","title":{"rendered":"Error script generating sinusoidal geometry"},"content":{"rendered":"

<p><p>Context: <\/strong>Need to generate a parmetric cylindrical geometry as seen below.<\/p><p><\/p><p>Problem: <\/strong>I generated the points needed for the surface but I don’t know how to join the dots to generate the surface.<\/p><p>Summary of code: <\/strong>Generates x,y an z coordinates of points along surface. Draws curves that links up the surface (this is the best I can do right now). Image below shows this.<br><\/p><p>Help wanted:<\/strong><\/p><p>1) Is there a method or function that can help me connect these points to form a surface?<\/p><p>2) Is there a better way of creating this geometry? E.g. instead of generating x,y and z locations for points and generating a surface you can use the equation feature in space-claim and revolve it (NB: Tried this method – was unable to include record equation in script and as a result unable to parametrise it)<\/p><p>Code so far:<br><\/strong><\/p><p>import math<\/p><p># Define parameters<br>D = 0.06  # Diameter of the cylinder (meters)<br>A = 0.0016  # Amplitude of sinusoidal wave (meters)<br>W = 0.02  # Wavelength of sinusoidal wave (meters)<br>length = 0.15  # Length of the cylinder (meters)<br>segments = 100  # Number of points used to make the structure<br>revolution_steps = 36  # Number of steps for full 360° revolution<\/p><p># Create a list to store the 3D coordinates<br>ptList = List[Point]()<\/p><p># Generate the 3D points by revolving the profile<br>for i in range(segments):<br>    x = i * (length \/ (segments – 1))<br>    y_base = (D \/ 2) + A * math.sin((2 * math.pi \/ W) * x)  # Radius at current x<\/p><p>    # Generate points along the revolution<br>    for j in range(revolution_steps):<br>        angle = (2 * math.pi \/ revolution_steps) * j  # Angle for revolution<br>        y = y_base * math.cos(angle)  # Rotate radius in the yz-plane<br>        z = y_base * math.sin(angle)<br>        <br>        # Create a point and add to the list<br>        ptList.Add(Point.Create(x, y, z))<\/p><p><br># Create the NURBS curve from the points<br>ncurve = NurbsCurve.CreateThroughPoints(False, ptList, 0.0001)<\/p><p># Create a curve segment from the NURBS curve<br>curveSegment = CurveSegment.Create(ncurve)<\/p><p># Create a design curve from the segment<br>designCurve = DesignCurve.Create(GetRootPart(), curveSegment)<\/p><\/p><\/p>\n","protected":false},"template":"","class_list":["post-407961","topic","type-topic","status-answered","hentry","topic-tag-moderator-spaceclaim","topic-tag-scripting","topic-tag-spaceclaim"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t