基于UIDL和APDL的ANSYS二次开发实例：ANSYS-SwiftComp GUI

banghuazhao

本文介绍一个基于UIDL和APDL的ANSYS二次开发实例：ANSYS-SwiftComp GUI。SwiftComp是AnalySwift公司的软件，可以用于复合材料仿真。
本文运用UIDL和APDL，在ANSYS里为SwiftComp创建了用户界面，如下图所示：（新创建的菜单由红色方框框出）


新创新的菜单有：Common SGs, Homogenization, Dehomogenization
它们在Main Menu里的显示由UIDL控制，对应的UIDL文件如下图所示：


新创建的功能有：1D SG: Fast Generate, Advanced Generate, Input File. 2D SG: Beam Section, Square Pack, Hexagonal Pack. 3D SG: Square Pack, Spherical Inclusion, Honeycomb. Homogenization: Beam Model, Plate/Shell Model, Solid Model. Dehomogenization: Beam Model, Plate/Shell Model, Solid Model. 它们功能的实现由APDL控制，对应的APDL文件如下图所示：


ANSYS-SwiftComp GUI可以用于多尺度结构仿真，对复合材料仿真特别合适。
这个二次开发所有的文件和说明，均可以在cdmHUB上免费下载（需要注册）：https://cdmhub.org/resources/1136
ANSYS-SwiftComp GUI的介绍和实例可以在我的领英上查看：https://www.linkedin.com/pulse/ansys-swiftcomp-gui-ansys-plugin-multiscale-modeling-banghua-zhao
下面通过一个实例介绍ANSYS-SwiftComp GUI：
Example Problem
ANSYS-SwiftComp GUI in principle can be used for multiscale constitutive modeling of any structures. We are using a simple example to demonstrate some of its features. As shown in Figure 2, an artificial, unidirectional fiber reinforced composite is made of graphite fiber and epoxy matrix. The fiber can be assumed to be transversely isotropic with E1 = 276 GPa, E2 = E3 = 19.5 GPa, G12 = G13 = 70 GPa, ν12 = ν13 = 0.28, ν23 = 0.70. The epoxy matrix can be assumed to be isotropic with E = 4.76 GPa and ν = 0.37. A laminate is made of a 0-degree ply on the bottom and a 90-degree ply on the top. Assuming the bottom layer contains 10 fibers and the top layer contains 50 fibers. The microstructure can be assumed as square packing with fiber volume fraction equal to 60%. Note the lines between unit cells in Figure 2 are artificial and they are used to show the unit cell. The length of the laminate is 25 mm, and the width is 5 mm, and the thickness of each layer is 0.5 mm. It is clamped at one end and free at the other end. It is subject to 100 Pa pressure on the top surface.


Figure 2 Example problem

ANSYS-SwiftComp GUI is used to solve this problem. There are two ways to attack this problem by MSG and its companion code SwiftComp™: Classical Plate Model and Classical Beam Model. For each way, the detail steps of modeling and analysis are given.

Classical Plate Model:

Using MSG, the original problem is decoupled, as shown in Figure 3, into a constructive modeling over a 3D SG (left) and a structural analysis of a plate (right).





Figure 3 3D SG and structural analysis for classical plate model

The detailed steps are given in the following:

First, create the 3D SG.

Then, use ANSYS-SwiftComp GUI to performed homogenization: Solution → Homogenization → Plate/Shell Model. Use default parameter. A screenshot is shown in Figure 4.



Figure 4 Homogenization for Plate/Shell Model

Next, create plate model and read result from homogenization to perform the plate analysis

Next, use ANSYS-SwiftComp GUI to performed dehomogenization: Solution → Dehomogenization → Plate/Shell Model. Input global response from the previous step. A screenshot is shown in Figure 5.




Figure 5 Dehomogenization for Plate/Shell Model

The contour plot for the nodal solution of Von Mises stress result is shown in Figure 6.



Figure 6 Contour Plot of Nodal Von Mises Stress

Classical Beam Model:

Using MSG, the original problem is decoupled, as shown in Figure 7, into a constructive modeling over the 3D SG (left) and a structural analysis of a beam (right).



Figure 7 3D SG and structural analysis for classical beam model

The detailed steps are given in the following:

First, create the 3D SG.

Then, use ANSYS-SwiftComp GUI to performed homogenization: Solution → Homogenization → Beam Model. Use default parameter. A screenshot is shown in Figure 8.



Figure 8 Homogenization for Beam Model

Next, create beam model and read the result from homogenization to perform global beam analysis.

Next, use ANSYS-SwiftComp GUI to performed dehomogenization: Solution → Dehomogenization → Beam Model. Input global response from the previous step. A screenshot is shown in Figure 9.



Figure 9 Dehomogenization for Beam Model

The contour plot for the nodal solution of Von Mises stress result is shown in Figure 10.




Figure 10 Contour Plot of Nodal Von Mises Stress

Results and Discussion

The comparisons of stress sigma11, sigma22 and sigma33 (along the path (12.75, 0.25, x3) in DNS) are shown from Figure 11 to Figure 13.

As we can see, both MSG plate and MSG beam agree very well with DNS.

Note, the current model for DNS has around 8 million DOFs, which needs 1 day for calculation with 16 CPUs. However, both MSG plate and MSG beam only take less than a minute for computation.


Figure 11 Comparison of sigma11 along path (12.75, 0.25, x3)


Figure 12 Comparison of sigma22 along path (12.75, 0.25, x3)


Figure 13 Comparison of sigma33 along path (12.75, 0.25, x3)



如果大家觉得这个二次开发有帮助，请给我的领英文章点个赞吧：https://www.linkedin.com/pulse/ansys-swiftcomp-gui-ansys-plugin-multiscale-modeling-banghua-zhao
大家有什么问题，欢迎和我讨论，我会及时予以回答，问题也可以提在这里：https://cdmhub.org/groups/yugroup/forum