On 02/04/2012 02:27 AM, Tom Aldcroft wrote:
On Fri, Feb 3, 2012 at 7:35 PM, Robert Cimrmancimr...@ntc.zcu.cz wrote:
Hi Tom,
On Fri, 3 Feb 2012, Tom Aldcroft wrote:
Hi -
I'm trying to specify a localized traction pressure which is a providing a radial outward force over a square patch in my adjustable x-ray optics model. A schematic picture is shown in the link below, where the dots represent nodes, and red dots with arrows are forces:
http://hea-www.harvard.edu/~aldcroft/tmp/force_plot.png
This is a 2-d plot, but in reality the nodes are on a 3-d thin shell cylinder, so the forces won't all lie in the same plane. I've looked at the traction pressure example but I am unable to extrapolate from that example to specify my problem.
It should be possible to define such a load - I will prepare an example showing that, with dw_surface_ltr and/or dw_point_load (with a tweak) terms.
So here it is. In plate_load.py, there is get_load_vectors() that defines the load towards the patch centre depending on normals. For dw_point_load (in this file) everything is defined in mesh nodes. For dw_surface_ltr (in plate_load2.py) the normals and load vectors are defined in surface quadrature points. I have used some ad-hoc values for the material parameters, etc. so the examples are by no means physical, but they should get you started. They use some new utility functions that were just added - use the latest git version to try them.
The mesh has been generated using
./script/blockgen.py -d[1,1,0.01] -s[21,21,3] -o plate.mesh
- if you try a different one, do not forget to change the region definitions accordingly.
Note that after we add the plate/shell elements, the load will be given in a slightly different way, depending on the model used.
Let us know how it goes.
This sounds great. I'll be presenting the results of these calculations at the San Diego SPIE meeting in August, so FWIW your efforts will get a viewgraph or two in my talk and reference in the conference paper.
That would be great, thanks!
As an aside, it's quite possible I'm going about this all wrong. Recall that I'm trying to simulate a piezoelectric patch mounted on a glass thin walled cylinder. I think that the radial force input represents the situation when the piezo is actuated (but maybe not). Certainly this can be modeled as a thermo-elastic problem, with two materials: the cylinder with a CTE of 0.0 and the piezo material with a non-zero CTE. Then change the temperature of the whole system and calculate the resultant displacements. I don't see any examples that seem to show the thermo-elastic equations.
This is certainly possible as well, but imho defining the right thermal expansion coefficient is difficult, no? Maybe it would be more feasible to model directly the piezo effect? Anyway, I think I have a thermo-elasticity example somewhere, I will send it to you...
There is also new examples/thermo_elasticity/thermo_elasticity.py, in case you would like to explore this way.
r.