This paper presents a modified method of discontinuous deformation analysis(DDA). In the presented method, open-close iteration may not be needed, small penetration is permitted among blocks, and springs are added bet...This paper presents a modified method of discontinuous deformation analysis(DDA). In the presented method, open-close iteration may not be needed, small penetration is permitted among blocks, and springs are added between contacting block pairs only when a penetration takes place. The three contact patterns(i.e. sliding, locking and opening) in original DDA method are not involved, and the recognition of these contact patterns and treatment of transformation among patterns are not required either,signi fi cantly saving the computing time. In a convex to concave contact, there are two candidate entrance edges which may cause uncertainty. In this case, we propose the angle bisector criterion to determine the entrance edge. The spring stiffness is much larger than Young's modulus in the original DDA, however we fi nd that the correct results can still be obtained when it is much smaller than Young's modulus. Finally,the penetrations by using penalty method and augmented Lagrangian method are compared. Penetration of the latter is 1/4 of the former. The range of spring stiffness for the latter is wider than the former,being 0.01-1 of the former. Both methods can lead to correct contact forces.展开更多
基金supported by CRSRI Open Research Program (No.CKWV2014206/KY)the National Basic Research Program of China (No.2011CB710603)
文摘This paper presents a modified method of discontinuous deformation analysis(DDA). In the presented method, open-close iteration may not be needed, small penetration is permitted among blocks, and springs are added between contacting block pairs only when a penetration takes place. The three contact patterns(i.e. sliding, locking and opening) in original DDA method are not involved, and the recognition of these contact patterns and treatment of transformation among patterns are not required either,signi fi cantly saving the computing time. In a convex to concave contact, there are two candidate entrance edges which may cause uncertainty. In this case, we propose the angle bisector criterion to determine the entrance edge. The spring stiffness is much larger than Young's modulus in the original DDA, however we fi nd that the correct results can still be obtained when it is much smaller than Young's modulus. Finally,the penetrations by using penalty method and augmented Lagrangian method are compared. Penetration of the latter is 1/4 of the former. The range of spring stiffness for the latter is wider than the former,being 0.01-1 of the former. Both methods can lead to correct contact forces.
