Hi everybody, i am Ravikumar doing M. Tech. project on multi scale analysis of composites.
I gave the input as nodal displacements and different material properties for glass and resin. I applied all the normal stresses and shear stresses on the micro scale model which is of hexagonal array structure. (i.e., one circular fiber at the centre and quarter fibers at the four corners of the model and the remaining portion is filled up by resin).
My aim is to findout all the engineering constants in the compliance matrix by knowing the stresses which i will get from Ansys and Strains (as i am giving as the input).
But, i dont know how can i get the average stress over an element. Because in the paper what i am refering, author mentioned as average stress.
I know stresses at each node and elemental stresses along the perticular direction. From any of these, can i get average stress on perticular element?
The paper states that for average stress, the following equation is used.
Avge. σij-- = (1/V)∫(σ)ij dV
What does the above equation means ?
I thought to multiply the stress of the element in perticular direction with volume of the element and to sum those values for all the elements and devide the result by the total volume of the structure (micro RUC)
Here, i am attaching the reference paper also. Can any one explain me what equation 5 exactly means?
Please help me and Thanks in advance
| Attachment | Size |
|---|---|
| Periodicity boundary conditions.pdf | 553.4 KB |
σij-- denotes
σij-- denotes average stress tensor of the RVE(unit cell). For the calculation of the average stress, please see http://www.imechanica.org/node/9431
In reply to σij-- denotes by Lianhua Ma
Dear Lianhua Ma, Thank
Dear Lianhua Ma,
Thank you very much for your help. I proceed in the way as you told by taking the reactional forces and devide those by the area, it worked well for finding out the normal stresses.
Can you tell me how to solve it for shear stress calculation?
I am applying positive displacement on the top and right planes (strain*lenghts of planes respectively) of RUC (Representative Unit Cell) and negative dispalcement on the left and bottom planes of RUC for finding out the average shear stress on XY plane. (as YZ plane in reference paper. Fig8)
Is it the way to apply the boundary conditions for that?
After solving, i am getting the same deformation mode as in reference as in Fig9.
I tried with multiplying the elemental shear stress (Sxy) with volume of the elements and sum all these values and finally devided the result by the volume of RUC. But, i am getting higher value. I think, it may not be correct.
Can you tell me how to calculate the shear stress?
Once again thank you for your help.......
In reply to Dear Lianhua Ma, Thank by Ravikumar. G
Hi Ravikumar, This
Hi Ravikumar,
This approach mentioned above is suitable for both normal and shear deformation. Please check the prescribed displacement on every edge of your RVE.
For the calcultion of shear stress in x-y plane:
Assuming a, b and c are the edge lengths of your RVE along x，y and z direction, the applied positive displacement on the top edge should be 1/4*b, whereas the displacement prescribed on the right edge should be 1/4*a. Similarly, the displacement applied on the left and bottom edges should be -1/4*a ,-1/4*b, respectively.
Hope this helps
Lianhua
In reply to Hi Ravikumar, This by Lianhua Ma
Hello Lianhua Ma, I am
Hello Lianhua Ma,
I am applying the shear load as you mentioned. But, how can i get the average shear stress over the plane.OR how can i get the average shear stress on XY plane by knowing reaction forces of X and Y on the edges?
Thank you....
In reply to Hello Lianhua Ma, I am by Ravikumar. G
As noted in your post.
As noted in your post. Average shear stress over the domain of RVE= "multiplying the elemental shear stress (Sxy) with volume of the elements and sum
all these values and finally devided the result by the volume of RUC."
In reply to As noted in your post. by Lianhua Ma
Yes....Exactly i am doing
Yes....Exactly i am doing in the same way..... But my results are not matching. I am getting higher values than the author what he mentioned in the paper.
I dont know may be any other method is available to solve for that. If you have any idea please let me know.
Thnak you......
In reply to Yes....Exactly i am doing by Ravikumar. G
You can try another method
You can try another method to solve the x-y shear stress.
Firstly extract the reaction force of every node in a face.
Secondly, sum them together.
And finally divide the result by the area of the face.
In reply to You can try another method by Lianhua Ma
If you use VAMUCH , a
If you use VAMUCH , a general-purpose micromechanics analysis code, you don't have to worry about how to apply loads and boundary conditions. What you need as inputs are the geometry of the unit cell and the constitutive properties of each phases. The 6x6 stiffness and compliance matrices are generated within one analysis. An interface is also developed for ANSYS users to take advantage of its powerful pre/post processing capabilities.