Skip to content

Ansys mech #105

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
Nick-Hemenway wants to merge 2 commits into
base: develop
Choose a base branch
from
Open

Ansys mech #105

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
0f69daa
first commit to branch. Add AnsysMech class and begin
Nick-Hemenway Sep 20, 2019
3287de4
add documentation for using aas
Nick-Hemenway Jan 3, 2020
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
50 changes: 50 additions & 0 deletions examples/AnsysMechExample.m
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,50 @@
clear

%% Create Cross Section

stator1 = CrossSectOuterRotorStator( ...
'name', 'stator1', ...
'dim_alpha_st', DimDegree(30), ...
'dim_alpha_so', DimDegree((30 / 2) * 0.25), ...
'dim_r_si', DimMillimeter(15), ...
'dim_d_sy', DimMillimeter(7.5), ...
'dim_d_st', DimMillimeter(7.5), ...
'dim_d_sp', DimMillimeter(5), ...
'dim_d_so', DimMillimeter(3), ...
'dim_w_st', DimMillimeter(7.5), ...
'dim_r_st', DimMillimeter(0), ...
'dim_r_sf', DimMillimeter(0), ...
'dim_r_sb', DimMillimeter(0), ...
'dim_Q', 8, ...
'location', Location2D( ...
'anchor_xy', DimMillimeter([0,0]), ...
'theta', DimDegree(0).toRadians() ...
) ...
);

%% Instantiate AnsysMech Tool
scratch_path = "C:\Users\nhemenway\Desktop\APDL_Learning\Ansys_Scratch";

toolMech = AnsysMech();
toolMech.open(scratch_path, 'test_job');
toolMech.reset();
toolMech.preprocesser();

%% Draw Cross Section and Save to IGES File

stator1.draw(toolMech);

%after drawing the cross section we can save it to an iges file that can
%then be opened and viewed with any CAD software
iges_path = 'C:\Users\nhemenway\Desktop\APDL_Learning\stator_test';
toolMech.save_IGES(iges_path);










213 changes: 213 additions & 0 deletions tools/ansysMech/AnsysMech.m
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,213 @@
classdef AnsysMech < ToolBase & DrawerBase
%AnsysMech Encapsulation for the Ansys Mechanical FEA software
% TODO: add more description
% TODO: add more description
% TODO: add more description
% TODO: add more description

properties (GetAccess = 'public', SetAccess = 'private')
iCoMapdlUnit; %Ansys server object
end

methods
function obj = AnsysMech(varargin)
obj = obj.createProps(nargin,varargin);
obj.validateProps();
end

function obj = open(obj, scratch_path, job_name)

%scratch_path is folder where ansys writes to as well as places
%server key

%This code starts up the server from matlab... it is not
%working yet and the server must be started from the mechanical
%apdl launcher with the -aas parameter specified under the
%additional parameters section

exec_path = '"C:\Program Files\ANSYS Inc\v190\ansys\bin\winx64\MAPDL.exe"';
job_flag = sprintf('-j %s', job_name);
misc_flags = ' -p aa_r -np 2 -lch -s read -l en-us -aas -t -d win32';
dir_flag = sprintf('-dir "%s"', scratch_path);

cmd_str = sprintf('%s %s %s %s', exec_path, dir_flag, job_flag, misc_flags);
system(cmd_str); %start ansys in as a server (aaS) mode

orb = initialize_orb(); %this function comes from ansys matlab toolbox
load_ansys_aas; %aas stands for ansys as a server
key_str = sprintf('%s\\aaS_MapdlId.txt', scratch_path);
obj.iCoMapdlUnit = actmapdlserver(orb, key_str);

end

function close(obj)
%this may still not work properly yet
obj.iCoMapdlUnit.executeCommandToString('/EXIT,NOSAVE');
end

function reset(obj)
obj.iCoMapdlUnit.executeCommandToString('FINISH');
obj.iCoMapdlUnit.executeCommandToString('/CLEAR,NOSTART');
end

function preprocesser(obj)
obj.iCoMapdlUnit.executeCommandToString('FINISH');
obj.iCoMapdlUnit.executeCommandToString('/PREP7');
end

function [kp_num] = create_keypoint(obj, x,y,z)

%creates a keypoint at x,y,z and returns the id of the keypoint
%created as an integer. keypoints are numbered automatically

cmd_str = sprintf('K,0,%f,%f,%f', x,y,z); %create APDL command string

out_str = obj.iCoMapdlUnit.executeCommandToString(cmd_str); %returns java.lang.String object
out_str = string(out_str);%conver to matlab string

%extract keypoint number out of string and convert to int type
t = strsplit(strip(out_str));
kp_num = uint64(str2double(t(end)));

end

function [tokenLine] = drawLine(obj, startxy, endxy)

%DRAWLINE Draw a line in the current Ansys APDL document.
% drawLine([start_x, _y], [end_x, _y]) draws a line

kp_start = obj.create_keypoint(startxy(1), startxy(2), 0);
kp_end = obj.create_keypoint(endxy(1), endxy(2), 0);

kp_nums = [kp_start, kp_end];

line_cmd = sprintf('LSTR,%d,%d',kp_nums(1), kp_nums(2));
line_output = obj.iCoMapdlUnit.executeCommandToString(line_cmd);
line_output = string(line_output);

% tokenLine = [kp_nums, line_output];
tokenLine = kp_nums(1);

end

function [tokenArc] = drawArc(obj, centerxy, startxy, endxy)
%DRAWARC Draw an arc in the current Ansys APDL document.
% drawarc([center_x,_y], [start_x, _y], [end_x, _y])
% draws an arc

p_center = obj.create_keypoint(centerxy(1), centerxy(2), 0);
p_zero = obj.create_keypoint(startxy(1), startxy(2), 0);
p_axis = obj.create_keypoint(centerxy(1), centerxy(2), 1);

kp_nums = [p_center, p_zero, p_axis];

v1 = [startxy - centerxy, 0];
v2 = [endxy - centerxy, 0];

c = cross(v1,v2);
deg = sign(c(3))*atan2d(norm(c),dot(v1,v2));

if deg <0
deg = deg + 360;
end

rad = norm(v1);

arc_cmd = sprintf('CIRCLE,%d,%f,%d,%d,%f', p_center, rad, p_axis, p_zero, deg);
arc_output = obj.iCoMapdlUnit.executeCommandToString(arc_cmd);
arc_output = string(arc_output);

% tokenArc = [kp_nums, arc_output];
tokenArc = kp_nums(1);

end


function [save_output] = save_IGES(obj, path)

save_cmd = sprintf("IGESOUT,'%s','IGES',,1", path);
save_output = obj.iCoMapdlUnit.executeCommandToString(save_cmd);
save_output = string(save_output);

end

function select(obj)
%SELECT Selects something from canvas (?)
% select()

% TODO:
% Implement this...
%
% This will need to take in arguments, or maybe
% CrossSect objects which then store internally all their
% lines and surfaces that need to be selected
end

function new = revolve(obj, name, material, center, axis, angle, token)
%REVOLVE Revolve a cross-section along an arc
%new = revolve(obj, name, material, center, axis, angle)
% name - name of the newly extruded component
% center - x,y coordinate of center point of rotation
% axis - x,y coordinate on the axis of ration (negative reverses
% direction) (0, -1) to rotate clockwise about the y axis
% angle - Angle of rotation (dimAngular)

end

function new = prepareSection(obj, csToken)

%TO DO: implement function if necessary

end

function setDefaultLengthUnit(obj, userUnit, makeAppDefault)
%SETDEFAULTLENGTHUNIT Set the default unit for length.
% setDefaultLengthUnit(userUnit, makeAppDefault)
% Sets the units for length.

% userUnit can be one of these options:
% 'kilometers'
% 'meters'
% 'centimeters'
% 'millimeters'
% 'microns'
% 'miles'
% 'yards'
% 'feet'
% 'inches'
%
% This is a wrapper for Document::setDefaultLengthUnit

%TO DO: implement if necessary

end

function viewAll(obj)
%TO DO: implement if necessary
end


function setVisibility(obj, visibility)
%TO DO: implement if necessary
end
end

methods(Access = protected)
function validateProps(obj)
%VALIDATE_PROPS Validate the properties of this component

% Use the superclass method to validate the properties
validateProps@ToolBase(obj);
validateProps@DrawerBase(obj);
end

function obj = createProps(obj, len, args)
%CREATE_PROPS Add support for value pair constructor

validateattributes(len, {'numeric'}, {'even'});
for i = 1:2:len
obj.(args{i}) = args{i+1};
end
end
end
end
Binary file added tools/ansysMech/AnsysMechExample.docx
View file
Open in desktop
Binary file not shown.